Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A104019
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A104019 Years in the Gregorian calendar for which Easter falls on the 25th day of the month. +0
1
1663, 1666, 1674, 1731, 1734, 1742, 1883, 1886, 1894, 1943, 1951, 2035, 2038, 2046, 2103, 2187, 2190, 2198, 2255, 2258, 2266, 2323, 2326, 2334, 2407, 2410, 2418, 2491, 2559, 2570, 2573, 2581, 2627, 2630, 2638, 2779, 2782, 2790, 2874, 2877, 2885, 2931 (list; graph; listen)
OFFSET

1,1

COMMENT

The starting point for the sequence is explained by the fact that the Gregorian calendar was only introduced in 1582.

The complete Easter cycle lasts 5700000 years. In this cycle, Mar 25 occurs 110200 times and Apr 25 occurs 42000 times for a total of 152200 times. This reduces to 761 occurrences every 28500 years (~2.67%). - Hans Havermann (pxp(AT)rogers.com), Jan 27 2008

LINKS

Author?, Calculation of Easter.

M. Montes, Frequency of the Date of Easter over one complete Gregorian Easter Cycle.

FORMULA

The formula is based on the algorithm of Oudin (1940) taken from the link.

MATHEMATICA

(* first do *) Needs["Miscellaneous`Calendar`"] (* then *) Select[ Range[1582, 2941], EasterSunday[ # ][[3]] == 25 &] (from Robert G. Wilson v (rgwv(AT)rgwv.com), Apr 06 2005)

PROGRAM

(PARI) edate(yr1, yr2, day) = { local(flag=1, d, y, y2, ct, dt); for(d=day, day, ct=0; for(y=yr1, yr2, dt=oudin(y); if(eval(mid(dt, 4, 2))==d, if(flag, y2=y; flag=0); ct++; \ print(ct" "dt" "y-y2); print1(y", "); if(y2<>y, y2=y); ); ); \ print1(ct", "); ) } oudin(y) = \This is based on the algorithm of Oudin (1940) { local(c, n, k, i1, i2, i3, a1, a2, m, d, l, dt, dat=""); c=floor(y/100); n=y-19*floor(y/19); k=floor((c-17)/25); i1=c-floor(c/4)-floor((c-k)/3)+19*n+15; i2=i1-30*floor(i1/30); i3=i2-floor(i2/28)*(1-floor(i2/28)*floor(29/(i2+1))*floor((21-n)/11)); a1=y+floor(y/4)+i3+2-c+floor(c/4); a2=a1-7*floor(a1/7); l=i3-a2; m=3+floor((l+40)/44); d=l+28-31*floor(m/4); dat = concat(dat, right(Str(m+100), 2)); dat = concat(dat, " "); dat = concat(dat, right(Str(d+100), 2)); dat = concat(dat, " "); dat = concat(dat, Str(y)); return(dat); }

CROSSREFS

Cf. A104034.

Sequence in context: A154075 A054809 A163273 this_sequence A054810 A068744 A054811

Adjacent sequences: A104016 A104017 A104018 this_sequence A104020 A104021 A104022

KEYWORD

nonn

AUTHOR

Cino Hilliard (hillcino368(AT)gmail.com), Mar 31 2005

EXTENSIONS

More terms from Robert G. Wilson v (rgwv(AT)rgwv.com), Apr 06 2005.

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified November 25 14:49 EST 2009. Contains 167514 sequences.


AT&T Labs Research