Va a ser más sencillo de lo que sigue para el uso de dos pasos a la vez, así que podemos empezar por calcular la relación de $a_{n+2} = \frac{a_n(n+2) + n}{a_n + n}$, y dejando $b_n = a_n / \sqrt n$ obtenemos $b_{n+2} = \frac {b_n (1 + 2/n) + 1/\sqrt n}{b_n\sqrt{1/n+2/n^2} + \sqrt{1+2/n}}$.
Vamos a denotar $\begin{bmatrix} a & b \\ c & d \end{bmatrix} x = \frac {ax+b}{cx+d}$.
Si tenemos una relación $b_{n+2} = \begin{bmatrix} 1+A & B \\ C & 1+D \end{bmatrix} b_n$ donde cada una de las $A,B,C,D$ $O(1/\sqrt n)$ y nos fijamos entonces en algún $c_n = (b_n-k)\sqrt n$, obtenemos la secuencia $c_n$ el mismo tipo de relación (así, con dominante término que todavía se comienzan $I_2$) con
$$A' = (1+kC)\sqrt{1+2/n}-1, \\B' = (B-kD+kA-k^2C)\sqrt{n+2}, \\
C' = C/\sqrt n, \\ D' = D+kC$$.
Ahora mismo, $C$ es todavía tan grande como todos los demás, para hacer de $B'$ no de un orden de magnitud más grande que el resto, tenemos que resolver un grado $2$ ecuación en $k$, lo que da $k= \pm 1$.
Si definimos $c_n = (b_n-1)\sqrt n$, obtenemos $A' = -n^{-1/2} + O(n^{-1}), D' = +n^{-1/2} + O(n^{-1}), C' = O(n^{-1}), B' = O(n^{-1/2})$.
A partir de ahora $C'$ seguirá cayendo, la dominante en términos de $A'$ $D'$ no será capaz de cambiar y de hacer la próxima $B$ a una $O(n^{-1/2})$ plazo, se tendrá que recoger $k = $ el término constante de $B/(A-D)$.
Lo que esto nos dice es que hay exactamente dos asintótica de los desarrollos a $a_n$ cuyos términos de error en cada nivel obedecer recurrentes en la relación de la forma $r_{n+2} = [I_2 + O(n^{-1/2})] r_n$.
Ya que todo esto se puede hacer de manera algebraica, los dos de alimentación de la serie se $\sqrt n$ los tiempos de las dos soluciones en $\Bbb Q[[n^{-1/2}]]$ a la ecuación original en $a_n/\sqrt n$ (y de ir de uno a otro por el cambio en el signo de $\sqrt n$).
Ahora podemos ver el $u_n = ((((a_n/\sqrt n - 1)\sqrt n - \frac 12)\sqrt n + \frac 18)\sqrt n + \frac 18)\sqrt n$.
Satisface $u_{n+2} = \begin{bmatrix}1 - n^{-1/2} + O(n^{-1}) && \frac 7 {64}n^{-1} + O(n^{-3/2}) \\ O(n^{-5/2}) && 1 + n^{-1/2} + O(n^{-1}) \end{bmatrix} u_n$,
(y cada una constante en los $O$ puede ser, efectivamente, calculada)
Esto significa que si usted tiene $u_n$ lo suficientemente cerca como para $0$,$u_{n+2} = (1-2n^{-1/2})u_n + O(n^{-1})$. Por lo que debe significar que para cada uno de los obligados $\delta$ hay un $n_0$ tal que $\forall n > n_0, |u_n < \delta | \implies |u_{n+2} < \delta|$.
Si efectivamente calcular $n_0$, por ejemplo, $\delta = 1$, esto permite que usted seguramente comprobar que una secuencia $u_n$ se queda delimitada por la búsqueda de un plazo $u_n$ $n$ bastante grande que está por debajo de la $\delta$. Si la asintótica de desarrollo es válida, de un número finito de cálculo será capaz de demostrarlo.
Me encantaría ver una prueba de que el desarrollo asintótico es válido para cada valor de inicio de $a_n$ (a excepción del valor que se va a quedar fuera infinitamente a partir de la fórmula asintótica y en lugar de obedecer a otros asintótica de desarrollo)