LMST=lmst(mjdat,lambda)

function LMST=lmst(mjdat,lambda)

% function LMST=lmst(mjdat,lambda)
%
% local mean sidereal time
% algorithms by Montenbruck/Pfleger: "Astronomy with the Personal Computer"
%
% 05.05.04 M.Penzkofer

mjdat0=fix(mjdat);
UT = (mjdat-mjdat0)*24;
T  = (mjdat0-51544.5)/36525.0;
GMST = 6.697374558 + 1.0027379093*UT+(8640184.812866+(0.093104-6.2E-6*T)*T)*T/3600.0;
LMST = 24.0*FRAC( (GMST-lambda/15.0) / 24.0 );