Suppose we want to solve the two dimensional equation :
The exact solution of this equation is given by : e2tsin(x+y)cosh(x+y).
We could use two maps of isosurfaces: one for the calculated solution, another for the exact solution. We will also generate a plot of the grid nodes, with the initial conditions. The initial and the four boundary conditions are chosen to be the exact solution.