
function z=fun(y,t)
%this is function f(t,y)=y-t^2+ 1

z=y-t.^2 +1;
