function exact = f(t);
exact = 0.2*t*exp(3.0*t) - 0.04*exp(3.0*t) + 0.04*exp(-2.0*t);
