|
|
A012290
|
|
Expansion of e.g.f. arctan(sin(x)*exp(x)).
|
|
1
|
|
|
0, 1, 2, 0, -24, -140, -8, 8960, 92736, 12240, -14154208, -210017280, -50776704, 58549241920, 1133642065792, 435228385280, -500888609593344, -11981516394489600, -6690495105515008, 7684815276420464640
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
LINKS
|
|
|
FORMULA
|
a(n) = Sum_{m=0..n/2} ((4^(-m)*Sum_{j=m..n/2} (binomial(n+1,2*j+1)*((2*m+1)^(n-2*j-1)*sum(i=0..m+1/2, (2*i-2*m-1)^(2*j+1)*(-1)^(j+1-i)*binomial(2*m+1,i)))))). - Vladimir Kruchinin, Jun 30 2011
|
|
EXAMPLE
|
arctan(sin(x)*exp(x)) = x + (2/2!)*x^2 - (24/4!)*x^4 - (140/5!)*x^5 - (8/6!)*x^6 + ...
|
|
MATHEMATICA
|
CoefficientList[Series[ArcTan[Sin[x]*Exp[x]], {x, 0, 20}], x]*Range[0, 20]! (* Vaclav Kotesovec, Jan 02 2014 *)
|
|
PROG
|
(Maxima)
a(n):=sum((4^(-m)*sum(binomial(n+1, 2*j+1)*((2*m+1)^(n-2*j-1)*sum((2*i-2*m-1)^(2*j+1)*(-1)^(j+1-i)*binomial(2*m+1, i), i, 0, m+1/2)), j, m, (n)/2)), m, 0, (n)/2); /* Vladimir Kruchinin, Jun 30 2011 */
(PARI) x='x+O('x^30); concat([0], Vec(serlaplace(atan(sin(x)*exp(x))))) \\ G. C. Greubel, Oct 26 2018
(Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Arctan(Sin(x)*Exp(x)) )); [0] cat [Factorial(n)*b[n]: n in [1..m-1]]; // G. C. Greubel, Oct 26 2018
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
Patrick Demichel (patrick.demichel(AT)hp.com)
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|