Mein MATLAB Forum

function JD=date2jul (y,m,d,ut)

% Funktion zur Berechnung des Julianischen Datums
% bei gegebenen Datums im gregorianischen Kalender
% und Uhrzeit (Weltzeit UT)

% Eingabeparameter:
%   y = Jahr, m = Monat, d = Tag, ut = Uhrzeit
% Ausgabeparameter:
%   JD = julianisches Datum

if m >= 3
    f=y;
    g=m;
else
    f=y-1;
    g=m+12;
end

A=2-fix(f/100)+fix(f/400);
JD=fix(365.25*f)+fix(30.6001*(g+1))+d+A+1720994.5+ut/24;