|
|
A141627
|
|
E.g.f. satisfies: A(x) = exp(x*A(sin(x)*x)).
|
|
1
|
|
|
1, 1, 1, 7, 25, 101, 721, 8653, 66753, 574489, 7509601, 76972281, 938682625, 14560953109, 226747016497, 4226839991221, 73582415930881, 1133352458491185, 23260378180229569, 478118972286336241, 9648062520586144001, 239696533730926824301, 6090975028325867045521
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,4
|
|
COMMENTS
|
First negative term is a(94).
|
|
LINKS
|
Alois P. Heinz, Table of n, a(n) for n = 0..200
|
|
MAPLE
|
A:= proc(n) option remember; if n=0 then 1 else unapply(convert(
series(exp(x*A(n-1)(sin(x)*x)), x, n+1), polynom), x) fi
end:
a:= n-> coeff(A(n)(x), x, n)*n!:
seq(a(n), n=0..30);
|
|
MATHEMATICA
|
A[n_] := A[n] = If[n == 0, 1, Normal[Series[Exp[x*A[n-1][Sin[x]*x]], {x, 0, n+1}]] /. x -> #]&; a[n_] := Coefficient[A[n][x], x, n]*n!; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Feb 14 2014, after Maple *)
|
|
CROSSREFS
|
Sequence in context: A035509 A332942 A247173 * A289606 A102900 A155271
Adjacent sequences: A141624 A141625 A141626 * A141628 A141629 A141630
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
Alois P. Heinz, Aug 24 2008
|
|
STATUS
|
approved
|
|
|
|