OFFSET
1583,5
LINKS
MATHEMATICA
Module[{computusGauss=Function[{year}, Module[{a, b, c, d, e, f, g, h, i, j, k, month, day}, a=Mod[year, 19]; b=Quotient[year, 100]; c=Mod[year, 100]; d=Quotient[b, 4]; e=Mod[ b, 4]; f=Quotient[8 b+13, 25]; g=Mod[19 a+b-d-f+15, 30]; h=Quotient[c, 4]; i=Mod[ c, 4]; j=Quotient[a+11 g, 319]; k=Mod[2 e+2 h-i-g+j+32, 7]; month=Quotient[g-j+k+90, 25]; day=Mod[g-j+k+month+19, 32]; {month, day}]]}, Table[computusGauss[k], {k, start, end}]]]
If[#[[1]]==4, 0, #[[2]]]&/@easter[1583, 1583+120] (* Harvey P. Dale, Aug 05 2021 *)
PROG
(PARI) for(y=1583, 1664, a=lift(Mod(y, 19)); b=floor(y/100); c=lift(Mod(y, 100)); d=floor(b/4); e=lift(Mod(b, 4)); f=floor((b+8)/25); g=floor((b-f+1)/3); h=lift(Mod(19*a+b-d-g+15, 30)); i=floor(c/4); k=lift(Mod(c, 4)); l=lift(Mod(32+2*e+2*i-h-k, 7)); m=floor((a+11*h+22*l)/451); p=lift(Mod(h+l-7*m+114, 31)); if(floor((h+l-7*m+114)/31)==3, print1(p+1, ", "), print1(0, ", ")));
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Arkadiusz Wesolowski, Aug 16 2014
STATUS
approved