function y = sol_exact(t) %y=6-2*exp(-t)-3*exp(-0.5*t); y=exp(t-2/3*t.^2); end