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. 32 more words