function mjd_axis(mjd,dmjd);
% Ðàçìåòêà âðåìåííîé îñè ïî NJD ñ øàãîì dmjd ñóòîê
%
% Âûçîâ: mjd_axis(mjd,dmjd)
%
%  mjd - àðãóìåíò â ôóíêöèè plot(mjd,Y) â ñóòêàõ
% dmjd - øàã ïî âðåìåíè (â ñóòêàõ)
%

% 19.09.2003 - Khomutov

mjd1=fix(min(mjd));
mjd2=fix(max(mjd))+1;
mjd_tick=(mjd1:dmjd:mjd2)';

[year,mon,day]=mjdcdt(mjd_tick);

if year(1)<1900,
   error('Error in MJD (year<1900)')
end

mjd_label=[];
for k=1:length(mjd_tick),
   mjd_label=strvcat(mjd_label,sprintf('%02d/%02d',mon(k),day(k)));
end
set(gca,'Xlim',[mjd1 mjd2],'Xtick',mjd_tick,'Xticklabel',mjd_label)
