Ir al contenido principal

¿Como resolver una ecuación diferencial por el método de series de potencias? 2

La ecuación diferencial a resolver es:
\[y'+y=0\]
Proponemos una solución como en el punto anterior con sus correspondientes derivadas hasta primer orden:
\[y=\sum_{n=0}^{\infty}c_{n}x^{n}\]
\[y'=\sum_{n=1}^{\infty}c_{n}nx^{n-1}\]
Reemplazamos en nuestra ecuación diferencial y tenemos:
\[\sum_{n=1}^{\infty}c_{n}nx^{n-1}+\sum_{n=0}^{\infty}c_{n}x^{n}=0\]
Realizamos un cambio de indices para la primera sumatoria con $k=n-1$, y por lo tanto $n=k+1$ y para la segunda sumatoria con $k=n$, reemplazamos en nuestra ecuación diferencial y nos quedara:
\[\sum_{k=0}^{\infty}c_{k+1}(k+1)x^{k}+\sum_{k=0}^{\infty}c_{k}x^{k}=0\]
Que podemos escribir de la forma:
\[\sum_{k=0}^{\infty}[c_{k+1}(k+1)+c_{k}]x^{k}=0\]
Luego:
\[\sum_{k=0}^{\infty}\underbrace{[c_{k+1}(k+1)+c_{k}]}_{=0}\underbrace{x^{k}}_{\neq0}=0\]
Obtenemos una ecuación característica de la cual vamos a sacar los valores de la ecuación diferencial:
\[c_{k+1}(k+1)+c_{k}=0\]
\[c_{k+1}(k+1)=-c_{k}\]
\[c_{k+1}=-\frac{c_{k}}{(k+1)} \]
Para $k=0$:
\[c_{0+1}=-\frac{c_{0}}{(0+1)}\]
\[c_{1}=-\frac{c_{0}}{(1)}=-c_{0} \]
Para $k=1$:
\[c_{1+1}=-\frac{c_{1}}{(1+1)}\]
\[c_{2}=-\frac{c_{1}}{(2)}=-\frac{(-c_{0})}{2}=\frac{c_{0}}{2}\]
Para $k=2$:
\[c_{2+1}=-\frac{c_{2}}{(2+1)}\]
\[c_{3}=-\frac{c_{2}}{(3)}=-\frac{c_{0}}{3 \cdot 2}\]
Para $k=3$:
\[c_{3+1}=-\frac{c_{3}}{(3+1)}\]
\[c_{4}=-\frac{c_{3}}{(4)}=-\frac{-c_{0}}{4 \cdot 3 \cdot 2}=\frac{c_{0}}{4 \cdot 3 \cdot 2}\]
Para $k=4$:
\[c_{4+1}=-\frac{c_{4}}{(4+1)}\]
\[c_{5}=-\frac{c_{4}}{(5)}=-\frac{c_{0}}{5 \cdot 4 \cdot 3 \cdot 2}\]
Luego la solución en series de potencias queda de la forma:
\[y=c_{0}+c_{1}x+c_{2}x^{2}+c_{3}x^{3}+c_{4}x^{4}+c_{5}x^{5}+...\]
Reemplazamos las constantes y tenemos:
\[y=c_{0}-c_{0}x+c_{0}\frac{x^{2}}{2}-c_{0}\frac{x^{3}}{3 \cdot 2}+c_{0}\frac{x^{4}}{4 \cdot 3 \cdot 2}-c_{0}\frac{x^{5}}{5 \cdot 4 \cdot 3 \cdot 2}+...\]
Factorizando $c_{0}$ y después es fácil obtener el término n-esimo:
\[y=c_{0}\left(1-x+\frac{x^{2}}{2}-\frac{x^{3}}{3 \cdot 2}+\frac{x^{4}}{4 \cdot 3 \cdot 2}-\frac{x^{5}}{5 \cdot 4 \cdot 3 \cdot 2}+...(-1)^{n}\frac{x^{n}}{n!}\right)\]
Y si nos damos cuenta el desarrollo en series de potencias de la función $e^{x}$ es de la forma:
\[e^{x}=1+x+\frac{x^{2}}{2}+\frac{x^{3}}{3 \cdot 2}+\frac{x^{4}}{4 \cdot 3 \cdot 2}+\frac{x^{5}}{5 \cdot 4 \cdot 3 \cdot 2}+...\frac{x^{n}}{n!}\]
Luego si hacemos el desarrollo en serie de potencias de la función $e^{-x}$ tendremos que cambiar el signo en las potencias impares:
\[e^{-x}=1-x+\frac{x^{2}}{2}-\frac{x^{3}}{3 \cdot 2}+\frac{x^{4}}{4 \cdot 3 \cdot 2}-\frac{x^{5}}{5 \cdot 4 \cdot 3 \cdot 2}+...(-1)^{n}\frac{x^{n}}{n!}\]
Que corresponde al término entre paréntesis que acompaña a la constante $c_{0}$, por lo tanto la solución de la ecuación diferencial es:
\[y=c_{0}e^{-x}\]

Etiquetas:

Comentarios

Publicar un comentario

Entradas populares de este blog

Problema de probabilidad 1

A un mono se le dan 9 bloques: 3 en forma de cuadrados ,3 como rectángulos y 3 como triángulos. si saca tres de cada clase en orden, es decir, tres triángulos, luego la misma cantidad de cuadrados y así sucesivamente. ¿ sospecharía usted que el mono asoció figuras que tengan forma idéntica?. calcule la probabilidad de este evento. Como tenemos 3 figuras, debemos tener en cuenta todas las posibles formas de organizarlas, Luego una forma de representar este orden es de la siguiente manera: Tomamos $C$ como Cuadrado, $R$ como Rectángulo, y $T$ como Triángulo \[CRT \quad RTC \quad TCR\] \[CTR \quad RCT \quad TRC\] Que corresponden a $6$ formas de organizar las tres figuras, este procedimiento se puede representar fácilmente con la operación factorial, para nuestro caso queda definido como: \[3!=3 \times 2 \times 1 = 6\] Luego por cada vez que el mono saca tres bloques, estos vendrán organizados de $3!$ formas, y como son 3 tandas de veces que el mono va a sacar cantidades
Publicar un comentario
Leer más

¿Cómo resolver la ecuación diferencial de un oscilador armónico simple con gravedad?

En este caso vamos a solucionar la ecuación diferencial para un oscilador armónico unidimensional que solo tiene el movimiento en el eje y, y además está en presencia de la gravedad, su ecuación diferencial es la siguiente: \[\ddot{y}+\omega^{2}y=g\] Esta es una ecuación diferencial no homogénea y la vamos a resolver por el método de variación de parámetros. Para este método necesitamos resolver la ecuación diferencial homogenea: \[\ddot{y}+\omega^{2}y=0\] Que ya la hemos resuelto en este enlace  aunque resuelta con la coordenada $x$, pero eso no va a importar, porque va a ser la misma solución, solo que cambiamos de $x$ a $y$: \[y(t)=c_{1}cos(\omega t)+c_{2}sen(\omega t)\] Para resolver la ecuación diferencial debemos hallar el wronskiano: \[W(f(t),g(t))=\begin{vmatrix}f(t) & g(t)\\f'(t) & g'(t)\end{vmatrix}\] Identificamos las funciones $f(t)=c_{1}cos(\omega t)$ y $g(t)=c_{2}sen(\omega t)$, hallamos las derivadas $f'(t)=-c_{1}\omega sen(\omega t)$ y $g'(t)=c_{
Publicar un comentario
Leer más

¿Cómo realizar una sustitución para integrar una función en términos de arcoseno?

La integral que queremos resolver es la siguiente: \[\sqrt{\frac{m}{2E}}\int \frac{dx}{\sqrt{1-\frac{1}{2}\frac{k}{E}x^{2}}}\] Realizamos un cambio de variable de la siguiente forma: \[\frac{k}{2E}x^{2}=y^{2}\] Para poder hallar la derivada sacamos raíz cuadrada a ambos términos: \[\sqrt{\frac{k}{2E}}x=y\] Sacamos las derivadas y obtenemos: \[\sqrt{\frac{k}{2E}}dx=dy\] Luego pasando el factor que multiplica el $dx$ al lado del dy: \[dx=\sqrt{\frac{2E}{k}}dy\] Por tanto la integral queda de la forma: \[\sqrt{\frac{m}{2E}}\sqrt{\frac{2E}{k}}\int\frac{dy}{\sqrt{1-y^{2}}}\] \[\sqrt{\frac{m}{k}}\int\frac{dy}{\sqrt{1-y^{2}}}\] Que corresponde a una integral de arcoseno, por lo tanto la respuesta nos queda como: \[\sqrt{\frac{m}{k}}\int\frac{dy}{\sqrt{1-y^{2}}}=\sqrt{\frac{m}{k}}arcsen(y)+C\] Deshaciendo el cambio de variable, obtenemos finalmente el resultado de nuestra integral: \[\sqrt{\frac{m}{2E}}\int \frac{dx}{\sqrt{1-\frac{1}{2}\frac{k}{E}x^{2}}}=\sqrt{\frac{m}{k}}arcse
Publicar un comentario
Leer más