A230114 E.g.f. 1/(1 - sin(8*x))^(1/8).

1, 1, 9, 89, 1521, 32401, 869049, 27608489, 1019581281, 42824944801, 2017329504489, 105299243488889, 6032850630082641, 376363074361201201, 25396689469918450329, 1843101478742259481289, 143145930384321475601601, 11846611289341729822881601, 1040750126963789832859930569

Offset

0,3

Comments

Generally, for e.g.f. 1/(1-sin(p*x))^(1/p) is a(n) ~ n! * 2^(n+3/p) * p^n / (Gamma(2/p) * n^(1-2/p) * Pi^(n+2/p)). - Vaclav Kotesovec, Jan 03 2014

Links

Table of n, a(n) for n=0..18.

Formula

E.g.f. A(x) satisfies:

(1) A(x) = (cos(4*x) - sin(4*x))^(-1/4).

(2) A(x)^4/A(-x)^4 = 1/cos(8*x) + tan(8*x).

(3) A(x) = exp( Integral A(x)^4/A(-x)^4 dx ).

O.g.f.: 1/G(0) where G(k) = 1 - (8*k+1)*x - 8*(k+1)*(4*k+1)*x^2/G(k+1) [continued fraction formula from A144015 due to Sergei N. Gladkovskii].

a(n) ~ n! * 2^(4*n+3/8) / (Gamma(1/4) * n^(3/4) * Pi^(n+1/4)). - Vaclav Kotesovec, Jan 03 2014

Example

E.g.f.: A(x) = 1 + x + 9*x^2/2! + 89*x^3/3! + 1521*x^4/4! + 32401*x^5/5! + ...
where A(x)^4 = 1 + 4*x + 48*x^2/2! + 704*x^3/3! + 14592*x^4/4! + 369664*x^5/5! + ...
and 1/A(x)^4 = 1 - 4*x - 16*x^2/2! + 64*x^3/3! + 256*x^4/4! - 1024*x^5/5! + ...
which illustrates 1/A(x)^4 = cos(4*x) - sin(4*x).
O.g.f.: 1/(1-x - 8*1*1*x^2/(1-9*x - 8*2*5*x^2/(1-17*x - 8*3*9*x^2/(1-25*x - 8*4*13*x^2/(1-33*x - 8*5*17*x^2/(1-...)))))), a continued fraction.

Mathematica

CoefficientList[Series[1/(1-Sin[8*x])^(1/8), {x, 0, 20}], x]* Range[0, 20]! (* Vaclav Kotesovec, Jan 03 2014 *)

Prog

(PARI) {a(n)=local(X=x+x*O(x^n)); n!*polcoeff((cos(4*X)-sin(4*X))^(-1/4), n)}
for(n=0, 20, print1(a(n), ", "))
(PARI) {a(n)=local(A=1+x+x*O(x^n)); for(i=0, n, A=exp(intformal(A^4/subst(A^4, x, -x)))); n!*polcoeff(A, n)}
for(n=0, 20, print1(a(n), ", "))

Crossrefs

Cf. A001586 (p=2), A007788 (p=3), A144015 (p=4), A230134 (p=5), A227544 (p=6), A235128 (p=7).

Sequence in context: A094935 A258388 A296618 * A187090 A078248 A186943

Adjacent sequences: A230111 A230112 A230113 * A230115 A230116 A230117

Keyword

nonn

Author

Paul D. Hanna, Dec 20 2013

Status

approved