function fty = ty4dy(t,y)
%
% derivative of test function ty4 for a stiff ODE
%
fty =  10 *exp(5*t).* (y-t);


