numerical solution of poisson equation