4.4.2: Raíces de la Ecuación Indicial

EN ESTA SECCIÓN CONSIDERAREMOS los tipos de soluciones que se pueden obtener de la ecuación diferencial,

\[y^{\prime \prime}(x)+a(x) y^{\prime}(x)+b(x) y(x)=0 \nonumber \]

cuando\(x=0\) es un punto singular regular. En este caso, asumimos que\(x a(x)\) y\(x^{2} b(x)\) tenemos las expansiones convergentes de la serie Maclaurin

\[ \begin{aligned} x a(x) &=a_{0}+a_{1} x+a_{2} x^{2}+\ldots \\ x^{2} b(x) &=b_{0}+b_{1} x+b_{2} x^{2}+\ldots \end{aligned} \label{4.40} \]

Usando el Método Frobenius, asumimos\(y(x)\) y sus derivados son de la forma

\[\begin{aligned} y(x) &=\sum_{n=0}^{\infty} c_{n} x^{n+r} \\ y^{\prime}(x) &=\sum_{n=0}^{\infty} c_{n}(n+r) x^{n+r-1} \\ y^{\prime \prime}(x) &=\sum_{n=0}^{\infty} c_{n}(n+r)(n+r-1) x^{n+r-2} \end{aligned} \label{4.41} \]

Insertando estas series en la ecuación diferencial, obtenemos

\[\sum_{n=0}^{\infty} c_{n}\left[(n+r)(n+r-1)+(n+r) x a(x)+x^{2} b(x)\right] x^{n+r-2}=0 \nonumber \]

Usando las expansiones para\(x a(x)\) y\(x^{2} b(x)\), observamos que la potencia más baja de\(x\) es\(n+r-2\) cuándo\(n=0\). El coeficiente para el\(n=0\) término debe desaparecer:

\[c_{0}\left[r(r-1)+a_{0} r+b_{0}\right]=0\nonumber \]

Suponiendo que\(c_{0} \neq 0\), tenemos la ecuación indicial

\[r(r-1)+a_{0} r+b_{0}=0\nonumber \]

Las raíces de la ecuación indicial determinan el tipo de comportamiento de la solución. Esto equivale a considerar tres casos diferentes. Que las raíces de la ecuación sean\(r_{1}\) y\(r_{2}\). Entonces,

i. Raíces distintas con\(r_{1}-r_{2} \neq\) número entero.

En este caso tenemos dos soluciones linealmente independientes,

\[\begin{aligned} &y_{1}(x)=|x|^{r_{1}} \sum_{n=0}^{\infty} c_{n} x^{n}, \quad c_{0}=1 \\ &y_{2}(x)=|x|^{r_{2}} \sum_{n=0}^{\infty} d_{n} x^{n}, \quad d_{0}=1 \end{aligned} \nonumber \]

ii. Raíces iguales:\(r_{1}=r_{2}=r\).

La forma para\(y_{1}(x)\) es la misma, pero se necesita usar el Método de Reducción del Orden para buscar la segunda solución linealmente independiente.

\[\begin{aligned} &y_{1}(x)=|x|^{r} \sum_{n=0}^{\infty} c_{n} x^{n}, \quad c_{0}=1 \\ &y_{2}(x)=|x|^{r} \sum_{n=0}^{\infty} d_{n} x^{n}+y_{1}(x) \ln |x|, \quad d_{0}=1 \end{aligned} \nonumber \]

iii. Raíces distintas con entero\(r_{1}-r_{2}=\) positivo.

Al igual que en el último caso, se necesita encontrar una segunda solución linealmente independiente.

\[\begin{aligned} &y_{1}(x)=|x|^{r_{1}} \sum_{n=0}^{\infty} c_{n} x^{n}, \quad c_{0}=1 \\ &y_{2}(x)=|x|^{r_{2}} \sum_{n=0}^{\infty} d_{n} x^{n}+\alpha y_{1}(x) \ln |x|, \quad d_{0}=1 \end{aligned} \nonumber \]

La constante\(\alpha\) puede determinarse posteriormente y en algunos casos puede desaparecer.

Para soluciones cercanas a puntos singulares regulares\(x=x_{0}\),, uno tiene un conjunto similar de casos pero para expansiones de la forma\(y_{1}(x)=\left|x-x_{0}\right|^{r_{1}} \sum_{n=0}^{\infty} c_{n}\left(x-x_{0}\right)^{n} .\)

Ejemplo\(\PageIndex{1}\)

\(x^{2} y^{\prime \prime}+x y^{\prime}+\left(x^{2}-\dfrac{1}{4}\right) y=0\).

En este ejemplo\(x=0\) se encuentra un punto singular. Tenemos\(a(x)=1 / x\) y\(b(x)=\left(x^{2}-1 / 4\right) / x^{2}\). Por lo tanto,

\[\begin{aligned} x a(x) &=1 \\ x^{2} b(x) &=x^{2}-\dfrac{1}{4} \end{aligned} \nonumber \]

Entonces,\(a_{0}=1\) y\(b_{0}=-\dfrac{1}{4} .\) La ecuación indicial se convierte

\[r(r-1)+r-\dfrac{1}{4}=0 \nonumber \]

Simplificando, tenemos\(0=r^{2}-\dfrac{1}{4}\), o\(r=\pm \dfrac{1}{2}\).

Para\(r=+\dfrac{1}{2}\), insertamos la serie\(y(x)=\sqrt{x} \sum_{n=0}^{\infty} c_{n} x^{n}\) en la ecuación diferencial, recogemos términos similares mediante la reindexación y encontramos

\[ \begin{aligned} 0 &=\sum_{n=0}^{\infty} c_{n}\left[\left(n+\dfrac{1}{2}\right)\left(n-\dfrac{1}{2}\right)+\left(n+\dfrac{1}{2}\right)+x^{2}-\dfrac{1}{4}\right] x^{n-\dfrac{3}{2}} \\ &=\sum_{n=0}^{\infty} c_{n}\left[n^{2}+n+x^{2}\right] x^{n-\dfrac{3}{2}} \\ &=\left[\sum_{n=0}^{\infty} n(n+1) c_{n} x^{n}+\sum_{n=0}^{\infty} c_{n} x^{n+2}\right] x^{-\dfrac{3}{2}} \\ &=\left[\sum_{k=0}^{\infty} k(k+1) c_{k} x^{k}+\sum_{n=2}^{\infty} c_{k-2} x^{k}\right] x^{-\dfrac{3}{2}} \\ &=2 c_{1} x+\sum_{n=2}^{\infty}\left[k(k+1) c_{k}+c_{k-2}\right] x^{k-\dfrac{3}{2}} \end{aligned} \label{4.42} \]

Esto da\(c_{1}=0\) y

\[c_{k}=-\dfrac{1}{k(k+1)} c_{k-2}, \quad k \geq 2 \nonumber \]

Iterando, tenemos\(c_{k}=0\) para\(k\) impar y

\[\begin{array}{ll} k=2: & c_{2}=-\dfrac{1}{3 !} c_{0} . \\ k=4: & c_{4}=-\dfrac{1}{20} c_{2}=\dfrac{1}{5 !} c_{0} . \\ k= 6: & c_{6}=-\dfrac{1}{42} c_{4}=\dfrac{1}{7 !} c_{0} . \end{array} \nonumber \]

Esto da

\[\begin{aligned} y_{1}(x) &=\sqrt{x} \sum_{n=0}^{\infty} c_{n} x^{n} \\ &=\sqrt{x}\left(c_{0}-\dfrac{1}{3 !} c_{0} x^{2}+\dfrac{1}{5 !} c_{0} x^{4}-\ldots\right) \\ &=\dfrac{c_{0}}{\sqrt{x}}\left(x-\dfrac{1}{3 !} x^{3}+\dfrac{1}{5 !} x^{5}-\ldots\right)=\dfrac{c_{0}}{\sqrt{x}} \sin x \end{aligned} \nonumber \]

De igual manera, para\(r=-\dfrac{1}{2}\), se obtiene la segunda solución,\(y_{2}(x)=\)\(\dfrac{d_{0}}{\sqrt{x}} \cos x .\) Setting\(c+0=1\) y\(d_{0}=1\), dar las dos soluciones linealmente independientes. Esta ecuación diferencial es la ecuación de Bessel de orden medio y las soluciones son funciones de Bessel de orden medio:

\[\begin{aligned} J_{\dfrac{1}{2}}(x) &=\sqrt{\dfrac{2}{\pi x}} \sin x, x>0 \\ J_{-\dfrac{1}{2}}(x) &=\sqrt{\dfrac{2}{\pi x}} \cos x, x>0 \end{aligned} \nonumber \]

Ejemplo\(\PageIndex{2}\)

\(x^{2} y^{\prime \prime}+3 x y^{\prime}+(1-2 x) y=0, x>0\)Para este problema\(x a(x)=3\) y\(x^{2} b(x)=1-2 x .\) Así, la ecuación indicial es

\[0=r(r-1)+3 r+1=(r+1)^{2} \nonumber \]

Este es un caso con dos raíces iguales,\(r=-1\). Una solución de la forma

\[y(x)=\sum_{n=0}^{\infty} c_{n} x^{n+r}\nonumber \]

sólo dará como resultado una solución. Utilizaremos esta forma de solución para llegar a una segunda solución linealmente independiente.

Todavía no vamos\(r=-1\) a insertar en el formulario de solución. Escribiendo la ecuación diferencial en forma de operador\(L[y]=0\),, tenemos

\[\begin{aligned} L[y]=& x^{2} y^{\prime \prime}+3 x y^{\prime}+(1-2 x) y \\ =& x^{2} \sum_{n=0}^{\infty} c_{n}(n+r)(n+r-1) x^{n+r-2}+3 x \sum_{n=0}^{\infty} c_{n}(n+r) x^{n+r-1} \\ &+(1-2 x) \sum_{n=0}^{\infty} c_{n} x^{n+r} \\ =& \sum_{n=0}^{\infty}[(n+r)(n+r-1)+2(n+r)+1] c_{n} x^{n+r}-\sum_{n=0}^{\infty} 2 c_{n} x^{n+r+1} \\ =& \sum_{n=0}^{\infty}(n+r+1)^{2} c_{n} x^{n+r}-\sum_{n=1}^{\infty} 2 c_{n-1} x^{n+r} \end{aligned} \nonumber \]

Estableciendo los coeficientes de términos similares iguales a cero para\(n \geq 1\), tenemos

\[c_{n}=\dfrac{2}{(n+r+1)^{2}} a_{n-1}, \quad n \geq 1.\nonumber \]

Iterando, encontramos

\[c_{n}=\dfrac{2^{n}}{[(r+2)(r+3) \cdots(r+n+1)]^{2}} c_{0}, \quad n \geq 1 \nonumber \]

Ajuste\(c_{0}=1\), tenemos la expresión

\[y(x, r)=x^{r}+\sum_{n=1}^{\infty} \dfrac{2^{n}}{[(r+2)(r+3) \cdots(r+n+1)]^{2}} x^{n+r}.\nonumber \]

Esta no es una solución de la ecuación diferencial porque no usamos la raíz\(r=-1\). En cambio, tenemos

\[L[y(x, c)]=(r+1)^{2} x^{r} \nonumber \]

del\(n=0\) término. Si\(r=-1\), entonces\(y(x,-1)\) es una solución de la ecuación diferencial. A saber,

\[y_{1}(x)=x^{-1}+\sum_{n=1}^{\infty} c_{n}(-1) x^{n-1}\nonumber \]

Ahora considere lo que sucede si diferenciamos Ecuación\(\PageIndex{4}\) con respecto a\(r\):

\[\begin{aligned} \dfrac{\partial}{\partial r} L[y(x, r)] &=\left[\dfrac{\partial y(x, r)}{\partial r}\right] \\ &=2(r+1) x^{r}+(r+1)^{2} x^{r} \ln x, \quad x>0 \end{aligned} \nonumber \]

Por lo tanto, también\(\dfrac{\partial y(x, r)}{\partial r}\) es una solución a la ecuación diferencial cuando\(r=-1\).

Desde

\[y(x, r)=x^{r}+\sum_{n=1}^{\infty} c_{n}(r) x^{n+r} \nonumber \]

tenemos

\[ \begin{aligned} \dfrac{\partial y(x, r)}{\partial r} &=x^{r} \ln x+\sum_{n=1}^{\infty} c_{n}(r) x^{n+r} \ln x+\sum_{n=1}^{\infty} c_{n}^{\prime}(r) x^{n+r} \\ &=y(x, r) \ln x+\sum_{n=1}^{\infty} c_{n}^{\prime}(r) x^{n+r} \end{aligned} \label{4.44} \]

Por lo tanto, la segunda solución viene dada por

\[y_{2}(x)=\left.\dfrac{\partial y(x, r)}{\partial r}\right|_{r=-1}=y(x,-1) \ln x+\sum_{n=1}^{\infty} c_{n}^{\prime}(-1) x^{n-1} \nonumber \]

Para determinar las soluciones, necesitamos evaluar\(c_{n}(-1)\) y\(c_{n}^{\prime}(-1) .\) Recordar eso (configuración\(c_{0}=1\))

\[c_{n}(r)=\dfrac{2^{n}}{[(r+2)(r+3) \cdots(r+n+1)]^{2}}, \quad n \geq 1 \nonumber \]

Por lo tanto,

\[ \begin{aligned} c_{n}(-1) &=\dfrac{2^{n}}{[(1)(2) \cdots(n)]^{2}}, \quad n \geq 1 \\ &=\dfrac{2^{n}}{(n !)^{2}} \end{aligned} \label{4.45} \]

A continuación calculamos\(c_{n}^{\prime}(-1)\). Esto se puede hacer usando diferenciación logarítmica. Consideramos

\[\begin{aligned} \ln c_{n}(r) &=\ln \left(\dfrac{2^{n}}{[(r+2)(r+3) \cdots(r+n+1)]^{2}}\right) \\ &=\ln 2^{n}-2 \ln (r+2)-2 \ln (r+3) \cdots-2 \ln (r+n+1) \end{aligned} \nonumber \]

Diferenciando con respecto\(r\) y evaluando en\(r=-1\), tenemos

\[ \begin{aligned} \dfrac{c_{n}^{\prime}(r)}{c_{n}(r)} &=-2\left(\dfrac{1}{r+2}+\dfrac{1}{r+2}+\cdots+\dfrac{1}{r+3}\right) \\ c_{n}^{\prime}(r) &=-2 c_{n}(r)\left(\dfrac{1}{r+2}+\dfrac{1}{r+3}+\cdots+\dfrac{1}{r+n+1}\right) . \\ c_{n}^{\prime}(-1) &=-2 c_{n}(-1)\left(\dfrac{1}{1}+\dfrac{1}{2}+\cdots+\dfrac{1}{n}\right) \\ &=-\dfrac{2^{n+1}}{(n !)^{2}} H_{n} \end{aligned} \label{4.46} \]

donde hemos definido

\[H_{n}=\dfrac{1}{1}+\dfrac{1}{2}+\cdots+\dfrac{1}{n}=\sum_{k=1}^{n} \dfrac{1}{k} \nonumber \]

Esto da la segunda solución como

\[y_{2}(x)=y_{1}(x) \ln x-\sum_{n=1}^{\infty} \dfrac{2^{n+1}}{(n !)^{2}} H_{n} x^{n-1} \nonumber \]

Ejemplo\(\PageIndex{3}\)

\(x^{2} y^{\prime \prime}+x y^{\prime}+\left(x^{2}-1\right) y=0 .\)

Esta ecuación es similar al último ejemplo, pero es la ecuación de Bessel de orden uno. La ecuación indicial viene dada por

\[0=r(r-1)+r-1=r^{2}-1 \nonumber \]

Las raíces son\(r_{1}=1, r_{2}=-1 .\) En este caso las raíces difieren en un entero,\(r_{1}-r_{2}=2\).

La primera solución se puede obtener usando

\[ \begin{aligned} y(x) &=\sum_{n=0}^{\infty} c_{n} x^{n+1} \\ y^{\prime}(x) &=\sum_{n=0}^{\infty} c_{n}(n+1) x^{n} \\ y^{\prime \prime}(x) &=\sum_{n=0}^{\infty} c_{n}(n+1)(n) x^{n-1} \end{aligned} \label{4.47} \]

Insertando estas series en la ecuación diferencial, tenemos

\[ \begin{aligned} 0 &=x^{2} y^{\prime \prime}+x y^{\prime}+\left(x^{2}-1\right) y \\ &=x^{2} \sum_{n=0}^{\infty} c_{n} n(n+1) x^{n-1}+x \sum_{n=0}^{\infty} c_{n}(n+1) x^{n} \\ &+\left(x^{2}-1\right) \sum_{n=0}^{\infty} c_{n} x^{n+1} \\ &=\sum_{n=0}^{\infty} c_{n} n(n+1) x^{n+1}+\sum_{n=0}^{\infty} c_{n}(n+1) x^{n+1} \\ & \sum_{n=0}^{\infty} c_{n} x^{n+3}-\sum_{n=0}^{\infty} c_{n} x^{n+1} \\ &=\sum_{n=0}^{\infty}[n(n+1)+(n+1)-1] c_{n} x^{n+1}+\sum_{n=0}^{\infty} c_{n} x^{n+3} \\ &=\sum_{n=0}^{\infty}\left[(n+1)^{2}-1\right] c_{n} x^{n+1}+\sum_{n=0}^{\infty} c_{n} x^{n+3} \end{aligned} \label{4.48} \]

Reindexamos la última suma con\(k=n+2\), o\(n=k-2\), para obtener

\[\sum_{n=2}^{\infty}\left(\left[(k+1)^{2}-1\right] c_{k}+c_{k-2}\right) x^{k+1}+3 c_{1} x^{2}=0 \nonumber \]

Así,\(c_{1}=0\) y

\[c_{k}=-\dfrac{1}{(k+1)^{2}-1} c_{k-2}=-\dfrac{1}{k(k+2)} c_{k-2}, \quad k \geq 2 \nonumber \]

Ya que\(c_{1}=0\), todos\(c_{k}{ }^{\prime} \mathrm{s}\) desaparecen por impar\(k .\) Para\(k=2 n\), tenemos

\[c_{2 n}=-\dfrac{1}{2 n(2 n+2)} c_{2 n-2}=-\dfrac{1}{4 n(n+1)} c_{2(n-1)}, \quad n \geq 1 .\nonumber \]

\[\begin{array}{ll} n=1: & c_{2}=-\dfrac{1}{4(1)(2)} c_{0} \\ n= & 2: \quad c_{4}=-\dfrac{1}{4(2)(3)} c_{2}=\dfrac{1}{4^{2} 2 ! 3 !} c_{0} \\ n= 3: & c_{6}=-\dfrac{1}{4(3)(4)} c_{4}=\dfrac{1}{4^{3} 3 ! 4 !} c_{0} . \end{array} \nonumber \]

Continuando, esto da

\[c_{2 n}=\dfrac{(-1)^{n}}{4^{n} n !(n+1) !} c_{0}\nonumber \]

y la primera solución es

\[y_{1}(x)=c_{0} \sum_{n=0}^{\infty} \dfrac{(-1)^{n}}{4^{n} n !(n+1) !} x^{2 n+1}\nonumber \]

Ahora buscamos una segunda solución linealmente independiente de la forma

\[y_{2}(x)=\sum_{n=0}^{\infty} d_{n} x^{n-1}+\alpha y_{1}(x) \ln x, \quad x>0 \nonumber \]

La función y sus derivadas están dadas por

\[\begin{aligned} y_{2}(x)=& \sum_{n=0}^{\infty} d_{n} x^{n-1}+\alpha y_{1}(x) \ln x \\ y_{2}^{\prime}(x)=& \sum_{n=0}^{\infty}(n-1) d_{n} x^{n-2}+\alpha\left[y_{1}^{\prime}(x) \ln x+y_{1}(x) x^{-1}\right] \\ y_{2}^{\prime \prime}(x)=& \sum_{n=0}^{\infty}(n-1) d_{n}(n-2) x^{n-3} \\ &+\alpha\left[y_{1}^{\prime \prime}(x) \ln x+2 y_{1}^{\prime}(x) x^{-1}-y_{1}(x) x^{-2}\right] \end{aligned} \nonumber \]

Insertando estas series en la ecuación diferencial, tenemos

\[\begin{aligned} 0=& x^{2} y_{2}^{\prime \prime}+x y_{2}^{\prime}+\left(x^{2}-1\right) y_{2} \\ =& x^{2} \sum_{n=0}^{\infty}(n-1)(n-2) d_{n} x^{n-3} \\ &+\alpha x^{2}\left[y_{1}^{\prime \prime}(x) \ln x+2 y_{1}^{\prime}(x) x^{-1}-y_{1}(x) x^{-2}\right] \\ &+x \sum_{n=0}^{\infty}(n-1) d_{n} x^{n-2}+\alpha x\left[y_{1}^{\prime}(x) \ln x+y_{1}(x) x^{-1}\right] \\ &+\left(x^{2}-1\right)\left[\sum_{n=0}^{\infty} d_{n} x^{n-1}+\alpha y_{1}(x) \ln x\right] \\ =& \sum_{n=0}^{\infty}[(n-1)(n-2)+(n-1)-1] d_{n} x^{n-1}+\sum_{n=0}^{\infty} d_{n} x^{n+1} \\ &+\alpha\left[x^{2} y_{1}^{\prime \prime}(x)+x y_{1}^{\prime}(x)+\left(x^{2}-1\right) y_{1}(x)\right] \ln x \\ &+\alpha\left[2 x y_{1}^{\prime}(x)-y_{1}(x)\right]+\alpha y_{1}(x) \\ =& \sum_{n=0}^{\infty} n(n-2) d_{n} x^{n-1}+\sum_{n=0}^{\infty} d_{n} x^{n+1}+2 \alpha x y_{1}^{\prime}(x) \\ & \sum_{n=0}^{\infty} n(n-2) d_{n} x^{n-1}+\sum_{n=0}^{\infty} d_{n} x^{n+1}+2 \alpha \sum_{n=0}^{\infty} \dfrac{(-1)^{n}(2 n+1)}{4^{n}(n+1) ! n !} x^{2 n+1} \end{aligned} \nonumber \]

Ahora tratamos de combinar como poderes de\(x\). Primero, combinamos los términos que involucran\(d_{n}^{\prime}\) s,

\[-d_{1}+d_{0} x+\sum_{k=3}^{\infty}\left[k(k-2) d_{k}+d_{k-2}\right] x^{k-1}=-2 \alpha \sum_{n=0}^{\infty} \dfrac{(-1)^{n}(2 n+1)}{4^{n} n !(n+1) !} x^{2 n+1} \nonumber \]

Como no hay potencias pares en el lado derecho de la ecuación, encontramos\(d_{1}=0\), y\(k(k-2) d_{k}+d_{k-2}=0, k \geq 3\) e\(k\) impar. Por lo tanto, todos\(d_{k}\) los impares desaparecen.

A continuación, establecemos\(k-1=2 n+1\), o\(k=2 n+2\), en los términos restantes, obtener

\[d_{0} x+\sum_{n=1}^{\infty}\left[(2 n+2)(2 n) d_{2 n+2}+d_{2 n}\right] x^{2 n+1}=-2 \alpha x-2 \alpha \sum_{n=1}^{\infty} \dfrac{(-1)^{n}(2 n+1)}{4^{n}(n+1) ! n !} x^{2 n+1} \nonumber \]

Así,\(d_{0}=-2 \alpha .\) elegimos\(\alpha=-\dfrac{1}{2}\), haciendo que\(d_{0}=1 .\) los términos restantes satisfagan la relación

\[(2 n+2)(2 n) d_{2 n+2}+d_{2 n}=\dfrac{(-1)^{n}(2 n+1)}{4^{n} n !(n+1) !}, \quad n \geq 1\nonumber \]

O

\[d_{2 n+2}=\dfrac{d_{2 n}}{4(n+1)(n)}+\dfrac{(-1)^{n}(2 n+1)}{(n+1)(n) 4^{n+1}(n+1) ! n !}, \quad n \geq 1 \nonumber \]

\[ \begin{aligned} & d_{4}=-\dfrac{1}{4(2)(1)} d_{2}-\dfrac{3}{(2)(1) 4^{2} 2 ! 1 !} \\ & =-\dfrac{1}{4^{2} 2 ! 1 !}\left(4 d_{2}+\dfrac{3}{2}\right) \\ & d_{6}=-\dfrac{1}{4(3)(2)} d_{4}+\dfrac{5}{(3)(2) 4^{3} 3 ! 2 !} \\ & =-\dfrac{1}{4(3)(2)}\left(-\dfrac{1}{4^{2} 2 ! 1 !}\left(4 d_{2}+\dfrac{3}{2}\right)\right)+\dfrac{5}{(3)(2) 4^{3} 3 ! 2 !} \\ & =\dfrac{1}{4^{3} 3 ! 2 !}\left(4 d_{2}+\dfrac{3}{2}+\dfrac{5}{6}\right) \\ & d_{8}=-\dfrac{1}{4(4)(3)} d_{6}-\dfrac{7}{(4)(3) 4^{4} 4 ! 3 !} \\ & =-\dfrac{1}{4^{4} 4 ! 3 !}\left(4 d_{2}+\dfrac{3}{2}+\dfrac{5}{6}\right)-\dfrac{7}{(4)(3) 4^{4} 4 ! 3 !} \\ & =-\dfrac{1}{4^{4} 4 ! 3 !}\left(4 d_{2}+\dfrac{3}{2}+\dfrac{5}{6}+\dfrac{7}{12}\right) . \end{aligned} \label{4.49} \]

Al elegir\(4 d_{2}=1\), los coeficientes toman una forma interesante. A saber,

\[ \begin{aligned} 1+\dfrac{3}{2} &=1+\dfrac{1}{2}+1 \\ 1+\dfrac{3}{2}+\dfrac{5}{6} &=1+\dfrac{1}{2} \dfrac{1}{3}+1+\dfrac{1}{2} \\ 1+\dfrac{3}{2}+\dfrac{5}{6}+\dfrac{7}{12} &=1+\dfrac{1}{2} \dfrac{1}{3}+\dfrac{1}{4}+1+\dfrac{1}{2}+\dfrac{1}{3} \end{aligned} \label{4.50} \]

Definir las sumas parciales de la serie armónica,

\[H_{n}=\sum_{k=1}^{n} \dfrac{1}{k}, \quad H_{0}=0 \nonumber \]

estos coeficientes se convierten\(H_{n}+H_{n-1}\) y los coeficientes en la expansión son

\[d_{2 n}=\dfrac{(-1)^{n-1}\left(H_{n}+H_{n-1}\right)}{4^{n} n !(n-1) !}, n=1,2, \ldots \nonumber \]

Podemos verificar esto computando\(d_{10}\):

\[\begin{aligned} d_{10} &=-\dfrac{1}{4(4)(3)} d_{8}+\dfrac{9}{(5)(4) 4^{5} 5 ! 4 !} \\ &=\dfrac{1}{4^{5} 5 ! 4 !}\left(4 d_{2}+\dfrac{3}{2}+\dfrac{5}{6}+\dfrac{7}{12}\right)+\dfrac{9}{(5)(4) 4^{5} 5 ! 4 !} \\ &=\dfrac{1}{4^{5} 5 ! 4 !}\left(4 d_{2}+\dfrac{3}{2}+\dfrac{5}{6}+\dfrac{7}{12}+\dfrac{9}{20}\right) \\ &=\dfrac{1}{4^{5} 5 ! 4 !}\left(1+\dfrac{1}{2}+\dfrac{1}{3}+\dfrac{1}{4}+\dfrac{1}{5}+1+\dfrac{1}{2}+\dfrac{1}{3}+\dfrac{1}{4}\right) \\ &=\dfrac{1}{4^{55} 5 ! 4 !}\left(H_{5}+H_{4}\right) \end{aligned} \nonumber \]

Esto da la segunda solución como

\(y_{2}(x)=\sum_{n=1}^{\infty} \dfrac{(-1)^{n-1}\left(H_{n}+H_{n-1}\right)}{4^{n} n !(n-1) !} x^{2 n-1}-\dfrac{1}{2} y_{1}(x) \ln x+x^{-1}\)