FSailEarthSun:

--------------------------------------------------------------------------
   Sail dynamics in the Earth-Sun three-body system.
   Uses SunV1 for ephemeris.

   Designed to be called by ode113. 

   Since version 7.
--------------------------------------------------------------------------
   Form:
   varargout = FSailEarthSun( t, x, flag, d )
--------------------------------------------------------------------------

   ------
   Inputs
   ------
   t                      Time (sec)
   x             (:,1)    The state vector [r;v]
   d             (1,1)    Data structure
                          .mass  (1,1) Mass
                          .force (3,1) External force
                          .earth (1,1)
                                      .mu (1.1) Gravitational potential
                                      .nZ (1.1) Number of zonal harmonics
                                      .nT (1.1) Number of tesseral harmonics
                                      .s  (:,:) S terms
                                      .c  (:,:) C terms
                                      .j  (:)   m = 0 terms
                                      .a  (1.1) Planet radius
                          .muSun (1,1) Sun gravitational parameter
                          .jDStart (1,1) Initial Julian date
                          .accel (1,1) Sail characteristic acceleration (km/s2)
                          .theta (1,1) Sail pitch angle

   -------
   Outputs
   -------
   vargout          (:,:)     Depends on flag

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

Children:

Math: Linear/Mag
SC: Ephem/SunV1