A141626 E.g.f. satisfies: A(x) = exp(x*A(sin(x^2))).

1, 1, 1, 7, 25, 121, 841, 9031, 74257, 779185, 9465841, 98157511, 1372854121, 19401758377, 301086307705, 5668190911111, 101297827616161, 1892733268852321, 38807453736318817, 751339911589120135

Offset

0,4

Links

Alois P. Heinz, Table of n, a(n) for n = 0..100

Maple

A:= proc(n) option remember; if n=0 then 1 else unapply (convert (series ( exp (x*A(n-1)(sin(x^2))), x, n+1), polynom), x) fi end; a:= n-> coeff (A(n)(x), x, n)*n!: seq (a(n), n=0..21);

Mathematica

A[n_] := A[n] = If[n == 0, 1, Normal[Series[Exp[x*A[n-1][Sin[x^2]]], {x, 0, n+1}]] /. x -> #]&; a[n_] := Coefficient[A[n][x], x, n]*n!; Table[a[n], {n, 0, 21}] (* Jean-François Alcover, Feb 14 2014, after Maple *)

Crossrefs

Sequence in context: A199893 A129791 A216688 * A138292 A138738 A057573

Adjacent sequences: A141623 A141624 A141625 * A141627 A141628 A141629

Keyword

nonn

Author

Alois P. Heinz, Aug 24 2008

Status

approved