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

