login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A290879 E.g.f. L(x) of aerated sequence satisfies: L(x) = Integral 1 / sqrt( cosh(2*L(x)) ) dx. 5

%I #16 Apr 02 2019 11:48:56

%S 1,-2,44,-2840,367760,-79719200,26016555200,-11921650083200,

%T 7300922254496000,-5759532173685056000,5688335131502291840000,

%U -6875441991877610827520000,9983390897443366347676160000,-17148197258942716354314368000000,34392433372153876998446324480000000,-79646680171456811546338888517120000000,210930739548407111241584046599398400000000

%N E.g.f. L(x) of aerated sequence satisfies: L(x) = Integral 1 / sqrt( cosh(2*L(x)) ) dx.

%H Robert Israel, <a href="/A290879/b290879.txt">Table of n, a(n) for n = 1..200</a>

%F E.g.f.: L(x) = Series_Reversion( Integral sqrt( cosh(2*x) ) dx ).

%F Let C(x) and S(x) be the e.g.f.s of A290880 and A290881, respectively, then e.g.f. L(x) satisfies:

%F (1) L(x) = log(C(x) + S(x)),

%F (2) cosh(2*L(x)) = C(x)^2 + S(x)^2,

%F (3) cosh(L(x)) = C(x) and sinh(L(x)) = S(x),

%F where C(x)^2 - S(x)^2 = 1 and C'(x)^2 + S'(x)^2 = 1.

%e E.g.f.: L(x) = x - 2*x^3/3! + 44*x^5/5! - 2840*x^7/7! + 367760*x^9/9! - 79719200*x^11/11! + 26016555200*x^13/13! - 11921650083200*x^15/15! + 7300922254496000*x^17/17! - 5759532173685056000*x^19/19! + 5688335131502291840000*x^21/21! +...

%e where

%e 1/sqrt( cosh(2*L(x)) ) = 1 - 2*x^2/2! + 44*x^4/4! - 2840*x^6/6! + 367760*x^8/8! - 79719200*x^10/10! + 26016555200*x^12/12! - 11921650083200*x^14/14! +...

%p Q:= series(Int(sqrt(cosh(2*t)),t),t,100):

%p S:= series(RootOf(Q-y,t),y,100):

%p seq(coeff(S,y,j)*j!,j=1..100,2); # _Robert Israel_, Aug 13 2017

%t terms = 17; m = 2 terms; cc = CoefficientList[InverseSeries[Integrate[Sqrt[ Cosh[2 x]] + O[x]^m, x], x], x]; DeleteCases[cc * Range[0, m-1]!, 0] (* _Jean-François Alcover_, Apr 02 2019 *)

%o (PARI) {a(n) = my(L=x); for(i=1,n, L = intformal( 1/sqrt(cosh(2*L + O(x^(2*n+2)))) )); (2*n-1)!*polcoeff(L,2*n-1)}

%o for(n=1,20, print1(a(n),", "))

%o (PARI) {a(n) = my(L=x); L = serreverse( intformal( sqrt(cosh(2*x + O(x^(2*n+2)))) )); (2*n-1)!*polcoeff(L,2*n-1)}

%o for(n=1,20, print1(a(n),", "))

%Y Cf. A290880, A290881, A290882, A290883, A153302.

%K sign

%O 1,2

%A _Paul D. Hanna_, Aug 13 2017

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 04:13 EDT 2024. Contains 371235 sequences. (Running on oeis4.)