OFFSET
0,3
COMMENTS
Radius of convergence |x| < r, where r = cos(r) = 0.739085133215160... (A003957).
The continued fraction converges in the whole complex plane, cut along |z|=infinity.
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..427
S. N. Gladkovskii, Analysis Of The Continued Fractions (in Russian).
FORMULA
E.g.f.: E(x)=1/(cos(x) - x) = (1-x^2/((x-1)*G(0) + x^3))/(1-x); G(k)= 2*(2*k+1)*(k+1) - x^2 + 2*x^2*(2*k+1)*(k+1)/G(k+1); (continued fraction Euler's kind, 1-step). - Sergei N. Gladkovskii, Jan 08 2012
E.g.f.: 1/(G(0) - x) where G(k) = 1 - x^2/((4*k+1)*(4*k+2) - x^2*(4*k+1)*(4*k+2)/(x^2 - 4*(k+1)*(4*k+3)/G(k+1) )); - Sergei N. Gladkovskii, Dec 16 2012
a(n) ~ n!/((sin(r)+1)*r^(n+1)), where r = 0.73908513321516... is the root of the equation r = cos(r). - Vaclav Kotesovec, Jun 27 2013
EXAMPLE
E.g.f.: A(x) = 1 + x + 3*x^2/2! + 12*x^3/3! + 65*x^4/4! + 440*x^5/5! +...
where 1/A(x) = 1 - x - x^2/2! + x^4/4! - x^6/6! +...
Special values.
A(Pi/5) = 10/(5*(sqrt(5)+1)/2 - 2*Pi) = 5.534081362740...
A(Pi/6) = 6/(3*sqrt(3) - Pi) = 2.920333635550...
A(Pi/8) = 8/(4*sqrt(2+sqrt(2)) - Pi) = 1.882599403781...
A(Pi/10) = 10/(5*sqrt(10+sqrt(20))/2 - Pi) = 1.5701119741529...
A(Pi/12) = 12/(6*sqrt(2+sqrt(3)) - Pi) = 1.4201994774470...
MATHEMATICA
CoefficientList[Series[1/(Cos[x]-x), {x, 0, 20}], x]* Range[0, 20]! (* Vaclav Kotesovec, Jun 27 2013 *)
PROG
(PARI) {a(n)=n!*polcoeff(1/(cos(x+x*O(x^n))-x), n)}
(PARI) x='x+O('x^30); Vec(serlaplace(1/(cos(x) - x))) \\ G. C. Greubel, Jul 10 2018
(Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(1/(Cos(x) - x))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Jul 10 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Nov 15 2011
STATUS
approved