login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A221706 Kochanski approximates to sqrt(2) starting with R_0=3, S_0=2. 0
2, 4, 4, 15, 17, 77, 101, 119, 143, 250, 362, 1401, 31168, 88629, 184654, 259251, 298769, 4196069, 38538873, 616984562, 1975413034, 5345718056, 27843871196, 54516286512, 334398528973, 445879679625, 495957494385, 2450869042060, 2629541150528, 4088114099884 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,1

LINKS

Table of n, a(n) for n=0..29.

Henryk Fuks, Adam Adamandy Kochanski's approximations of pi: reconstruction of the algorithm, arXiv preprint arXiv:1111.1739, 2011. Math. Intelligencer, Vol. 34 (No. 4), 2012, pp. 40-45.

FORMULA

Definitions 1 and 2 of Fuks (2011) give formulas.

PROG

(PARI)

galpha(alpha, R, S) = {floor((alpha - floor(alpha))/(R - alpha*S)); }

fuks() = { n = 29; default(realprecision, 200); alpha = sqrt(2); R = 3; S = 2; x = galpha(alpha, R, S); print1(x, ", "); for (i=1, n, R = R*(x+1) + floor(alpha); S = S*(x+1) + 1; x = galpha(alpha, R, S); print1(x, ", "); ); }

\\ Michel Marcus, Feb 07 2013

CROSSREFS

Cf. A191642.

Sequence in context: A263023 A193848 A218974 * A145891 A077815 A218075

Adjacent sequences:  A221703 A221704 A221705 * A221707 A221708 A221709

KEYWORD

nonn

AUTHOR

N. J. A. Sloane, Jan 23 2013

EXTENSIONS

More terms from Michel Marcus, Feb 07 2013

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 18 22:08 EDT 2019. Contains 322237 sequences. (Running on oeis4.)