Delta V required to change the inclination
[dv,anv] = CL_man_inclination(ai,ei,inci,pomi,incf [,posman,mu])
Computes the velocity increment needed to change the inclination.
The maneuver consists in making the orbit rotate around the line of nodes, thus changing the velocity direction at one of the nodes.
The output argument dv is the velocity increment vector expressed in spherical coordinates in the QSW local frame: [lambda; phi; |dv|], where lambda is the in-plane angle (+%pi: towards planet and +%pi/2: ~along velocity), phi is the out-of-plane angle, positive towards the angular momentum vector (the angular momentum vector is perpendicular to the orbit plane and oriented according to right hand convention), dv is the norm of the velocity increment.
anv is the true anomaly at the maneuver position.
The optional flag posman can be used to define the position of the maneuver (1->ascending node, -1->descending node, 0->minimal |DV|)
Semi major axis [m] (1xN)
Eccentricity (1xN)
Initial inclination [rad] (1xN)
Argument of periapsis [rad] (1xN)
Final inclination [rad] (1xN)
(optional) Flag specifying the position of maneuver (1->ascending node, -1->descending node, 0->minimal |dv|). Default is ascending node) (1xN)
(optional) Gravitational constant. [m^3/s^2] (default value is %CL_mu)
velocity increment vector in spherical coordinates in the QSW frame [lambda;phi;|dv|] [rad,rad,m/s] (3xN)
true anomaly at maneuver position [rad] (1xN)
CNES - DCT/SB