function y = f( t,u)
    %y=3+exp(-t)-0.5*u;
y=(1-4/3*t)*u;
end