How to determine the general term of the $(u_n)$ sequence with $u_0 \in \mathbb R_+$ and $u_{n+1}=\frac{u_n^2+1}{1+2u_n}$ ?

Source : les dattes à Dattier

  • $\begingroup$ With the trick you can do the same think with : $u_{n+1}=u_{n}^2-2$ $\endgroup$
    – Dattier
    Jan 26, 2017 at 10:51

3 Answers 3


Part 1: An interesting phenomenon occurs for $u_0=2$:

$$u_1=1, u_2=\dfrac{2}{3}, u_3=\dfrac{13}{21}, u_4=\dfrac{610}{987}, u_5=\dfrac{1346269}{2178309}...$$

$$\text{i.e.,} \ \ u_1=\dfrac{F_1}{F_2}, u_2=\dfrac{F_3}{F_4}, u_3=\dfrac{F_7}{F_8}, u_4=\dfrac{F_{15}}{F_{16}}, u_5=\dfrac{F_{31}}{F_{32}}...$$

where the numerators and denominators are consecutive elements of Fibonacci sequence ($F_n$). See(https://oeis.org/A000045) with general form:


It is not so surprising because we know the connection between the Golden Ratio and Fibonacci sequence, with striking jumps (classical of the quick convergence of Newton's method ; see below). Relationship:


is a consequence of identity $F_{2n-1}=F_{n}^2+F_{n-1}^2$ see this.

Part 2 (new Edit): Let us have a look at the general case now. We are going to transform the nice formula obtained by Francesco Alem., taking into account the fact that it has a narrow domain of validity (we must have $-1<\frac{2 u_0+1}{\sqrt{5}}<1$, which means $u_0<\dfrac{\sqrt{5}-1}{2}=\Phi-1.$)

Let us write this formula under the form:

$$\dfrac{2u_n+1}{\sqrt{5}}=\coth \left(2^n \coth ^{-1}\left(\frac{2 u_0+1}{\sqrt{5}}\right)\right)$$


$$\tag{*}v_n=\coth \left(2^n \coth ^{-1}(v_0)\right)$$

by setting $$v_n:=\frac{2u_n+1}{\sqrt{5}}.$$

Let us transform $(*)$ by using the following formulas:

$$\coth(x)=\dfrac{e^{2x}+1}{e^{2x}-1} \ \ \text{and} \ \ \coth^{-1}(x)=\frac12\ln\left(\dfrac{1+x}{1-x}\right)$$

(see (https://en.wikipedia.org/wiki/Inverse_hyperbolic_function)). Setting


relationship (*) becomes

$$v_n=\left(exp(2 \ln(a^{2^{n-1}})+1)\right)/\left(exp(2 \ln(a^{2^{n-1}})-1)\right)$$


giving a very tractable explicit formula for $u_n$.

An important consequence of (**) is that we can in this way assert the convergence of sequence $v_n$, thus of sequence $u_n.$

The fact that $|a|>1$ implies that $v_n\to1$. As a consequence, $u_n\to \frac{\sqrt{5}-1}{2}=\Phi-1.$

Remark : We have to look for another recurrence formula, complementary to (*) for the values of $v_0$ that are larger than $\Phi-1.$

Part 3: An interpretation of sequence $u_{n+1}=\frac{u_n^2+1}{1+2u_n}$.

It could come from Newton's method applied to the roots of $f(x)=x^2+x-1$:


that may converge towards one of the roots of $f$, i.e., $-\Phi$ or $\Phi-1$ (where $\Phi$ denotes the golden ratio: $\Phi\approx 1.618...$, i.e.) under the condition of being in one of the basins of attraction.

As $u_0>0$, convergence is necessarily towards the positive root $\Phi-1.$

"may converge" $\rightarrow$ means that a further analysis is necessary. The condition of convergence/divergence of Newton's method are not that evident...

  • $\begingroup$ Sorry, misread the question. I have some doubts about the existence of a tractable expression for the general term... Do you have a track ? $\endgroup$
    – Jean Marie
    Jan 26, 2017 at 8:48
  • $\begingroup$ Most of the problem I propose are based on simple tricks (easily understood) but hard to find, If you do not already know them. $\endgroup$
    – Dattier
    Jan 26, 2017 at 8:55
  • $\begingroup$ See my edit. Am I heading in the "right" direction ? $\endgroup$
    – Jean Marie
    Jan 26, 2017 at 9:22
  • $\begingroup$ In the general case, I think is difficult to have regularity. $\endgroup$
    – Dattier
    Jan 26, 2017 at 9:30

1/The trick :

Let $t,u,v$ functions with $t\circ v=v\circ u$, then : $t^n\circ v=v\circ u^n$.

Here, we have choosen, $t,u,v$ of the type : $$t(x)=\frac{x^2+1}{2x+1}, u(x)=x^2 \text{ and } v(x)=\frac{ax+b}{cx+d}$$

2/For the solution proposed by Francesco, the form is :

$$h_n=f\circ g^n\circ f^{-1}$$

with $$f(x)=\frac{\sqrt 5}{2}\coth(x)-\frac{1}{2}$$ $$g(x)=2x$$

We have : $h_n\circ h_m=h_{n+m}$

and a conjecture: $$h_1(x)=\frac{1+x^2}{2x+1}$$


we go to show : $h_1(x)=\frac{1+x^2}{2x+1}$

we have : $$(1):\coth(2x)=\frac{\coth(x)^2+1}{2\coth(x)}$$ More : $$f=w \circ \coth$$ with $$w(x)=\frac{\sqrt{5}}{2}x-\frac{1}{2}$$ and $$w^{-1}(x)=\frac{(2x+1)}{\sqrt{5}}$$ so

$$h_1=w \circ \coth \circ g \circ \coth^{-1} \circ w^{-1} $$

$$\text{ with } w(x)=\frac{\sqrt{5}}{2}x-\frac{1}{2} \text{ and } g(x)=2x$$

with using of $(1)$ then : $$ \coth \circ g \circ \coth^{-1}(x)=\frac{x^2+1}{2x} $$ so $$h_1(x)=w\left(\frac{w^{-1}(x)^2+1}{2w^{-1}(x)}\right)=\frac{\sqrt{5}}{2}\left(\frac{(\frac{2x+1}{\sqrt{5}})^2+1}{2(\frac{2x+1}{\sqrt{5}})}\right)-\frac{1}{2}$$







about the case $a_{n+1}=a_n^2-2$

This case is interested because, for the previous case, resolved by Francesco and Jean-Marie, we could make changes of invertible variables, here not.

We use, a function $v$, no invertible, so $u(x)=x^2$ and $t(x)=x^2-2$ are not conjugate, but the trick continues to work and here


we have :

$$\text{ if } a_{n+1}=a_n^2-2, \text{ and }b \in \mathbb C \text{ with } a_0=b+\frac{1}{b} \text{ then } a_n=b^{2^{n}}+\frac{1}{b^{2^n}}$$

  • $\begingroup$ Waooo... I would have never thought. Thank you for this nice problem. Besides, have you seen the answer I have made to Francesco Alem ? I am almost certain that there is some kind of orthogonal polynomials family behind that. $\endgroup$
    – Jean Marie
    Jan 26, 2017 at 10:15
  • 1
    $\begingroup$ The same door can have several keys, and the same key can open several doors. $\endgroup$
    – Dattier
    Jan 26, 2017 at 10:18
  • $\begingroup$ $f$ and $h$ are conjugated through $g$. $\endgroup$
    – Jean Marie
    Jan 26, 2017 at 10:35
  • $\begingroup$ Show to $h_1(x)=\frac{1+x^2}{2x+1}$ resolves the problem. $\endgroup$
    – Dattier
    Jan 26, 2017 at 10:39
  • $\begingroup$ almost 6 months later, how are you ? $\endgroup$
    – Jean Marie
    Jun 18, 2017 at 18:52

Mathematica has solved the problem. $$ u_n=\frac{1}{2} \left(\sqrt{5} \coth \left(2^n \coth ^{-1}\left(\frac{2 u_0+1}{\sqrt{5}}\right)\right)-1\right) $$ There is nothing intuitive about this solution imho.

Mathematica command:


  • 2
    $\begingroup$ How do you write the Mathematica request giving this ? $\endgroup$
    – Jean Marie
    Jan 26, 2017 at 9:23
  • $\begingroup$ added just now. $\endgroup$ Jan 26, 2017 at 9:28
  • $\begingroup$ There is a problem with solution, There should be terms at the power $2^n$ $\endgroup$
    – Dattier
    Jan 26, 2017 at 9:42
  • $\begingroup$ No problem with the solution. check again. $\endgroup$ Jan 26, 2017 at 9:46
  • $\begingroup$ I do not have this software, I trust you $\endgroup$
    – Dattier
    Jan 26, 2017 at 9:49

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .