2D Poisson equation

Solving the 2D Poisson equation iteratively, using the 5-point finite difference stencil

The 2D Poisson equation is solved in an iterative manner (number of iterations is to be specified) on a square 2x2 domain using the standard 5-point stencil. Homogenous neumann boundary conditions have been used.