pmoon.m

来自「球面天文学（岁差、章动）和天体力学（行星星历表）的原始计算公式、算法和程序。ma」· M 代码 · 共 16 行

function an=pmoon(TDT)
% geocentric ecliptic polar coordinates of moon at TDT (J2000)
% input: TDT
% return:geocentric ecliptic polar coordinates of moon (J2000)
header
pm=g2plan(TDT);
s=g1plan(TDT);
pm(2)=s;

ecrm=pol2rec(pm);  % rectangular ecliptic coordinates

ecrm=precess(ecrm,TDT,1); % precess to J2000

pm=rec2pol(ecrm);

an=pm;