function y = solution(t)
y = 4/1.3*(exp(0.8*t)-exp(-0.5*t))+2*exp(-0.5*t);