function ff = f(t,y)
%
  ff = t^2*y;