
function zs = fun(t,y)
%-------------------
zs = 1 - t + 4*y;

