pearth.m

球面天文学（岁差、章动）和天体力学（行星星历表）的原始计算公式、算法和程序。matlab编译通过。主要包括：公历/儒略历转换为儒略日

function an=pearth(TDT)
% heliocentric ecliptic polar coordinates of earth at TDT (J2000)
% input: TDT
% return:heliocentric ecliptic polar coordinates of earth (J2000)

pemb=g3plan(TDT);
pm=pmoon(TDT);

rect_emb=pol2rec(pemb);
rect_m=pol2rec(pm);
rect_earth=rect_emb-(1/82.3007)*rect_m;

an=rec2pol(rect_earth);