|
|
A219010
|
|
Numerators in a product expansion for sqrt(5).
|
|
4
|
|
|
|
OFFSET
|
0,1
|
|
COMMENTS
|
a(4) has 262 digits.
Iterating the algebraic identity sqrt(1 + 4/x) = (1 + 2*(x + 2)/(x^2 + 3*x + 1)) * sqrt(1 + 4/(x*(x^2 + 5*x + 5)^2)) produces the rapidly converging product expansion sqrt(1 + 4/x) = Product_{n = 0..inf} (1 + 2*a(n)/b(n)), where a(n) and b(n) are integer sequences when x is an integer: a(n) satisfies the recurrence a(n+1) = a(n)^5 - 5*a(n)^3 + 5*a(n) with the initial condition a(0) = x + 2 and b(n) satisfies the recurrence b(n+1) = 5/2*(b(n)^4 - b(n)^2)*sqrt(4*b(n) + 5) + b(n)^5 + 15/2*b(n)^4 - 25/2*b(n)^2 + 5 with the initial condition b(0) = x^2 + 3*x + 1.
The present case is when x = 1. The denominators b(n) are in A219011. See also A219012 (x = 2) and A219014 (x = 4).
For another product expansion for sqrt(5) see A002814.
|
|
LINKS
|
|
|
FORMULA
|
Let tau = (3 + sqrt(5))/2. Then a(n) = tau^(5^n) + 1/tau^(5^n).
a(n) = (1/2)*(1 + 5*(Fibonacci(5^n)*Fibonacci(5^n + 1) + (Fibonacci(5^n - 1))^2)).
Recurrence equation: a(n+1) = a(n)^5 - 5*a(n)^3 + 5*a(n) with initial condition a(0) = 3.
|
|
EXAMPLE
|
The first three terms of the product give 51 correct decimal places of sqrt(5): (1 + 2*3/5)*(1 + 2*123/15005)*(1 + 2*28143753123/792070839820228500005) = 2.23606 79774 99789 69640 91736 68731 27623 54406 18359 61152 5(4...).
|
|
MATHEMATICA
|
Table[(1 + 5*(Fibonacci[5^n] Fibonacci[5^n + 1] + Fibonacci[5^n - 1]^2))/2, {n, 0, 3}] (* or *)
|
|
PROG
|
(Maxima) A219010(n):=1/2*(1 + 5*(fib(5^n)*fib(5^n+1)+(fib(5^n - 1))^2))$
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|