function exForwardEuler(y0,tfinal,nsteps)
%
  [y,t] = forwardEuler(y0,tfinal,nsteps);
  ye = y0*exp(t.^3/3);
  plot(t,[y' ye']);
  axis([0 tfinal y0 y0*exp(tfinal^3/3)]);