login
a(n) = 4*n^2 + 3.
4

%I #47 Jul 11 2020 07:32:50

%S 3,7,19,39,67,103,147,199,259,327,403,487,579,679,787,903,1027,1159,

%T 1299,1447,1603,1767,1939,2119,2307,2503,2707,2919,3139,3367,3603,

%U 3847,4099,4359,4627,4903,5187,5479,5779,6087,6403,6727,7059,7399

%N a(n) = 4*n^2 + 3.

%C 2/a(n) = R(n)/r, n >= 0, with R(n) the n-th radius of the clockwise Pappus chain of the arbelos with semicircle radii r, r1 = 2r/3, r2 = r/3. See the MathWorld link for Pappus chain (there only the counterclockwise chain is shown). The counterclockwise chain companion has circle radii R(n)/r = 2/A114949(n), n>=0.

%C Binomial transform of (3, 4, 8, 0, 0, 0, 0, 0, 0, 0, ...). - _Philippe Deléham_, Mar 07 2013

%H Ivan Panchenko, <a href="/A222465/b222465.txt">Table of n, a(n) for n = 0..1000</a>

%H Kival Ngaokrajang, <a href="/A222465/a222465.pdf">Illustration of clockwise Pappus chain</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PappusChain.html">Pappus chain</a>

%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (3,-3,1).

%F a(n) = 4*n^2 + 3, n >= 0.

%F O.g.f.: (3 - 2*x + 7*x^2)/(1-x)^3.

%F a(n) = A016742(n) + 3. - _Omar E. Pol_, Mar 02 2013

%F a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n>2, a(0) = 3, a(1) = 7, a(2) = 19. - _Philippe Deléham_, Mar 05 2013

%F From _Amiram Eldar_, Jul 11 2020: (Start)

%F Sum_{n>=0} 1/a(n) = 1/6 + sqrt(3)*Pi*coth(sqrt(3)*Pi/2)/12.

%F Sum_{n>=0} (-1)^n/a(n) = 1/6 + sqrt(3)*Pi*cosech(sqrt(3)*Pi/2)/12. (End)

%e The dimensionless radii R(n)/r of the clockwise Pappus chain for the arbelos (r,r1,r2=r-r1) = r*(1,2/3,1/3) are [2/3, 2/7, 2/19, 2/39, 2/67, 2/103, 2/147, 2/199 ,...], for n >= 0. The circle for n=0 has radius r1=2/3 and center (2/3,0) with the origin at the left tip of the arbelos. The n=1 circle coincides with the one of the counterclockwise companion chain.

%p A222465(n):=n->4*n^2 + 3; seq(A222465(n), n=0..50); # _Wesley Ivan Hurt_, Feb 06 2014

%t Table[4 n^2 + 3, {n, 0, 50}] (* _Wesley Ivan Hurt_, Feb 06 2014 *)

%t Array[4 #^2 + 3 &, 44, 0] (* _Luiz Roberto Meier_, Jan 22 2015 *)

%o (PARI) a(n)=4*n^2+3 \\ _Charles R Greathouse IV_, Aug 20 2013

%Y Cf. A016742, A114949.

%K nonn,easy

%O 0,1

%A _Wolfdieter Lang_, Mar 01 2013