Алгоритм вычисления даты католической Пасхи

Данный алгоритм принадлежит неаполитанскому астроному Алоизию Лилию и немецкому монаху-иезуиту Кристоферу Клавию (конец XVI века).
Вычисление даты католической пасхи в год Y:

1. G=(Y mod 19) + 1 (G — так называемое «золотое число в метоновом» цикле — 19-летнем цикле полнолуний)
2. C=[Y/100] + 1 (если Y не кратен 100, то С — номер века)
3. X=[3C/4] – 12 (поправка на то, что три из четырех лет кратных 100 не високосный)
4. Z=[(8C+5)/25] – 5 (синхронизация с лунной орбитой, год не кратен лунному месяцу)
5. D=[5Y/4] – X – 10 (в марте день −D mod 7 будет воскресенье)
6. E=(10G+20+Z-X) MOD30 (эпакта — указывает на день наступления полнолуния)
7. ЕСЛИ (E = 24) ИЛИ (E = 25 И G > 11) ТО увеличить E на 1
8. N = 44 - E (N-е марта — день календарного полнолуния)
9. ЕСЛИ N < 21 ТО увеличить N на 30
10. N=N+7-((D+N) MOD7
11.ЕСЛИ N > 31 ТО дата Пасхи (N − 31) апреля ИНАЧЕ дата Пасхи N марта

Алгоритм вычисления даты католической Пасхи