function deriv = f(t,y);
deriv = 4*t*cos(y);
%deriv = 20.0*y;
