La naturaleza está escrita en lenguaje matemático.
– Galileo Galilei

Introducción

En las dos últimas entradas hemos desarrollado métodos de resolución de ecuaciones diferenciales lineales de segundo orden con coeficientes variables. El primer caso fue cuando $x_0=0$ es un punto ordinario y en el segundo caso cuando $x_0=0$ es un punto singular regular. En esta y la siguiente entrada aplicaremos estos métodos para resolver algunas ecuaciones diferenciales especiales, tan especiales que cada una de ellas tiene su propio nombre y son de bastante utilidad en otras áreas del conocimiento como la física e ingeniería.

A continuación presentamos las ecuaciones diferenciales que resolveremos:

• Ecuación de Hermite.

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

• Ecuación de Laguerre.

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

• Ecuación de Legendre.

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

• Ecuación de Bessel.

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

• Ecuación de Chebyshev.

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

• Ecuación Hipergeométrica de Gauss.

$$x(1-x){\displaystyle \frac{d^{2}y}{d{x}^2}}+[\gamma -(\alpha +\beta +1)x]{\displaystyle \frac{dy}{dx}}-\alpha \beta y=0$$

• Ecuación de Airy.

$${\displaystyle \frac{d^{2}y}{d{x}^2}}-xy=0$$

Algunos ejemplos en los que aparecen este tipo de ecuaciones son en el estudio de potenciales en campos conservativos y no conservativos, esfuerzos de torsión, distribución de temperaturas, propagación de calor, vibraciones de cuerdas y membranas, propagación de ondas sonoras, luminosas, de radio entre muchas otras aplicaciones.

Es importante aclarar que todas estas ecuaciones, y las soluciones de cada una, tienen importantes propiedades matemáticas que no serán expuestas en este curso, nuestro propósito es el de sólo dar con la solución aplicando los métodos ya mencionados. Sin embargo, estos resultados seguramente serán de bastante utilidad más adelante cuando en semestres posteriores se estudien con mayor detalle. Por supuesto, si en estos momentos se desea conocer más acerca de estas ecuaciones diferenciales se puede consultar bibliografía existente para cada una de ellas.

Comencemos con la ecuación de Hermite.

Ecuación de Hermite

La ecuación de Hermite es

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

Con $x\in \mathbb{R}$ y $\lambda$ una constante.

Esta ecuación diferencial es llamada así en honor al matemático francés Charles Hermite (1822 – 1901), quien realizó investigaciones sobre teoría de números, formas cuadráticas, teoría de invariantes, polinomios ortogonales y funciones elípticas entre otros. Varias entidades matemáticas se llaman hermitianas en su honor.

La ecuación de Hermite se encuentra en forma estándar lo que nos permite notar que el punto $x_0=0$ es un punto ordinario, esto nos indica que su solución es de la forma

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

Cuyas derivadas son

$$\begin{array}{}\text{(3)}& {\displaystyle \frac{dy}{dx}}=\sum _{n=1}^{\mathrm{\infty }}n{c}_n{x}^{n-1}y{\displaystyle \frac{d^{2}y}{d{x}^2}}=\sum _{n=2}^{\mathrm{\infty }}n(n-1)c_n{x}^{n-2}\end{array}$$

Sustituyamos en la ecuación de Hermite.

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

Introducimos la $x$ en la segunda serie.

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

En la primer serie hacemos la sustitución $k=n-2$ y en las otras dos hacemos $k=n$.

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

Extraemos el primer término de la primera y última serie para que todas comiencen en $k=1$.

$$\begin{array}{}\text{(4)}& 2c_2+\lambda c_0=0\end{array}$$

de donde,

$$c_2=-{\displaystyle \frac{\lambda }{2}}{c}_0$$

Ahora tenemos la ecuación

$$\begin{array}{}\text{(5)}& \sum _{k=1}^{\mathrm{\infty }}(k+2)(k+1)c_{k+2}{x}^k-2\sum _{k=1}^{\mathrm{\infty }}k{c}_k{x}^k+\lambda \sum _{k=1}^{\mathrm{\infty }}{c}_k{x}^k=0\end{array}$$

Ahora que todas las series comienzan con el mismo índice y tienen la misma potencia en la variable $x$, podemos juntar todo en una sola serie.

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

De donde necesariamente debe de ocurrir que

$$\begin{array}{}\text{(6)}& (k+2)(k+1)c_{k+2}-(2k-\lambda )c_k=0\end{array}$$

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

$$\begin{array}{}\text{(7)}& c_{k+2}={\displaystyle \frac{2k-\lambda }{(k+2)(k+1)}}{c}_k,k=0,1,2,3\cdots \end{array}$$

Determinemos los coeficientes. Ya vimos que para $k=0$,

$c_2=-{\displaystyle \frac{\lambda }{2!}}{c}_0$

$k=1$.

$$c_3={\displaystyle \frac{2(1)-\lambda }{(3)(2)}}{c}_1={\displaystyle \frac{2-\lambda }{3!}}{c}_1$$

$k=2$.

$$c_4={\displaystyle \frac{2(2)-\lambda }{(4)(3)}}{c}_2={\displaystyle \frac{4-\lambda }{(4)(3)}}(-{\displaystyle \frac{\lambda }{2}}{c}_0)=-{\displaystyle \frac{\lambda (4-\lambda )}{4!}}{c}_0$$

$k=3$.

$$c_5={\displaystyle \frac{2(3)-\lambda }{(5)(4)}}{c}_3={\displaystyle \frac{6-\lambda }{(5)(4)}}\left({\displaystyle \frac{2-\lambda }{(3)(2)}}{c}_1\right)={\displaystyle \frac{(6-\lambda )(2-\lambda )}{5!}}{c}_1$$

$k=4$.

$$c_6={\displaystyle \frac{2(4)-\lambda }{(6)(5)}}{c}_4={\displaystyle \frac{8-\lambda }{(6)(5)}}(-{\displaystyle \frac{\lambda (4-\lambda )}{4!}}{c}_0)=-{\displaystyle \frac{\lambda (4-\lambda )(8-\lambda )}{6!}}{c}_0$$

$k=5$.

$$c_7={\displaystyle \frac{2(5)-\lambda }{(7)(6)}}{c}_5={\displaystyle \frac{10-\lambda }{(7)(6)}}\left({\displaystyle \frac{(6-\lambda )(2-\lambda )}{5!}}{c}_1\right)={\displaystyle \frac{(2-\lambda )(6-\lambda )(10-\lambda )}{7!}}{c}_1$$

Etcétera, si tomamos como factores comunes a $C_1=c_0$ y $C_2=c_1$, entonces podemos escribir a la solución general de la ecuación de Hermite como

$$\begin{array}{rl}y(x)& =C_1[1-{\displaystyle \frac{\lambda }{2!}}{x}^2-{\displaystyle \frac{\lambda (4-\lambda )}{4!}}{x}^4-{\displaystyle \frac{\lambda (4-\lambda )(8-\lambda )}{6!}}{x}^6–\cdots ]\\ \text{(8)}& & +C_2[x+{\displaystyle \frac{(2-\lambda )}{3!}}{x}^3+{\displaystyle \frac{(2-\lambda )(6-\lambda )}{5!}}{x}^5+{\displaystyle \frac{(2-\lambda )(6-\lambda )(10-\lambda )}{7!}}+\cdots ]\end{array}$$

Un caso interesante ocurre cuando el parámetro $\lambda$ es positivo y es par, es decir de la forma $\lambda =2k$, en este caso la relación de recurrencia muestra que

$$c_{k+2}=c_{k+4}=\cdots =0$$

Notemos que si $\lambda =2k$ y además $k$ es par y se toma $C_2=0$, entonces la solución se reduce a un polinomio de grado $k$, lo mismo ocurre si $k$ es impar y se toma $C_1=0$, la solución se reduce a otro polinomio de grado $k$.

Con una adecuada elección de $C_1$ y $C_2$, de tal manera que el coeficiente de $x^k$ sea $2^k$, resultan los denominados polinomios de Hermite.

$$\begin{array}{rl}{H}_0(x)& =1\\ H_1(x)& =2x\\ H_2(x)& =4x^2-2\\ H_3(x)& =8x^3-12x\\ H_4(x)& =16x^4-48x^2+12\\ H_5(x)& =32x^5-160x^3+120x\\ ⋮\end{array}$$

Cada polinomio de Hermite es solución particular de la ecuación de Hermite con $\lambda =0,2,4,6\cdots$, respectivamente. En general, el $n$-ésimo polinomio de Hermite es solución particular de la ecuación de Hermite con $\lambda =2n$.

Los polinomios de Hermite aparecen en la resolución del problema del oscilador armónico unidimensional en Mecánica Cuántica.

Pasemos a resolver la ecuación de Laguerre.

Ecuación de Laguerre

La ecuación de Laguerre es

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

Con $\lambda$ una constante.

Los polinomios de Laguerre son una familia de polinomios ortogonales que surgen de examinar las soluciones de la ecuación diferencial ($\text{9}$). Edmond Nicolás Laguerre (1834 – 1886) fue un matemático francés conocido principalmente por la introducción de los polinomios que llevan su nombre.

Resolvamos la ecuación, para ello dividimos todo por $x$ para obtener la forma estándar.

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

Identificamos que

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

Es claro que ambas funciones no están definidas en $x=0$, de manera que este punto es un punto singular. Si definimos las funciones

$$p(x)=xP(x)yq(x)=x^{2}Q(x)$$

obtenemos que

$$p(x)=1-xyq(x)=\lambda x$$

Si calculamos los límites se obtiene lo siguiente.

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

Los límites existen, esto nos indica que el punto $x_0=0$ es un punto singular regular. La solución para este caso es de la forma

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

Las derivadas son

$$\begin{array}{}\text{(12)}& {\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}\end{array}$$

Sustituyamos en la ecuación de Laguerre.

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

Expandiendo y simplificando se tiene

$$\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}-\sum _{n=0}^{\mathrm{\infty }}(n+r)c_n{x}^{n+r}+\lambda \sum _{n=0}^{\mathrm{\infty }}{c}_n{x}^{n+r}=0$$

En las dos primeras series hacemos $k=n$ y en las dos últimas series hacemos $n=k-1$.

$$\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}-\sum _{k=1}^{\mathrm{\infty }}(k-1+r)c_{k-1}{x}^{k+r-1}+\lambda \sum _{k=1}^{\mathrm{\infty }}{c}_{k-1}{x}^{k+r-1}=0$$

Extraemos los términos para $k=0$ y así hacer que todas las series comiencen en $k=1$.

$$\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\\ r(r-1)+r& =0\end{array}$$

La ecuación indicial es

$$\begin{array}{}\text{(13)}& r^2=0\end{array}$$

de donde $r_1=r_2=r=0$. Como las raíces indiciales son iguales, la forma de las soluciones es

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

y

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

Continuemos con la ecuación que teníamos.

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

Ahora que todas inician en $k=1$ y tienen la misma potencia podemos agruparlas en una sola serie.

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

De donde es necesario que

$$\begin{array}{rl}(k+r)(k+r-1)c_k+(k+r)c_k-(k-1+r)c_{k-1}+\lambda c_{k-1}& =0\\ c_k[(k+r)(k+r-1)+(k+r)]+c_{k-1}[\lambda -(k-1+r)]& =0\end{array}$$

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

$$\begin{array}{}\text{(16)}& c_k={\displaystyle \frac{(k-1+r)-\lambda }{(k+r)(k+r-1)+(k+r)}}{c}_{k-1}\end{array}$$

De tarea moral muestra que la relación de recurrencia se puede reescribir como

$$\begin{array}{}\text{(17)}& c_k={\displaystyle \frac{(k+r)-(\lambda +1)}{(k+r{)}^2}}{c}_{k-1}\end{array}$$

Sabemos que la raíz indicial es $r=0$, entonces la relación de recurrencia se reduce a

$$\begin{array}{}\text{(18)}& c_k={\displaystyle \frac{k-(\lambda +1)}{k^2}}{c}_{k-1},k=1,2,3,\cdots \end{array}$$

Determinemos los coeficientes.

$k=1$.

$$c_1={\displaystyle \frac{1-(\lambda +1)}{1^2}}{c}_0=-\lambda c_0$$

$k=2$.

$$c_2={\displaystyle \frac{2-(\lambda +1)}{2^2}}{c}_1={\displaystyle \frac{1-\lambda }{4}}{c}_1={\displaystyle \frac{\lambda (\lambda -1)}{4}}{c}_0$$

$k=3$.

$$c_3={\displaystyle \frac{3-(\lambda +1)}{3^2}}{c}_2={\displaystyle \frac{2-\lambda }{9}}{c}_2=-{\displaystyle \frac{\lambda (\lambda -1)(\lambda -2)}{36}}{c}_0$$

Continuando es posible encontrar el patrón y establecer que

$$\begin{array}{}\text{(19)}& c_k=(-1{)}^k{\displaystyle \frac{\lambda (\lambda -1)(\lambda -2)\cdots (\lambda -k+1)}{(k!{)}^2}}{c}_0\end{array}$$

De tarea moral demuestra por inducción el resultado anterior.

Entonces la solución de la ecuación de Laguerre es

$$\begin{array}{}\text{(20)}& y(x)& =c_0(1-{\displaystyle \frac{\lambda }{(1!{)}^2}}x+{\displaystyle \frac{\lambda (\lambda -1)}{(2!{)}^2}}{x}^2-{\displaystyle \frac{\lambda (\lambda -1)(\lambda -2)}{(3!{)}^2}}{x}^3+\cdots +(-1{)}^k{\displaystyle \frac{\lambda (\lambda -1)(\lambda -2)\cdots (\lambda -k+1)}{(k!{)}^2}}{x}^k+\cdots )\end{array}$$

Recordemos que el método de Frobenius nos dice que existe una segunda solución de la forma

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

Obtener la segunda solución resulta ser una tarea muy complicada debido a la enorme cantidad de cálculos que se deben realizar, en el video correspondiente se hace notar esta dificultad, sin embargo la solución obtenida suele ser suficiente para trabajar y es la que se utiliza en las aplicaciones que aparecen principalmente en Física.

Observemos que si $\lambda \in {\mathbb{Z}}^{+}$, entonces la serie solución se hace finita, ya que cada coeficiente de la serie contiene un término $(\lambda -m)$ con $m\in {\mathbb{Z}}^{+}$ que se repite cada vez que aparece por primera vez, por ejemplo el término $(\lambda -2)$ aparece por primera vez en el coeficiente de $x^3$ y a partir de ahí aparece en el resto de coeficientes de la serie, de manera que si $\lambda =2$, entonces todos los coeficientes que contengan el término $(\lambda -2)$ se anularán y sólo nos quedará un polinomio de grado $n=2$. Estos polinomios resultantes son los llamados polinomios de Laguerre.

Para $\lambda =0,1,2,3,\cdots$ y con el valor adecuado de $c_0$ se obtienen los siguientes polinomios de Laguerre.

$$\begin{array}{rl}{L}_0(x)& =1\\ L_1(x)& =1-x\\ L_2(x)& =1-2x+{\displaystyle \frac{1}{2}}{x}^2\\ L_3(x)& =1-3x+{\displaystyle \frac{3}{2}}{x}^2-{\displaystyle \frac{1}{6}}{x}^3\\ ⋮\end{array}$$

En general, el $n$-ésimo polinomio de Laguerre será solución particular de la ecuación de Laguerre cuando $\lambda =n$.

Finalicemos esta entrada con la ecuación de Legendre.

Ecuación de Legendre

La ecuación de Legendre es

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

Con $\lambda$ una constante.

Esta ecuación lleva este nombre en honor al matemático francés Adrien – Marie Legendre (1752 – 1833). Legendre hizo importantes contribuciones a la estadística, la teoría de números, el álgebra abstracta y el análisis matemático.

Resolvamos la ecuación, dividimos todo por el coeficiente de la segunda derivada de $y$.

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

Identificamos que

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

Vemos que ambas funciones no están definidas en $x=1$ ni $x=-1$, pero si en en el punto $x_0=0$, de manera que dicho punto es un punto ordinario, entonces la forma de la solución de la ecuación de Legendre es

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

Con primera y segunda derivada dadas como

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

Sustituyamos en la ecuación de Legendre.

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

Expandiendo y simplificando, se tiene

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

En la primer serie hacemos $k=n-2$ y en el resto $k=n$.

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

Extraemos los términos para $k=0$ y $k=1$ y con ello lograr que todas las series comiencen en $k=2$.

Por un lado, para $k=0$,

$$2c_2+\lambda (\lambda +1)c_0=0$$

De donde

$$c_2=-{\displaystyle \frac{\lambda (\lambda +1)}{2}}{c}_0$$

Por otro lado, para $k=1$,

$$\begin{array}{rl}3(2)c_{3}x-2c_{1}x+\lambda (\lambda +1)c_{1}x& =0\\ [6c_3-2c_1+\lambda (\lambda +1)c_1]x& =0\\ 6c_3-2c_1+\lambda (\lambda +1)c_1& =0\end{array}$$

De donde

$$c_3={\displaystyle \frac{2-\lambda (\lambda +1)}{6}}{c}_1$$

Veremos más adelante que es conveniente escribir este resultado como

$$c_3=-{\displaystyle \frac{(\lambda -1)(\lambda +2)}{6}}{c}_1$$

Ahora tenemos la ecuación

$$\begin{array}{}\text{(23)}& \sum _{k=2}^{\mathrm{\infty }}(k+2)(k+1)c_{k+2}{x}^k-\sum _{k=2}^{\mathrm{\infty }}k(k-1)c_k{x}^k-2\sum _{k=2}^{\mathrm{\infty }}k{c}_k{x}^k+\lambda (\lambda +1)\sum _{k=2}^{\mathrm{\infty }}{c}_k{x}^k=0\end{array}$$

Juntemos todo en una sola serie.

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

De donde es necesario que

$$\begin{array}{}\text{(24)}& (k+2)(k+1)c_{k+2}-[k(k-1)+2k-\lambda (\lambda +1)]c_k=0\end{array}$$

Despejando a $c_{k+2}$ obtenemos la relación de recurrencia.

$$\begin{array}{}\text{(25)}& c_{k+2}={\displaystyle \frac{k(k-1)+2k-\lambda (\lambda +1)}{(k+2)(k+1)}}{c}_k,k=0,1,2,\cdots \end{array}$$

Es conveniente reescribir a la ecuación de recurrencia de la siguiente manera.

$$\begin{array}{}\text{(26)}& c_{k+2}=-{\displaystyle \frac{(\lambda -k)(\lambda +k+1)}{(k+2)(k+1)}}{c}_k,k=0,1,2,\cdots \end{array}$$

Determinemos los coeficientes. Ya vimos que para $k=0$,

$$c_2=-{\displaystyle \frac{\lambda (\lambda +1)}{2!}}{c}_0$$

y para $k=1$,

$$c_3=-{\displaystyle \frac{(\lambda -1)(\lambda +2)}{3!}}{c}_1$$

$k=2$.

$$c_4=-{\displaystyle \frac{(\lambda -2)(\lambda +3)}{(4)(3)}}{c}_2={\displaystyle \frac{(\lambda -2)\lambda (\lambda +1)(\lambda +3)}{4!}}{c}_0$$

$k=3$.

$$c_5=-{\displaystyle \frac{(\lambda -3)(\lambda +4)}{(5)(4)}}{c}_3={\displaystyle \frac{(\lambda -3)(\lambda -1)(\lambda +2)(\lambda +4)}{5!}}{c}_1$$

$k=4$.

$$c_6=-{\displaystyle \frac{(\lambda -4)(\lambda +5)}{(6)(5)}}{c}_4=-{\displaystyle \frac{(\lambda -4)(\lambda -2)\lambda (\lambda +1)(\lambda +3)(\lambda +5)}{6!}}{c}_0$$

$k=5$.

$$c_7=-{\displaystyle \frac{(\lambda -5)(\lambda +6)}{(7)(6)}}{c}_5=-{\displaystyle \frac{(\lambda -5)(\lambda -3)(\lambda -1)(\lambda +2)(\lambda +4)(\lambda +6)}{7!}}{c}_1$$

Etcétera, si tomamos como factores comunes a $C_1=c_0$ y $C_2=c_1$, entonces podemos escribir a la solución general de la ecuación de Legendre como

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

Donde,

$$\begin{array}{rl}{y}_1(x)& =1-{\displaystyle \frac{\lambda (\lambda +1)}{2!}}{x}^2+{\displaystyle \frac{(\lambda -2)\lambda (\lambda +1)(\lambda +3)}{4!}}{x}^4\\ \text{(28)}& & -{\displaystyle \frac{(\lambda -4)(\lambda -2)\lambda (\lambda +1)(\lambda +3)(\lambda +5)}{6!}}{x}^6+\cdots \end{array}$$

y

$$\begin{array}{rl}{y}_2(x)& =x-{\displaystyle \frac{(\lambda -1)(\lambda +2)}{3!}}{x}^3+{\displaystyle \frac{(\lambda -3)(\lambda -1)(\lambda +2)(\lambda +4)}{5!}}{x}^5\\ \text{(29)}& & -{\displaystyle \frac{(\lambda -5)(\lambda -3)(\lambda -1)(\lambda +2)(\lambda +4)(\lambda +6)}{7!}}{x}^7+\cdots \end{array}$$

Para $\lambda =0,1,2,3,\cdots$ y con el valor adecuado de $C_1$ y de $C_2$ se obtienen los conocidos polinomios de Legendre:

$$\begin{array}{rl}{P}_0(x)& =1\\ P_1(x)& =x\\ P_2(x)& ={\displaystyle \frac{1}{2}}(3x^2-1)\\ P_3(x)& ={\displaystyle \frac{1}{2}}(5x^3-3x)\\ P_4(x)& ={\displaystyle \frac{1}{8}}(35x^4-30x^2+3)\\ P_5(x)& ={\displaystyle \frac{1}{8}}(63x^5-70x^3+15x)\\ ⋮\end{array}$$

En general, el $n$-ésimo polinomio de Legendre será solución particular de la ecuación de Legendre cuando $\lambda =n$.

La ecuación de Legendre aparece con mucha frecuencia en problemas de Física, en particular en electromagnetismo en problemas de valor límite en esferas.

Los polinomios de Legendre aparecen cuando se resuelve la ecuación de Helmholtz (un tipo de ecuación en derivadas parciales) en coordenadas esféricas mediante el método de separación de variables.

Hasta aquí concluimos esta primer entrada sobre la resolución de algunas ecuaciones diferenciales especiales de segundo orden, en la siguiente entrada continuaremos resolviendo el resto de ecuaciones.

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. Los primeros 6 polinomios de Hermite son solución de la ecuación de Hermite para $\lambda =0,2,4,6,8,10$ respectivamente. Determinar el valor de las constantes $C_1$ y $C_2$, tal que se obtengan los primeros 6 polinomios de Hermite.

2. Resolver la siguiente ecuación de Hermite realizando todo el procedimiento del método.

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

3. Los primeros 4 polinomios de Laguerre son solución de la ecuación de Laguerre para $\lambda =0,1,2,3$ respectivamente. Determinar el valor del coeficiente $c_0$, tal que se obtengan los primeros 4 polinomios de Laguerre.

4. Resolver la siguiente ecuación de Laguerre realizando todo el procedimiento del método.

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

5. Los primeros 6 polinomios de Legendre son solución de la ecuación de Legendre para $\lambda =0,1,2,3,4,5$ respectivamente. Determinar el valor correspondiente de $C_1$ y $C_2$, tal que se obtengan los primeros 6 polinomios de Legendre.

6. Los puntos $x_0=1$ y $x_0=-1$ son puntos singulares de la ecuación de Legendre. Usando el método de Frobenius determinar la solución de la ecuación de Legendre con respecto al punto singular $x_0=1$.
   Hint: Usar el cambio de variable $t=x-x_0$ y la regla de la cadena.

Más adelante…

Hemos resuelto 3 de las 7 ecuaciones diferenciales especiales que deseamos resolver, en la siguiente entrada concluiremos con el resto de ecuaciones y así mismo estaremos concluyendo con la unidad 2 del curso.

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 singulares
• Siguiente entrada del curso: Ecuaciones de Bessel, Chebyshev e Hipergeométrica

• Video relacionado al tema: Ecuaciones de Hermite y Laguerre
• Video relacionado al tema: Ecuaciones de Bessel y Legendre

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»