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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A142970 Numerators of n-th approximants of a continued fraction for Pi-3. 5
0, 1, 6, 61, 660, 8901, 133266, 2303865, 43808040, 928665225, 21386693790, 537861526965, 14540730176700, 423407835413325, 13140639311294250, 434929825450371825, 15237733330856005200, 565064979900590948625, 22056613209702152061750, 905913636742121921038125 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

COMMENTS

The corresponding denominators are A001879(n), n >= 0.

Pi = 3 + 1^2/(6 + 3^2/(6 + 5^2/(6 + ... ))). See the J.-P. Delahaye reference. R. Rosenthal mentioned this continued fraction in an e-mail to the author Jul 16 2008.

For the approximants in lowest terms cf. the ones for 3*(Pi-3) given by A130411(n)/A130412(n) in lowest terms.

The above continued fraction for Pi is the particular case n = 0, x = 3 of a result of Ramanujan, previously given by Euler - see Berndt et al., Chapter 12, Entry 25. - Peter Bala, Feb 19 2015

REFERENCES

J.-P. Delahaye, Le fascinant nombre pi, Pour la Science, Paris 1997. In German: Pi - die Story, Birkhäuser, 1999 Basel, p. 87.

LINKS

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

B. C. Berndt, R. L. Lamphere, and B. M. Wilson Chapter 12 of Ramanujan's second notebook: Continued fractions, Rocky Mountain Journal of Mathematics, Volume 15, Number 2 (1985), 235-310

L. Euler, De fractionibus continuis observationes, The Euler Archive, Index Number 123, Section 67.

W. Lang, Approximants for Pi-3 and more.

L. J. Lange, An Elegant Continued Fraction for π, The American Mathematical Monthly, 106 (1999), 456-458.

FORMULA

a(n) = 6*a(n-1) + ((2*n-1)^2)*a(n-2), a(0)=0, a(1)=1.

E.g.f.: (-3*(1+x-sqrt(1-4*x^2))+ 2*(1+x)*arcsin(2*x))/(1-2*x)^(5/2) from the solution of the linear second order differential equation (1-4*x^2)*y''(x) - 2*(8*x+3)*y'(x) - 9*y(x)=0, obtained from the recurrence, with inputs y(0)=0 and y'(0)=1. A special solution is the e.g.f. of the denominators A001879: (1+x)/(1-2*x)^(5/2).

a(n) ~ (Pi-3) * 2^(n+3/2) * n^(n+2) / exp(n). - Vaclav Kotesovec, Oct 05 2013

EXAMPLE

Approximants a(n)/A001879(n) (not in lowest terms): [0/1]; [1/6]; [6/45]; [61/420]; [660/4725]; [8901/62370];..

Approximants in lowest terms: [0/1]; [1/6]; [2/15]; [61/420]; [44/315]; [989/6930]; ...

MATHEMATICA

RecurrenceTable[{a[0]==0, a[1]==1, a[n]==6 a[n-1] + (2 n-1)^2 a[n-2]}, a, {n, 30}] (* Vincenzo Librandi, Feb 20 2015 *)

PROG

(MAGMA) I:=[1, 6]; [0] cat [n le 2 select I[n] else 6*Self(n-1)+(2*n-1)^2*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Feb 20 2015

CROSSREFS

Cf. A254795.

Sequence in context: A053462 A160751 A236670 * A034659 A064088 A191803

Adjacent sequences:  A142967 A142968 A142969 * A142971 A142972 A142973

KEYWORD

nonn,easy,frac,cofr

AUTHOR

Wolfdieter Lang, Sep 15 2008

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 August 20 10:17 EDT 2019. Contains 326149 sequences. (Running on oeis4.)