FIG. 5. Solution error at t = 1 using Crank–Nicolson (stars) and TGA (circles), for the Neumann problem for (19) with fixed boundaries. Left-hand plot shows max norm, right-hand plot shows 1-norm. We see that both jj jj1 = O(h2) and jj jj1 = O(h2), indicating second-order accuracy.
...read more