function deriv = f(t,y);
deriv = t*exp(3.0*t) - 2.0*y;