login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A012303
Expansion of e.g.f. sech(sin(x)*sin(x)), even powers only.
1
1, 0, -12, 240, 4368, -942720, 42639168, 8638506240, -2287049555712, 38746394603520, 145600412857955328, -42974022434701332480, -7760313189299657306112, 11863589839842518338928640
OFFSET
0,3
LINKS
EXAMPLE
E.g.f. = 1 - 12*x^4/4! + 240*x^6/6! + 4368*x^8/8! + ...
MAPLE
seq(coeff(series(factorial(n)*sech(sin(x)*sin(x)), x, n+1), x, n), n = 0 .. 26, 2); # Muniru A Asiru, Oct 26 2018
MATHEMATICA
With[{nn = 100}, CoefficientList[Series[Sech[Sin[x]^2], {x, 0, nn}], x] Range[0, nn]!][[1 ;; ;; 2]] (* G. C. Greubel, Oct 25 2018 *)
PROG
(PARI) x='x+O('x^50); v=Vec(serlaplace(1/cosh(sin(x)*sin(x)))); vector(#v\2, n, v[2*n-1]) \\ G. C. Greubel, Oct 25 2018
(Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( 1/Cosh(Sin(x)^2) )); [Factorial(2*n)*b[2*n+1]: n in [0..Floor((m-2)/2)]]; // G. C. Greubel, Oct 25 2018
CROSSREFS
Sequence in context: A012544 A009052 A213449 * A119837 A012538 A006588
KEYWORD
sign
AUTHOR
Patrick Demichel (patrick.demichel(AT)hp.com)
STATUS
approved