function DemoCalLune(Annee)
% DEMOCALLUNE(ANNEE) affiche le calendrier des phases de la lune¨
% pour l'ANNEE.

% NL contient les dates des nouvelles lunes
nl=NouvLune(Annee,'long');

% nb de jours par mois
jpm = [31 28 31 30 31 30 31 31 30 31 30 31];
if ((rem(Annee,4) == 0 & rem(Annee,100) ~= 0) |  rem(Annee,400) == 0)
   jpm(2) = 29;	% Année bissextile
end

% matrice des phases lunaires
L=zeros(31,12);
h=zeros(12,1);

for i=1:12
	[j,m,a,h(i)]=caldat(nl(i));
	L(j:jpm(i),i)=(1:jpm(i)-j+1)';
	L(1:j-1,i)=(29-j+2:29)';
end

% représentation du calendrier lunaire
clf 
hold on
text(16,64,int2str(Annee),'color','w')

% titre = mois
moisT='Jan  Fév  Mar  Avr  Mai  Juin Juil Aoû  Sep  Oct  Nov  Déc  ';
for i=1:12
	text(2+3*(i-1),62,moisT(5*(i-1)+1:5*i),'color','w')
end

%jours
for i=1:31
	text(0,2*i-2.3,num2str(32-i),'color','w')
end

% phases de la lune
[li,co]=size(L);
for i=1:li
	for j=1:co
		if L(i,j)~=0
			DessineLune(L(i,j),31-i,j,h(j))
		end
	end
end

set(gcf,'color','b')
axis equal
axis off
hold off
