function DessineLune(ind,jour,mois,heure)
n=30;
pas=100;
t=0:pi/pas:2*pi;
xc=cos(t);
yc=sin(t);
t1=pi/2:pi/pas:3*pi/2;
x1=cos(t1);
y1=sin(t1);
t2=5*pi/2:-pi/pas:3*pi/2;
x2=cos(t2);
y2=sin(t2);
dx1=1/7;
dx2=2/15;

if ind<=15
   i=1-(ind-1)*dx1;
else
   i=1-(ind-15)*dx2;
end

h=(i^2-1)/(2*i+eps);
r=sqrt(h^2+1);
alfa=atan2(1,-h);
if h<=0
   t=-alfa:alfa/pas:alfa;
else
   t=2*pi-alfa:-alfa/pas:alfa;
end

x=r*cos(t)+h;
y=r*sin(t);

if ind==1
  	text(3*mois-0.7,2*jour,sprintf('%2i',round(heure)),'color','w')
elseif ind<15
   fill([x2 x]+3*mois,[y2 y]+2*jour,'y','edgecolor','y')
elseif ind==15
   fill(xc+3*mois,yc+2*jour,'y','edgecolor','y')
else
   fill([x x1]+3*mois,[y y1]+2*jour,'y','edgecolor','y')
end
plot(xc+3*mois,yc+2*jour,'color',[.5 .5 .5])


