SailClosedLoopNLEqns:

--------------------------------------------------------------------------
   Nonlinear sail RHS with area control.

   Since version 7.
--------------------------------------------------------------------------
   Form:
   xd = SailClosedLoopNLEqns(t,x,kappa0,rho0,z0,G,mode)
--------------------------------------------------------------------------

   ------
   Inputs
   ------
   t            (1,1)     Time
   x            (6,1)     The state vector [rho;z;rhoDot;zDot;theta;thetaDot]
   kappa0       (1,1)     Sail characteristic acceleration
   rho0         (1,1)     Target orbit amplitude
   z0           (1,1)     Target orbit displacements
   G            (4,1)     Control gains
   mode         (1,1)     Control mode, (1) for area and (2) for pitch

   -------
   Outputs
   -------
   xd           (6,1)     State derivative

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