% J.-B. ROUX, 3S, 1990
%
% Ce programme détermine le jour de Pâques 
% le plus tôt et le plus tard dans l'intervalle
% 1900 à 2099 inclus.
%
dpaques=[];
minimum=31;
maximum=0;
for Y=1900:2099
   N=Y-1900;
   A=rem(N,19);
   B=floor((7*A+1)/19);
   M=rem(11*A+4-B,29);
   Q=floor(N/4);
   W=rem(N+Q+31-M,7);
   dateP=25-M-W;
   dpaques=[dpaques,dateP];
   if dateP < minimum
      minimum = dateP;
      datemin = Y;
   end
   if dateP > maximum
      maximum = dateP;
      datemax = Y;
   end
end
disp(['Pâques est au plus tôt le ',int2str(31+minimum),...
     ' mars ',int2str(datemin)])
disp(['Pâques est au plus tard le ',int2str(maximum),...
     ' avril ',int2str(datemax)])

