Sin matemáticas, no hay nada que puedas hacer. Todo a tu alrededor
es matemáticas. Todo a tu alrededor son números.
– Shakuntala Devi

Introducción

Hemos comenzado con el estudio de las ecuaciones diferenciales lineales de segundo orden con coeficientes variables. Ya hemos aprendido cómo obtener soluciones con respecto a puntos ordinarios, ahora aprenderemos a obtener soluciones con respecto a puntos singulares.

En la entrada anterior vimos que para resolver ecuaciones de la forma

$$\begin{array}{}\text{(1)}& {\displaystyle \frac{d^{2}y}{d{x}^2}}+P(x){\displaystyle \frac{dy}{dx}}+Q(x)y=0\end{array}$$

se proponía una solución de la forma

$$\begin{array}{}\text{(2)}& y(x)=\sum _{n=0}^{\mathrm{\infty }}{c}_n(x-x_0{)}^n\end{array}$$

donde $x_0$ es un punto ordinario de la ecuación diferencial ($\text{1}$).

En ocasiones no se pueden encontrar soluciones como ($\text{2}$), así que se propone una solución de la forma

$$\begin{array}{}\text{(3)}& y(x)=(x-x_0{)}^r\sum _{n=0}^{\mathrm{\infty }}{c}_n(x-x_0{)}^n\end{array}$$

Donde $r$ es una constante. En realidad, la solución ($\text{3}$) es una generalización ya que si $r=0$ regresamos a la forma ($\text{2}$).

En esta entrada aprenderemos a resolver ecuaciones diferenciales en las que su solución es de la forma ($\text{3}$).

Puntos singulares

El que la solución de una ecuación diferencial sea de la forma ($\text{3}$) esta directamente relacionado con que el punto $x_0$ sea un punto singular y no un punto ordinario. En la entrada anterior definimos estos conceptos, sin embargo en esta entrada es necesario profundizar más acerca de los puntos singulares. Recordemos la definición de punto singular.

Nota: Las siguientes definiciones se basan en la forma estándar ($\text{1}$) de una ecuación diferencial lineal de segundo orden.

Lo nuevo ahora es que un punto singular puede ser clasificado como regular o irregular.

Para fines prácticos en conveniente definir los puntos singulares regulares e irregulares a través de un límite.

Realicemos algunos ejemplos.

Ejemplo: Clasificar los puntos singulares de la ecuación diferencial

$$x^3(x^2-9){\displaystyle \frac{d^{2}y}{d{x}^2}}+(x+3){\displaystyle \frac{dy}{dx}}+(x-3{)}^{3}y=0$$

Solución: El primer paso es escribir a la ecuación diferencial en la forma estándar ($\text{1}$), así que dividimos toda la ecuación por el coeficiente de la segunda derivada de $y$.

$${\displaystyle \frac{d^{2}y}{d{x}^2}}+{\displaystyle \frac{1}{x^3(x-3)}}{\displaystyle \frac{dy}{dx}}+{\displaystyle \frac{(x-3{)}^2}{x^3(x+3)}}y=0$$

Identificamos que

$P(x)={\displaystyle \frac{1}{x^3(x-3)}}yQ(x)={\displaystyle \frac{(x-3{)}^2}{x^3(x+3)}}$

Notamos que la función $P(x)$ no está definida en los puntos $x=3$ y $x=0$, mientras que la función $Q(x)$ no está definida en $x=-3$ y $x=0$, de manera que los puntos singulares son $x_0=3$, $x_0=0$ y $x_0=-3$. El resto de puntos en $\mathbb{R}$ son puntos ordinarios de la ecuación diferencial.

Para determinar si son regulares o irregulares definamos las nuevas funciones de acuerdo a ($\text{4}$) y observemos si dichas funciones son analíticas o no en el correspondiente punto singular.

• Caso 1: $x_0=3$.

Definamos las nuevas funciones.

$$p(x)=(x-3)P(x)={\displaystyle \frac{1}{x^3}}yq(x)=(x-3{)}^{2}Q(x)={\displaystyle \frac{(x-3{)}^4}{x^3(x+3)}}$$

Es claro que las nuevas funciones $p(x)$ y $q(x)$ si son analíticas en $x_0=3$, por lo que dicho punto es un punto singular regular. Usando la definición de límite, se tiene

$$\underset{x\to 3}{lim}p(x)=\underset{x\to 3}{lim}{\displaystyle \frac{1}{x^3}}={\displaystyle \frac{1}{9}}y\underset{x\to 3}{lim}q(x)=\underset{x\to 3}{lim}{\displaystyle \frac{(x-3{)}^4}{x^3(x+3)}}=0$$

Los límites existen, así que llegamos a la misma conclusión.

• Caso 2: $x_0=0$.

Definamos las nuevas funciones.

$$p(x)=xP(x)={\displaystyle \frac{1}{x^2(x-3)}}yq(x)=x^{2}Q(x)={\displaystyle \frac{(x-3{)}^2}{x(x+3)}}$$

En este caso las funciones $p(x)$ y $q(x)$ siguen sin estar definidas para $x=0$ lo que significa que no se pueden representar mediante una serie de potencias, es decir, no son analíticas en dicho punto. Veamos que ocurre con los limites. Por un lado,

$$\underset{x\to 0}{lim}p(x)=\underset{x\to 0}{lim}{\displaystyle \frac{1}{x^2(x-3)}}=-\mathrm{\infty }$$

Por otro lado,

$$\underset{x\to 0^{+}}{lim}q(x)=\underset{x\to 0^{+}}{lim}{\displaystyle \frac{(x-3{)}^2}{x(x+3)}}=\mathrm{\infty }y\underset{x\to 0^{-}}{lim}q(x)=\underset{x\to 0^{-}}{lim}{\displaystyle \frac{(x-3{)}^2}{x(x+3)}}=-\mathrm{\infty }$$

Vemos que el limite de $p(x)$ es divergente, mientras que el límite de $q(x)$ no existe en $x=0$.

En conclusión, $x_0=0$ es un punto singular irregular.

• Caso 3: $x_0=-3$.

Definamos las nuevas funciones.

$$p(x)=(x+3)P(x)={\displaystyle \frac{x+3}{x^3(x-3)}}yq(x)=(x+3{)}^{2}Q(x)={\displaystyle \frac{(x-3{)}^2(x+3)}{x^3}}$$

Las nuevas funciones son analíticas en $x_0=-3$, confirmemos que los límites existen.

$$\underset{x\to -3}{lim}p(x)=\underset{x\to -3}{lim}{\displaystyle \frac{x+3}{x^3(x-3)}}=0y\underset{x\to -3}{lim}q(x)=\underset{x\to -3}{lim}{\displaystyle \frac{(x-3{)}^2(x+3)}{x^3}}=0$$

En efecto, los limites existen, así que $x_0=-3$ es un punto singular regular.

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Realicemos un ejemplo más.

Ejemplo: Determinar el punto singular de la ecuación diferencial

$$(x+1{)}^2{\displaystyle \frac{d^{2}y}{d{x}^2}}+x{\displaystyle \frac{dy}{dx}}+x^{2}y=0$$

Solución: Escribimos a la ecuación diferencial en su forma estándar.

$${\displaystyle \frac{d^{2}y}{d{x}^2}}+{\displaystyle \frac{x}{(x+1{)}^2}}{\displaystyle \frac{dy}{dx}}+{\displaystyle \frac{x^2}{(x+1{)}^2}}y=0$$

Identificamos que

$$P(x)={\displaystyle \frac{x}{(x+1{)}^2}}yQ(x)={\displaystyle \frac{x^2}{(x+1{)}^2}}$$

Notamos que el único punto singular es $x_0=-1$. Definamos las funciones $p(x)$ y $q(x)$.

$$p(x)=(x+1)P(x)={\displaystyle \frac{x}{x+1}}yq(x)=(x+1{)}^{2}Q(x)=x^2$$

Aunque la función $q(x)$ si es analítica en $x_0=-1$, $p(x)$ no lo es. Por lo tanto, la ecuación diferencial no se puede desarrollar en potencias de $x+1$ y por definición $x_0=-1$ es un punto singular irregular.

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Solución a ecuaciones diferenciales

Ahora que sabemos identificar puntos singulares de una ecuación diferencial podemos resolverlas con respecto a dichos puntos proponiendo una solución de la forma ($\text{3}$). Ahora bien, debido a la complejidad de los cálculos, sólo estudiaremos el caso en el que el punto $x_0=0$ es un punto singular regular.

A continuación enunciamos el teorema que establece que ($\text{3}$) es una solución de la ecuación diferencial ($\text{1}$) con respecto al punto singular $x_0$.

Con este teorema podemos establecer lo siguiente:

• Si $x_0$ es un punto ordinario, entonces $r=0$ y ($\text{2}$) es la solución general.

• Si $x_0$ es un punto singular regular, entonces ($\text{6}$) dará una solución o la solución general.

• Si $x_0$ es un punto singular irregular, entonces pueden o no existir soluciones de la forma ($\text{6}$).

No demostraremos este teorema, pero será la base para resolver ecuaciones diferenciales.

La manera de resolver ecuaciones diferenciales con respecto a puntos singulares es bastante similar al caso de soluciones con respecto a puntos ordinarios, sin embargo en este caso, además de obtener una relación de recurrencia, obtendremos una ecuación cuadrática para $r$ que deberemos de resolver, a dicha ecuación se le conoce como ecuación indicial.

A continuación desarrollaremos el método de resolución que nos permitirá obtener la expresión general de la ecuación indicial, dicho método se conoce como método de Frobenius.

Método de Frobenius

Queremos resolver una ecuación diferencial en su forma estándar con respecto al punto singular regular $x_0=0$.

$${\displaystyle \frac{d^{2}y}{d{x}^2}}+P(x){\displaystyle \frac{dy}{dx}}+Q(x)y=0$$

Multipliquemos esta ecuación por $x^2$.

$$x^2{\displaystyle \frac{d^{2}y}{d{x}^2}}+x[xP(x)]{\displaystyle \frac{dy}{dx}}+[x^{2}Q(x)]y=0$$

Si usamos las definiciones ($\text{4}$) para $x_0=0$, entonces podemos escribir la ecuación anterior de la siguiente manera.

$$\begin{array}{}\text{(7)}& x^2{\displaystyle \frac{d^{2}y}{d{x}^2}}+xp(x){\displaystyle \frac{dy}{dx}}+q(x)y=0\end{array}$$

Con $p(x)$ y $q(x)$ son funciones analíticas en $x=0$, esto significa que se pueden representar mediante una serie de potencias con respecto a dicho punto, sean

$$\begin{array}{}\text{(8)}& p(x)=p_0+p_{1}x+p_2{x}^2+\cdots =\sum _{n=0}^{\mathrm{\infty }}{p}_n{x}^n\end{array}$$

y

$$\begin{array}{}\text{(9)}& q(x)=q_0+q_{1}x+q_2{x}^2+\cdots =\sum _{n=0}^{\mathrm{\infty }}{q}_n{x}^n\end{array}$$

dichas series. Una observación interesante es que si todos los coeficientes son cero excepto $p_0$ y $q_0$, entonces recuperamos la ecuación de Cauchy – Euler.

$$\begin{array}{}\text{(10)}& x^2{\displaystyle \frac{d^{2}y}{d{x}^2}}+p_{0}x{\displaystyle \frac{dy}{dx}}+q_{0}y=0\end{array}$$

El teorema anterior nos indica que la forma de la solución es

$$y(x)=\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r}$$

La primera y segunda derivada son

$${\displaystyle \frac{dy}{dx}}=\sum _{n=0}^{\mathrm{\infty }}(n+r)c_n{x}^{n+r-1}y{\displaystyle \frac{d^{2}y}{d{x}^2}}=\sum _{n=0}^{\mathrm{\infty }}(n+r)(n+r-1)c_n{x}^{n+r-2}$$

Sustituyamos todos estos resultados en la ecuación diferencial ($\text{7}$).

$$x^2\sum _{n=0}^{\mathrm{\infty }}(n+r)(n+r-1)c_n{x}^{n+r-2}+x\left[\sum _{n=0}^{\mathrm{\infty }}{p}_n{x}^n\right]\sum _{n=0}(n+r)c_n{x}^{n+r-1}+\left[\sum _{n=0}^{\mathrm{\infty }}{q}_n{x}^n\right]\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r}=0$$

Introducimos los términos $x^2$ y $x$ a las series de las derivadas de $y$.

$$\sum _{n=0}^{\mathrm{\infty }}(n+r)(n+r-1)c_n{x}^{n+r}+\left[\sum _{n=0}^{\mathrm{\infty }}{p}_n{x}^n\right]\sum _{n=0}(n+r)c_n{x}^{n+r}+\left[\sum _{n=0}^{\mathrm{\infty }}{q}_n{x}^n\right]\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r}=0$$

Tomemos los términos para $n=0$.

$$\begin{array}{rl}r(r-1)c_0{x}^r+p_{0}r{c}_0{x}^r+q_0{c}_0{x}^r& =0\\ c_0{x}^r[r(r-1)+p_{0}r+q_0]& =0\end{array}$$

Sabemos que $x^r\ne 0$ y el método nos obliga a considerar que siempre $c_0\ne 0$, entonces

$$r(r-1)+p_{0}r+q_0=0$$

o bien,

$$\begin{array}{}\text{(11)}& r^2+(p_0-1)r+q_0=0\end{array}$$

Esta relación corresponde a la ecuación indicial con raíces $r_1$ y $r_2$ reales. En todos los casos se le asigna a $r_1$ la raíz mayor, es decir, debe ocurrir que $r_1>r_2$, siempre y cuando no sean raíces repetidas. A las raíces $r_1$ y $r_2$ se les denomina raíces indiciales.

El siguiente paso en el método es continuar igualando cada término a cero a través de una relación de recurrencia y con ello determinar los coeficientes de la solución propuesta $y(x)$, todo de manera similar que en el método de la entrada anterior.

En el enunciado del teorema enfatizamos que hay al menos una solución, esto significa que no siempre puede obtenerse una segunda serie solución que junto con la primera serie forme la solución general de la ecuación diferencial. No lo demostraremos, pero a continuación se muestra la forma de ambas soluciones linealmente independientes de acuerdo a los casos que pueden ocurrir con las raíces indiciales.

De acuerdo a la ecuación indicial ($\text{11}$) se distinguen tres casos:

• Caso 1: $r_1-r_2\ne$ número entero.

En este caso las soluciones de la ecuación diferencial ($\text{1}$), son

$$\begin{array}{}\text{(12)}& y_1(x)=\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r_1},c_0\ne 0\end{array}$$

$$\begin{array}{}\text{(13)}& y_2(x)=\sum _{n=0}^{\mathrm{\infty }}{\hat{c}}_n{x}^{n+r_2},{\hat{c}}_0\ne 0\end{array}$$

• Caso 2: $r_1=r_2=r$.

En el caso en el que ambas raíces indiciales son iguales, las soluciones de la ecuación diferencial ($\text{1}$), son

$$\begin{array}{}\text{(14)}& y_1(x)=\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r},c_0\ne 0\end{array}$$

$$\begin{array}{}\text{(15)}& y_2(x)=y_1(x)\ln (x)+\sum _{n=0}^{\mathrm{\infty }}{\hat{c}}_n{x}^{n+r}\end{array}$$

• Caso 3: $r_1-r_2=$ entero positivo.

En este caso las soluciones de la ecuación diferencial ($\text{1}$), son

$$\begin{array}{}\text{(16)}& y_1(x)=\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r_1},c_0\ne 0\end{array}$$

$$\begin{array}{}\text{(17)}& y_2(x)=C{y}_1(x)\ln (x)+\sum _{n=0}^{\mathrm{\infty }}{\hat{c}}_n{x}^{n+r_2},{\hat{c}}_0\ne 0\end{array}$$

Donde $C$ es una constante que podría ser cero.

En todos los casos $y_1(x)$ y $y_2(x)$ son linealmente independientes. Por lo tanto, la solución general es

$$\begin{array}{}\text{(18)}& y(x)=C_1{y}_1(x)+C_2{y}_2(x)\end{array}$$

En los casos en los que el método de Frobenius no nos de una segunda solución es posible obtenerla con métodos que ya hemos estudiado antes. El primero de ellos es usar variación de parámetros, en este caso se propone la solución

$$y_2(x)=u(x)y_1(x)$$

y se sustituye, junto con las derivadas correspondientes, en la ecuación diferencial, esto nos permitirá obtener una ecuación diferencial para $u(x)$ que debemos resolver.

Otro método es usar directamente la forma de las soluciones $y_2(x)$ propuestas anteriormente para cada caso, calcular las derivadas correspondientes y sustituir en la ecuación diferencial.

Un tercer método se puede aplicar una vez que ya hemos determinado la primer solución $y_1(x)$ y es usando la expresión que deducimos en entradas anteriores.

$$\begin{array}{}\text{(19)}& y_2(x)=y_1(x)\int {\displaystyle \frac{e^{-\int P(x)dx}}{y_1^2(x)}}dx\end{array}$$

Una de las mejores maneras para comprender algo es a través de ejemplos y práctica, así que hemos decidido resolver tres ejemplos, uno para cada caso y así poder comprender del todo en qué consiste el método de Frobenius.

Cabe mencionar que a lo largo de esta entrada hemos dado las herramientas para trabajar, pero no se ha dado un fundamento formal de los resultados, para conocerlos se pueden revisar los videos del tema correspondiente en la sección de videos de este curso, en él se encontrarán los fundamentos de cómo es que se obtienen las soluciones linealmente independientes dadas para cada condición de las raíces indiciales.

Para concluir esta entrada realicemos los 3 ejemplos antes mencionados.

Solución cuando la diferencia de las raíces indiciales difiere de un número entero

Ejemplo: Resolver la ecuación diferencial

$$3x^2{\displaystyle \frac{d^{2}y}{d{x}^2}}-x{\displaystyle \frac{dy}{dx}}+(1-x)y=0$$

con respecto al punto singular $x_0=0$.

Solución: Dividimos la ecuación diferencial por el coeficiente de la segunda derivada de $y$ para obtener la forma estándar.

$${\displaystyle \frac{d^{2}y}{d{x}^2}}-{\displaystyle \frac{1}{3x}}{\displaystyle \frac{dy}{dx}}+{\displaystyle \frac{1-x}{3x^2}}y=0$$

Identificamos que

$$P(x)=-{\displaystyle \frac{1}{3x}}yQ(x)={\displaystyle \frac{1-x}{3x^2}}$$

Ninguna función está definida en $x=0$. Definimos las funciones $p(x)$ y $q(x)$ de acuerdo a ($\text{4}$).

$$p(x)=-{\displaystyle \frac{1}{3}}yq(x)={\displaystyle \frac{1-x}{3}}$$

Vemos que

$$\underset{x\to 0}{lim}p(x)=-{\displaystyle \frac{1}{3}}y\underset{x\to 0}{lim}q(x)={\displaystyle \frac{1}{3}}$$

Esto nos muestra que $p(x)$ y $q(x)$ son analíticas en $x=0$ y que dicho punto es un punto singular regular.

Obtendremos la ecuación indicial directamente de la expresión ($\text{11}$).

Vemos que

$$p(x)=\sum _{n=0}^{\mathrm{\infty }}{p}_n{x}^n=p_0+p_{1}x+p_2{x}^2+\cdots =-{\displaystyle \frac{1}{3}}$$

de donde,

$$p_0=-{\displaystyle \frac{1}{3}}$$

y $p_k=0$ $\mathrm{\forall }$ $k⩾1$ con $k\in \mathbb{N}$. Por otro lado

$$q(x)=\sum _{n=0}^{\mathrm{\infty }}{q}_n{x}^n=q_0+q_{1}x+q_2{x}^2+\cdots ={\displaystyle \frac{1}{3}}-{\displaystyle \frac{1}{3}}x$$

de donde,

$$q_0={\displaystyle \frac{1}{3}}y{q}_1=-{\displaystyle \frac{1}{3}}$$

y $q_k=0$ $\mathrm{\forall }$ $k⩾2$ con $k\in \mathbb{N}$.

Sustituimos $p_0$ y $q_0$ en la ecuación indicial ($\text{11}$).

$$r^2+(-{\displaystyle \frac{1}{3}}-1)r+{\displaystyle \frac{1}{3}}=r^2-{\displaystyle \frac{4}{3}}r+{\displaystyle \frac{1}{3}}=0$$

Resolviendo para $r$ se obtiene que las raíces son

$$r_1=1y{r}_2={\displaystyle \frac{1}{3}}$$

Notemos que

$$r_1-r_2={\displaystyle \frac{2}{3}}$$

es decir, la diferencia de las raíces indiciales difiere de un número entero, esto nos indica que estamos en condiciones del caso 1, en donde las soluciones están dadas por las funciones ($\text{12}$) y ($\text{13}$).

$$y_1(x)=\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+1},c_0\ne 0$$

y

$$y_2(x)=\sum _{n=0}^{\mathrm{\infty }}{\hat{c}}_n{x}^{n+1/3},{\hat{c}}_0\ne 0$$

Para continuar con el método de Frobenius consideremos la solución general

$$y(x)=\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r}$$

y sus derivadas

$${\displaystyle \frac{dy}{dx}}=\sum _{n=0}^{\mathrm{\infty }}(n+r)c_n{x}^{n+r-1}y{\displaystyle \frac{d^{2}y}{d{x}^2}}=\sum _{n=0}^{\mathrm{\infty }}(n+r)(n+r-1)c_n{x}^{n+r-2}$$

Una vez obtenida la relación de recurrencia ya se podrá sustituir los valores correspondientes de $r$. Sustituyamos en la ecuación diferencial.

$$3x^2[\sum _{n=0}^{\mathrm{\infty }}(n+r)(n+r-1)c_n{x}^{n+r-2}]-x[\sum _{n=0}^{\mathrm{\infty }}(n+r)c_n{x}^{n+r-1}]+(1-x)\left[\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r}\right]=0$$

$$3\sum _{n=0}^{\mathrm{\infty }}(n+r)(n+r-1)c_n{x}^{n+r}-\sum _{n=0}^{\mathrm{\infty }}(n+r)c_n{x}^{n+r}+\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r}-\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r+1}=0$$

En la última serie hacemos $k=n+1$ y en el resto $k=n$.

$$3\sum _{k=0}^{\mathrm{\infty }}(k+r)(k+r-1)c_k{x}^{k+r}-\sum _{k=0}^{\mathrm{\infty }}(k+r)c_k{x}^{k+r}+\sum _{k=0}^{\mathrm{\infty }}{c}_k{x}^{k+r}-\sum _{k=1}^{\mathrm{\infty }}{c}_{k-1}{x}^{k+r}=0$$

Para que todas las series comiencen en $k=1$ extraemos el primer término de las tres primeras series y la suma la igualamos a cero.

$$\begin{array}{rl}3r(r-1)c_0{x}^r-r{c}_0{x}^r+c_0{x}^r& =0\\ c_0{x}^r[3r(r-1)-r+1]& =0\end{array}$$

Como $x^r\ne 0$ y $c_0\ne 0$, entonces

$$\begin{array}{rl}3r(r-1)-r+1& =0\\ 3r^2-4r+1& =0\end{array}$$

Con este otro procedimiento podemos obtener la ecuación indicial. Ahora nos queda la ecuación

$$3\sum _{k=1}^{\mathrm{\infty }}(k+r)(k+r-1)c_k{x}^{k+r}-\sum _{k=1}^{\mathrm{\infty }}(k+r)c_k{x}^{k+r}+\sum _{k=1}^{\mathrm{\infty }}{c}_k{x}^{k+r}-\sum _{k=1}^{\mathrm{\infty }}{c}_{k-1}{x}^{k+r}=0$$

Podemos juntar todas las series en una sola.

$$\sum _{k=1}^{\mathrm{\infty }}[3(k+r)(k+r-1)c_k-(k+r)c_k+c_k-c_{k-1}]x^{k+r}=0$$

Para satisfacer la igualdad es necesario que

$$c_k[3(k+r)(k+r-1)-(k+r)+1]-c_{k-1}=0$$

Despejando a $c_k$ obtenemos la relación de recurrencia.

$$c_k={\displaystyle \frac{c_{k-1}}{3(k+r)(k+r-1)-(k+r)+1}},k=1,2,3,\cdots$$

Hay que determinar los coeficientes para cada valor de las raíces indiciales. Para el valor de la primer raíz indicial $r=1$, la relación de recurrencia es

$$c_k={\displaystyle \frac{c_{k-1}}{k(3k+2)}},k=1,2,3,\cdots$$

Determinemos los coeficientes para este caso.

$k=1$.

$$c_1={\displaystyle \frac{c_0}{1(3(1)+2)}}={\displaystyle \frac{c_0}{5}}$$

$k=2$.

$$c_2={\displaystyle \frac{c_1}{2(3(2)+2)}}={\displaystyle \frac{c_1}{16}}={\displaystyle \frac{c_0}{80}}$$

$k=3$.

$$c_3={\displaystyle \frac{c_2}{3(3(3)+2)}}={\displaystyle \frac{c_2}{33}}={\displaystyle \frac{c_0}{2640}}$$

$k=4$.

$$c_4={\displaystyle \frac{c_3}{4(3(4)+2)}}={\displaystyle \frac{c_3}{56}}={\displaystyle \frac{c_0}{147840}}$$

Etcétera, entonces la primer solución es de la forma

$$\begin{array}{rl}{y}_1(x)& =x^1(c_0+c_{1}x+c_2{x}^2+c_3{x}^3+c_4{x}^4+\cdots )\\ & =x(c_0+{\displaystyle \frac{c_0}{5}}x+{\displaystyle \frac{c_0}{80}}{x}^2+{\displaystyle \frac{c_0}{2640}}{x}^3+{\displaystyle \frac{c_0}{147840}}{x}^4+\cdots )\\ & =c_{0}x(1+{\displaystyle \frac{x}{5}}+{\displaystyle \frac{x^2}{80}}+{\displaystyle \frac{x^3}{2640}}+{\displaystyle \frac{x^4}{147840}}+\cdots )\end{array}$$

Por otro lado, para $r={\displaystyle \frac{1}{3}}$ la relación de recurrencia es

$${\hat{c}}_k={\displaystyle \frac{{\hat{c}}_{k-1}}{k(3k-2)}},k=1,2,3,\cdots$$

Usamos la notación ${\hat{c}}_k$ sólo para hacer referencia de que son los coeficientes de la segunda solución, pero se obtiene de la misma relación de recurrencia obtenida por el método, sólo que ahora usando $r={\displaystyle \frac{1}{3}}$.

Determinemos los coeficientes para este caso.

$k=1$.

$${\hat{c}}_1={\displaystyle \frac{{\hat{c}}_0}{1(3(1)-2)}}={\hat{c}}_0$$

$k=2$.

$${\hat{c}}_2={\displaystyle \frac{{\hat{c}}_1}{2(3(2)-2)}}={\displaystyle \frac{{\hat{c}}_0}{8}}$$

$k=3$.

$${\hat{c}}_3={\displaystyle \frac{{\hat{c}}_2}{3(3(3)-2)}}={\displaystyle \frac{{\hat{c}}_2}{21}}={\displaystyle \frac{{\hat{c}}_0}{168}}$$

$k=4$.

$${\hat{c}}_4={\displaystyle \frac{{\hat{c}}_3}{4(3(4)-2)}}={\displaystyle \frac{{\hat{c}}_3}{40}}={\displaystyle \frac{{\hat{c}}_0}{6720}}$$

Etcétera, entonces la segunda solución es de la forma

$$\begin{array}{rl}{y}_2(x)& =x^{1/3}({\hat{c}}_0+{\hat{c}}_{1}x+{\hat{c}}_2{x}^2+{\hat{c}}_3{x}^3+{\hat{c}}_4{x}^4+\cdots )\\ & =x^{1/3}({\hat{c}}_0+{\hat{c}}_{0}x+{\displaystyle \frac{{\hat{c}}_0}{8}}{x}^2+{\displaystyle \frac{{\hat{c}}_0}{168}}{x}^3+{\displaystyle \frac{{\hat{c}}_0}{6720}}{x}^4+\cdots )\\ & ={\hat{c}}_0{x}^{1/3}(1+x+{\displaystyle \frac{x^2}{8}}+{\displaystyle \frac{x^3}{168}}+{\displaystyle \frac{x^4}{6720}}+\cdots )\end{array}$$

Si definimos $C_1=c_0$ y $C_2={\hat{c}}_0$, entonces la solución general de la ecuación diferencial es

$$y(x)=C_{1}x(1+{\displaystyle \frac{x}{5}}+{\displaystyle \frac{x^2}{80}}+{\displaystyle \frac{x^3}{2640}}+{\displaystyle \frac{x^4}{147840}}+\cdots )+C_2{x}^{1/3}(1+x+{\displaystyle \frac{x^2}{8}}+{\displaystyle \frac{x^3}{168}}+{\displaystyle \frac{x^4}{6720}}+\cdots )$$

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Con este ejemplo podemos aclarar algunas cosas.

La primera de ellas es que desarrollando el método mismo obtendremos la ecuación indicial, así que no necesariamente debemos sustituir en la ecuación ($\text{11}$), sin embargo sustituir en la ecuación ($\text{11}$) nos permitirá, desde un inicio, conocer las raíces indiciales y con ello podremos determinar la forma de la segunda solución según sea el caso.

Otra cosa importante es que se pueden calcular los coeficientes que se deseen, en el ejemplo sólo calculamos los primeros $5$ coeficientes, es decir hasta $k=4$, pero se puede continuar, lo interesante de continuar es que en algunas ocasiones es posible determinar una relación que generaliza la forma de los coeficientes y con ello formar una serie que incluso puede converger a una función conocida. Los siguientes ejercicios son un ejemplo de esto.

También hay que mencionar que en este ejemplo el método de Frobenius sí nos proporcionó la segunda solución usando la relación de recurrencia, esto no ocurrirá en algunos otros casos, como el que sigue a continuación, en estos casos será necesario aplicar algunos de los métodos que ya mencionamos antes.

Solución cuando las raíces indiciales son repetidas

Ejemplo: Resolver la ecuación diferencial

$$x^2{\displaystyle \frac{d^{2}y}{d{x}^2}}+(x^2-x){\displaystyle \frac{dy}{dx}}+y=0$$

con respecto al punto singular $x_0=0$.

Solución: Escribimos la ecuación en su forma estándar.

$${\displaystyle \frac{d^{2}y}{d{x}^2}}+{\displaystyle \frac{x-1}{x}}{\displaystyle \frac{dy}{dx}}+{\displaystyle \frac{1}{x^2}}y=0$$

Identificamos que

$$P(x)={\displaystyle \frac{x-1}{x}}yQ(x)={\displaystyle \frac{1}{x^2}}$$

Mientras que las funciones $p(x)$ y $q(x)$ están dadas por

$$p(x)=x-1yq(x)=1$$

Como los límites existen

$$\underset{x\to 0}{lim}p(x)=-1y\underset{x\to 0}{lim}q(x)=1$$

entonces $x=0$ es un punto singular regular. En esta ocasión vamos a obtener las raíces indiciales directamente de la expresión resultante para $k=0$. Sustituyamos las funciones correspondientes en la ecuación diferencial.

$$x^2[\sum _{n=0}^{\mathrm{\infty }}(n+r)(n+r-1)c_n{x}^{n+r-2}]+(x^2-x)[\sum _{n=0}^{\mathrm{\infty }}(n+r)c_n{x}^{n+r-1}]+\left[\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r}\right]=0$$

$$\sum _{n=0}^{\mathrm{\infty }}(n+r)(n+r-1)c_n{x}^{n+r}+\sum _{n=0}^{\mathrm{\infty }}(n+r)c_n{x}^{n+r+1}-\sum _{n=0}^{\mathrm{\infty }}(n+r)c_n{x}^{n+r}+\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r}=0$$

En la segunda serie hacemos $k=n+1$ y en el resto $k=n$.

$$\sum _{k=0}^{\mathrm{\infty }}(k+r)(k+r-1)c_k{x}^{k+r}+\sum _{k=1}^{\mathrm{\infty }}(k-1+r)c_{k-1}{x}^{k+r}-\sum _{k=0}^{\mathrm{\infty }}(k+r)c_k{x}^{k+r}+\sum _{k=0}^{\mathrm{\infty }}{c}_k{x}^{k+r}=0$$

Extraemos el primer término de las series que comienzan con $k=0$ para que todas comiencen con $k=1$ y la suma de dichos términos la igualamos a cero.

$$\begin{array}{rl}r(r-1)c_0{x}^r-r{c}_0{x}^r+c_0{x}^r& =0\\ c_0{x}^r[r(r-1)-r+1]& =0\end{array}$$

como $x^r\ne 0$ y $c_0\ne 0$, entonces

$$\begin{array}{rl}r(r-1)-r+1& =0\\ r^2-2r+1& =0\end{array}$$

Hemos obtenido la ecuación indicial. Resolviendo para $r$ se obtiene que

$$r_1=r_2=1$$

Las raíces indiciales son iguales, de manera que estamos en condiciones del caso 2 en el que las soluciones son de la forma ($\text{14}$) y ($\text{15}$).

$$y_1(x)=\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+1},c_0\ne 0$$

y

$$y_2(x)=\ln (x)\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+1}+\sum _{n=0}^{\mathrm{\infty }}{\hat{c}}_n{x}^{n+1}$$

Ahora tenemos la ecuación en la que todas las series tienen la misma potencia y comienzan con el mismo índice.

$$\sum _{k=1}^{\mathrm{\infty }}(k+r)(k+r-1)c_k{x}^{k+r}+\sum _{k=1}^{\mathrm{\infty }}(k-1+r)c_{k-1}{x}^{k+r}-\sum _{k=1}^{\mathrm{\infty }}(k+r)c_k{x}^{k+r}+\sum _{k=1}^{\mathrm{\infty }}{c}_k{x}^{k+r}=0$$

Juntamos todo en una sola serie.

$$\sum _{k=1}^{\mathrm{\infty }}[(k+r)(k+r-1)c_k+(k-1+r)c_{k-1}-(k+r)c_k+c_k]x^{k+r}=0$$

de donde

$$c_k[(k+r)(k+r-1)-(k+r)+1]+c_{k-1}(k-1+r)=0$$

despejando a $c_k$ se obtiene la relación de recurrencia.

$$c_k={\displaystyle \frac{c_{k-1}(k-1+r)}{(k+r)-1-(k+r)(k+r-1)}}={\displaystyle \frac{c_{k-1}}{1-k-r}},k=1,2,3,\cdots$$

Cómo $r=1$, entonces la relación de recurrencia es

$$c_k=-{\displaystyle \frac{c_{k-1}}{k}},k=1,2,3,\cdots$$

Determinemos los coeficientes.

$k=1$.

$$c_1=-{\displaystyle \frac{c_0}{1}}=-c_0$$

$k=2$.

$$c_2=-{\displaystyle \frac{c_1}{2}}={\displaystyle \frac{c_0}{2}}$$

$k=3$.

$$c_3=-{\displaystyle \frac{c_2}{3}}=-{\displaystyle \frac{c_0}{6}}$$

$k=4$.

$$c_4=-{\displaystyle \frac{c_3}{4}}={\displaystyle \frac{c_0}{24}}$$

$k=5$.

$$c_5=-{\displaystyle \frac{c_4}{5}}=-{\displaystyle \frac{c_0}{120}}$$

Etcétera, la primera solución es

$$\begin{array}{rl}{y}_1(x)& =x(c_0+c_{1}x+c_2{x}^2+c_3{x}^3+c_4{x}^4+c_5{x}^5+\cdots )\\ & =x(c_0-c_{0}x+{\displaystyle \frac{c_0}{2}}{x}^2-{\displaystyle \frac{c_0}{6}}{x}^3+{\displaystyle \frac{c_0}{24}}{x}^4-{\displaystyle \frac{c_0}{120}}{x}^5+\cdots )\\ & =c_{0}x(1-x+{\displaystyle \frac{x^2}{2}}-{\displaystyle \frac{x^3}{3!}}+{\displaystyle \frac{x^4}{4!}}-{\displaystyle \frac{x^5}{5!}}+\cdots )\end{array}$$

Lo que mencionamos antes, la solución va teniendo forma de una serie que conocemos, pues sabemos que

$$e^{-x}=\sum _{n=0}^{\mathrm{\infty }}{\displaystyle \frac{(-x{)}^n}{n!}}=1-x+{\displaystyle \frac{x^2}{2}}-{\displaystyle \frac{x^3}{3!}}+{\displaystyle \frac{x^4}{4!}}-{\displaystyle \frac{x^5}{5!}}+\cdots$$

Entonces,

$$y_1(x)=c_{0}x{e}^{-x}$$

Consideremos que $c_0=1\ne 0$, así la primer solución de la ecuación diferencial es

$$y_1(x)=x{e}^{-x}$$

Notemos que el método ya no nos ofrece una segunda solución. Para obtener la segunda solución se pueden usar los tres métodos antes mencionados. Uno de ellos es usando variación de parámetros. Un segundo método puede ser por derivación de la solución propuesta

$$y_2(x)=y_1\ln (x)+\sum _{n=0}^{\mathrm{\infty }}{\hat{c}}_n{x}^{n+1}$$

las derivadas son

$${\displaystyle \frac{d{y}_2}{dx}}={\displaystyle \frac{y_1}{x}}+\ln (x){\displaystyle \frac{d{y}_1}{dx}}+\sum _{n=0}^{\mathrm{\infty }}(n+1){\hat{c}}_n{x}^n$$

y

$${\displaystyle \frac{d^2{y}_2}{d{x}^2}}=-{\displaystyle \frac{y_1}{x^2}}+{\displaystyle \frac{2}{x}}{\displaystyle \frac{d{y}_1}{dx}}+\ln (x){\displaystyle \frac{d^2{y}_1}{d{x}^2}}+\sum _{n=0}^{\mathrm{\infty }}(n+1)n{\hat{c}}_n{x}^{n-1}$$

Se sustituyen estos resultados en la ecuación diferencial y se procede igual que antes con la diferencia de que ahora no obtendremos una ecuación indicial, pero sí una relación de recurrencia para obtener los coeficientes ${\hat{c}}_k$. ¡Seguro este método es un camino largo!.

Un tercer método es aplicar directamente la formula ($\text{19}$). Debido a que este es el camino menos largo, obtendremos la segunda solución por este método.

Recordemos que

$$P(x)={\displaystyle \frac{x-1}{x}}$$

y que la primer solución es

$$y_1(x)=x{e}^{-x}$$

Notemos que

$$-\int P(x)dx=-\int {\displaystyle \frac{x-1}{x}}dx=\int ({\displaystyle \frac{1}{x}}-1)dx=\ln (x)-x$$

Sustituimos en ($\text{19}$).

$$\begin{array}{rl}{y}_2(x)& =x{e}^{-x}\int {\displaystyle \frac{e^{\ln (x)-x}}{(x{e}^{-x}{)}^2}}dx\\ & =x{e}^{-x}\int {\displaystyle \frac{x{e}^{-x}}{x^2{e}^{-2x}}}dx\\ & =x{e}^{-x}\int {\displaystyle \frac{e^x}{x}}dx\end{array}$$

La integral resultante es conocida como integral exponencial $Ei(x)$ y corresponde a una función especial definida en el plano complejo. Para nuestro caso es conveniente escribir a la exponencial como serie e integrar término a término.

$$\begin{array}{rl}{y}_2(x)& =x{e}^{-x}\int {\displaystyle \frac{1}{x}}(1+x+{\displaystyle \frac{x^2}{2}}+{\displaystyle \frac{x^3}{3!}}+{\displaystyle \frac{x^4}{4!}}+\cdots )dx\\ & =x{e}^{-x}\int ({\displaystyle \frac{1}{x}}+1+{\displaystyle \frac{x}{2}}+{\displaystyle \frac{x^2}{3!}}+{\displaystyle \frac{x^3}{4!}}+\cdots )dx\\ & =x{e}^{-x}[\ln (x)+x+{\displaystyle \frac{x^2}{2(2!)}}+{\displaystyle \frac{x^3}{3(3!)}}+{\displaystyle \frac{x^4}{4(4!)}}+\cdots ]\\ & =x{e}^{-x}\ln (x)+x{e}^{-x}\sum _{n=1}^{\mathrm{\infty }}{\displaystyle \frac{x^n}{n(n!)}}\end{array}$$

Vemos que

$$\begin{array}{rl}x{e}^{-x}\sum _{n=1}^{\mathrm{\infty }}{\displaystyle \frac{x^n}{n(n!)}}& =x(1-x+{\displaystyle \frac{x^2}{2!}}-{\displaystyle \frac{x^3}{3!}}+{\displaystyle \frac{x^4}{4!}}-\cdots )(x+{\displaystyle \frac{x^2}{2(2!)}}+{\displaystyle \frac{x^3}{3(3!)}}+{\displaystyle \frac{x^4}{4(4!)}}+\cdots )\\ & =(x-x^2+{\displaystyle \frac{x^3}{2}}-{\displaystyle \frac{x^4}{6}}+{\displaystyle \frac{x^5}{24}}-\cdots )(x+{\displaystyle \frac{x^2}{4}}+{\displaystyle \frac{x^3}{18}}+{\displaystyle \frac{x^4}{96}}+\cdots )\\ & =x^2+({\displaystyle \frac{x^3}{4}}-x^3)+({\displaystyle \frac{x^4}{18}}-{\displaystyle \frac{x^4}{4}}+{\displaystyle \frac{x^4}{2}})+({\displaystyle \frac{x^5}{96}}-{\displaystyle \frac{x^5}{18}}+{\displaystyle \frac{x^5}{8}}-{\displaystyle \frac{x^5}{6}})+\cdots \\ & =x^2-{\displaystyle \frac{3}{4}}{x}^3+{\displaystyle \frac{11}{36}}{x}^4-25x^5+\cdots \end{array}$$

Entonces la segunda solución es

$$y_2(x)=x{e}^{-x}\ln (x)+x^2-{\displaystyle \frac{3}{4}}{x}^3+{\displaystyle \frac{11}{36}}{x}^4-25x^5+\cdots$$

Por lo tanto, la solución general de la ecuación diferencial es

$$y(x)=C_{1}x{e}^{-x}+C_2(x{e}^{-x}\ln (x)+x^2-{\displaystyle \frac{3}{4}}{x}^3+{\displaystyle \frac{11}{36}}{x}^4-25x^5+\cdots )$$

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Solución cuando la diferencia de las raíces indiciales es un número entero positivo

Ejemplo: Resolver la ecuación diferencial

$$x{\displaystyle \frac{d^{2}y}{d{x}^2}}-{\displaystyle \frac{dy}{dx}}+4x^{3}y=0$$

con respecto al punto singular $x_0=0$.

Solución: Dividimos toda la ecuación por $x$ para obtener la forma estándar.

$${\displaystyle \frac{d^{2}y}{d{x}^2}}-{\displaystyle \frac{1}{x}}{\displaystyle \frac{dy}{dx}}+4x^{2}y=0$$

Identificamos que

$$P(x)=-{\displaystyle \frac{1}{x}}yQ(x)=4x^2$$

Es claro que $x=0$ es un punto ordinario de $Q(x)$, sin embargo es un punto singular regular de $P(x)$, pues

$$\underset{x\to 0}{lim}xP(x)=\underset{x\to 0}{lim}-1=-1$$

Sustituimos las funciones correspondientes en la ecuación diferencial.

$$x[\sum _{n=0}^{\mathrm{\infty }}(n+r)(n+r-1)c_n{x}^{n+r-2}]-[\sum _{n=0}^{\mathrm{\infty }}(n+r)c_n{x}^{n+r-1}]+4x^3\left[\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r}\right]=0$$

$$\sum _{n=0}^{\mathrm{\infty }}(n+r)(n+r-1)c_n{x}^{n+r-1}-\sum _{n=0}^{\mathrm{\infty }}(n+r)c_n{x}^{n+r-1}+4\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r+3}=0$$

Con el propósito de que en la tercer serie $x$ tenga la misma potencia que las dos primeras, hacemos $k=n+4$ y en el resto $k=n$.

$$\sum _{k=0}^{\mathrm{\infty }}(k+r)(k+r-1)c_k{x}^{k+r-1}-\sum _{k=0}^{\mathrm{\infty }}(k+r)c_k{x}^{k+r-1}+4\sum _{k=4}^{\mathrm{\infty }}{c}_{k-4}{x}^{k+r-1}=0$$

Para $k=0$, se tiene

$$\begin{array}{rl}r(r-1)c_0{x}^{r-1}-r{c}_0{x}^{r-1}& =0\\ c_0{x}^{r-1}[r(r-1)-r]& =0\end{array}$$

de donde se obtiene la ecuación indicial

$$r^2-2r=0$$

cuyas raíces son

$$r_1=2y{r}_2=0$$

Como

$$r_1-r_2=2$$

Es decir, la diferencia es un número entero, entonces estamos en condiciones del caso 3 y por tanto las soluciones son de la forma ($\text{16}$) y ($\text{17}$).

$$y_1(x)=\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+2},c_0\ne 0$$

y

$$y_2(x)=C\ln (x)\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+2}+\sum _{n=0}^{\mathrm{\infty }}{\hat{c}}_n{x}^n,{\hat{c}}_0\ne 0$$

Recordemos que $C$ puede ser cero.

Necesitamos que todas las series comiencen en $k=4$ para poder obtener la relación de recurrencia. Extraemos los términos para $k=1$, $k=2$ y $k=3$ y cada suma correspondiente la igualamos a cero.

$k=1$.

$$\begin{array}{rl}(1+r)(r)c_1{x}^r-(1+r)c_1{x}^r& =0\\ x^r[(1+r)(r)-(1+r)]c_1& =0\end{array}$$

Debido a que

$$(1+r)(r)-(1+r)\ne 0$$

de acuerdo a los valores de las raíces indiciales, entonces necesariamente $c_1=0$.

$k=2$.

$$\begin{array}{rl}(2+r)(1+r)c_2{x}^{r+1}-(2+r)c_2{x}^{r+1}& =0\\ x^{r+1}[(2+r)(1+r)-(2+r)]c_2& =0\end{array}$$

de donde necesariamente $c_2=0$.

$k=3$.

$$\begin{array}{rl}(3+r)(2+r)c_3{x}^{r+2}-(3+r)c_3{x}^{r+2}& =0\\ x^{r+2}[(3+r)(2+r)-(3+r)]c_3& =0\end{array}$$

Igualmente obtenemos que $c_3=0$.

Ahora tenemos la ecuación

$$\sum _{k=4}^{\mathrm{\infty }}(k+r)(k+r-1)c_k{x}^{k+r-1}-\sum _{k=4}^{\mathrm{\infty }}(k+r)c_k{x}^{k+r-1}+4\sum _{k=4}^{\mathrm{\infty }}{c}_{k-4}{x}^{k+r-1}=0$$

La reescribimos en una sola serie.

$$\sum _{k=4}^{\mathrm{\infty }}[(k+r)(k+r-1)c_k-(k+r)c_k+4c_{k-4}]x^{k+r-1}=0$$

De donde,

$$c_k[(k+r)(k+r-1)-(k+r)]+4c_{k-4}=0$$

Despejando $c_k$ obtenemos la relación de recurrencia.

$$c_k={\displaystyle \frac{4c_{k-4}}{(k+r)-(k+r)(k+r-1)}},k=4,5,6,\cdots$$

Para el caso en el que $r=2$ la relación de recurrencia es

$$c_k=-{\displaystyle \frac{4c_{k-4}}{k(k+2)}},k=4,5,6,\cdots$$

Determinemos los coeficientes.

$k=4$.

$$c_4=-{\displaystyle \frac{4c_0}{4(4+2)}}=-{\displaystyle \frac{4c_0}{24}}=-{\displaystyle \frac{c_0}{6}}$$

Para $k=5$, $k=6$ y $k=7$ obtendremos que $c_5=0$, $c_6=0$ y $c_7=0$ respectivamente.

$k=8$.

$$c_8=-{\displaystyle \frac{4c_4}{8(8+2)}}=-{\displaystyle \frac{4c_4}{80}}=-{\displaystyle \frac{c_4}{20}}={\displaystyle \frac{c_0}{120}}$$

De la misma manera $c_9=c_{10}=c_{11}=0$.

$k=12$.

$$c_{12}=-{\displaystyle \frac{4c_8}{12(12+2)}}=-{\displaystyle \frac{4c_8}{168}}=-{\displaystyle \frac{c_8}{42}}=-{\displaystyle \frac{c_0}{5040}}$$

Etcétera, entonces

$$\begin{array}{rl}{y}_1(x)& =x^2(c_0-{\displaystyle \frac{c_0}{6}}{x}^4+{\displaystyle \frac{c_0}{120}}{x}^8-{\displaystyle \frac{c_0}{5040}}{x}^{12}+\cdots )\\ & =c_0(x^2-{\displaystyle \frac{x^6}{3!}}+{\displaystyle \frac{x^{10}}{5!}}-{\displaystyle \frac{x^{14}}{7!}}+\cdots )\end{array}$$

Sabemos que

$$\sin (x)=x-{\displaystyle \frac{x^3}{3!}}+{\displaystyle \frac{x^5}{5!}}-{\displaystyle \frac{x^7}{7!}}+\cdots =\sum _{n=0}^{\mathrm{\infty }}{\displaystyle \frac{(-1{)}^n{x}^{2n+1}}{(2n+1)!}}$$

Entonces la primer solución es

$$y_1(x)=c_0\sin (x^2)$$

Para obtener la segunda solución $y_2$ podemos probar con la relación de recurrencia que obtuvimos o por alguno de los métodos que ya conocemos.

Consideremos la relación de recurrencia obtenida

$$c_k={\displaystyle \frac{4c_{k-4}}{(k+r)-(k+r)(k+r-1)}},k=4,5,6,\cdots$$

Usemos la notación ${\hat{c}}_k$ y el valor de la segunda raíz indicial $r=0$, en este caso la relación de recurrencia es

$${\hat{c}}_k=-{\displaystyle \frac{4c_{k-4}}{k(k-2)}},k=4,5,6,\cdots$$

Los mismos coeficientes que fueron cero en el caso anterior serán cero en este caso, así que sólo consideraremos que $k=4,8,12,\cdots$. Determinemos los coeficientes.

$k=4$.

$${\hat{c}}_4=-{\displaystyle \frac{4{\hat{c}}_0}{4(4-2)}}=-{\displaystyle \frac{4{\hat{c}}_0}{8}}=-{\displaystyle \frac{{\hat{c}}_0}{2}}$$

$k=8$.

$${\hat{c}}_8=-{\displaystyle \frac{4{\hat{c}}_4}{8(8-2)}}=-{\displaystyle \frac{4{\hat{c}}_4}{48}}=-{\displaystyle \frac{{\hat{c}}_4}{12}}={\displaystyle \frac{{\hat{c}}_0}{24}}$$

$k=12$.

$${\hat{c}}_{12}=-{\displaystyle \frac{4{\hat{c}}_8}{12(12-2)}}=-{\displaystyle \frac{4{\hat{c}}_8}{120}}=-{\displaystyle \frac{{\hat{c}}_8}{30}}=-{\displaystyle \frac{{\hat{c}}_0}{720}}$$

Etcétera, entonces

$$\begin{array}{rl}y& ={\hat{c}}_0-{\displaystyle \frac{{\hat{c}}_0}{2}}{x}^4+{\displaystyle \frac{{\hat{c}}_0}{24}}{x}^8-{\displaystyle \frac{{\hat{c}}_0}{720}}{x}^{12}+\cdots \\ & ={\hat{c}}_0(1-{\displaystyle \frac{x^4}{2!}}+{\displaystyle \frac{x^8}{4!}}-{\displaystyle \frac{x^{12}}{6!}}+\cdots )\end{array}$$

Sabemos que

$$\cos (x)=1-{\displaystyle \frac{x^2}{2!}}+{\displaystyle \frac{x^4}{4!}}-{\displaystyle \frac{x^6}{6!}}+\cdots =\sum _{n=0}^{\mathrm{\infty }}{\displaystyle \frac{(-1{)}^n{x}^{2n}}{(2n)!}}$$

Entonces la segunda solución es

$$y_2(x)={\hat{c}}_0\cos (x^2)$$

Vemos que el método no nos indica la existencia de la función $\ln (x)$ y nosotros esperamos una solución de la forma

$$y_2(x)=C\ln (x)\sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+2}+\sum _{n=0}^{\mathrm{\infty }}{\hat{c}}_n{x}^n,{\hat{c}}_0\ne 0$$

Entonces podemos concluir que $C=0$, así

$y_2(x)={\hat{c}}_0\cos (x^2)$

Veamos que se obtiene usando la fórmula ($\text{19}$). Recordemos que

$$P(x)=-{\displaystyle \frac{1}{x}}$$

y consideremos que $c_0=1$, tal que

$$y_1(x)=\sin (x^2)$$

Vemos que

$$-\int P(x)dx=\int {\displaystyle \frac{dx}{x}}=\ln (x)$$

Sustituyamos en ($\text{19}$).

$$y_2(x)=\sin (x^2)\int {\displaystyle \frac{e^{\ln (x)}}{(\sin (x^2){)}^2}}dx=\sin (x^2)\int {\displaystyle \frac{x}{(\sin (x^2){)}^2}}dx$$

Resolviendo la integral se obtiene que

$$\int {\displaystyle \frac{x}{(\sin (x^2){)}^2}}dx=-{\displaystyle \frac{1}{2}}\cot (x^2)$$

Entonces,

$$y_2(x)=-{\displaystyle \frac{1}{2}}\sin (x^2)\left({\displaystyle \frac{\cos (x^2)}{\sin (x^2)}}\right)=-{\displaystyle \frac{1}{2}}\cos (x^2)={\hat{c}}_0\cos (x^2)$$

Este método nos indica que efectivamente $C=0$. Si $C_1=c_0$ y $C_2={\hat{c}}_0$, entonces la solución general de la ecuación diferencial es

$$y(x)=C_1\sin (x^2)+C_2\cos (x^2)$$

$\mathrm{<img\; draggable="false"\; role="img"\; class="emoji"\; alt="◻"\; src="https://s.w.org/images/core/emoji/14.0.0/svg/25fb.svg">}$

Hemos concluido con esta entrada.

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

1. Usar el método de Frobenius para obtener la solución general de las siguientes ecuaciones diferenciales en el punto singular $x_0=0$. Verificar que dicho punto es singular.

• $2x{\displaystyle \frac{d^{2}y}{d{x}^2}}+(x+1){\displaystyle \frac{dy}{dx}}+3y=0$.

• $x^2{\displaystyle \frac{d^{2}y}{d{x}^2}}-{\displaystyle \frac{1}{6}}x{\displaystyle \frac{dy}{dx}}+{\displaystyle \frac{1}{3}}y=0$.

• $x^2{\displaystyle \frac{d^{2}y}{d{x}^2}}+6x{\displaystyle \frac{dy}{dx}}+(6-x^2)y=0$.

• $2x^2{\displaystyle \frac{d^{2}y}{d{x}^2}}-x^2{\displaystyle \frac{dy}{dx}}-(x+4)y=0$.

• $x{\displaystyle \frac{d^{2}y}{d{x}^2}}+(x-1){\displaystyle \frac{dy}{dx}}+({\displaystyle \frac{1}{x}}-1)y=0$.

• $(x^2-x){\displaystyle \frac{d^{2}y}{d{x}^2}}+(3x-1){\displaystyle \frac{dy}{dx}}+y=0$.

Más adelante…

Ahora que sabemos resolver ecuaciones diferenciales lineales de segundo orden con coeficientes variables con respecto a puntos ordinarios y puntos singulares, en las siguientes entradas resolveremos algunas ecuaciones diferenciales especiales cuya utilidad es de suma importancia en otras áreas del conocimiento como la física, biología e ingeniería entre otras.

Entradas relacionadas

• Página principal del curso: Ecuaciones Diferenciales I
• Entrada anterior del curso: Ecuaciones lineales de segundo orden con coeficientes variables – Soluciones en series de potencias respecto a puntos ordinarios
• Siguiente entrada del curso: Ecuaciones del Hermite, Laguerre y Legendre

• Video relacionado al tema: Ecuaciones lineales de segundo orden con coeficientes variables. Solución por series cerca de un punto singular regular

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»