( Avertissement: Les divisions portent sur des entiers dont le reste n'est pas pris en compte )
G = Année mod 19
Calendrier Julien
I = (19 * G + 15) mod 30
J = (année + (année / 4) + I) mod 7
Calendrier Grégorien
C = année / 100
H = (C - C/4 - (8 * C + 13)/25 + (19 * G) + 15) mod 30
I = H - (H/28)*(1 - (H/28)*(29/(H + 1))*((21 - G)/11))
J = (année + année/4 + I + 2 - C + C/4) mod 7
Ensuite pour les deux calendriers
L = I - J
MoisPaques = 3 + (L + 40)/44
JourPaques = L + 28 - 31*(MoisPaques / 4)
Cet algorithme est basé sur l'algorithme de Oudin (1940) cité dans "supplément explicatif de l'almanach astronomique" aux éditions P. Kenneth Seidelmann.
Pour ceux qui veulent creuser un peu plus
G est le nombre d'or - 1
H est 23 - Epacte (modulo 30)
I est le nombre de jours à partir du 21 mars pour atteindre la pleine lune Pascale
J est le jour de la semaine pour la pleine lune Pascale (0=dimanche)
L est le nombre de jours du 21 mars jusqu'au dimanche de ou avant la pleine lune Pascale (un nombre entre -6 et 28)