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}\]

Comentarios

Entradas populares de este blog

¿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_{...

Demostración ecuación diferencial lineal homogénea de segundo orden con Wronskiano

En esta ocasión, vengo a mostrarles la solución del problema 47 de ecuaciones diferenciales de segundo orden homogeneas del libro de  Ecuaciones diferenciales con aplicaciones 2da edición,  del autor  Dennis G. Zill. El problema es el siguiente: Sean $y_{1}$ y $y_{2}$ dos soluciones de: \[a_{2}(x)y''+a_{1}(x)y'+a_{0}y=0\] (a) Sí $W(y_{1},y_{2})$ es el wronskiano de $y_{1}$ y $y_{2}$, demuestre que \[a_{2}(x)\frac{dW}{dx}+a_{1}(x)W=0\] Como $y_{1}$ y $y_{2}$ son soluciones de la ecuación diferencial homogénea de segundo orden, luego también deben ser soluciones de las ecuaciones diferenciales: \[a_{2}(x)y_{1}''+a_{1}(x)y_{1}'+a_{0}y_{1}=0\] \[a_{2}(x)y_{2}''+a_{1}(x)y_{2}'+a_{0}y_{2}=0\] El wronskiano es de la forma: \[W(y_{1},y_{2})=\begin{vmatrix}y_{1} & y_{2}\\y_{1}' & y_{2}'\end{vmatrix}\] \[W(y_{1},y_{2})=y_{1}y_{2}'-y_{2}y_{1}'\] La derivada del Wronskiano es: \[\frac{dW}{dx}=(y_{1}y_{2}')'-(y_{2}y_{1}')'\...

¿Cómo hallar el centro de masa de una lámina con integrales dobles? 2

 A petición de un suscriptor y seguidor voy a mostrar dos ejercicios para el cálculo del centro de masa con integrales dobles, comencemos: El segundo ejercicio es: 2. Calcular el centro de masa de una lámina representada por la región $R$ que se encuentra por encima del eje $x$ y entre las líneas $y=x$; $y=-x$, $x^{2}+y^{2}=4y$; $x^{2}+y^{2}=6y$; $y>0$; donde la densidad es $\sqrt{x^{2}+y^{2}}$ Primero vamos a representar el área, al cuál vamos a hallar el centro de masa (voy a utilizar Geogebra). Este problema resulta más simple si lo resolvemos mediante coordenadas polares, así que la densidad $\rho(x,y)$ en coordenadas polares de acuerdo a las reglas de transformación $x=rcos\theta$ y $y=rsen\theta$ es: \[\rho(x,y)=\rho(rcos\theta,rsen\theta)=\sqrt{x^{2}+y^{2}}=\sqrt{(rcos\theta)^{2}+(rsen\theta)^{2}}=\sqrt{r^{2}}=r\] Como vimos en el primer ejercicio , las integrales que tiene cada punto del centro de masa se pueden representar en coordenadas polares, escogemos este tipo de ...