function st=time_s(MJD,UT,lon);
%Вычисление местного звездного времени
%
%Вызов:
%        st=time_s(MJD,UT,lon);
%
% MJD - модифицированна? юлианска? дата наблюдени?;
%  UT - Всемирное врем? (в часах) момента наблюдени?;
%          UT = LT-7 (летом)
%          UT = LT-6 (зимой)
% lon - долгота (в часах; положительна к востоку !!!);
%
%  st - местное звездное врем? (в часах).

% Хомутов С.Ю.    28.08.98
% Источник:  "Астрономический ежегодник на 1992 год", с.619
% 06.06.2002 - замена имени ф-ции на time_s из-за того, что times - встроенна?

myu1 = 1.0027379093;   % Посто?нна? перевода солнечного времени
                       % в звездное

[year,mon,day] = mjdcdt(MJD);
S0 = times0(year,mon,day);   % Гринвичское звездное врем? в 0:00 UT

st = S0 + UT*myu1 + lon;
st=mod(st,24);
