OFFSET
0,1
COMMENTS
This number is the inverse of the number whose decimal expansion is A136319. It is also the latter number minus phi. - Geoffrey Caveney, Apr 18 2014
From A.H.M. Smeets, Apr 01 2026: (Start)
Repeat s = s + phi; s=1/s. The initial value of s is irrelevant, as long as s != -phi.
Similar constants defined by continued fractions 1/(r+1/(r+1/(r + ...))): A010527 (real part for r=i), A086773 (r=Pi), A086774 (r=exp(1)), A101263 (r=sqrt(2)), A188485-1 (r=1/2), A340616 (r=2*sqrt(2)), A394765 (r=(sqrt(7+2*sqrt(2))-sqrt(2)-1)/2).
In general, for the metallic number M = [N; N, N, N, ...], the number defined by the continued fraction [0; M, M, M, ...] is the smallest positive real root of x^4 + N*x^3 - 3*x^2 - N*x + 1 = 0. (End)
LINKS
FORMULA
Repeat s = s + phi; s=1/s. The initial value of s is irrelevant, as long as it isn't -phi.
Equals (sqrt(5+phi)-phi)/2, where phi = A001622. - R. J. Mathar, Sep 15 2012
Equals 1/A136319. - Geoffrey Caveney, Apr 18 2014
Smallest positive real root of x^4 + x^3 - 3*x^2 - x + 1 = 0. - A.H.M. Smeets, Mar 29 2026
EXAMPLE
1
------
phi + 1
------
phi + 1
-------
phi + etc
Equals 0.4772599964740196445422298845...
MATHEMATICA
First[RealDigits[(Sqrt[5 + GoldenRatio] - GoldenRatio)/2, 10, 100]] (* Paolo Xausa, Apr 07 2026 *)
PROG
(PARI) ?\p 2000 ?f(n) = phi=(sqrt(5)+1)/2; s=0; for(x=1, n, s=s+phi; s=1/s); print(s)
(PARI) polrootsreal(Pol([1, 1, -3, -1, 1]))[3] \\ A.H.M. Smeets, Mar 29 2026
(PARI) solve(x = 0.01, 1, log(x + (1 + sqrt(5))/2) + log(x)) \\ A.H.M. Smeets, Apr 01 2026
CROSSREFS
KEYWORD
AUTHOR
Cino Hilliard, Aug 03 2003
EXTENSIONS
Name clarified by Sean A. Irvine, Mar 31 2026
STATUS
approved
