OFFSET
0,5
FORMULA
E.g.f. A(x) satisfies: A'(x) = A(x)^2 - exp(-x).
From Vaclav Kotesovec, Jun 11 2020: (Start)
E.g.f.: exp(-x/2) * (BesselI(1, 2*exp(-x/2)) * (BesselK(0, 2) + BesselK(1, 2)) + (BesselI(0, 2) - BesselI(1, 2)) * BesselK(1, 2*exp(-x/2))) / ((BesselI(1, 2) - BesselI(0, 2)) * BesselK(0, 2*exp(-x/2)) + BesselI(0, 2*exp(-x/2)) * (BesselK(0, 2) + BesselK(1, 2))).
a(n) ~ n! / r^(n+1), where r = 1.4982609322383959128764444062824740935658895762... is the real root of the equation (BesselI(0, 2) - BesselI(1, 2)) * BesselK(0, 2*exp(-r/2)) = (BesselK(0, 2) + BesselK(1, 2)) * BesselI(0, 2*exp(-r/2)). (End)
MATHEMATICA
a[n_] := a[n] = (-1)^n + Sum[Binomial[n - 1, k] a[k] a[n - k - 1], {k, 0, n - 1}]; Table[a[n], {n, 0, 24}]
terms = 24; A[_] = 1; Do[A[x_] = Normal[Integrate[A[x]^2 - Exp[-x], x] + O[x]^(terms + 1)], terms]; CoefficientList[A[x], x] Range[0, terms]!
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jun 11 2020
STATUS
approved