Las matemáticas son un lugar donde puedes hacer
cosas que no puedes hacer en el mundo real.
– Marcus du Sautoy

Introducción

Con lo que hemos estudiado en las dos últimas entradas somos capaces de resolver ecuaciones diferenciales lineales de segundo orden homogéneas y no homogéneas con coeficientes constantes, es decir, ecuaciones de la forma

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = g(x) \label{1} \tag{1}$$

Con $a, b, c$ constantes y $g(x) = 0$ en el caso homogéneo o $g(x) \neq 0$ en el caso no homogéneo, en éste último caso aún estamos limitados a la forma que puede tener la función $g$, pues sabemos resolver las ecuaciones diferenciales en el caso en el que la función $g$ es una constante, una función polinomial, una función exponencial, funciones seno o coseno, o una combinación entre ellas. La pregunta ahora es, ¿cómo resolver este tipo de ecuaciones para cualquier tipo de función $g(x)$?.

En esta entrada desarrollaremos un método que nos permite obtener la solución general independientemente de la forma que tenga la función $g(x)$. A dicho método se le conoce como variación de parámetros.

El nombre de este método resulta familiar. En la unidad anterior desarrollamos éste método para el caso de las ecuaciones diferenciales lineales de primer orden como método alterno al método por factor integrante. Lo que haremos en esta entrada es una adaptación del método de variación de parámetros para el caso en el que las ecuaciones diferenciales son de orden superior, en particular, de segundo orden.

Variación de parámetros

Consideremos la ecuación diferencial

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

Si $a_{2}(x) \neq 0$ para toda $x$ en el intervalo $\delta$ en el que está definida la solución, entonces podemos definir las funciones

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

de manera que la ecuación (\ref{2}) la podemos escribir en su forma estándar como

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

En el caso de primer orden se hizo la suposición de que la solución particular era de la forma

$$y_{p}(x) = k(x)y_{1}(x) = k(x) e^{-\int{P(x)dx}}$$

Manteniendo esta idea, en el caso de segundo orden se busca una solución de la forma

$$y_{p}(x) = k_{1}(x)y_{1}(x) + k_{2}(x)y_{2}(x) \label{4} \tag{4}$$

Donde $y_{1}$ y $y_{2}$ forman un conjunto fundamental de soluciones en $\delta$ de la ecuación homogénea asociada a (\ref{3}). Determinemos la primera y segunda derivada de $y_{p}$ para sustituir los resultados en la ecuación diferencial (\ref{3}).

$$\dfrac{dy_{p}}{dx} = k_{1}\dfrac{dy_{1}}{dx} + y_{1}\dfrac{dk_{1}}{dx} + k_{2}\dfrac{dy_{2}}{dx} + y_{2}\dfrac{dk_{2}}{dx} \label{5} \tag{5}$$

y

$$\dfrac{d^{2}y_{p}}{dx^{2}} = k_{1}\dfrac{d^{2}y_{1}}{dx^{2}} + 2\dfrac{dy_{1}}{dx}\dfrac{dk_{1}}{dx} + y_{1}\dfrac{d^{2}k_{1}}{dx^{2}} + k_{2}\dfrac{d^{2}y_{2}}{dx^{2}} + 2\dfrac{dy_{2}}{dx}\dfrac{dk_{2}}{dx} + y_{2}\dfrac{d^{2}k_{2}}{dx^{2}} \label{6} \tag{6}$$

Sustituyendo en (\ref{3}) y reorganizando los términos obtenemos lo siguiente.

\begin{align*}
k_{1} \left[ \dfrac{d^{2}y_{1}}{dx^{2}} + P\dfrac{dy_{1}}{dx} + Q y_{1} \right] + y_{1} \dfrac{d^{2}k_{1}}{dx^{2}} + \dfrac{dk_{1}}{dx} \dfrac{dy_{1}}{dx} + k_{2} \left[ \dfrac{d^{2}y_{2}}{dx^{2}} + P\dfrac{dy_{2}}{dx} + Q y_{2} \right] \\
+ y_{2} \dfrac{d^{2}k_{2}}{dx^{2}} + \dfrac{dk_{2}}{dx} \dfrac{dy_{2}}{dx} + P \left[ y_{1}\dfrac{dk_{1}}{dx} + y_{2}\dfrac{dk_{2}}{dx} \right] + \dfrac{dy_{1}}{dx} \dfrac{dk_{1}}{dx} + \dfrac{dy_{2}}{dx} \dfrac{dk_{2}}{dx} = g(x)
\end{align*}

Como $y_{1}$ y $y_{2}$ son soluciones de la ecuación homogénea asociada, entonces

$$\dfrac{d^{2}y_{1}}{dx^{2}} + P\dfrac{dy_{1}}{dx} + Q y_{1} = 0 \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y_{2}}{dx^{2}} + P\dfrac{dy_{2}}{dx} + Q y_{2} = 0$$

y además notamos que

$$\dfrac{d}{dx} \left[ y_{1}\dfrac{dk_{1}}{dx} \right] = y_{1} \dfrac{d^{2}k_{1}}{dx^{2}} + \dfrac{dk_{1}}{dx} \dfrac{dy_{1}}{dx} \label{7} \tag{7}$$

y

$$\dfrac{d}{dx} \left[ y_{2}\dfrac{dk_{2}}{dx} \right] = y_{2} \dfrac{d^{2}k_{2}}{dx^{2}} + \dfrac{dk_{2}}{dx} \dfrac{dy_{2}}{dx} \label{8} \tag{8}$$

Considerando lo anterior la ecuación diferencial queda como

$$\dfrac{d}{dx} \left[ y_{1}\dfrac{dk_{1}}{dx} + y_{2}\dfrac{dk_{2}}{dx} \right] + P \left[ y_{1}\dfrac{dk_{1}}{dx} + y_{2}\dfrac{dk_{2}}{dx} \right] + \dfrac{dy_{1}}{dx} \dfrac{dk_{1}}{dx} + \dfrac{dy_{2}}{dx} \dfrac{dk_{2}}{dx} = g(x) \label{9} \tag{9}$$

Nuestro propósito es determinar a las funciones $k_{1}(x)$ y $k_{2}(x)$ de (\ref{4}), esto implica que debemos formar un sistema con dos ecuaciones que debemos resolver para obtener dichas funciones. De acuerdo al resultado obtenido vamos a establecer la restricción de que las funciones $k_{1}$ y $k_{2}$ satisfacen la relación

$$y_{1}\dfrac{dk_{1}}{dx} + y_{2}\dfrac{dk_{2}}{dx} = 0 \label{10} \tag{10}$$

Considerando esto la ecuación se reduce a

$$\dfrac{dy_{1}}{dx} \dfrac{dk_{1}}{dx} + \dfrac{dy_{2}}{dx} \dfrac{dk_{2}}{dx} = g(x) \label{11} \tag{11}$$

Las ecuaciones (\ref{10}) y (\ref{11}) corresponden al sistema de dos ecuaciones que debemos resolver.

Como podemos notar, es un sistema para determinar las derivadas de las funciones $k_{1}$ y $k_{2}$ y no las funciones mismas, esto implica que una vez que determinemos a las derivadas será necesario hacer una integración a cada una de ellas. Resolvamos el sistema.

Multipliquemos la ecuación (\ref{10}) por $\dfrac{dy_{2}}{dx}$ y la ecuación (\ref{11}) por $y_{2}$.

$$y_{1}\dfrac{dk_{1}}{dx} \dfrac{dy_{2}}{dx} + y_{2}\dfrac{dk_{2}}{dx} \dfrac{dy_{2}}{dx} = 0 \label{12} \tag{12}$$

$$y_{2} \dfrac{dy_{1}}{dx} \dfrac{dk_{1}}{dx} + y_{2} \dfrac{dy_{2}}{dx} \dfrac{dk_{2}}{dx} = y_{2} g(x) \label{13} \tag{13}$$

Si a la ecuación (\ref{12}) le restamos la ecuación (\ref{13}) obtenemos lo siguiente.

$$\dfrac{dk_{1}}{dx} \left( y_{1} \dfrac{dy_{2}}{dx} -y_{2} \dfrac{dy_{1}}{dx} \right) = -y_{2}(x)g(x) \label{14} \tag{14}$$

Recordemos que el Wronskiano esta definido como

$$W(y_{1}, y_{2}) = y_{1} \dfrac{dy_{2}}{dx} -y_{2} \dfrac{dy_{1}}{dx} \label{15} \tag{15}$$

Entonces la ecuación (\ref{14}) la podemos escribir como

$$\dfrac{dk_{1}}{dx} \left[ W(y_{1}, y_{2}) \right] = -y_{2}(x)g(x) \label{16} \tag{16}$$

Como $y_{1}$ y $y_{2}$ forman un conjunto fundamental de soluciones de la ecuación homogénea asociada, entonces

$$W(y_{1}, y_{2}) \neq 0$$

Así, de la ecuación (\ref{16}) obtenemos que

$$\dfrac{dk_{1}}{dx} = -\dfrac{y_{2}(x)g(x)}{W(y_{1}, y_{2})} \label{17} \tag{17}$$

Hemos encontrado el valor de la derivada de la función $k_{1}(x)$, integrando obtenemos finalmente que

$$k_{1}(x) = -\int{ \dfrac{y_{2}(x)g(x)}{W(y_{1}, y_{2})} dx} \label{18} \tag{18}$$

En un proceso totalmente análogo, si multiplicamos a la ecuación (\ref{10}) por $\dfrac{dy_{1}}{dx}$ y a la ecuación (\ref{11}) por $y_{1}$ y realizamos los mismos pasos obtendremos la ecuación para la derivada de la función $k_{2}(x)$.

$$\dfrac{dk_{2}}{dx} = \dfrac{y_{1}(x)g(x)}{W(y_{1}, y_{2})} \label{19} \tag{19}$$

Integrando obtendremos la función que buscamos

$$k_{2}(x) = \int{ \dfrac{y_{1}(x)g(x)}{W(y_{1}, y_{2})} dx} \label{20} \tag{20}$$

Sustituyendo los resultados (\ref{18}) y (\ref{20}) en la solución particular (\ref{4}) obtenemos finalmente la solución que buscábamos

$$y_{p}(x) = -y_{1}(x) \int{ \dfrac{y_{2}(x)g(x)}{W(y_{1}, y_{2})} dx} + y_{2}(x) \int{ \dfrac{y_{1}(x)g(x)}{W(y_{1}, y_{2})} dx} \label{21} \tag{21}$$

El resultado (\ref{21}) corresponde a la solución particular de la ecuación diferencial (\ref{2}) que, a diferencial del método de coeficientes indeterminados, se aplica para cualquier función $g(x)$, aunque cabe mencionar que si la función $g$ es muy compleja, entonces nos resultará, en algunas ocasiones, complicado resolver las integrales involucradas.

A lo largo del curso hemos motivado a no memorizar las formulas y en su lugar desarrollar el procedimiento del método correspondiente, sin embargo, en esta ocasión se trata de un método muy largo y complicado para usarse cada vez que se intente resolver una ecuación diferencial, por lo que se recomienda seguir los siguientes pasos.

• Primero se determina la solución complementaria $$y_{c} = c_{1}y_{1} + c_{2}y_{2}$$ de la ecuación diferencial homogénea asociada, esto nos permitirá determinar el conjunto fundamental de soluciones $\{y_{1}, y_{2}\}$.

• Una vez conocido el conjunto fundamental de soluciones se procede a calcular el Wronskiano $W(y_{1}, y_{2})$.

• Posteriormente se divide la ecuación diferencial por $a_{2}$ para escribir la ecuación es su forma estándar (\ref{3}) y así obtener la forma de la función $g(x)$.

• Se sustituyen los valores correspondientes en (\ref{18}) y (\ref{20}) para obtener las funciones $k_{1}$ y $k_{2}$ respectivamente.

• Finalmente se sustituyen los resultados en la solución particular $$y_{p} = k_{1}y_{1} + k_{2}y_{2}$$ y posteriormente en la solución general $$y = y_{c} + y_{p}$$

Cuando se calculan las integrales indefinidas (\ref{18}) y (\ref{20}) no es necesario considerar las constantes de integración. Para mostrar esto consideremos las constantes $c_{3}$ y $c_{4}$, tales que

\begin{align*}
y(x) &= y_{c}(x) + y_{p}(x) \\
&= c_{1}y_{1}(x) + c_{2}y_{2}(x) + \left[ k_{1}(x) + c_{3} \right] y_{1}(x) + \left[ k_{2}(x) + c_{4} \right] y_{2}(x) \\
&= \left[ c_{1} + c_{3} \right] y_{1}(x) + \left[ c_{2} + c_{4} \right] y_{2}(x) + k_{1}(x)y_{1}(x) + k_{2}(x)y_{2}(x) \\
&= C_{1}y_{1}(x) + C_{2}y_{2}(x) + k_{1}(x)y_{1}(x) + k_{2}(x)y_{2}(x)
\end{align*}

Es decir, las constantes de la solución complementaria contienen todas las constantes que puedan aparecer en el método.

Realicemos algunos ejemplos.

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

$$3 \dfrac{d^{2}y}{dx^{2}} -6 \dfrac{dy}{dx} + 6y = e^{x} \sec(x)$$

Solución: El primer paso es obtener la solución complementaria. La ecuación auxiliar es

$$3k^{2} -6k + 6 = 0$$

De donde $k_{1} = 1 + i$ y $k_{2} = 1 -i$, identificamos que $\alpha = \beta = 1$, entonces la forma de la solución complementaria es

$$y_{c}(x) = c_{1}e^{x} \cos(x) + c_{2}e^{x} \sin(x)$$

El conjunto fundamental de soluciones esta conformado por las funciones

$$y_{1}(x) = e^{x} \cos(x) \hspace{1cm} y \hspace{1cm} y_{2}(x) = e^{x} \sin(x)$$

La derivada de ambas soluciones son

$$\dfrac{dy_{1}}{dx} = e^{x} \cos(x) -e^{x} \sin(x) \hspace{1cm} y \hspace{1cm} \dfrac{dy_{2}}{dx} = e^{x} \sin(x) + e^{x} \cos(x)$$

Usando estos resultados calculemos el Wronskiano.

\begin{align*}
W &= \begin{vmatrix}
e^{x} \cos(x) & e^{x} \sin(x) \\
e^{x} \cos(x) -e^{x} \sin(x) & e^{x} \sin(x) + e^{x} \cos(x)
\end{vmatrix} \\
&= e^{2x} \cos(x) \sin(x) + e^{2x} \cos^{2}(x) -e^{2x} \cos(x) \sin(x) + e^{2x} \sin^{2}(x) \\
&= e^{2x}
\end{align*}

El Wronskiano es

$$W(x) = e^{2x}$$

¡Cuidado!, como en la ecuación diferencial la segunda derivada tiene un coeficiente, debemos dividir toda la ecuación por dicho coeficiente para obtener la forma estándar y así la función $g(x)$. La ecuación en su forma estándar es

$$\dfrac{d^{2}y}{dx^{2}} -2 \dfrac{dy}{dx} + 2y = \dfrac{e^{x} \sec(x)}{3}$$

En este caso la función $g$ es

$$g(x) = \dfrac{e^{x} \sec(x)}{3}$$

Ahora que ya conocemos los valores que necesitábamos, recurrimos a las ecuaciones (\ref{18}) y (\ref{20}) para obtener las funciones que buscamos.

Para la función $k_{1}(x)$, tenemos lo siguiente.

\begin{align*}
k_{1}(x) &= -\int{ \dfrac{y_{2}(x)g(x)}{W(y_{1}, y_{2})} dx} \\
&= -\int{\dfrac{\left( e^{x} \sin(x) \right) \left( \dfrac{e^{x} \sec(x)}{3} \right)}{e^{2x}} dx} \\
&= -\int{\dfrac{e^{2x} \sin(x) \sec(x)}{3e^{2x}} dx} \\
&= -\dfrac{1}{3} \int{\tan(x) dx} \\
&= \dfrac{1}{3} \ln|\cos(x)|
\end{align*}

La integral de la tangente es común. Por tanto, la función $k_{1}$ es

$$k_{1}(x) = \dfrac{1}{3} \ln|\cos(x)|$$

Para el caso de la función $k_{2}(x)$, tenemos lo siguiente.

\begin{align*}
k_{2}(x) &= \int{ \dfrac{y_{1}(x)g(x)}{W(y_{1}, y_{2})} dx} \\
&= \int{\dfrac{\left( e^{x} \cos(x) \right) \left( \dfrac{e^{x} \sec(x)}{3} \right)}{e^{2x}} dx} \\
&= \int{\dfrac{e^{2x} \cos(x) \sec(x)}{3e^{2x}} dx} \\
&= \dfrac{1}{3} \int{dx} \\
&= \dfrac{1}{3}x
\end{align*}

La función $k_{2}$ es

$$k_{2}(x) = \dfrac{1}{3}x$$

Ya podemos establecer que la solución particular, de acuerdo a (\ref{4}), es

$$y_{p}(x) = \dfrac{1}{3} \ln|\cos(x)| \left[ e^{x} \cos(x) \right] + \dfrac{1}{3}x \left[ e^{x} \sin(x) \right]$$

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

$$3 \dfrac{d^{2}y}{dx^{2}} -6 \dfrac{dy}{dx} + 6y = e^{x} \sec(x)$$

es

$$y(x) = c_{1}e^{x} \cos(x) + c_{2}e^{x} \sin(x) + \dfrac{1}{3}e^{x} \cos(x) \ln|\cos(x)| + \dfrac{1}{3}x e^{x} \sin(x)$$

$\square$

Con este ejemplo encontramos un buen momento para reflexionar y darnos cuenta de que ya hemos avanzado mucho, tan sólo observa el tipo de ecuación que acabamos de resolver y no sólo eso, observa que tan compleja es la solución general.

¡Sigamos adelante!.

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

$$\dfrac{d^{2}y}{dx^{2}} + 2\dfrac{dy}{dx} + y = \dfrac{e^{-x}}{x}$$

Solución: Como la ecuación ya está es su forma estándar la función $g$ es

$$g(x) = \dfrac{e^{-x}}{x}$$

Determinemos la solución complementaria, la ecuación auxiliar es

$$k^{2} + 2k + 1 = 0$$

De donde $k_{1} = k_{2} = -1$, la multiplicidad de la solución nos indica que la forma de la solución complementaria es

$$y_{c}(x) = c_{1}e^{-x} + c_{2}xe^{-x}$$

El conjunto fundamental de soluciones esta conformado por las funciones

$$y_{1}(x) = e^{-x} \hspace{1cm} y \hspace{1cm} y_{2}(x) = xe^{-x}$$

Usando estas soluciones y sus derivadas calculemos el Wronskiano.

\begin{align*}
W &= \begin{vmatrix}
e^{-x} & xe^{-x} \\
-e^{-x} & e^{-x} -xe^{-x}
\end{vmatrix} \\
&= e^{-2x} -xe^{2x} + xe^{-2x} \\
&= e^{-2x}
\end{align*}

El Wronskiano es

$$W(x) = e^{-2x}$$

Sustituyamos estos resultados directamente en la ecuación (\ref{21}).

\begin{align*}
y_{p}(x) &= -e^{-x} \int {\dfrac{ \left( xe^{-x} \right) \left( \dfrac{e^{-x}}{x} \right) }{e^{-2x}}dx} + xe^{-x} \int {\dfrac{ \left( e^{-x} \right) \left( \dfrac{e^{-x}}{x} \right) }{e^{-2x}}dx} \\
&= -e^{-x}\int dx+xe^{-x}\int \dfrac{dx}{x} \\
&= -xe^{-x} + xe^{-x}\ln(x)
\end{align*}

La solución particular es

$$y_{p}(x) = -xe^{-x} + xe^{-x}\ln(x)$$

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

$$y(x) = c_{1}e^{-x} + c_{2}xe^{-x} -xe^{-x} + xe^{-x}\ln(x)$$

$\square$

Un ejemplo más.

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

$$\dfrac{d^{2}y}{dx^{2}} -y = 4x^{3}e^{x}$$

Solución: La función $g$ es

$$g(x) = 4x^{3}e^{x}$$

y la ecuación auxiliar es

$$k^{2} -1 = 0$$

De donde $k_{1} = 1$ y $k_{2} = -1$. Entonces, la solución complementaria es

$$y_{c}(x) = c_{1}e^{x} + c_{2}e^{-x}$$

El conjunto fundamental de soluciones esta conformado por

$$y_{1}(x) = e^{x} \hspace{1cm} y \hspace{1cm} y_{2}(x) = e^{-x}$$

Usando estas funciones y sus derivadas calculemos el Wronskiano.

\begin{align*}
W &= \begin{vmatrix}
e^{x} & e^{-x} \\
e^{x} & -e^{-x} \end{vmatrix} = -2
\end{align*}

El Wronskiano es $W = -2$. Sustituyendo estos resultados directamente en la ecuación (\ref{21}), obtenemos lo siguiente.

\begin{align*}
y_{p}(x) &= -e^{x} \int {\dfrac{(e^{-x})(4x^{3}e^{x})}{-2} dx} + e^{-x} \int {\dfrac{(e^{x})(4x^{3}e^{x})}{-2} dx} \\
&= 2e^{x} \int {x^{3} dx} -2e^{-x} \int {x^{3}e^{2x} dx} \\
&= \dfrac{1}{2}x^{4}e^{x} -2e^{-x} \int {x^{3}e^{2x} dx}
\end{align*}

La integral que nos falta se puede resolver por partes tomando $u = x^{2}$ y $v^{\prime} = e^{2x}$. Resolviendo la integral obtendremos lo siguiente.

$$\int {x^{3}e^{2x} dx} = \dfrac{1}{2}e^{2x}x^{3} -\dfrac{3}{4}e^{2x}x^{2} + \dfrac{3}{4}e^{2x}x -\dfrac{3}{8}e^{2x}$$

Sustituyendo en la solución particular tenemos

\begin{align*}
y_{p}(x) &= \dfrac{1}{2}x^{4}e^{x} -2e^{-x} \left( \dfrac{1}{2}e^{2x}x^{3} -\dfrac{3}{4}e^{2x}x^{2} + \dfrac{3}{4}e^{2x}x -\dfrac{3}{8}e^{2x} \right) \\
&= \dfrac{1}{2}x^{4}e^{x} -x^{3}e^{x} + \dfrac{3}{2}x^{2}e^{x} -\dfrac{3}{2}xe^{x} + \dfrac{3}{4}e^{x}
\end{align*}

Finalmente obtenemos como solución particular a la función

$$y_{p}(x) = e^{x} \left( \dfrac{1}{2}x^{4} -x^{3} + \dfrac{3}{2}x^{2} -\dfrac{3}{2}x + \dfrac{3}{4} \right)$$

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

$$y(x) = c_{1}e^{x} + c_{2}e^{-x} + e^{x} \left( \dfrac{1}{2}x^{4} -x^{3} + \dfrac{3}{2}x^{2} -\dfrac{3}{2}x + \dfrac{3}{4} \right)$$

Este resultado es válido, sin embargo se puede simplificar más, ya que se puede reescribir a la solución como

$$y(x) = e^{x} \left( c_{1} + \dfrac{3}{4} \right) + c_{2}e^{-x} + e^{x} \left( \dfrac{1}{2}x^{4} -x^{3} + \dfrac{3}{2}x^{2} -\dfrac{3}{2}x \right)$$

y definir la constante $C_{1} = c_{1} + \dfrac{3}{4}$ para finalmente escribir la solución como

$$y(x) = C_{1} e^{x} + c_{2}e^{-x} + e^{x} \left( \dfrac{1}{2}x^{4} -x^{3} + \dfrac{3}{2}x^{2} -\dfrac{3}{2}x \right)$$

$\square$

Variación de parámetros en ecuaciones de orden superior

Este método se puede generalizar a ecuaciones de orden $n$ aunque, por su puesto, los cálculos se vuelven más extensos.

A continuación mostraremos el panorama general para ecuaciones diferenciales de orden $n$ y mostraremos los resultados para el caso $n = 3$ que nos mostrará la forma en que aumenta la complejidad de los cálculos.

La ecuación de orden $n$ es su forma estándar es

$$\dfrac{d^{n}y}{dx^{n}} + P_{n -1}(x) \dfrac{d^{n -1}y}{dx^{n -1}} + \cdots + P_{1}(x) \dfrac{dy}{dx} + P_{0}(x) y = g(x) \label{22} \tag{22}$$

Si la solución complementaria de (\ref{22}) es

$$y_{c}(x) = c_{1}y_{1}(x) + c_{2}y_{2}(x) + \cdots + c_{n}y_{n}(x) \label{23} \tag{23}$$

Entonces una solución particular debe ser

$$y_{p}(x) = k_{1}(x)y_{1}(x) + k_{2}(x)y_{2}(x) + \cdots + k_{n}(x)y_{n}(x) \label{24} \tag{24}$$

Análogo a las ecuaciones (\ref{10}) y (\ref{11}), las derivadas $\dfrac{dk_{i}}{dx} = k^{\prime}_{i}$ con $i = 1, 2, \cdots, n$ se determinan del sistema de $n$ ecuaciones

\begin{align*}
y_{1}k^{\prime}_{1} + y_{2}k^{\prime}_{2} + \cdots + y_{n}k^{\prime}_{n} &= 0 \\
y^{\prime}_{1}k^{\prime}_{1} + y^{\prime}_{2}k^{\prime}_{2} + \cdots + y^{\prime}_{n}k^{\prime}_{n} &= 0 \\
\vdots \\
y^{(n -1)}_{1}k^{\prime}_{1} + y^{(n -1)}_{2}k^{\prime}_{2} + \cdots + y^{(n-1)}_{n}k^{\prime}_{n} &= g(x) \label{25} \tag{25}
\end{align*}

Al igual que el caso de segundo orden, las primeras $n -1$ ecuaciones del sistema son suposiciones que se hacen para simplificar la ecuación resultante de sustituir la solución (\ref{24}) en la ecuación (\ref{22}).

Usando la regla de Cramer para resolver el sistema se obtiene que

$$\dfrac{dk_{i}}{dx} = \dfrac{W_{i}}{W}; \hspace{1cm} i = 1, 2 , \cdots, n \label{26} \tag{26}$$

Donde $W$ es el Wronskiano del conjunto fundamental $\{ y_{1}(x), y_{2}(x), \cdots, y_{n}(x) \}$ y $W_{i}$ es el determinante que se obtiene de remplazar la $i$-ésima columna del Wronskiano por la columna formada por el lado derecho de (\ref{25}), es decir, la columna que consta de $(0, 0, \cdots, g(x))$.

Para que quede más claro lo anterior, en el caso $n = 3$ las $\dfrac{dk_{i}}{dx}$, $i = 1, 2, 3$ quedan como

$$\dfrac{dk_{1}}{dx} = \dfrac{W_{1}}{W}, \hspace{1cm} \dfrac{dk_{2}}{dx} = \dfrac{W_{2}}{W}, \hspace{1cm} \dfrac{dk_{3}}{dx} = \dfrac{W_{3}}{W} \label{27} \tag{27}$$

Donde

$$W = \begin{vmatrix}
y_{1} & y_{2} & y_{3} \\
y^{\prime}_{1} & y^{\prime}_{2} & y^{\prime}_{3} \\
y^{\prime \prime}_{1} & y^{\prime \prime}_{2} & y^{\prime \prime}_{3}
\end{vmatrix} \label{28} \tag{28}$$

y

\begin{align*}
W_{1} = \begin{vmatrix}
0 & y_{2} & y_{3} \\
0 & y^{\prime}_{2} & y^{\prime}_{3} \\
g(x) & y^{\prime \prime}_{2} & y^{\prime \prime}_{3}
\end{vmatrix}, \hspace{1cm} W_{2} = \begin{vmatrix}
y_{1} & 0 & y_{3} \\
y^{\prime}_{1} & 0 & y^{\prime}_{3} \\
y^{\prime \prime}_{1} & g(x) & y^{\prime \prime}_{3}
\end{vmatrix}, \hspace{1cm} W_{3} = \begin{vmatrix}
y_{1} & y_{2} & 0 \\
y^{\prime}_{1} & y^{\prime}_{2} & 0 \\
y^{\prime \prime}_{1} & y^{\prime \prime}_{2} & g(x)
\end{vmatrix}
\end{align*}

Habrá que integrar las ecuaciones de (\ref{27}) para obtener las funciones $k_{i}$, $i = 1, 2, 3$ y así obtener la solución particular

$$y_{p}(x) = k_{1}(x)y_{1}(x) + k_{2}(x)y_{2}(x) + k_{3}(x)y_{3}(x) \label{29} \tag{29}$$

Notemos que usando esta notación, los resultados (\ref{17}) y (\ref{19}) del caso $n = 2$ se pueden escribir como

$$\dfrac{dk_{1}}{dx} = \dfrac{W_{1}}{W} = -\dfrac{y_{2}g(x)}{W} \hspace{1cm} y \hspace{1cm} \dfrac{dk_{2}}{dx} = \dfrac{W_{2}}{W} = \dfrac{y_{1}g(x)}{W} \label{30} \tag{30}$$

Donde

\begin{align*} W = \begin{vmatrix}
y_{1} & y_{2} \\
y^{\prime}_{1} & y^{\prime}_{2}
\end{vmatrix}, \hspace{1cm} W_{1} = \begin{vmatrix}
0 & y_{2} \\
g(x) & y^{\prime}_{2}
\end{vmatrix} \hspace{1cm} y \hspace{1cm} W_{2} = \begin{vmatrix}
y_{1} & 0 \\
y^{\prime}_{1} & g(x)
\end{vmatrix} \label{31} \tag{31}
\end{align*}

Realicemos un ejemplo con una ecuación de orden $3$.

Ejemplo: Resolver la ecuación diferencial de tercer orden

$$\dfrac{d^{3}y}{dx^{3}} + \dfrac{dy}{dx} = \tan(x)$$

Solución: La función $g$ es

$$g(x) = \tan(x)$$

y la ecuación auxiliar es

$$k^{3} + k = k(k^{2} + 1) = 0$$

De donde $k_{1} = 0$, $k_{2} = i$ y $k_{3} = -i$. Dos raíces son complejas conjugadas con $\alpha = 0$ y $\beta = 1$. La primer raíz nos indica que la forma de una solución es

$$y_{1}(x) = e^{k_{1}x} = 1$$

mientras que las dos raíces restantes nos indican dos solución de la forma

$$y_{2}(x) = e^{\alpha x} \cos(\beta x) = \cos(x)$$

y

$$y_{3}(x) = e^{\alpha x} \sin(\beta x) = \sin(x)$$

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

$$y_{c}(x) = c_{1} + c_{2} \cos(x) + c_{3} \sin(x)$$

Como vimos, el conjunto fundamental de soluciones es $\{ y_{1} = 1, y_{2} = \cos(x) , y_{3} = \sin(x) \}$, las primeras y segundas derivadas correspondientes son

$$\dfrac{dy_{1}}{dx} = 0, \hspace{1cm} \dfrac{dy_{2}}{dx} = -\sin(x), \hspace{1cm} \dfrac{dy_{3}}{dx} = \cos(x)$$

y

$$\dfrac{d^{2}y_{1}}{dx^{2}} = 0, \hspace{1cm} \dfrac{d^{2}y_{2}}{dx^{2}} = -\cos(x), \hspace{1cm} \dfrac{d^{2}y_{3}}{dx^{2}} = -\sin(x)$$

Ahora calculemos los determinantes correspondientes, el primero de ellos es el Wronskiano

$$W = \begin{vmatrix}
1 & \cos(x) & \sin(x) \\
0 & -\sin(x) & \cos(x) \\
0 & -\cos(x) & -\sin(x)
\end{vmatrix} = \sin^{2}(x) + \cos^{2}(x) = 1$$

Para el resto de determinantes, tenemos

$$W_{1} = \begin{vmatrix}
0 & \cos(x) & \sin(x) \\
0 & -\sin(x) & \cos(x) \\
\tan(x) & -\cos(x) & -\sin(x)
\end{vmatrix} = \tan(x) \left[ \cos^{2}(x) + \sin^{2}(x) \right] = \tan(x)$$

$$W_{2} = \begin{vmatrix}
1 & 0 & \sin(x) \\
0 & 0 & \cos(x) \\
0 & \tan(x) & -\sin(x)
\end{vmatrix} = -\cos(x) \tan(x) = -\sin(x)$$

y

$$W_{3} = \begin{vmatrix}
1 & \cos(x) & 0 \\
0 & -\sin(x) & 0 \\
0 & -\cos(x) & \tan(x)
\end{vmatrix} = -\sin(x) \tan(x) = -\dfrac{\sin^{2}(x)}{\cos(x)}$$

Sustituyendo estos resultados en (\ref{27}), obtenemos

$$\dfrac{dk_{1}}{dx} = \tan(x), \hspace{1cm} \dfrac{dk_{2}}{dx} = -\sin(x), \hspace{1cm} \dfrac{dk_{3}}{dx} = -\dfrac{\sin^{2}(x)}{\cos(x)}$$

Procedemos a integrar cada ecuación (sin considerar constantes) para obtener las funciones que buscamos.

La primer integral es común,

$$k_{1}(x) = \int{\tan(x) dx} = -\ln|\cos(x)|$$

La segunda integral es directa

$$k_{2}(x) = -\int{\sin(x) dx} = \cos(x)$$

Mientras que para la tercer integral si se requiere de un mayor cálculo.

\begin{align*}
k_{3}(x) &= -\int{\dfrac{\sin^{2}(x)}{\cos(x)} dx} \\
&= -\int{\dfrac{1 -\cos^{2}(x)}{\cos(x)} dx} \\
&= -\int{\dfrac{1}{\cos(x)} dx} + \int{\cos(x) dx} \\
\end{align*}

Por un lado

$$-\int{\dfrac{1}{\cos(x)} dx} = -\int{\sec(x) dx} = -\ln|\tan(x) + \sec(x)|$$

por otro lado,

$$\int{\cos(x) dx} = \sin(x)$$

entonces

$$k_{3}(x) = -\ln|\tan(x) + \sec(x)| + \sin(x)$$

Ahora que conocemos las funciones incógnita concluimos que la solución particular de la ecuación diferencial es

\begin{align*}
y_{p}(x) &= -\ln|\cos(x)|(1) + \cos(x) [\cos(x)] + \left[ -\ln|\tan(x) + \sec(x)| + \sin(x) \right] (\sin(x)) \\
&= -\ln|\cos(x)| + \cos^{2}(x) -\sin(x) \ln|\tan(x) + \sec(x)| + \sin^{2}(x) \\
&= -\ln|\cos(x)| -\sin(x) \ln|\tan(x) + \sec(x)| \\
\end{align*}

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

$$y(x) = c_{1} + c_{2} \cos(x) + c_{3} \sin(x) -\ln|\cos(x)| -\sin(x) \ln|\tan(x) + \sec(x)|$$

$\square$

Como podemos notar, los cálculos se hacen más extensos, sin embargo los pasos a seguir son los mismos para cualquier orden.

El método de variación de parámetros, a diferencia del método de coeficientes indeterminados, tiene la ventaja de que siempre produce una solución de la ecuación diferencial independientemente de la forma de la función $g(x)$, siempre y cuando se pueda resolver la ecuación homogénea asociada. Además, el método de variación de parámetros es aplicable a ecuaciones diferenciales lineales con coeficientes variables.

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. Obtener la solución general de las siguientes ecuaciones diferenciales.

• $\dfrac{d^{2}y}{dx^{2}} -9\dfrac{dy}{dx} = 18x^{2}e^{9x}$

• $\dfrac{d^{2}y}{dx^{2}} + 9y = 18e^{x} \sin(x)$

• $4 \dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + y = e^{x/2}\sqrt{1 -x^{2}}$

2. Resolver las siguientes ecuaciones diferenciales para las condiciones iniciales dadas.

• $\dfrac{d^{2}y}{dx^{2}} -3\dfrac{dy}{dx} = 12e^{4x}(x + 1); \hspace{1cm} y(0) = 0, \hspace{0.6cm} y^{\prime}(0) = 4$

• $\dfrac{d^{2}y}{dx^{2}} + \dfrac{dy}{dx} -6y = 10e^{x} \sin(x); \hspace{1cm} y(0) = \dfrac{2}{17}, \hspace{0.6cm} y^{\prime}(0) = 0$

3. Obtener la solución general de las siguientes ecuaciones diferenciales de tercer orden. Simplificar la forma de la solución redefiniendo las constantes.

• $\dfrac{d^{3}y}{dx^{3}} + 4 \dfrac{dy}{dx} = \sec(2x)$

• $\dfrac{d^{3}y}{dx^{3}} + \dfrac{d^{2}y}{dx^{2}} -\dfrac{dy}{dx} -y = 2e^{-x}$

Más adelante…

Hemos concluido con el estudio de las ecuaciones diferenciales lineales de orden superior con coeficientes constantes.

Lo que sigue es estudiar este mismo tipo de ecuaciones, pero en el caso en el que los coeficientes no son constantes, es decir, son coeficientes variables. Estas ecuaciones suelen ser mucho más difícil de resolver, sin embargo existe un tipo de ecuación especial, conocida como ecuación de Cauchy – Euler, que contiene coeficientes variables, pero que su método de resolución es bastante similar a lo que hemos desarrollado en el caso de coeficientes constantes, pues su resolución involucra resolver una ecuación auxiliar. En la siguiente entrada estudiaremos dicha ecuación.

Entradas relacionadas

• Página principal del curso: Ecuaciones Diferenciales I
• Entrada anterior del curso: Ecuaciones lineales no homogéneas de segundo orden – Método de coeficientes indeterminados
• Siguiente entrada del curso: Ecuación de Cauchy – Euler

• Video relacionado al tema: Ecuaciones lineales no homogéneas de segundo orden. Solución por variación de parámetros

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»