login
A243001
a(n) = date of Easter in March of n, in the Gregorian calendar after A.D. 1582. a(n) = 0 if Easter occurs in April.
1
0, 0, 0, 0, 29, 0, 0, 0, 0, 29, 0, 0, 26, 0, 0, 22, 0, 0, 0, 0, 30, 0, 0, 26, 0, 0, 0, 0, 0, 0, 0, 30, 0, 0, 26, 0, 31, 0, 0, 27, 0, 0, 30, 0, 0, 0, 0, 31, 0, 0, 27, 0, 0, 23, 0, 0, 0, 0, 31, 0, 0, 27, 0, 0, 0, 0, 0, 0, 0, 31, 0, 0, 28, 0, 0, 0, 0, 28, 0, 0, 25, 0
OFFSET
1583,5
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
Cf. A243002.
Sequence in context: A366269 A123269 A067158 * A092995 A055785 A040847
KEYWORD
nonn,easy
AUTHOR
STATUS
approved