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”).

A012297
Expansion of e.g.f. arctan(sin(x)*sin(x)), even powers only.
1
2, -8, -208, 13312, 161792, -138042368, 9993992192, 4147254525952, -1386315007787008, -171033573395529728, 284350624150760456192, -45336566464274820497408, -87136608072231186709086208
OFFSET
0,1
LINKS
EXAMPLE
E.g.f. = 2*x^2/2! - 8*x^4/4! - 208*x^6/6! + 13312*x^8/8! + ...
MATHEMATICA
With[{nn = 40}, CoefficientList[Series[ArcTan[Sin[x]^2], {x, 0, nn}], x] Range[0, nn]!][[1 ;; ;; 2]] (* G. C. Greubel, Oct 26 2018 *)
PROG
(PARI) x='x+O('x^40); v=Vec(serlaplace(atan(sin(x)^2))); vector(#v\2, n, v[2*n-1]) \\ G. C. Greubel, Oct 26 2018
(Magma) m:=40; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Arctan(Sin(x)^2) )); [Factorial(2*n+2)*b[2*n+1]: n in [0..Floor((m-4)/2)]]; // G. C. Greubel, Oct 26 2018
CROSSREFS
Sequence in context: A012535 A012300 A009801 * A000893 A030177 A009611
KEYWORD
sign
AUTHOR
Patrick Demichel (patrick.demichel(AT)hp.com)
STATUS
approved