|
%I #26 Jan 06 2014 19:44:37
%S 1,1,5,29,265,3001,42125,696149,13296145,287706481,6959431445,
%T 186061833869,5448382252825,173418192216361,5961442393047965,
%U 220112963745653189,8687730877758518305,365023930617143804641,16266420334783460443685,766297734521812843642109
%N E.g.f. 1/(1 - sin(4*x))^(1/4).
%C Row sums of A186492 - _Peter Bala_, Feb 22 2011.
%F E.g.f. A(x) satisfies:
%F (1) A(x) = (cos(2*x) - sin(2*x))^(-1/2).
%F (2) A(x)^2/A(-x)^2 = 1/cos(4*x) + tan(4*x).
%F (3) A(x) = exp( Integral A(x)^2/A(-x)^2 dx).
%F (4) A'(x) = A(x)^3/A(-x)^2 with A(0) = 1.
%F (5) A(x) = 1/sqrt(1 - 2*Series_Reversion( Integral 1/sqrt(1+4*x-4*x^2) dx )).
%F G.f.: 1/G(0) where G(k) = 1 - x*(4*k+1) - 4*x^2*(k+1)*(2*k+1)/G(k+1); (continued fraction). - _Sergei N. Gladkovskii_, Jan 11 2013.
%F a(n) ~ 2^(3*n+5/4)*n^n/(exp(n)*Pi^(n+1/2)). - _Vaclav Kotesovec_, Jun 26 2013
%e E.g.f.: A(x) = 1 + x + 5*x^2/2! + 29*x^3/3! + 265*x^4/4! + 3001*x^5/5! +...
%e log(A(x)) = x + 4*x^2/2! + 16*x^3/3! + 128*x^4/4! + 1280*x^5/5! +...
%e A(x)^2/A(-x)^2 = 1 + 4*x + 16*x^2/2! + 128*x^3/3! +...+ 4^n*A000111(n)*x^n/n! +...
%e O.g.f.: 1/(1-x - 4*1*1*x^2/(1-5*x - 4*2*3*x^2/(1-9*x - 4*3*5*x^2/(1-13*x - 4*4*7*x^2/(1-17*x - 4*5*9*x^2/(1-...)))))) [continued fraction by Sergei Gladkovskii].
%t CoefficientList[Series[1/(1-Sin[4*x])^(1/4), {x, 0, 20}], x]* Range[0, 20]! (* _Vaclav Kotesovec_, Jun 26 2013 *)
%o (PARI) {a(n)=local(X=x+x*O(x^n)); n!*polcoeff((cos(2*X)-sin(2*X))^(-1/2), n)}
%o for(n=0, 20, print1(a(n), ", "))
%o (PARI) {a(n)=local(A=1+x+x*O(x^n));for(i=0,n,A=exp(intformal(A^2/subst(A^2,x,-x))));n!*polcoeff(A,n)}
%o for(n=0, 20, print1(a(n), ", "))
%o (PARI) /* From A'(x) = A(x)^3 / A(-x)^2: */
%o {a(n)=local(A=1); for(i=0, n, A=1+intformal(A^3/subst(A, x, -x)^2 +x*O(x^n) )); n!*polcoeff(A, n)}
%o for(n=0, 20, print1(a(n), ", "))
%o (PARI) /* 1/sqrt(1-2*Series_Reversion(Integral 1/sqrt(1+4*x-4*x^2) dx)): */
%o {a(n)=local(A=1);A=1/sqrt(1-2*serreverse(intformal(1/sqrt(1+4*x-4*x^2 +x*O(x^n)))));n!*polcoeff(A, n)}
%o for(n=0,20,print1(a(n),", "))
%Y Cf. A000111, A001586, A007788, A186492, A230134, A227544, A230114, A235329.
%K nonn
%O 0,3
%A _Paul D. Hanna_, Sep 09 2008
| 