I have to create a charge distribution of circular shape. I already created a square charge distribution and it look like this. How can I change this square distribution to circular?
My code for the square distribution:
Nmax = 120; Niter = 10
V = zeros((Nmax, Nmax+100) , float)
for k in range(0, 40): V[40+k,40] = 2400.0
for k in range(0, 40): V[40,k+40] = 2400.0
for k in range(0, 40): V[80,k+40] = 2400.0
for k in range(0, 40): V[40+k,80] = 2400.0