%I #15 Aug 24 2017 04:58:34
%S 1,3,33,819,34209,2189187,198717057,24278289651,3842052205761,
%T 764478393601923,186805627856569953,54994092004290217779,
%U 19197418751181422089569,7840711973025043515377667,3704137338316764145483007937,2004220869541285849551954747891
%N E.g.f.: S(x) satisfies: S(x) = Integral sqrt( (1 + S(x)^2)*(1 + 2*S(x)^2) ) dx.
%H Vaclav Kotesovec, <a href="/A289695/b289695.txt">Table of n, a(n) for n = 1..228</a>
%F E.g.f.: S(x) = Series_Reversion( Integral 1/sqrt( (1 + x^2)*(1 + 2*x^2) ) dx ).
%F E.g.f.: S(x) = sinh( Series_Reversion( Integral 1/sqrt( cosh(2*x) ) dx ) ).
%F E.g.f.: S(x) = Integral sqrt(1 + S(i*x)^2) / (1 + 2*S(i*x)^2) dx, where i^2 = -1.
%F Let C(x) be the e.g.f. of A193543, then
%F (1) S(x) = sqrt(C(x)^2 - 1),
%F (2) S(x) = Integral C(x) * sqrt(C(x)^2 + S(x)^2) dx,
%F (3) C(x) = 1 + Integral S(x) * sqrt(C(x)^2 + S(x)^2) dx,
%F (4) C(x) + S(x) = exp( Integral sqrt(1 + 2*S(x)^2) ) dx ).
%e E.g.f.: S(x) = x + 3*x^3/3! + 33*x^5/5! + 819*x^7/7! + 34209*x^9/9! + 2189187*x^11/11! + 198717057*x^13/13! + 24278289651*x^15/15! + 3842052205761*x^17/17! + 764478393601923*x^19/19! + 186805627856569953*x^21/21! +...
%e where
%e sqrt((1 + S(x)^2)*(1 + 2*S(x)^2)) = 1 + 3*x^2/2! + 33*x^4/4! + 819*x^6/6! + 34209*x^8/8! +...
%e RELATED SERIES.
%e C(x) = sqrt(1 + S(x)^2) equals the e.g.f. of A193543, and begins
%e C(x) = 1 + x^2/2! + 9*x^4/4! + 153*x^6/6! + 4977*x^8/8! + 261009*x^10/10! + 20039481*x^12/12! + 2121958377*x^14/14! + 296297348193*x^16/16! +...
%o (PARI) {a(n) = my(S=1); S = sinh( serreverse( intformal( 1/sqrt( cosh(2*x +O(x^(2*n+1))) ) ) ) ); (2*n-1)!*polcoeff(S,2*n-1)}
%o for(n=1,20,print1(a(n),", "))
%Y Cf. A193543.
%K nonn
%O 1,2
%A _Paul D. Hanna_, Aug 14 2017
|