The OEIS is supported by the many generous donors to the OEIS Foundation.

 Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 59th year, we have over 358,000 sequences, and we’ve crossed 10,300 citations (which often say “discovered thanks to the OEIS”). Other ways to Give
 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A248833 The curvature of touching circles inscribed in a special way in the larger segment of circle of radius 1/6 divided by a chord of length sqrt(8/75). 2
 10, 25, 160, 1225, 9610, 75625, 595360, 4687225, 36902410, 290532025, 2287353760, 18008298025, 141779030410, 1116223945225, 8788012531360, 69187876305625, 544714997913610, 4288532107003225, 33763541858112160, 265819802757894025, 2092794880205040010, 16476539238882426025 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,1 COMMENTS Refer to comment of A240926. Consider a circle C of radius 1/6 (in some length units) with a chord of length sqrt(8/75). This has been chosen such that the larger sagitta has length 1/5. The input, besides the circle C is the circle C_0 with radius R_0 = 1/10, touching the chord and circle C. The following sequence of circles C_n with radii R_n, n >= 1, is obtained from the condition that C_n touches i) the circle C, ii) the chord and iii) the circle C_(n-1). The circle curvatures C_n = 1/R_n, n >= 0, are conjectured to be a(n). If one considers the curvature of touching circles inscribed in the smaller segment (sagitta length 2/15), the sequence would be A248834 See an illustration given in the link. LINKS Vincenzo Librandi, Table of n, a(n) for n = 0..300 Kival Ngaokrajang, Illustration of initial terms Eric Weisstein's World of Mathematics,Sagitta Index entries for linear recurrences with constant coefficients, signature (9,-9,1). FORMULA From Colin Barker, Oct 15 2014: (Start) a(n) = 9*a(n-1) - 9*a(n-2) + a(n-3). G.f.: -5*(5*x^2-13*x+2) / ((x-1)*(x^2-8*x+1)). (End) a(n) = 5*(2+(4-sqrt(15))^n+(4+sqrt(15))^n)/2. - Colin Barker, Mar 03 2016 MATHEMATICA CoefficientList[Series[- 5 (5 x^2 - 13 x + 2)/((x - 1) (x^2 - 8 x + 1)), {x, 0, 30}], x] (* Vincenzo Librandi, Oct 29 2014 *) LinearRecurrence[{9, -9, 1}, {10, 25, 160}, 30] (* G. C. Greubel, Dec 20 2017 *) PROG (PARI) { r=0.6; print1(round(6/r), ", "); r1=r; for (n=1, 40, if (n<=1, ab=2-r, ab=sqrt(ac^2+r^2)); ac=sqrt(ab^2-r^2); if (n<=1, z=0, z=(Pi/2)-atan(ac/r)+asin((r1-r)/(r1+r)); r1=r); b=acos(r/ab)-z; r=r*(1-cos(b))/(1+cos(b)); print1(round(6/r), ", "); ) } (PARI) Vec(-5*(5*x^2-13*x+2)/((x-1)*(x^2-8*x+1)) + O(x^100)) \\ Colin Barker, Oct 15 2014 (Magma) I:=[10, 25, 160]; [n le 3 select I[n] else 9*Self(n-1)-9*Self(n-2)+Self(n-3): n in [1..30]]; // Vincenzo Librandi, Oct 29 2014 CROSSREFS Cf. A240926, A078986, A097315, A247512, A247335, A247512, A248834. Sequence in context: A251194 A071289 A268303 * A220039 A219377 A156183 Adjacent sequences: A248830 A248831 A248832 * A248834 A248835 A248836 KEYWORD nonn,easy AUTHOR Kival Ngaokrajang, Oct 15 2014 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified December 4 21:40 EST 2022. Contains 358570 sequences. (Running on oeis4.)