Consider the sequence defined by
\begin{align}
c\_0 &{}= 1 \\
c\_n &{}= 2\,n\,c\_{n-1}-\frac{1}{2}\sum\_{m=1}^{n-1}c\_m\,c\_{n-m}.
\end{align}
How can you prove that it has the following asymptotics (strongly suggested by numerics)
$$
c\_n\sim \frac{2}{\pi}\Gamma(n)\,2^n\ \ \ ?
$$
To make the question more motivated: the associated divergent series is the asymptotic expansion of a simple ratio of modified Bessel functions
$$
R(z) = \frac{K\_0(-\frac{1}{4z})}{K\_1(-\frac{1}{4z})}\sim \sum\_{n\ge 1}c\_n z^n = 1+2z+6z^2+24z^3+\dotsb.
$$
Indeed, as pointed out in the [comment](https://mathoverflow.net/questions/418658/asymptotics-of-a-quadratic-recursion#comment1075111_418658) by Richard Stanley, one has
$$
4z^2 R'+4z R+1-R^2 = 0.
$$
So the question boils down to what can be said about the large order behaviour of the terms of the asymptotic expansion of a particular known function. If $R$ were analytic one could have used Darboux theorems to relate the answer to the type of nearest singularity. For a divergent non-Borel summable function I don't know whether there are general results. By the way, having the Borel transform in closed form could help because then one could deform the integration contour off the positive real axis to get information, in the spirit of Nevanlinna theorems.



Asymptotics of a quadratic recursion
**TL;DR:** I have a proof of your conjectured asymptotic formula, modulo the correctness of a certain alternative description of your $c\_n$ sequence.




---


I tried to complete Iosif Pinelis's elegant analysis by finding a way to derive the value of the constant $2/\pi$ (denoted $a$ in Iosif's answer) from the quadratic recurrence relation. The problem with that relation is that the behavior it implies for $c\_n$ for large $n$ seems to depend in a very sensitive way on the initial values of the sequence, so I concluded that this approach has little chance of working.


Fortunately, I've now discovered another, *linear* recurrence relation for the same sequence $c\_n$ that has better behavior and gives your claimed asymptotics without much effort. The relation is:


$$
c\_0=1, \qquad c\_n = g\_n + \sum\_{k=1}^n h\_k c\_{n-k} \quad (n\ge 1), \qquad (\*)
$$
where I define
\begin{align}
g\_n &= \frac{((2n)!)^2}{2^{3n}(n!)^3},
\\
h\_n &= \frac{((2n)!)^2}{2^{3n}(n!)^3}\cdot \frac{2n-1}{2n+1}.
\end{align}
(**Edit:** this corrects a small typo from the earlier version. As you pointed out in a comment, $g\_n$ can also be expressed as $\frac{2^n \Gamma(n+1/2)^2}{\pi \Gamma(n)}$.)


I haven't verified rigorously that this relation is equivalent to your quadratic relation, but numerically it gives the correct sequence 1, 2, 6, 24, 126, 864, 7596, ..., and I believe this should be straightforward to prove. The reasoning that led me to it involves your description of the sequence as coming from an asymptotic expansion for a ratio of two Bessel functions. I started with the relation
$$
K\_0(-1/4z) = K\_1(-1/4z) \times \sum\_{n} c\_n z^n, 
$$
and, expanding both Bessel functions in a power series (actually not quite a traditional power series because of some nasty-looking transcendental terms, but those can be factored out), massaged this into a linear system of equations satisfied by the $c\_n$'s, which gave me the recurrence after a bit of additional guesswork.


To rigorously prove the relation, one can work with this Bessel function picture and do the analysis more carefully, or one can try to prove directly that the linear recurrence is equivalent to the quadratic recurrence without any reference to Bessel functions. I suspect this is doable through an inductive argument, probably involving formulating and proving some auxiliary hypergeometric summation identities.


Finally, if we assume that $(\*)$ is correct, we can prove your claim that
$$
c\_n \sim \frac{2}{\pi} 2^n (n-1)!.
$$
Observe that
$$
\frac{c\_n}{2^n (n-1)!} = \frac{g\_n}{2^n (n-1)!} + \frac{h\_n}{2^n (n-1)!}
+ \frac{1}{2^n (n-1)!} \sum\_{k=1}^{n-1} h\_k c\_{n-k}
$$
Using Stirling's formula, you can check that each of the first two terms in this expression converges to $1/\pi$. The third term (the normalized sum) can be easily shown to be $O(1/n)$ (the first summand is $O(1/n)$, and the remaining summands are $O(1/n^2)$ and there are $n-2$ of them). Here I am using the fact that the sequence $c\_n/2^n (n-1)!$ is bounded, as shown in Iosif Pinelis's answer (and as can probably also be shown from the linear recurrence without much effort).


