getdate.m

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

function [J,TDT]=getdate
% get julian day and TDT from input
% return: Julian and TDT

y=input('input year	(2000):  ');
if isempty(y)
   y=2000;
end
m=input('input month (1)	:  ');
if isempty(m)
   m=1;
end
d=input('input day	(1)	:  ');
if isempty(d)
   d=1;
end

%Julian day
J=cal2julian(y,m,d);
%TDT
TDT=J+deltat(J)/86400.;