Electrostatics through simulation

Here are the results of solving the poisson equation in 2D boundary: I wrote this code for our class which will go opensource book soon-involving the introduction of python in physics teaching My solver uses fft(fast fourier trnasform) hence its pretty much fast. Enjoy the Pics.Write to me if you find problems,bugs,or dont agree on a thing : 1.Arbitrary charge distribution:

2.Dipole Charge distribution

3.Quardapole Charge distribution.

4. Unstable twi dipole side by side.