CL_mod_equinoxesEquation — Equation of the equinoxes (IERS 1996)
[equi[,equip]]=CL_mod_equinoxesEquation(dPsi,eps0,F5[,dPsip,epsOp,F5p])
longitude (see CL_mod_nutationAngles) [rad](1xN)
(optional with eps0p and F5p) longitude first time derivative (see CL_mod_nutationAngles) [rad/s] (1xN)
mean obliquity (see CL_mod_meanObliquity) [rad] (1xN)
(optional with dPsip and F5p) mean obliquity first time derivative (see CL_mod_meanObliquity) [rad/s] (1xN)
Mean longitude of Moon's ascending node (see CL_mod_nutationArg) [rad] (1xN)
(optional with eps0p and dPsip) Mean longitude of Moon's ascending node first time derivative (see CL_mod_nutationArg) [rad/s] (1xN)
equinoxes equation (1xN)
(optional, needs dPsip, eps0p and F5p) first time derivative of equi (1xN)
1 IERS Conventions (1996), Dennis D. McCarthy
2 Explanatory Supplement to the Astronomical Almanac, Seidelman (1992)
jj_tai = [19500:1:20500]; ss_tai = 43200*ones(jj_tai); [NUT,F,NUTP,FP,NUTPP,FPP]=CL_mod_nutationAngles(jj_tai,ss_tai,"s"); //calculation of nutation angles (dPsi,dEps) [eps0,eps0p,eps0pp]=CL_mod_meanObliquity(jj_tai,ss_tai,"s"); //calculation of mean obliquity [equi,equip]=CL_mod_equinoxesEquation(NUT(1,:),eps0,F(5,:),NUTP(1,:),eps0p,FP(5,:)); //equinoxes equation plot(jj_tai'+ss_tai'/86400,equi');