OFFSET
1,1
COMMENTS
The previous definition was "a(n) = normalized inverse radius of the inscribed circle that is tangent to the left circle of the symmetric arbelos and the n-th and (n-1)-st circles in the Pappus chain".
See links section for image of these circles, via Wolfram MathWorld (there an asymmetric arbelos is shown).
The Rothman-Fukagawa article has another picture of the circles, based on a Japanese 1788 sangaku problem. - N. J. A. Sloane, Jan 02 2020
LINKS
N. J. A. Sloane, Table of n, a(n) for n = 1..1000
Brady Haran and Simon Pampena, Epic Circles, Numberphile video (2014).
Tony Rothman and Hidetoshi Fukagawa, Japanese temple geometry, Scientific American, Vol. 278, No. 5, May 1998, pp. 85-91.
Eric Weisstein's World of Mathematics, Image of inscribed circles (in red).
Eric Weisstein's World of Mathematics, Pappus Chain.
Wikipedia, Pappus chain.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).
FORMULA
a(n) = 4*n^2 - 4*n + 15.
From Colin Barker, May 14 2014: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
G.f.: -x*(15*x^2 - 22*x + 15)/(x-1)^3. (End)
From Descartes three circle theorem:
a(n) = 2 + c(n) + c(n-1) + 2*sqrt(2*(c(n) + c(n-1) + c(n)*c(n-1)), with c(n) = A059100(n) = n^2 + 2, n >= 1, which produces 4*n^2 - 4*n + 15. - Wolfdieter Lang, Jul 01 2015
From Elmo R. Oliveira, Nov 17 2024: (Start)
E.g.f.: exp(x)*(4*x^2 + 15) - 15.
a(n) = A060747(n)^2 + 14. (End)
EXAMPLE
For n = 1, the radius of the outermost circle divided by the radius of a circle drawn tangent to all three of the initial inner circle, the opposite inner circle (the 0th circle in the chain), and the 1st circle in the chain is 15.
For n = 2, the radius of the outermost circle divided by the radius of a circle drawn tangent to all three of the initial inner circle, the 1st circle in the chain, and the 2nd circle in the chain is 23.
MAPLE
MATHEMATICA
Table[4 n^2 - 4 n + 15, {n, 50}] (* Wesley Ivan Hurt, May 13 2014 *)
LinearRecurrence[{3, -3, 1}, {15, 23, 39}, 50] (* Harvey P. Dale, Feb 22 2023 *)
PROG
(Magma) [4*n^2 - 4*n + 15: n in [1..50]]; // Wesley Ivan Hurt, May 13 2014
(PARI) a(n) = 4*n^2 - 4*n + 15 \\ Charles R Greathouse IV, May 14 2014
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Aaron David Fairbanks, May 13 2014
EXTENSIONS
More terms from Wesley Ivan Hurt, May 13 2014
More terms and links from Robert G. Wilson v, May 13 2014
Edited: Name reformulated (with consent of the author). - Wolfdieter Lang, Jul 01 2015
Edited by N. J. A. Sloane, Jan 02 2020, simplifying the definition and adding a reference to the fact that this sequence arose in a sangaku problem from 1788 in a temple in Tokyo Prefecture.
STATUS
approved