function [a_alpha,a_delta]=ap_pos(m_alpha,m_delta,mu_alpha,mu_delta,T0,T)
%Âû÷èñëåíèå ïð?ìîãî âîñõîæäåíè? è ñêëîíåíè? çâåçäû ñ ó÷åòîì ïðåöåññèè
%è ñîáñòâåííîãî äâèæåíè?
%
%Âûçîâ:
%   [a_alpha,a_delta]=a_posit(m_alpha,m_delta,mu_alpha,mu_delta,T0,T);
%
%   m_alpha, m_delta  - ñðåäíèå ìåñòà íà ýïîõó T0;
%   mu_alpha,mu_delta - âåêîâîå ñîáñòâåííîå äâèæåíèå;
%                  T0 - ýïîõà êàòàëîãà (â äîë?õ ñòîëåòè? îò J2000.0);
%                   T - ýïîõà íàáëþäåíèé (â äîë?õ ñòîëåòè? îò J2000.0);
%
%     a_alpha,a_delta - èñïðàâëåííûå çà ïðåöåññèþ è ñîáñòâåííîå äâèæåíèå
%                       ìåñòà
%
% Âåëè÷èíû äë? ïð?ìîãî âîñõîæäåíè? (alpha) ä.á. â ÷àñàõ,
% äë? ñêëîíåíè? (delta) - â ãðàäóñàõ
%
% Èñòî÷íèê: Æ. Ìååñ "Àñòðîíîìè÷åñêèå ôîðìóëû äë? êàëüêóë?òîðîâ" - Ì: Ìèð,
%           1988. ñòð.54-55

T = 1+T;         %ïåðåõîä ê ýïîõàì îò J1900.0
T0 = 1+T0;

m = 3.07234 + 0.00186*T0; % â ñåê. âðåìåíè
n = 20.0468 - 0.0085*T0;  % â óãë. ñåê.

% Ãîäè÷íîå èçìåíåíèå êîîðäèíàò èç-çà ïðåöåññèè

da = m + (n/15)*sin(d2r(m_alpha*15)).*tan(d2r(m_delta));
dd = n*cos(d2r(m_alpha*15));

% Ãîäè÷íîå èçìåíåíèå êîîðäèíàò ñ ó÷åòîì ñîáñòâåííîãî äâèæåíè?

da = da+mu_alpha/100;
dd = dd+mu_delta/100;

% Èçìåíåíèå êîîðäèíàò çà âðåì? ìåæäó ýïîõàìè T0 è T

da = da*(T-T0)*100;
dd = dd*(T-T0)*100;

% Èñïðàâëåííûå ìåñòà çâåçä

a_alpha = m_alpha+da/3600;
a_delta = m_delta+dd/3600;
