function eps0=epsil0(YEAR,MONTH,DAY);
%Вычисление среднего наклона эклиптики к экватору
%
% Вызов:      eps0=epsil0(YEAR,MONTH,DAY);
%
%     YEAR - год;
%    MONTH - мес?ц;
%      DAY - день;
%
%       eps0 - средний наклон эклиптики

% Хомутов С.Ю.   07.10.98
% Источник: "Астрономический ежегодник СССР на 1992 год". с.637

MJD=cdtmjd(YEAR,MONTH,fix(DAY)); % модифиц. юлианска? дата
T=epoha(MJD);                    % Эпоха в дол?х столети?
                                 % от J2000
T=T+rem(DAY,1)/36525;            % учет долей суток

eps0 = 84381.448 - 46.8150*T - 0.00059*T.^2 + 0.001813*T.^3;

eps0=rem(eps0/3600,360);         % Приводим к диапазону (0-360) град.
k=find(eps0<0);
eps0(k)=eps0(k)+360;
