hel2geo.m

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

  function an=hel2geo(p,e,eps)
  % reduce heliocentric to geocentric coordinates
  % input: ecliptic polar
  % return : equatorial polar
  
  ecre=pol2rec(e);
  ecrp=pol2rec(p);
  v=ecrp-ecre;
     
  eqr=ec2eq(v,eps,1);
  
  an=rec2pol(eqr);