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 smaller sagitta has length 2/15. The input, besides the circle C, is the circle C_0 with radius R_0 = 1/15, touching the chord and circle C. The following sequence of circles C_n with radii R_n, n >= 1, is obtained from the conditions that C_n touches (i) the circle C, (ii) the chord and (iii) the circle C_(n-1). The curvature of the n-th circle is C_n = 1/R_n, n >= 0, and its numerator is conjectured to be a(n). The denominator is A000244 for n > 0. If one considers the curvature of touching circles inscribed in the larger segment (sagitta length 1/5), the sequence would be A248833. See an illustration given in the link.
LINKS
Kival Ngaokrajang, Illustration of initial terms
FORMULA
Conjecture: a(n) = 17*a(n-1) - 51*a(n-2) + 27*a(n-3) for n > 3. - Colin Barker, Oct 15 2014
Empirical g.f.: 5*(54*x^3-117*x^2+46*x-3) / ((3*x-1)*(9*x^2-14*x+1)). - Colin Barker, Oct 15 2014
PROG
(PARI)
{
r=0.4; print1(round(6/r), ", "); r1=r; dn=1;
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)*dn), ", ");
dn=dn*3
)
}
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Kival Ngaokrajang, Oct 15 2014
STATUS
approved