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