Archivo de la etiqueta: convergencia de series de potencias

Cálculo Diferencial e Integral II: Series de potencia

Por Miguel Ángel Rodríguez García

Introducción

En la sección anterior vimos el criterio de la convergencia absoluta para las series alternantes, en esta sección veremos las series de potencia, que, como bien dice el nombre, son series polinómicas, veamos la siguiente definición.

Series de potencia

Definición. Una serie de potencia es la serie de la siguiente forma:

$$\sum_{n=0}^{\infty}c_{n}x^{n}=c_{0}+c_{1}x+c_{2}x^{2}+…..$$

A la serie anterior, se le dice que es una serie de potencias alrededor de $x=0$, mientras que, la series de potencias alrededor de $x=a$ se le conoce como series de potencias centradas en $a$, y es de la siguiente forma:

$$\sum_{n=0}^{\infty}c_{n}(x-a)^{n}=c_{0}+c_{1}(x-a)+c_{2}(x-a)^{2}+…..$$

Donde $c_{n}$ son coeficientes en ambos casos.

Un ejemplo de estas series son las series geométricas que ya hemos visto, al hacer los n-coeficientes $c_{n}$ igual a 1:

$$\sum_{n=0}^{\infty}x^{n}=1+x+x^{2}+x^{3}…..$$

Veamos el siguiente teorema de convergencia llamado el teorema de Abel para las series de potencias.

Teorema de Abel:

Sea la siguiente serie: $\sum_{n=0}^{\infty}c_{n}(x-a)^{n}$, entonces se cumple una y solo una de las siguientes afirmaciones:

$a)$ La serie converge solo cuando $x=a$.

$b)$ Existe un número positivo $R$ tal que la serie converge $|x-a|<R$ y diverge si $|x-a|>R$.

$c)$ La serie converge para toda $x$.

La demostración de este teorema es extensa, por lo que sería más conveniente analizarla que demostrarla.

Al número $R$ se le llama el radio de convergencia de la serie, notemos que la serie converge en el intervalo $(a-R, \space a+R)$, si $R=0$ tenemos el primer caso $a)$, es decir, el intervalo consta de un solo punto $x=a$, si $R \to \infty$ entonces tenemos el caso $c)$, es decir, el intervalo de convergencia en este caso es $(-\infty, \infty)$, en el intervalo $b)$ se tiene 4 casos posibles:

$$(a-R, \space a+R)$$

$$[a-R, \space a+R]$$

$$(a-R, \space a+R]$$

$$[a-R, \space a+R)$$

Es decir, la serie puede diverger en ambos extremos o solo un extremo, al igual que la convergencia de la serie.

El teorema de Cauchy-Hadamard nos permite conocer la convergencia de la serie de potencias:

Teorema de Cauchy-Hadamard:

Consideremos la serie de potencias $\sum_{n=0}^{\infty}c_{n}(x-a)^{n}$ y consideremos a $A$ como:

$$A=\lim_{n \to \infty}\sqrt[n]{|c_{n}|}$$

Entonces la serie de potencias converge si el radio de convergencia $R$ se define como:

$$R=\frac{1}{A}$$

De este teorema podemos concluir lo siguiente, dependiendo del valor de $A$ podemos decir que si:

  • $$A=0 \Rightarrow R \to \infty$$
  • $$A \to \infty \Rightarrow R=0$$
  • $$0<A<\infty \Rightarrow R=\frac{1}{A}$$

Demostración:

Sin perdida de generalidad podemos suponer que $a=0$. Supongamos que $|x|<R$, entonces:

$$|c_{n}x^{n}|\leq |c_{n}|R^{n} \tag{1}$$

Ahora, como:

$$\lim_{n \to \infty}\sqrt[n]{|c_{n}|}=A= \frac{1}{R}$$

Para casi todos los índices de $n$, ya que:

$$ \sqrt[n]{|c_{n}|} \leq \frac{1}{R} \Rightarrow |c_{n}| \leq R^{-n}$$

Por lo que en $(1)$:

$$|c_{n}x^{n}|\leq |c_{n}|R^{n} \leq R^{n} R^{-n}=1 $$

Lo cual vemos que es una serie absolutamente convergente, por el criterio de la absoluta convergencia:

$$ \sum_{n=0}^{\infty}c_{n}(x-a)^{n} \space converge $$

Para el caso cuando $|x|>R $, de la misma manera anterior obtendremos que:

$$ |c_{n}x^{n}| \geq 1 $$

Lo que significa que no puede convergir a cero, lo que significa que la serie diverge.

$$ \sum_{n=0}^{\infty}c_{n}(x-a)^{n} \space diverge $$

$\square$

El teorema nos dice que podemos usar el criterio de la raíz, también podemos usar el criterio de la razón.

Veamos unos ejemplos.

Ejemplos

  • $$\sum_{n=0}^{\infty}(n!)x^{n}$$

En esta serie notamos que $c_{n}=n!$, entonces calculamos al valor $A$ como sigue:

$$A=\lim_{n \to \infty}\frac{|a_{n}+1|}{|a_{n}|}=\lim_{n \to \infty}\frac{|(n+1)!|}{|n!|}=\lim_{n \to \infty}(n+1) \rightarrow \infty \Rightarrow R=0$$

Por lo que el radio de convergencia es $R=0$ y la serie solo converge cuando $x=0$ según el teorema de Abel.

  • $$\sum_{n=0}^{\infty}\frac{n^{2n+1}}{2^{n^{2}+1}}x^{n}$$

Vemos en este caso que $c_{n}=\frac{n^{2n+1}}{2^{n^{2}+1}}$, utilizamos el criterio de la raíz como sigue:

$$A=\lim_{n \to \infty}\sqrt[n]{|c_{n}|}=\lim_{n \to \infty}\sqrt[n]{\frac{n^{2n+1}}{2^{n^{2}+1}}}=\lim_{n \to \infty}\frac{n^{2+1/n}}{2^{n+1/n}}=0 \Rightarrow R\rightarrow \infty$$

Por lo que el intervalo de convergencia es: $R \space \epsilon \space (-\infty, \infty)$ y la serie es convergente para cualquier valor de $x$.

Tarea moral

Los siguientes ejercicios no son para evaluación, pero son ejercicios para que practiques lo aprendido que te ayudaran en el desarrollo del entendimiento del tema, por lo que te invitamos a resolver los siguientes ejercicios propuestos relacionados con el tema visto.

Diga si la siguientes series convergen o diverge.

Los siguientes ejercicios te ayudarán a repasar los conceptos vistos en esta entrada.

  1. $$\sum_{n=0}^{\infty}\frac{n^{3}}{4^{n}}x^{n}$$
  2. $$\sum_{n=1}^{\infty}\frac{x^{n}}{n}$$
  3. $$\sum_{n=1}^{\infty}\frac{(2x)^{n}}{n^{2}}$$
  4. $$\sum_{n=1}^{\infty}(-1)^{n-1}\frac{x^{n}}{n}$$
  5. $$\sum_{n=1}^{\infty}\frac{(-3)^{n}x^{n}}{\sqrt{n+1}}$$

Más adelante…

En esta sección vimos las series de potencias y dos teoremas importantes para la convergencia de estas series que son el teorema de Abel y el teorema de Cauchy-Hadamard, en la siguiente sección veremos los polinomios de Taylor y de Mclaurin que están relacionados con estas series de potencias.

Entradas relacionadas

Ecuaciones Diferenciales I: Ecuaciones lineales de segundo orden con coeficientes variables – Soluciones en series de potencias respecto a puntos ordinarios

Por Omar González Franco

El mundo de las matemáticas no es un lugar aburrido en el que estar.
Es un lugar extraordinario; merece la pena pasar el tiempo allí.
– Marcus du Sautoy

Introducción

Hasta este punto de la unidad dos hemos desarrollado distintos métodos para resolver ecuaciones diferenciales lineales de orden superior, en particular de segundo orden con coeficientes constantes a excepción de la ecuación de Cauchy – Euler.

Para finalizar con la segunda unidad es el turno de estudiar las ecuaciones diferenciales lineales de segundo orden con coeficientes variables. Estas ecuaciones suelen ser mucho más complicadas de resolver ya que no se resuelven en términos de funciones elementales, sino que tienen forma de serie de potencias infinitas.

Nos parece adecuado comenzar esta entrada con un repaso sobre series de potencias, posteriormente veremos su utilidad en los métodos de resolución de las ecuaciones diferenciales antes mencionadas, así mismo, introduciremos algunos conceptos nuevos relacionados con el tipo de solución que tienen estas ecuaciones diferenciales.

Series de potencias

Algunas propiedades y conceptos importantes que debemos recordar son los siguientes.

Toda serie de potencias tiene un intervalo de convergencia.

Si $R> 0$, entonces la serie de potencias (\ref{1}) converge para $|x -a| < R$ y diverge para $|x -a| > R$.

Si la serie converge sólo en su centro $a$, entonces $R = 0$.

Si la serie converge para toda $x$, entonces se escribe $R = \infty$.

Una serie de potencias podría converger o no en los puntos extremos $a -R$ y $a + R$ de este intervalo.

El radio de convergencia también se puede determinar con las siguientes expresiones.

$$R = \left( \lim_{n \to \infty} \sqrt[n]{|c_{n}|} \right)^{-1} \hspace{1cm} o \hspace{1cm} R = \lim_{n \to \infty} \left| \dfrac{c_{n}}{c_{n + 1}} \right| \label{4} \tag{4}$$

Realicemos un ejemplo.

Ejemplo: Hallar el radio de convergencia y el intervalo de convergencia de la serie de potencias

$$\sum_{n = 1}^{\infty} \dfrac{n^{2}}{2^{n}} \left( x -1 \right)^{n}$$

Solución: Para determinar el radio de convergencia utilicemos la segunda expresión de (\ref{4}). De la serie de potencias identificamos que

$$c_{n} = \dfrac{n^{2}}{2^{n}} \hspace{1cm} y \hspace{1cm} c_{n + 1} = \dfrac{(n + 1)^{2}}{2^{n + 1}}$$

Calculemos el límite.

$$R = \lim_{n \to \infty} \left| \dfrac{c_{n}}{c_{n + 1}} \right| = \lim_{n \to \infty} \left| \dfrac{\dfrac{n^{2}}{2^{n}}}{\dfrac{(n + 1)^{2}}{2^{n + 1}}} \right| = 2 \lim_{n \to \infty} \left| \dfrac{n^{2}}{(n + 1)^{2}} \right|$$

Sabemos que

$$\lim_{n \to \infty} \left| \dfrac{n^{2}}{(n + 1)^{2}} \right| = 1$$

Por lo tanto, el radio de convergencia es $R = 2$.

Para determinar el intervalo de convergencia utilicemos la expresión (\ref{5}).

\begin{align*}
\lim_{n \to \infty} \left| \dfrac{c_{n + 1}(x -a)^{n + 1}}{c_{n}(x -a)^{n}} \right| &= \lim_{n \to \infty} \left| \dfrac{\dfrac{(n + 1)^{2}}{2^{n + 1}}(x -1)^{n + 1}}{\dfrac{n^{2}}{2^{n}}(x -1)^{n}} \right| \\
&= |x -1| \lim_{n \to \infty} \dfrac{2^{n}(n + 1)^{2}}{2^{n + 1}n^{2}} \\
&= \dfrac{1}{2} |x -1| \lim_{n \to \infty} \dfrac{n^{2} + 2n + 1}{n^{2}} \\
&= L
\end{align*}

Es claro que

$$\lim_{n \to \infty} \dfrac{n^{2} + 2n + 1}{n^{2}} = 1$$

Entonces,

$$\dfrac{1}{2} |x -1| = L$$

La condición de convergencia nos indica que $L < 1$, considerando esto tenemos que

\begin{align*}
\dfrac{1}{2} |x -1| &< 1 \\
|x -1| &< 2 \\
-2 < x -1 &< 2 \\
-1 < x &< 3
\end{align*}

Por lo tanto, el intervalo de convergencia es $I = (-1, 3)$.

Notemos que la mitad de la longitud del intervalo de convergencia efectivamente corresponde al valor del radio de convergencia obtenido.

$$R = \dfrac{3 -(-1)}{2} = \dfrac{4}{2} = 2$$

$\square$

Series de potencias como funciones

Nota: La convergencia en un extremo se podría perder por derivación o ganar por integración. Algo similar ocurre con los índices de una serie, supongamos que

$$y = \sum_{n = 0}^{\infty}c_{n}x^{n}$$

es una serie de potencias en $x$, las primeras dos derivadas están dadas como

$$\dfrac{dy}{dx} = \sum_{n = 0}^{\infty} n x^{n -1} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} = \sum_{n = 0}^{\infty}n(n -1)x^{n -2}$$

Sin embargo, notemos que el primer término en la primera derivada y los dos primeros términos de la segunda derivada son cero, entonces los podemos omitir y correr el índice para escribir

$$\dfrac{dy}{dx} = \sum_{n = 1}^{\infty} n x^{n -1} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} = \sum_{n = 2}^{\infty}n(n -1)x^{n -2}\label{7} \tag{7}$$

Un concepto de bastante importancia y utilidad en las próximas entradas es el siguiente.

Podemos hacer operaciones con series de potencias, a continuación se muestran algunas de ellas.

  • Suma: Dos series de potencias pueden sumarse término a término.

Sean

$$f(x) = \sum_{n = 0}^{\infty}c_{n}(x -a)^{n} \hspace{1cm} y \hspace{1cm} g(x) = \sum_{n = 0}^{\infty}b_{n}(x -a)^{n}$$

dos series de potencias con radio de convergencia $R> 0$, entonces

$$f(x) + g(x) = \sum_{n = 0}^{\infty}(c_{n} + b_{n})(x -a)^{n} \label{9} \tag{9}$$

Para toda $|x -a| < R$.

  • Producto: Dos series de potencias pueden multiplicarse término a término (cada término de la primera por cada término de la segunda).

Sean

$$f(x) = \sum_{n = 0}^{\infty}c_{n}(x -a)^{n} \hspace{1cm} y \hspace{1cm} g(x) = \sum_{n = 0}^{\infty}b_{n}(x -a)^{n}$$

dos series de potencias con radio de convergencia $R> 0$, entonces

$$f(x)g(x) = \sum_{n = 0}^{\infty}(c_{0}b_{n} + c_{1}b_{n -1} + \cdots + c_{n}b_{0})(x -a)^{n} \label{10} \tag{10}$$

Para toda $|x -a| < R$.

  • Derivación: Una serie de potencias puede derivarse término a término.

Sea

$$f(x) = \sum_{n = 0}^{\infty}c_{n}(x -a)^{n}$$

una serie de potencias convergente para $|x -a| < R$ con $R> 0$. La derivada de la serie $f$ es

$$F(x) = \dfrac{df}{dx} = \sum_{n = 1}^{\infty}nc_{n}(x -a)^{n -1} \label{11} \tag{11}$$

y también es convergente y tiene el mismo radio de convergencia que $f(x)$.

  • Integración: Una serie de potencias puede integrarse término a término.

Sea

$$f(x) = \sum_{n = 0}^{\infty}c_{n}(x -a)^{n}$$

una serie de potencias convergente para $|x -a| < R$ con $R> 0$. La integral de la serie $f$ es

$$F(x) = \int_{0}^{x}f(t)dt = \sum_{n = 0}^{\infty}\dfrac{c_{n}}{n + 1}(x -a)^{n + 1} \label{12} \tag{12}$$

y tiene a $R$ como radio de convergencia.

A lo largo de ésta y las siguientes entradas será de suma importancia y utilidad simplificar la suma de dos o más series de potencias, cada una expresada en notación de suma, en una sola expresión de suma, muchas veces esto implica que se deba hacer un cambio en el índice de la suma.

Para poder sumar dos series en necesario que ambos índices de las sumas comiencen con el mismo número y las potencias de $x$ sean las mismas y estén en fase. Por ejemplo, consideremos las siguientes dos series

$$f(x) = \sum_{n = 1}^{\infty} = \dfrac{n}{n+2}x^{n + 1} \hspace{1cm} y \hspace{1cm} g(x) = \sum_{n = 1}^{\infty} = \dfrac{1}{n^{2} + 1}x^{n + 1}$$

Como ambas series comienzan con el mismo número $n = 1$ y en ambas la potencia de $x$ es la misma $n + 1$, entonces podemos combinar ambas series en una sola de acuerdo a la expresión (\ref{9})

\begin{align*}
f(x) + g(x) &= \sum_{n = 1}^{\infty} \left[ \dfrac{n}{n+2} + \dfrac{1}{n^{2} + 1} \right]x^{n + 1} \\
&= \sum_{n = 1}^{\infty} \dfrac{n^{3} + 2n + 2}{n^{3} + 2n^{2} + n + 2}x^{n + 1}
\end{align*}

¿Pero que ocurre si no comienzan con el mismo número y/o las potencias de $x$ no coinciden?. En estos casos será necesario hacer un cambio en el índice de la suma y por tanto en la potencia de $x$. A continuación se muestra un ejemplo en el que describimos la forma de hacerlo.

Ejemplo: Reescribir la expresión

$$f(x) = \sum_{n = 1}^{\infty}2nc_{n}x^{n -1} + \sum_{n = 0}^{\infty}6c_{n}x^{n + 1}$$

como una sola serie de potencias cuyo término general tenga $x^{k}$.

Solución: Notemos que la potencia de $x$ en la primer serie para $n = 1$ es $x^{0}$, mientras que en la segunda serie para $n = 0$ es $x^{1}$, como ambas potencias son distintas decimos que no están en fase, para corregir esto y hacer que estén en fase extraemos el primer término de la primer serie.

$$\sum_{n = 1}^{\infty}2nc_{n}x^{n -1} = 2c_{1} + \sum_{n = 2}^{\infty}2nc_{n}x^{n -1}$$

Así, la potencia de $x$ para $n = 2$ es $x^{1}$. Con esto hemos logrado que ambas series estén en fase a pesar de que tengan distintas potencias en $x$ y comiencen con distintos números para $n$.

Procedemos a hacer el cambio de índice, para ello se toman como guía los exponentes de $x$. Para la primer serie tomamos $k = n -1$, de donde $n = k + 1$. Si $n = 2$, entonces $k = 1$ con esto podemos escribir a la primer serie de la siguiente manera.

$$\sum_{n = 1}^{\infty}2nc_{n}x^{n -1} = 2c_{1} + \sum_{k = 1}^{\infty}2(k + 1)c_{k + 1}x^{k}$$

Para la segunda serie tomamos $k = n + 1$, de donde $n = k -1$, si $n = 0$, entonces $k = 1$, así la segunda serie se puede escribir de la siguiente manera.

$$\sum_{n = 0}^{\infty}6c_{n}x^{n + 1} = \sum_{k = 1}^{\infty}6c_{k -1}x^{k}$$

Ahora podemos escribir

$$\sum_{n = 1}^{\infty}2nc_{n}x^{n -1} + \sum_{n = 0}^{\infty}6c_{n}x^{n + 1} = 2c_{1} + \sum_{k = 1}^{\infty}2(k + 1)c_{k + 1}x^{k} + \sum_{k = 1}^{\infty}6c_{k -1}x^{k}$$

Observemos que ambas series ya comienzan con el mismo número $k = 1$ y la potencia de $x$ es $k$ para ambas, entonces ya podemos combinar las series en una sola, de tal manera que

$$f(x) = 2c_{1} + \sum_{k = 1}^{\infty} \left[ 2(k + 1)c_{k + 1} + 6c_{k -1} \right] x^{k}$$

$\square$

En el caso de una sola serie es mucho mas sencillo pues basta tomar a $k$ como la potencia de $x$ y evaluar el valor del primer número en la serie, por ejemplo para la serie

$$\sum_{n = 1}^{\infty}nc_{n}x^{n + 2}$$

Si queremos que el termino $x$ tenga potencia $k$ hacemos $k = n + 2$, de donde $n = k -2$, la serie comienza en $n = 1$, sustituyendo en $k$ obtenemos que $k = 3$, por lo tanto la serie en términos del índice $k$ se puede escribir de la siguiente manera.

$$\sum_{n = 1}^{\infty}nc_{n}x^{n + 1} = \sum_{k = 3}^{\infty}(k -2)c_{k -2}x^{k}$$

Puedes desglosar ambas sumas para convencerte de la igualdad.

Hasta aquí concluimos nuestro repaso de series de potencias, es momento de aplicarlo en la resolución de ecuaciones diferenciales.

Soluciones en series de potencias de ecuaciones diferenciales

Las ecuaciones diferenciales lineales de segundo orden con coeficientes variables tienen la forma

$$a_{2}(x)\dfrac{d^{2}y}{dx^{2}} + a_{1}(x) \dfrac{dy}{dx} + a_{0}(x)y = g(x) \label{13} \tag{13}$$

Comenzaremos por considerar que $g(x) = 0$.

$$a_{2}(x)\dfrac{d^{2}y}{dx^{2}} + a_{1}(x) \dfrac{dy}{dx} + a_{0}(x)y = 0 \label{14} \tag{14}$$

Si dividimos la ecuación por $a_{2}(x) \neq 0$ y definimos

$$P(x) = \dfrac{a_{1}(x)}{a_{2}(x)} \hspace{1cm} y \hspace{1cm} Q(x) = \dfrac{a_{0}(x)}{a_{2}(x)}$$

podemos escribir la ecuación (\ref{14}) en su forma estándar.

$$\dfrac{d^{2}y}{dx^{2}} + P(x) \dfrac{dy}{dx} + Q(x)y = 0 \label{15} \tag{15}$$

En base a la ecuación estándar (\ref{15}) establecemos las siguientes definiciones.

De acuerdo a estas definiciones notamos que un punto singular $x_{0}$ es un punto no ordinario.

Realicemos un ejemplo.

Ejemplo: Hallar los puntos ordinarios y singulares de la ecuación diferencial

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

Solución: El primer paso es escribir a la ecuación diferencial en su forma estándar, para ello dividimos toda la ecuación por el coeficiente de la segunda derivada suponiendo que es distinto de cero.

\begin{align*}
\dfrac{d^{2}y}{dx^{2}} + \dfrac{x^{3}(x^{2} -1)}{x^{2}(x -1)} \dfrac{dy}{dx} + \dfrac{x}{x^{2}(x -1)}y &= 0 \\
\dfrac{d^{2}y}{dx^{2}} + x(x + 1) \dfrac{dy}{dx} + \dfrac{1}{x(x -1)}y &= 0
\end{align*}

Identificamos que

$$P(x) = x(x + 1) \hspace{1cm} y \hspace{1cm} Q(x) = \dfrac{1}{x(x -1)}$$

Para el caso de la función $P(x)$ notamos que es analítica para toda $x \in \mathbb{R}$, mientras que la función $Q(x)$ no está definida en $x = 0$ ni $x = 1$, es decir, no es analítica en dichos puntos.

Por lo tanto, los puntos ordinarios de la ecuación diferencial son todas las $x \in \mathbb{R}$ excepto $x = 0$ y $x = 1$, éstos puntos corresponde a los puntos singulares de la ecuación.

$\square$

Una observación interesante es que la ecuación de Cauchy-Euler

$$ax^{2} \dfrac{d^{2}y}{dx^{2}} + bx \dfrac{dy}{dx} + cy = 0 \label{16} \tag{16}$$

en su forma estándar

$$\dfrac{d^{2}y}{dx^{2}} + \dfrac{b}{ax} \dfrac{dy}{dx} + \dfrac{c}{ax^{2}}y = 0 \label{17} \tag{17}$$

nos muestra que las funciones

$$P(x) = \dfrac{b}{ax} \hspace{1cm} y \hspace{1cm} Q(x) = \dfrac{c}{ax^{2}}$$

no están definidas en $x = 0$, por tanto $x = 0$ es un punto singular y todos los demás puntos (reales o complejos) son puntos ordinarios, es por ello que toda la teoría realizada en la entrada correspondiente fue para $x > 0$.

De acuerdo al título de esta entrada, nos enfocaremos en soluciones respecto a puntos ordinarios, sin embargo, cabe mencionar que en la siguiente entrada estudiaremos soluciones respecto a puntos singulares y será necesario hacer una distinción entre dos tipos de puntos singulares que definiremos como punto singular regular y punto singular irregular. Estos conceptos los revisaremos en la siguiente entrada.

Como ya hemos mencionando, las soluciones de la ecuación diferencial (\ref{15}) son soluciones en forma de series de potencias. Si una ecuación diferencial es analítica en un punto $x_{0}$, entonces su solución también lo es en $x_{0}$, y como dicha solución será una función desarrollable en series de potencias, podemos suponer que, en forma general, tendrá la siguiente forma.

$$y(x) = \sum_{n = 0}^{\infty}c_{n}(x -x_{0})^{n} \label{18} \tag{18}$$

donde $c_{n}$ cambia para cada función específica.

A continuación enunciamos el teorema que establece la existencia y forma de las soluciones de (\ref{15}).

Una solución en serie converge, por lo menos, en un intervalo definido por $|x -x_{0}| < R$, donde $R$ es la distancia desde $x_{0}$ al punto singular más cercano, es decir, es el valor mínimo o límite inferior del radio de convergencia de las soluciones en serie de la ecuación diferencial respecto a $x_{0}$.

La demostración a este teorema suele ser bastante larga pero intuitiva. En esta ocasión no lo demostraremos y en su lugar desarrollaremos varios ejemplos que ilustran el resultado. Sin embargo, en la sección de videos de este mismo curso se puede encontrar con todo detalle la demostración de este teorema, además del método para hallar el radio de convergencia de la solución en serie de potencias cerca de un punto ordinario.

Método de resolución

Si bien, en la demostración del teorema de existencia y forma de la solución en series de potencias se describe el método de resolución, nosotros vamos a describirlo de manera breve y realizaremos algunos ejemplos para que quede bastante claro.

Recordemos que el método de coeficientes indeterminados desarrollado para ecuaciones diferenciales lineales no homogéneas de segundo orden con coeficientes constantes ya involucraba soluciones en forma de series de potencias y lo que hacíamos al final del método era igualar los coeficientes de ambos lados de la ecuación para satisfacer la igualdad, la diferencia ahora es que el lado derecho de la ecuación es cero y no una función $g(x)$, sin embargo el procedimiento es bastante similar.

Debido a que se trata de un método bastante laborioso, por simplicidad encontraremos soluciones en series de potencias sólo con respecto al punto ordinario $x_{0} = 0$. Así, las soluciones serán de la forma

$$y(x) = \sum_{n = 0}^{\infty}c_{n}x^{n} \label{19} \tag{19}$$

La descripción del método se muestra a continuación:

  • El método de resolución implica considerar la solución (\ref{19}) y su primera y segunda derivada (\ref{7}) para sustituirlas en la ecuación diferencial (\ref{14}).

$$a_{2}(x) \left[ \sum_{n = 2}^{\infty}n(n -1)c_{n}x^{n -2} \right] + a_{1}(x) \left[ \sum_{n = 1}^{\infty}nc_{n}x^{n -1} \right] + a_{0}(x) \left[ \sum_{n = 0}^{\infty}c_{n}x^{n} \right] = 0$$

  • El siguiente paso es reescribir toda la ecuación en una sola serie lo que, en la mayoría de los casos, requerirá de hacer cambios de índices para que se tenga la misma potencia de $x$.
  • Como el resultado será idénticamente cero será necesario que el coeficiente de cada potencia de $x$ se iguale a cero. Como veremos más adelante, esto nos generará una ecuación general para los coeficientes de $y(x)$, dicha expresión se conoce como relación de recurrencia.
  • La tarea final será usar la relación de recurrencia para obtener el valor de los coeficientes $c_{n}$ de (\ref{19}) y con ello la forma de la solución de la ecuación diferencial en cuestión.

Es importante aclarar que la sola suposición de la solución (\ref{19}) conduce a dos conjuntos de coeficientes, de manera que se tendrán dos series de potencias distintas $y_{1}$ y $y_{2}$, ambas desarrolladas respecto al punto ordinario $x_{0}$. Se puede demostrar que la solución general de la ecuación diferencial (\ref{14}) es

$$y(x) = C_{1}y_{1}(x) + C_{2}y_{2}(x) \label{20} \tag{20}$$

en donde $C_{1} = c_{0}$ y $C_{2} = c_{1}$, es decir, los primeros coeficientes de la serie (\ref{19}).

Este método no solo es aplicable a ecuaciones de la forma (\ref{14}), sino que se puede aplicar a distintas ecuaciones que satisfagan las propiedades necesarias descritas a lo largo de la entrada.

Para comprender el método resolvamos una ecuación bastante sencilla de primer orden y veamos que resultado obtenemos.

Ejemplo: Determinar la solución de la ecuación diferencial

$$\dfrac{dy}{dx} -y = 0$$

usando series de potencias respecto al punto ordinario $x_{0} = 0$.

Solución: La solución debe ser de la forma

$$y = \sum_{n = 0}^{\infty}c_{n}x^{n}$$

La derivada de esta función es

$$\dfrac{dy}{dx} = \sum_{n = 1}^{\infty}nc_{n}x^{n -1}$$

Sustituimos en la ecuación diferencial.

$$\sum_{n = 1}^{\infty}nc_{n}x^{n -1} -\sum_{n = 0}^{\infty}c_{n}x^{n} = 0$$

Hay que reescribir esta ecuación en una sola serie en la que la potencia de $x$ sea $k$.

Guiándonos en los exponentes de $x$, en la primer serie tomamos $k = n -1$, de donde $n = k + 1$, si la serie comienza en $n = 1$, entonces $k = 1 -1 = 0$. En el caso de la segunda serie basta hacer $k = n$, entonces tenemos que

$$\sum_{k = 0}^{\infty}(k + 1)c_{k + 1}x^{k} -\sum_{k = 0}^{\infty}c_{k}x^{k} = 0$$

Ahora si podemos unir las series en una sola de acuerdo a (\ref{9})

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

Como $x^{k}\neq 0$ por ser la solución propuesta, entonces necesariamente

$$(k + 1)c_{k + 1} -c_{k} = 0$$

Como $k$ es un número entero que comienza en cero hacía infinito, entonces $k$ no puede ser negativo, lo que significa que no hay valor de $k$, tal que $k + 1 = 0$, es así que podemos despejar a $c_{k + 1}$ de la expresión anterior sin problema.

$$c_{k + 1} = \dfrac{c_{k}}{k + 1}, \hspace{1cm} k = 0, 1, 2, 3, \cdots$$

Ésta última expresión corresponde a la relación de recurrencia, de la que se obtiene cada una de las constantes para cada uno de los términos de la serie solución.

Comencemos con $k = 0$.

$$c_{1} = \dfrac{c_{0}}{0 + 1} =c_{0}$$

Para $k = 1$, tenemos

$$c_{2} = \dfrac{c_{1}}{1 + 1} = \dfrac{c_{0}}{2}$$

$k = 2$.

$$c_{3} = \dfrac{c_{2}}{2 + 1} = \dfrac{c_{0}}{6}$$

$k = 3$.

$$c_{4} = \dfrac{c_{3}}{3 + 1} = \dfrac{c_{0}}{24}$$

Etcétera, entonces la solución va teniendo la siguiente forma.

\begin{align*}
y(x) &= c_{0} + c_{1}x + c_{2}x^{2} + c_{3}x^{3} + c_{4}x^{4} + \cdots \\
&= c_{0} + c_{0}x + \dfrac{c_{0}}{2}x^{2} + \dfrac{c_{0}}{6}x^{3} + \dfrac{c_{0}}{24}x^{4} + \cdots \\
&= c_{0} \left[1 + x + \dfrac{x^{2}}{2} + \dfrac{x^{3}}{6} + \dfrac{x^{4}}{24} + \cdots \right] \\
&= c_{0} \left[ 1 + x + \dfrac{x^{2}}{2!} + \dfrac{x^{3}}{3!} + \dfrac{x^{4}}{4!} + \cdots \right]
\end{align*}

En algunas ocasiones las series de potencias resultan ser series conocidas, como lo es en este caso, pues sabemos que

$$e^{x} = \sum_{n = 0}^{\infty}\dfrac{x^{n}}{n!} = 1 + x + \dfrac{x^{2}}{2!} + \dfrac{x^{3}}{3!} + \dfrac{x^{4}}{4!} + \cdots$$

Por lo tanto, si definimos $c = c_{0}$, la solución de la ecuación diferencial es

$$y(x) = ce^{x}$$

Para asegurarnos del resultado se puede sustituir en la ecuación diferencial y ver que la satisface, o bien, podemos usar separación de variables para resolver la ecuación y verificar el resultado.

\begin{align*}
\dfrac{dy}{dx} -y &= 0 \\
\dfrac{dy}{dx} &= y \\
\dfrac{1}{y} \dfrac{dy}{dx} &= 1 \\
\int{\dfrac{dy}{y}} &= \int{dx} \\
\ln(y) &= x + k \\
y &= e^{x + k} \\
y &= e^{k}e^{x} \\
y(x) &= ce^{x}
\end{align*}

¡Verificado!. Interesante ¿no?.

$\square$

Con este ejemplo se espera que se comprenda la noción del método, como se puede notar es un proceso largo a pesar de ser una ecuación muy simple. Concluiremos esta entrada resolviendo dos ecuaciones diferenciales de las que si estamos interesados en resolver, es decir, de la forma (\ref{14}).

Ejemplo: Resolver la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} + xy = 0$$

respecto al punto ordinario $x_{0} = 0$.

Solución: Debido a que no hay puntos singulares, el teorema garantiza dos soluciones en serie de potencias centradas en $x_{0} = 0$, convergentes para $|x|< \infty$.

Consideremos la solución

$$y = \sum_{n = 0}^{\infty}c_{n}x^{n}$$

y su segunda derivada

$$\dfrac{d^{2}y}{dx^{2}} = \sum_{n = 2}^{\infty}n(n -1)c_{n}x^{n -2}$$

Sustituyamos en la ecuación diferencial.

\begin{align*}
\dfrac{d^{2}y}{dx^{2}} + xy &= \left[ \sum_{n = 2}^{\infty}c_{n}n(n -1)x^{n -2} \right] + x \left[ \sum_{n = 0}^{\infty}c_{n}x^{n} \right] \\
&= \sum_{n = 2}^{\infty}c_{n}n(n -1)x^{n -2} + \sum_{n = 0}^{\infty}c_{n}x^{n + 1}
\end{align*}

Para que practiques muestra que

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

Por lo tanto,

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

Para que esta igualdad se cumpla es necesario que el coeficiente de cada potencia de $x$ se iguale a cero. Para el caso de la potencia $k = 0$ tenemos que $2c_{2} = 0$, de donde $c_{2} = 0$, para el resto de potencias formamos la relación de recurrencia.

$$(k + 1)(k + 2)c_{k + 2} + c_{k -1} = 0, \hspace{1cm} k = 1, 2, 3, \cdots$$

Esta expresión determina los coeficientes $c_{k}$ que buscamos. Como $(k + 1)(k + 2) \neq 0$ para los valores de $k$, podemos escribir $c_{k + 2}$ en términos de $c_{k -1}$.

$$c_{k + 2} = -\dfrac{c_{k -1}}{(k + 1)(k + 2)}, \hspace{1cm} k = 1, 2, 3, \cdots$$

Esta relación genera coeficientes consecutivos de la solución propuesta una vez que $k$ toma los valores enteros sucesivos indicados.

Comencemos con $k = 1$.

$$c_{3} = -\dfrac{c_{0}}{2\cdot 3}$$

Para $k = 2$, se tiene

$$c_{4} = -\dfrac{c_{1}}{3 \cdot 4}$$

Para $k = 3$ hacemos uso de que $c_{2} = 0$.

$$c_{5} = -\dfrac{c_{2}}{4 \cdot 5} = 0$$

A partir de $k = 4$ hacemos uso de los valores previos.

$$c_{6} = -\dfrac{c_{3}}{5 \cdot 6} = -\left( -\dfrac{c_{0}}{2\cdot 3} \right) \dfrac{1}{5 \cdot 6} = \dfrac{1}{2 \cdot 3 \cdot 5 \cdot 6} c_{0}$$

$k = 5$.

$$c_{7} = -\dfrac{c_{4}}{6 \cdot 7}=-\left( -\dfrac{c_{1}}{3 \cdot 4} \right) \dfrac{1}{6 \cdot 7} = \dfrac{1}{3 \cdot 4 \cdot 6 \cdot 7}c_{1}$$

Para $k = 6$ recordamos que $c_{5} = 0$.

$$c_{8} = -\dfrac{c_{5}}{7 \cdot 8} = 0$$

$k = 7$.

$$c_{9} = -\dfrac{c_{6}}{8 \cdot 9} = \dfrac{1}{2 \cdot 3 \cdot 5 \cdot 6 \cdot 8 \cdot 9}c_{0}$$

$k = 8$.

$$c_{10} = -\dfrac{c_{7}}{9 \cdot 10} = \dfrac{1}{3 \cdot 4 \cdot 6 \cdot 7 \cdot 9 \cdot 10}c_{1}$$

$k = 9$.

$$c_{11} = -\dfrac{c_{8}}{10 \cdot 11} = 0$$

Podemos hacer estos cálculos para la $k$ que deseemos, el objetivo es intentar determinar que tipo de serie numérica es la que se logra formar. En este caso nos detendremos hasta $k = 9$, con ello hemos logrado obtener los primeros $11$ coeficientes de la solución que buscamos (recordemos que $c_{0}$ y $c_{1}$ tienen valores arbitrarios).

\begin{align*}
y(x) &= c_{0} + c_{1} x + c_{2}x^{2} + c_{3}x^{3} +c_{4}x^{4} + c_{5}x^{5} + c_{6}x^{6} \\
&+ c_{7}x^{7} + c_{8}x^{8} + c_{9}x^{9} + c_{10}x^{10} + c_{11}x^{11} + \cdots
\end{align*}

Sustituyamos los coeficientes obtenidos.

\begin{align*}
y(x) &= c_{0} + c_{1}x + 0 -\dfrac{c_{0}}{2 \cdot 3}x^{3} -\dfrac{c_{1}}{3 \cdot 4}x^{4} + 0 + \dfrac{c_{0}}{2 \cdot 3 \cdot 5 \cdot 6}x^{6} + \dfrac{c_{1}}{3 \cdot 4 \cdot 6 \cdot 7}x^{7} + 0 \\
&-\dfrac{c_{0}}{2 \cdot 3 \cdot 5 \cdot 6 \cdot 8 \cdot 9}x^{9} -\dfrac{c_{1}}{3 \cdot 4 \cdot 6 \cdot 7 \cdot 9 \cdot 10}x^{10} + 0 + \cdots
\end{align*}

Para obtener la solución general

$$y(x) = c_{0}y_{1}(x) + c_{1}y_{2}(x)$$

agrupemos los términos que contienen $c_{0}$ y por otro lado los que tienen $c_{1}$.

\begin{align*}
y(x) &= c_{0} \left[ 1 -\dfrac{1}{2 \cdot 3}x^{3} + \dfrac{1}{2 \cdot 3 \cdot 5 \cdot 6}x^{6} -\dfrac{1}{2 \cdot 3 \cdot 5 \cdot 6 \cdot 8 \cdot 9}x^{9} + \cdots \right] \\
&+ c_{1} \left[ x -\dfrac{1}{3 \cdot 4}x^{4} + \dfrac{1}{3 \cdot 4 \cdot 6 \cdot 7}x^{7} -\dfrac{1}{3 \cdot 4 \cdot 6 \cdot 7 \cdot 9 \cdot 10}x^{10} + \cdots \right]
\end{align*}

Por lo tanto,

\begin{align*}
y_{1}(x) &= 1 -\dfrac{1}{2 \cdot 3}x^{3} + \dfrac{1}{2 \cdot 3 \cdot 5 \cdot 6}x^{6} -\dfrac{1}{2 \cdot 3 \cdot 5 \cdot 6 \cdot 8 \cdot 9}x^{9} + \cdots \\
&= 1 + \sum_{k = 1}^{\infty}\dfrac{(-1)^{k}}{2 \cdot 3 \cdots (3k -1)(3k)}
\end{align*}

y

\begin{align*}
y_{2}(x) &= x -\dfrac{1}{3 \cdot 4}x^{4} + \dfrac{1}{3 \cdot 4 \cdot 6 \cdot 7}x^{7} -\dfrac{1}{3 \cdot 4 \cdot 6 \cdot 7 \cdot 9 \cdot 10}x^{10} + \cdots \\
&= x + \sum_{k = 1}^{\infty}\dfrac{(-1)^{k}}{3 \cdot 4 \cdots (3k)(3k + 1)}x^{3k + 1}
\end{align*}

Con esto hemos concluido el ejercicio. Los coeficientes $c_{0}$ y $c_{1}$ quedan completamente indeterminados de manera que se pueden elegir de forma arbitraria.

Por el teorema de existencia y forma de la solución también se puede deducir que las series que forman a $y_{1}$ y $y_{2}$ convergen para $|x|< \infty$.

$\square$

Como dato interesante, la ecuación diferencial que acabamos de resolver es una forma de lo que se conoce como ecuación de Airy y se encuentra en el estudio de la difracción de la luz, la difracción de ondas de radio alrededor de la superficie de la tierra, la aerodinámica y la deflexión de una columna vertical delgada uniforme que se curva bajo su propio peso.

Realicemos un ejemplo más en el que los coeficientes de la ecuación no sean polinomios, esto nos permitirá poner en práctica la multiplicación de dos series de potencias.

Ejemplo: Resolver la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} + \cos (x) y = 0$$

respecto al punto ordinario $x_{0} = 0$.

Solución: Se puede comprobar que la función coseno es analítica en $x = 0$, esto verifica que efectivamente $x_{0} = 0$ es un punto ordinario. De hecho, al ser analítica en $x = 0$ su serie de Maclaurin es

\begin{align*}
\cos (x) &= 1 -\dfrac{x^2}{2!} + \dfrac{x^4}{4!} -\cdots + \dfrac{(-1)^kx^{2k}}{(2k)!} + \cdots \\
&= \sum_{n = 0}^{\infty} \dfrac{(-1)^nx^{2n}}{(2n)!}
\end{align*}

Resolvamos la ecuación. Consideremos la solución

$$y = \sum_{n = 0}^{\infty}c_{n}x^{n}$$

y su segunda derivada

$$\dfrac{d^{2}y}{dx^{2}} = \sum_{n = 2}^{\infty}n(n -1)c_{n}x^{n -2}$$

Sustituyamos en la ecuación diferencial.

$$\dfrac{d^{2}y}{dx^{2}} + \cos (x) y = \sum_{n = 2}^{\infty}n(n -1)c_{n}x^{n -2} + \left[ \sum_{n = 0}^{\infty} \dfrac{(-1)^nx^{2n}}{(2n)!} \right] \sum_{n = 0}^{\infty}c_{n}x^{n} = 0$$

En este caso no intentaremos reescribir la ecuación en una sola serie ya que puede ser más complicado al tratarse de un producto de series, en su lugar vamos a determinar el valor de los coeficientes de cada $x^{k}$, $k = 0, 1, 2, 3, \cdots$, realizando las operaciones correspondientes, para ello desglosemos las sumas para los primeros términos. Por un lado

$$\sum_{n = 2}^{\infty}n(n -1)c_{n}x^{n -2} = 2c_{2} + 6c_{3}x + 12c_{4}x^{2} + 20c_{5}x^{3} + \cdots$$

Por otro lado,

$$\left[ \sum_{n = 0}^{\infty} \dfrac{(-1)^nx^{2n}}{(2n)!} \right] \sum_{n = 0}^{\infty}c_{n}x^{n} = \left( 1 -\dfrac{x^{2}}{2!} + \dfrac{x^{4}}{4!} + \cdots \right) \left( c_{0} + c_{1}x + c_{2}x^{2} + c_{3}x^{3} + \cdots \right)$$

Si se hacen las cuentas correspondientes podremos obtener los coeficientes de cada $x^{k}$, $k = 0, 1, 2, 3, \cdots$.

Hasta $k = 3$ se obtiene lo siguiente.

$$(2c_{2} + c_{0}) + (6c_{3} + c_{1})x +\left( 12c_{4} + c_{2} -\dfrac{1}{2}c_{0} \right) x^{2} + \left( 20c_{5} + c_{3} -\dfrac{1}{2}c_{1} \right)x^{3} + \cdots = 0$$

Igualamos cada coeficiente a cero.

\begin{align*}
2c_{2} + c_{0} &= 0 \\
6c_{3} + c_{1} &= 0 \\
12c_{4} + c_{2} -\dfrac{1}{2}c_{0} &= 0 \\
20c_{5} + c_{3} -\dfrac{1}{2}c_{1} &= 0 \\
&\vdots
\end{align*}

etcétera. Esto nos da como resultados

\begin{align*}
c_{2} &= -\dfrac{1}{2}c_{0} \\
c_{3} &= -\dfrac{1}{6}c_{1} \\
c_{4} &= \dfrac{1}{12}c_{0} \\
c_{5} &= \dfrac{1}{30}c_{1} \\
&\vdots
\end{align*}

En este caso no se obtuvo una relación de recurrencia, pero $c_{0}$ y $c_{1}$ siguen siendo coeficientes indeterminados que pueden tomar valores arbitrarios. Sustituyendo los valores determinados en la solución propuesta se obtiene

\begin{align*}
y(x) &= c_{0} + c_{1}x -\dfrac{c_{0}}{2}x^{2} -\dfrac{c_{1}}{6}x^{3} + \dfrac{c_{0}}{12}x^{4} + \dfrac{c_{1}}{30}x^{5} + \cdots \\
&= c_{0} \left[ 1 -\dfrac{1}{2}x^{2} + \dfrac{1}{12}x^{4} + \cdots \right] + c_{1}\left[ x -\dfrac{1}{6}x^{3} + \dfrac{1}{30}x^{5} + \cdots \right]
\end{align*}

Recordando que la solución general es

$$y(x) = c_{0}y_{1}(x) + c_{1}y_{2}(x)$$

entonces,

$$y_{1}(x) = 1 -\dfrac{1}{2}x^{2} + \dfrac{1}{12}x^{4} + \cdots$$

y

$$y_{2}(x) = x -\dfrac{1}{6}x^{3} + \dfrac{1}{30}x^{5} + \cdots$$

Ambas series de potencias convergen para $|x| < \infty$.

$\square$

Con esto concluimos 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. Determinar el radio de convergencia y el intervalo de convergencia de las siguientes series de potencias.
  • $\sum_{n = 1}^{\infty}\dfrac{2^{n}}{n}x^{n}$
  • $\sum_{n = 1}^{\infty}\dfrac{n}{n + 2}x^{n}$
  • $\sum_{n = 1}^{\infty}\dfrac{(x -1)^{n}}{n!}$
  1. Reescribir la siguiente expresión como una sola serie de potencias cuyo término general tenga $x^{k}$.
  • $\sum_{n = 2}^{\infty}n(n -1)c_{n}x^{n} + 2 \sum_{n = 2}^{\infty}n(n -1)c_{n}x^{n -2} + 3 \sum_{n = 1}^{\infty}nc_{n}x^{n}$
  1. Comprobar por sustitución directa que la siguiente serie de potencias es una solución particular de la ecuación diferencial dada.
  • $y(x) = \sum_{n = 0}^{\infty}\dfrac{(-1)^{n}}{2^{2n}(n!)^{2}}x^{2n}, \hspace{1cm} x\dfrac{d^{2}y}{dx^{2}} + \dfrac{dy}{dx} + xy = 0$
  1. Encontrar la solución general en series de potencias de las siguientes ecuaciones diferenciales respecto al punto ordinario $x_{0} = 0$.
  • $\dfrac{d^{2}y}{dx^{2}} + x^{2} \dfrac{dy}{dx} + xy = 0$
  • $\dfrac{d^{2}y}{dx^{2}} + \sin (x) y = 0$
  1. Usar el método de series de potencias para resolver el siguiente problema con valores iniciales.
  • $(x + 1) \dfrac{d^{2}y}{dx^{2}} -(2 -x) \dfrac{dy}{dx} + y = 0, \hspace{1cm} y(0) = 2, \hspace{0.5cm} y^{\prime}(0) = -1$

Más adelante…

En esta entrada aprendimos a resolver ecuaciones diferenciales de segundo orden con coeficientes variables respecto al punto ordinario $x_{0} = 0$.

En la siguiente entrada resolveremos ecuaciones diferenciales del mismo tipo, pero ahora con respecto a puntos singulares. El método de resolución es conocido como Método de Frobenius.

Entradas relacionadas

Ecuaciones Diferenciales I – Videos: Ecuaciones lineales de segundo orden con coeficientes variables. Solución por series de potencias cerca de un punto ordinario

Por Eduardo Vera Rosales

Introducción

A lo largo de las entradas anteriores que forman parte de la segunda unidad hemos estudiado a detalle ecuaciones lineales de segundo orden con coeficientes constantes, es decir, ecuaciones de la forma $$a\frac{d^{2}y}{dt^{2}}+b\frac{dy}{dt}+cy=g(t), \,\,\,\,\, a\neq 0$$ y hemos desarrollado diversos métodos para resolverlas. Es momento de revisar ecuaciones lineales de segundo orden, pero ahora con coeficientes variables, es decir, del tipo $$a_{0}(t)\frac{d^{2}y}{dt^{2}}+a_{1}(t)\frac{dy}{dt}+a_{2}(t)y=g(t).$$

Hallar soluciones para este tipo de ecuaciones no resulta tan sencillo como para el caso con coeficientes constantes, y en ocasiones no podremos encontrar soluciones en términos de funciones elementales como polinomios, exponenciales, trigonométricas, etc., por lo que una manera de hallar soluciones es suponiendo que la solución puede escribirse como una serie de potencias alrededor de un punto dado.

Estudiaremos entonces soluciones por series de potencias en dos tipos de puntos: cuando los coeficientes tienen desarrollo en series de Taylor alrededor del punto dado, y cuando lo anterior no ocurre. En particular, en esta entrada revisaremos el primer caso. Definiremos los conceptos de puntos ordinarios y singulares, y demostraremos la existencia de soluciones en series de potencias cerca de un punto ordinario,.

¡Manos a la obra!

Soluciones en series de potencias cerca de un punto ordinario

En el primer video ofrecemos la definición de puntos ordinarios y puntos singulares, y probamos la existencia de soluciones en series de potencias cerca de un punto ordinario, a la ecuación diferencial $$a_{0}(t)\frac{d^{2}y}{dt^{2}}+a_{1}(t)\frac{dy}{dt}+a_{2}(t)y=0.$$ La solución encontrada será, además, la solución general a la ecuación diferencial.

Radio de convergencia de la solución en serie de potencias cerca de un punto ordinario

En el segundo video de la entrada encontramos el radio de convergencia para la solución en serie de potencias cerca de un punto ordinario.

Ejemplos

En el último video de la entrada resolvemos un par de ejemplos de ecuaciones diferenciales con coeficientes variables, con el método desarrollado a lo largo de esta misma entrada.

Tarea moral

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

  • ¿Qué sucede si suponemos que $a_{0}=0$ en la demostración del primer video?
  • ¿Qué pasa si suponemos que $c=1$ en la demostración del primer video?
  • Prueba que las series de potencias que aparecen en la solución general a la ecuación diferencial $$\frac{d^{2}y}{dt^{2}}-ty=0$$ son soluciones particulares a la misma ecuación, y que estas son linealmente independientes. Por tanto, la solución general efectivamente lo es para la ecuación diferencial.
  • Encuentra la solución general a la ecuación $$\frac{d^{2}y}{dt^{2}}-y=0$$ usando series de potencias alrededor de $t_{0}=0$.
  • Encuentra la solución al problema de valor inicial $$\frac{d^{2}y}{dt^{2}}-ty=0$$ $$y(1)=0; \,\,\,\,\, \frac{dy}{dt}(1)=2$$ calculando una solución por serie de potencias alrededor de $t_{0}=1$.

Más adelante

Terminamos de estudiar las soluciones cerca de un punto ordinario. Lo siguientes será revisar el caso cuando el punto en cuestión no es un punto ordinario, es decir, es un punto singular de nuestra ecuación diferencial $$a_{0}(t)\frac{d^{2}y}{dt^{2}}+a_{1}(t)\frac{dy}{dt}+a_{2}(t)y=0.$$

Pero antes analizaremos un caso particular sencillo de resolver: la ecuación de Euler que tiene la forma $$t^{2}\frac{d^{2}y}{dt^{2}}+\alpha t\frac{dy}{dt}+\beta y=0.$$

A partir de la solución para esta ecuación podremos generalizar más adelante el método a una clase más general de ecuaciones diferenciales con puntos singulares.

¡No se lo pierdan!

Entradas relacionadas