function[z] = fun(t,y)
z = 5*exp(5*t)*((y-t)^2) + 1;
end

