function S0=times0(YEAR,MONTH,DAY);
%Вычисление среднего гринвичского звездного времени в 0 час. UT
%
% Вызов:      S0=times0(YEAR,MONTH,DAY) или S0=times0(MJD);
%
%     YEAR - год;
%    MONTH - мес?ц;
%      DAY - день;
%      MJD - модифицир. юлианска? дата.
%
%       S0 - среднее гринвичское звездное врем?
%            в 0 час. UT (в часах).

% Хомутов С.Ю.   27.08.98
% Источник: "Астрономический ежегодник СССР на 1992 год". с.615
%     06.06.2002 - добавлена возможность ввода MJD

if nargin==1,     % вводитс? только MJD
   MJD=fix(YEAR);
else				   % вводитс? календарна? дата
   MJD=cdtmjd(YEAR,MONTH,fix(DAY));
end
T=epoha(MJD);     % Эпоха в дол?х столети? от J2000

S0 = 6*3600 +41*60 +50.54841 +...
     8640184.812866*T + 0.093104*T.^2 - 6.2e-6*T.^3; % в секундах

S0=mod(S0/3600,24);  % Приводим к диапазону (0-24) час.
