function y = sol_de_1(t)

y = 1+0.5*( exp(-4*t)-exp(-2*t) );