First of all, we will prove that $a_n \to 1$.
Step 1. $(a_n)_{n \geqslant 0}$ is bounded.
Proof. We are looking for $m,M$ such that $m \leqslant a_n \leqslant M$ for every $n \geqslant 0$. If it's true for $n,n+1$ then : $$\frac{1}{M}\leqslant \frac{2}{a_n+a_{n+1}} \leqslant \frac{1}{m}$$
thus we need $mM=1$. Then we just chose $m$ wisely so that $m \leqslant a_0,a_1$.
Since the sequence is bounded, we set $\ell$ the inferior limit and $L$ the superior one.
Step 2. $\ell L=1$
Proof. Let $\phi$ be an extraction such that $a_{\phi (n)} \to \ell$ and by successive extractions we can assume that $a_{\phi (n)-1)} \to \ell _1 \in [\ell, L]$ and $a_{\phi (n)-2} \to \ell _2 \in [\ell, L]$. We obtain : $$\ell = \frac{2}{\ell _1 + \ell _2} \geqslant \frac{1}{L}$$ and then $\ell L \geqslant 1$. By replacing $\ell$ with $L$ we have $\ell L=1$
Step 3. $\ell = L=1$
Proof. It is the same idea, we set $\phi$ an extraction such that $a_{\phi (n)} \to L$, $a_{\phi (n)-1} \to \ell _1 \in [\ell , L]$, $a_{\phi (n)-2} \to \ell _2 \in [\ell , L]$, $a_{\phi (n)-3} \to \ell _3 \in [\ell , L]$. Then we have : $$\ell _1 + \ell _2 = \frac{2}{L}=2\ell $$ and $$\ell _2+ \ell _3 = \frac{2}{\ell _1}$$
Since $\ell _1, \ell_2 \geqslant \ell$ we obtain $\ell _1 = \ell _2 = \ell$ and then $\ell_2 + \ell _3= \frac{2}{\ell}=2L$ then $\ell_2=\ell_3 =L$ and we conclude $\ell =L$ : the sequence converges to $1$.
Then let $N \geqslant 0$ such that for $n \geqslant N$ we have $|a_n+a_{n+1}| \geqslant 2- \varepsilon $
For your question.
let $U_n=(a_n,a_{n+1})$ so that $U_{n+1}=f(U_n)$ where $f(x,y)= \left(y,\frac{2}{x+y} \right)$. We can see that :
$$A:= \mathrm{d} f(1,1)= \begin{pmatrix}
0 & 1 \\ -\frac{1}{2} & -\frac{1}{2} \end{pmatrix} $$
Fact. Maple or direct computation shows that $\displaystyle \rho (A) := \max \left \{ |\lambda|, \lambda \in \operatorname{sp} ( A) \right \} < 1 $
A classical lemma states that we can find a norm $\| \cdot \|_A$ such that : $$\| \mathrm{d}f (1,1) \|_A \leqslant \rho (A) + \mu < 1 $$ for a wisely chosen $\mu >0$. Then we set a ball $B:=B((1,1), \varepsilon$ such that $\|df(x)\|_A \leqslant \rho_0 <1 $ for every $x \in B$ (because $x \mapsto \mathrm{d}f(x)$ is continuous) and by the mean value theorem we obtain : \begin{align} \left \|\begin{pmatrix}
a_n \\ a_{n+1} \end{pmatrix} - \begin{pmatrix}
1 \\ 1 \end{pmatrix} \right \|_A & = \left \| f^{(n)} \left( \begin{pmatrix}
a_n \\ a_{n+1} \end{pmatrix} \right) - f^{(n)} \left( \begin{pmatrix}
1 \\ 1 \end{pmatrix}\right) \right \|_A \\
& \leqslant \left( \sup_{x \in B} \|\mathrm{d}f (x)\|_A \right)^{n-N} \left \|\begin{pmatrix}
a_N \\ a_{N+1} \end{pmatrix} - \begin{pmatrix}
1 \\ 1 \end{pmatrix} \right \|_A
\end{align}
for $n \geqslant N$ where $N$ is chosen such that $a_n \in B$ for $n \geqslant N$.By norm equivalence in finite dimensional spaces there is a constant $C_1$ such that : \begin{align}
|a_n-1| & \leqslant \max \left \{ |a_n-1|,|a_{n+1}-1| \right \} \\
& \leqslant C_1 \left \| \begin{pmatrix}
a_n \\ a_{n+1} \end{pmatrix} - \begin{pmatrix}
1 \\ 1 \end{pmatrix} \right \|_A
\end{align}
Thus we have $|a_n-1| \leqslant c\lambda ^n $ for $n \geqslant N$ with : $$c= \sup_{x \in B} \|\mathrm{d}f (x)\|_A <1$$ and $$C=C_1 c^{-N} \left \|\begin{pmatrix}
a_N \\ a_{N+1} \end{pmatrix} - \begin{pmatrix}
1 \\ 1 \end{pmatrix} \right \|_A$$
Since we can let $C$ grow we can assume that the inequality holds for every $n \geqslant 0$.
Let $x_n := a_n - 1 \to 0$, and $X_n = (x_n,x_{n+1})$. We set $e_n := \|X_n\|$
The previous stuff proves that $\displaystyle e_n = \cal{O} \left( a^{ \frac{n}{2}}\right)$ for every $a> \dfrac{1}{2}$.
A simple computation and the fact that $x_nx_{n+1} \leqslant \frac{x_n^2 + x_{n+1}^2}{2}$ shows that : $$X_{n+1} = AX_n + {\cal{O}} \left( \| X_n \|^2 \right)$$
We set : $h_n:=2^{-n/2}e_n$.
We can find a norm that we will still write $\| \cdot \|$ so that the associate matrix norm satisfies $\|A\| = \dfrac{1}{\sqrt{2}}$. A simple computation shows that : $$h_{n+1} \leqslant h_n + {\cal O} \left( \frac{h_n^2}{2^{n/2}} \right)$$ (I've used that $\|A X_n\| \leqslant \|A\| \|X_n\| \leqslant \dfrac{\|X_n\|}{\sqrt{2}}$ )
Thus $ \displaystyle h_{n+1} - h_n \leqslant {\cal O} \left( (1- \varepsilon)^n \right) $ for a wisely chosen $a$. Then $\sum h_{n+1} - h_n$ converges so that $h_n$ converge and we have $e_n \sim C 2^{-n/2}$.
I will investigate if we can go further, but it seems that we can, just setting $Y_n:=A^{n}X_n$ and study this.
I hope it is clear, and please feel free to correct my so weak english.
- @math110 can you prove that $(1-a_n)(1-a_{n+1})<0$ for $n$ large enough ?
- Can someone gives the precise asymptotic behavior of $a_n-1$ ?
– M.LTA Oct 13 '15 at 16:22