TwoDEscape:

Path: Interplanetary/InterplanetaryOrbit

% Simulate a low thrust departure from a planet
 Treat the problem as planar. Creates a plot of the simulated spiral 
 and a time history plot of the delta-V, fuel consumed, and orbit
 eccentricity. You can add radiation belts to any planet, though not all
 have radiation belts. The spiral into orbit is the same.
--------------------------------------------------------------------------
   Form:
   [mX, xP, t] = TwoDEscape( mu, rI, uE, thrust, m0, mF, rB, planet )
--------------------------------------------------------------------------

  Inputs
  ------
  mu     (1,1) Gravitational parameter
  rI     (1,1) Initial radius (km)
  uE     (1,1) Exhaust velocity (m/s)
  thrust (1,1) Thrust (N)
  m0     (1,1) Dry mass  (kg)
  mF     (1,1) Fuel mass  (kg)
  rB     (1,2) Radiation belt start/stop
  planet (1,:) Planet name

  Outputs
  -------
  mX     (1,1) Fuel consumed (kg)
  xP     (6,n) State vector [r;v;m;DV]
  t      (1,n) Time vector corresponding to xP (s)

--------------------------------------------------------------------------