2
$\begingroup$

Let $a, b , c, d$ be reals. Define the sequence $(x_n)$ as:

$$x_0 = a,\,\, x_1 = b$$ $$x_n = \left(1 - \frac{b^2}{n^2}\right)x_{n-1} + \frac{1}{n-1}\sum_{k=0}^{n-2}\binom{2n+1}{2k+1}^{-1} (x_{k+1}-x_k)(c\, x_{n-k-1}- d\, x_{n-k-2}),\,\,\, n \geq2.$$

I'm actually studying the asymptotic series of a solution of some nonlinear differential equation, where $(-1)^n(2n+1)!x_n$ represent the coefficients of such series. I want to prove that $(x_n)$ is convergent. Here are two examples for different values of $(a, b , c, d).$

enter image description here enter image description here

It seems (after several numerical tests) that the sequence is bounded and monotone from specific $n_0.$ The boundness of the sequence imply that $$\sum_{k=0}^{n-2}\binom{2n+1}{2k+1}^{-1} (x_{k+1}-x_k)(c\, x_{n-k-1}- d\, x_{n-k-2})$$ is bounded and the term with the sum goes to zero.

Thank you for any hint

$\endgroup$
4
  • 2
    $\begingroup$ Can you say a few words about what motivated you to investigate this sequence? Are there any restrictions on $a,b,c,d$? $\endgroup$ Sep 10, 2021 at 10:23
  • $\begingroup$ @LiviuNicolaescu I edited the question. There are no restrictions for $a, b , c ,d.$ $\endgroup$
    – A. PI
    Sep 10, 2021 at 10:34
  • 1
    $\begingroup$ Could you include explicitly the nonlinear differential equation that lead to this recurrence? $\endgroup$ Sep 10, 2021 at 11:05
  • $\begingroup$ The setting will be then very complicated. I want to prove the convergence of such sequence independently. $\endgroup$
    – A. PI
    Sep 10, 2021 at 11:19

2 Answers 2

3
$\begingroup$

That is, indeed, more bark than bite, but the bark is somewhat louder than mlk presented it.

First of all, let us define $X_k=\max(|x_k|,|x_{k-1}|)$ and notice that ${2n+1\choose 2k+1}\ge 2^{\min(k,n-k)}n$ for $1\le k\le n-1$, provided that $n$ is not too small. With all this, we have the recursive inequality \begin{align} X_n &\le X_{n-1}+\frac{C}{n^2}\sum_{k=1}^{n-1}2^{-\min(k,n-k)}X_kX_{n-k} \\ &\le X_{n-1}+\frac{2C}{n^2}\sum_{k=1}^{n-1}2^{-k}X_kX_{n-k} \label{1}\tag{$\ast$} \end{align} for large enough $n$. Our task will be to show that any such sequence $X_n$ is bounded. We'll do it in several steps. We shall also consider the majorant $Y_n$ of $X_n$ that is the same as $X_n$ until the recurrent inequality \eqref{1} kicks in and $2C/n$ becomes less than $\frac 12$ and beyond this point is defined by the recurrence $$ Y_n=Y_{n-1}+\frac 1{2n}\sum_{k=1}^{n-1} 2^{-k}Y_kY_{n-k} \label{2}\tag{$\ast\ast$} $$

We'll start with showing that $Y_n\le A^n$ for all $n\ge 1$ with some $A<+\infty$

Indeed, just take any $A>2$ for which the inequality holds for $Y_n$ before \eqref{2} kicks in and then use the trivial induction on $n$.

This makes the function $$ f(z)=\sum_{n=0}^\infty Y_nz^n $$ well-defined and analytic in some neighborhood of $0$. Note now that it satisfies the differential equation $$ z\frac d{dz}[(1-z)f(z)]=\frac 12[f(z/2)-Y_0]f(z)+P(z)\label{3}\tag{${\ast\ast\ast}$} $$ where $P(z)$ is some polynomial. This equation allows one to analytically extend $f$ from any open disk to a twice larger one (treating $f(z/2)-Y_0$ as an already known function and thus viewing \eqref{3} as a linear differential equation for $f(z)$ and solving it by the standard formula) until you hit the singularity at $1$. For that last step (extending to the whole unit disk), we have the linear differential equation of the sort $$ f'(z)=\frac{Q(z)}{1-z}f(z)+\frac{R(z)}{1-z} $$ where $Q$ and $R$ are analytic up to the boundary ($f(z/2)$ is incorporated in $Q$), so the standard solution formula shows that $f$ can grow at most like $(1-|z|)^p$ for some large $p$ near the boundary, which, by Cauchy formulae, implies that $Y_n=O(n^p)$ and so is $X_n$.

Now it is time to use the full strength of \eqref{1}. Without loss of generality, we may assume that our initial $p$ is half-integer. Then, if $p\ge 3/2$, we can say that $$ X_n-X_{n-1}\le \frac{C'}{n^2}\sum_{k=1}^{n-1}2^{-k}k^p(n-k)^p=O(n^{p-2})\,. $$ which implies that $X_n=O(n^{p-1})$. Once we reach $O(n^{1/2})$ in this way, the next conclusion is that $X_n-X_{n-1}=O(n^{-3/2})$ and after that we conclude that $X_n$ is actually bounded.

But once it is bounded, the bark loudness reduces to the level described by mlk except there is one term in the sum where the coefficient is actually of order $n^{-2}$ rather than $n^{-3}$ (look at $k=0$). However, since it is just one, it is harmless.

The end.

$\endgroup$
0
$\begingroup$

This looks like one of these examples that is more bark than bite. If we separate the first term of the sum and do a bit of reordering we get $$x_n-x_{n-1} = \frac{b^2}{n^2} x_{n-1} + \frac{b-a}{(n-1)(2n+1)} (c x_{n-1}-d x_{n-2}) + \frac{1}{n-1}\sum_{k=1}^{n-2}\binom{2n+1}{2k+1}^{-1} (x_{k+1}-x_k)(c\, x_{n-k-1}- d\, x_{n-k-2})$$ Now the sum has less than $n-1$ terms and the binomial coefficient will always be at least of order $n^3$. So a rather brutal estimate would be (for some constants $c_1,c_2>0$) $$||x_n|-|x_{n-1}|| \leq |x_n-x_{n-1}| < \frac{c_1}{n^2}(|x_{n-1}|+|x_{n-2}|) + \frac{c_2}{n^3} (\max_{k<n} |x_k|)^2. \quad (*) $$

A simple comparison then yields that $|x_n| \leq a_n$ where $a_0=|a|,a_1=|b|$ and $$ a_n -a_{n-1} = \frac{c_1}{n^2} (a_{n-1}+a_{n-2}) + \frac{c_2}{n^3} a_{n-1}^2, $$ which looks like it should be bounded by some discrete Gronwall-argument. But if $a_n$ is bounded, then so is $|x_n|$ and then $(*)$ implies that $x_n$ is a Cauchy-sequence.

$\endgroup$
4
  • $\begingroup$ Thank you for this answer, but how will you use the "discrete-Gronwall" argument ! The last term in the r.h.s has a square. This was the mistake that killed the answer here math.stackexchange.com/questions/4221679/… f $\endgroup$
    – A. PI
    Sep 10, 2021 at 13:53
  • 3
    $\begingroup$ Alas, I do not think this will work. Indeed, the difference inequality $a_n -a_{n-1}\le\frac{c_1}{n^2} (a_{n-1}+a_{n-2}) + \frac{c_2}{n^3} a_{n-1}^2$ will hold if e.g. $a_n=(n+1)^2$, $c_1>0$, and $c_2=3$. $\endgroup$ Sep 10, 2021 at 14:21
  • $\begingroup$ @IosifPinelis True, I underestimated the square there. So only the "bounded -> convergent" part still holds. One can probably extract more terms to get a higher exponent $n^a$ in front of that term, but it might not be enough. Now I wonder if there might not be an unbounded sequence for large intial values. $\endgroup$
    – mlk
    Sep 10, 2021 at 18:18
  • 2
    $\begingroup$ @mlk : The problem is that $x_n$ may initially grow very fast. E.g., for $a = 1, b = 10, c = 100, d = 0;$ we have $(x_0,\dots,x_{10})\approx(1, 10, 1560, 106655, 4.23279*10^6, 1.10927*10^8, 2.07656*10^9, 2.93052*10^{10}, 3.24375*10^{11}, 2.90352*10^{12}, 2.15372*10^{13})$. However, I think at some point a subtle transition to a steady regime occurs, and I don't know how to capture this transition. $\endgroup$ Sep 10, 2021 at 18:33

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.