function fty = ty4(t,y)
%
% test function for a stiff ODE
%
fty =  5 *exp(5*t).* ((y-t).^2)+ 1;