Poisson's equation with scipy.sparse: two plates in a grounded box

Afterwards you can build the five-point Laplacian as a scipy.sparse matrix with kron, fix boundary values, solve with spsolve, and check it on a known case.