|
|
A191301
|
|
Expansion of exp(x*arcsin(x)) = sum_{n>=0} a(n)*x^(2n)/product(k, 0<k<2n).
|
|
2
|
|
|
1, 1, 4, 49, 1303, 63513, 5044601, 598488981, 99463522845, 22073876512113, 6308016788410641, 2256148067062888845, 987271755178677563541, 518851042593331302909225, 322466959923499314299129625, 233940983258826325064978994501, 195913817323275641425583595611805
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
Expansion of exp(sqrt(x)*arcsin(sqrt(x))) = sum_{n>=0} a(n)*x^(n)/product(k, 0<k<2n). - Daniel Forgues, Mar 05 2013
|
|
LINKS
|
|
|
FORMULA
|
a(n) = (2*n-1)!*sum(m=1..n, sum(k=0..n, (-1)^k*C((2*n-m-2)/2,n-m-k)* sum(i=0..2*k, (2^i*stirling1(i+m,m)*C(2*k+m-1,i+m-1))/(i+m)!))), n>0, a(0)=1.
a(n) ~ exp(Pi/2) * 2^(2*n-1) * n^(2*n-2) / exp(2*n). - Vaclav Kotesovec, Apr 05 2016
|
|
EXAMPLE
|
exp(x*arcsin(x)) =
1 + 1/1!*x^2 + 4/3!*x^4 + 49/5!*x^6 + 1303/7!*x^8 + 63513/9!*x^10 + ...
exp(sqrt(x)*arcsin(sqrt(x))) =
1 + 1/1!*x^1 + 4/3!*x^2 + 49/5!*x^3 + 1303/7!*x^4 + 63513/9!*x^5 + ...
|
|
MATHEMATICA
|
a[n_] := (2*n-1)!* Sum[ Sum[ (Sum[ (2^i*StirlingS1[i+m, m]* Binomial[2*k+m-1, i+m-1])/(i+m)!, {i, 0, 2*k}])*(-1)^k*Binomial[(2*n-m-2)/2, n-m-k], {k, 0, n}], {m, 1, n}]; a[0] = 1; Table[a[n], {n, 0, 16}] (* Jean-François Alcover, Mar 04 2013, translated from Maxima *)
|
|
PROG
|
(Maxima)
a(n):=(2*n-1)!*sum(sum((sum((2^i*stirling1(i+m, m)*binomial(2*k+m-1, i+m-1))/ (i+m)!, i, 0, 2*k))*(-1)^k*binomial((2*n-m-2)/2, n-m-k), k, 0, n), m, 1, n);
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|