OFFSET
1,1
COMMENTS
a(n) gives twice the inverse radius of the circles touching the large Arbelos (2/3,1/3) circle (radius 1) and the n-th and (n-1)-th circles of the counterclockwise Pappus chain.
For twice the curvatures (inverse radii) of the counterclockwise Pappus chain of the (2/3,1/3) arbelos see A114949, also for the MathWorld link to Pappus chain.
For the small curvatures touching the left circle of the (2/3,1/3) arbelos and the n-th and (n-1)-st circles of the counterclockwise Pappus chain see A259555.
The curvatures of the circles can be computed with Descartes' three (actually 5) circle theorem. See A259555 for links to Descartes' theorem.
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..1000
Kival Ngaokrajang, Illustration of initial terms.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).
FORMULA
a(n) = 4*n^2 - 4*n + 19, n >= 1.
O.g.f.: x*(19 - 30*x + 19*x^2)/(1 - x)^3.
From Elmo R. Oliveira, Nov 17 2024: (Start)
E.g.f.: exp(x)*(4*x^2 + 19) - 19.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 3. (End)
MATHEMATICA
Array[4 #^2 - 4 # + 19 &, 40] (* Michael De Vlieger, Jul 02 2015 *)
LinearRecurrence[{3, -3, 1}, {19, 27, 43}, 40] (* Harvey P. Dale, May 17 2016 *)
PROG
(PARI) vector(60, n, 4*n^2 - 4*n + 19) \\ Michel Marcus, Jul 03 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang and Kival Ngaokrajang, Jul 02 2015
EXTENSIONS
Terms a(37) and beyond from Andrew Howroyd, May 01 2020
STATUS
approved