OFFSET
1,1
COMMENTS
See A335828 for details.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..238
MATHEMATICA
c[1] = 1/2; c[n_] := c[n] = (2*Sum[(n - k)*(2*n - 2*k - 1)*c[n - k]*c[k], {k, 1, n - 1}] - Sum[(n - m)*(2*n - 2*m - 1)*c[n - m]*c[m - k]*c[k], {m, 2, n - 1}, {k, 1, m - 1}])/(n*(2*n - 1)); Denominator @ Array[c, 16]
(* or *)
Quiet[-Denominator @ CoefficientList[AsymptoticDSolveValue[{y[x]*y'[x]^2 == 2*(1-y[x]), y[0] == 1}, y[x], {x, 0, 25}], x][[3;; -1;; 2]]] (* requires Mathematica 11.3+ *)
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Amiram Eldar, Jun 25 2020
STATUS
approved