function y = true_sol(t)
y = 2*exp(-t) + t - 1;