0,3

Vaclav Kotesovec, Table of n, a(n) for n = 0..470

E.g.f. A(x) satisfies: A''(x) = exp(x) + A(x) * A'(x).

From Vaclav Kotesovec, Jun 11 2020: (Start)

E.g.f.: (BesselY(0, sqrt(2))*(BesselJ(1, sqrt(2)*exp(x/2)) - sqrt(2)*exp(x/2)*BesselJ(0, sqrt(2)*exp(x/2))) + BesselJ(0, sqrt(2))*(sqrt(2)*exp(x/2)*BesselY(0, sqrt(2)*exp(x/2)) - BesselY(1, sqrt(2)*exp(x/2)))) / (BesselJ(1, sqrt(2)*exp(x/2))*BesselY(0, sqrt(2)) - BesselJ(0, sqrt(2))*BesselY(1, sqrt(2)*exp(x/2))).

a(n) ~ 2 * n! / r^(n+1), where r = 1.35169030867903432729790416904526340210784862703704233748118252928787... is the smallest real root of the equation BesselY(0, sqrt(2))*BesselJ(1, sqrt(2)*exp(r/2)) = BesselJ(0, sqrt(2))*BesselY(1, sqrt(2)*exp(r/2)). (End)

a[n_] := a[n] = 1 + Sum[Binomial[n - 2, k - 1] a[k] a[n - k - 1], {k, 1, n - 1}]; Table[a[n], {n, 0, 23}]

terms = 23; A[_] = 0; Do[A[x_] = Normal[Integrate[Integrate[Exp[x] + A[x] D[A[x], x], x], x] + O[x]^(terms + 1)], terms]; CoefficientList[A[x], x] Range[0, terms]!

Cf. A000111, A006677, A052886, A054687, A080635, A332776.

Sequence in context: A242790 A013044 A110577 * A012948 A013182 A013170

Adjacent sequences: A335438 A335439 A335440 * A335442 A335443 A335444

Ilya Gutkovskiy, Jun 10 2020

approved