Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #12 Aug 05 2021 12:57:22
%S 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,
%T 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,
%U 27,0,0,0,0,0,0,0,31,0,0,28,0,0,0,0,28,0,0,25,0
%N 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.
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Easter">Easter</a>
%H <a href="/index/Ca#calendar">Index entries for sequences related to calendars</a>
%t 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}]]]
%t If[#[[1]]==4,0,#[[2]]]&/@easter[1583,1583+120] (* _Harvey P. Dale_, Aug 05 2021 *)
%o (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, ", ")));
%Y Cf. A243002.
%K nonn,easy
%O 1583,5
%A _Arkadiusz Wesolowski_, Aug 16 2014