Reference solutions correct to within a tolerance of 2e-14.
Computed with mono-adaptive cG(3) and a fixed time step of
size k = 0.00005.
