NKPlanetSim:

% Simulate a non-Keplerian orbit in a cylindrical planet-centered frame.
 The sail enables an orbit about a point in the anti-sun direction.
 The orbit period is chosen to be the same as a Kepler orbit of radius r.

 The initial orbit is unstable. A quadratic regulator in pitch is the added
 to control the orbit.

   Since version 7.
  ------------------------------------------------------------------------
   Reference: Colin R. McInnes, "Solar Sailing: Technology, Dynamics and Mission  
              Applications", Springer-Praxis, London, 1999, Fig. 5.18 on p. 207
  ------------------------------------------------------------------------
  See also NonKeplerianPlanet., QCR, S2Damp, Constant, InformDlg, Plot2D, 
  Plot3D, TimeLabl, Accel, FPlanetControl, FRotatingPlanet, VarEqSailPlanet, 
  Cyl2Cart
  ------------------------------------------------------------------------
%
-------------------------------------------------------------------------------
    Copyright (c) 2007 Princeton Satellite Systems, Inc.
    All rights reserved.
-------------------------------------------------------------------------------

Children:

Common: Control/QCR
Common: Control/S2Damp
Common: Database/Constant
Common: Graphics/InformDlg
Common: Graphics/Plot2D
Common: Graphics/Plot3D
Common: Graphics/TimeLabl
Sail: Analysis/NonKeplerianPlanet
Sail: Control/FPlanetControl
Sail: OrbitDynamics/FRotatingPlanet
Sail: OrbitDynamics/VarEqSailPlanet
Sail: Utilities/Cyl2Cart