角度轉弧度：

function hd=dms2rad(JD);%角度化弧度函式
%JD=99.3527;
  w=abs(JD);
  d=fix(w);
  df=(w-d)*100;
  f=fix(df);
  m=(df-f)*100;
  ddd=d+f/60+m/3600;
  hd=ddd*pi/180;
  if JD<0;
      hd=-hd;
  end
  return

弧度轉角度：

function JD=rad2dms(hd)
%弧度化角度函式
   w=abs(hd);
   df=w*180/pi;
   d=fix(df);
   mf=(df-d)*60;
   m=fix(mf);
   sf=(mf-m)*60;
   JD=d+m/100+sf/10000;
   if hd<0;
       JD=-JD;
   end
return