rec2pol.m

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

function an=rec2pol(rec)
% transfer from rectangular to polar  coordinates
r=sqrt(dot(rec,rec));
r0=sqrt(rec(1)*rec(1)+rec(2)*rec(2));
PHI=atan2(rec(3),r0);
TH=atan2(rec(2),rec(1));
if TH<0
   TH=TH+2*pi;
end


an=[TH,PHI,r];