RK45:

-------------------------------------------------------------------------------
   Fourth/fifth order Runge-Kutta. 
   Called function is of the form:

   Fun(x,t,p1,p2,...)

   Accepts optional arguments that are passed through to Fun.
   Time is also optional.

   This function will integrate Fun from the current t to t + hMax.
 
-------------------------------------------------------------------------------
   Form:
   [x, hLast] = RK45( Fun, x, h, hMax, hMin, tol, t, varargin )
-------------------------------------------------------------------------------

   ------
   Inputs
   ------
   Fun                Function    Fun(x,t,p1,p2...)
   x                  State (column vector)
   h                  Independent variable step
   hMax               Maximum step size
   hMin               Minimum step size
   tol                Tolerance on error
   t                  Current time
   varargin           Optional arguments

   -------
   Outputs
   -------
   x                  Updated state
   hLast              Independent variable step

-------------------------------------------------------------------------------
	Reference: Cash, J.R., A. H. Karp, "A Variable Order Runge-Kutta Method for 
              Initial Value Problems with Rapidly Varying Right-Hand Sides,"
              ACM Trans. on Math. Soft., Vol. 16, No.3, Sept. 1990, pp 201-222.
-------------------------------------------------------------------------------

Children: