Numerical approximation to Poisson’s equation over the square [a,b]x[a,b] with Dirichlet boundary conditions. Uses a uniform mesh with (n+2)x(n+2) total points (i.e, n x n interior grid points).