Las matemáticas puras son, en su forma, la poesía de las ideas lógicas.
– Albert Einstein

Introducción

Más adelante en esta unidad estudiaremos las ecuaciones diferenciales de orden superior con coeficientes variables, éstas ecuaciones suelen ser mucho más difícil de resolver ya que no se resuelven en términos de funciones elementales, una estrategia usual es suponer una solución en forma de series infinitas y proceder de manera similar al método de coeficientes indeterminados. Sin embargo, existe una ecuación diferencial de coeficientes variables que es una excepción, pues su solución general siempre se puede expresar en términos de potencias de $x$, senos, cosenos y funciones logarítmicas, dicha ecuación es conocida como ecuación de Cauchy – Euler y dedicaremos esta entrada a estudiarla, así como su método de resolución.

Decidimos estudiar esta ecuación en este momento debido a que el método de resolución es bastante similar al de las ecuaciones con coeficientes constantes en los que se debe resolver una ecuación auxiliar.

Ecuación de Cauchy – Euler

Enseguida nos damos cuenta de que los coeficientes

$$b_{n}(x) = a_{n}x^{n}, \hspace{0.5cm} b_{n -1}(x) = a_{n -1}x^{n -1}, \hspace{0.5cm} \cdots, \hspace{0.5cm} b_{1}(x) = a_{1}x^{1}, \hspace{0.5cm} b_{0}(x) = a_{0}x^{0}$$

son dependientes de $x$, es decir, son coeficientes variables, además la característica importante de esta ecuación es que el grado $k = n, n -1, \cdots, 1, 0$ de los coeficientes monomiales $x^{k}$ coincide con el orden $k$ de la derivación $\dfrac{d^{k}y}{dx^{k}}$.

Como se ha hecho a lo largo de la unidad, desarrollaremos con todo detalle el método de resolución de la ecuación de Cauchy – Euler para el caso de segundo orden, recordando que es posible extender el método a cualquier orden $n$ siguiendo el mismo razonamiento.

Iniciaremos nuestro análisis con un estudio detallado de las formas de las soluciones generales de la ecuación homogénea de segundo orden

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

con $a$, $b$ y $c$ constantes. Para resolver la ecuación no homogénea

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

con $g(x) \neq 0$, basta aplicar el método de variación de parámetros (o de coeficientes indeterminados) una vez que se ha determinado la función complementaria $y_{c}$, es decir, la solución general de la ecuación homogénea (\ref{2}).

Una consideración importante es que el coeficiente $ax^{2}$ de $\dfrac{d^{2}y}{dx^{2}}$ es cero en $x = 0$, para garantizar los resultados fundamentales del teorema de existencia y unicidad y sean aplicables a la ecuación de Cauchy – Euler debemos encontrar soluciones generales definidas en el intervalo $\delta = (0, \infty)$. Las soluciones en el intervalo $(-\infty, 0)$ se obtienen al sustituir $t = -x$ en la ecuación diferencial.

Método de resolución

En el caso de las ecuaciones diferenciales con coeficientes constantes propusimos como solución una función de la forma

$$y(x) = e^{kx}$$

De manera similar, en este caso se prueba una solución de la forma

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

Donde $k$ es un valor que se debe determinar. Al sustituir $x^{k}$, cada término de una ecuación de Cauchy – Euler se convierte en un polinomio en $k$ veces $x^{k}$, puesto que

\begin{align*}
a_{n}x^{n}\dfrac{d^{n}y}{dx^{n}} &= a_{n}x^{n} \left[ k(k -1)(k -2) \cdots (k -n + 1)x^{k -n} \right] \\
&= \left[ a_{n} k(k -1)(k -2) \cdots (k -n + 1) \right] x^{k}
\end{align*}

Por ejemplo, cuando sustituimos $y = x^{k}$ y las respectivas derivadas en la ecuación de segundo orden (\ref{2}), se obtiene

\begin{align*}
ax^{2} \dfrac{d^{2}y}{dx^{2}} + bx \dfrac{dy}{dx} + cy &= ax^{2}[k(k -1)x^{k -2}] + bx[kx^{k -1}] + cx^{k} \\
&= ak(k -1)x^{k} + bkx^{k} + cx^{k} \\
&= \left[ ak(k -1) + bk + c \right] x^{k} \label{4} \tag{4}
\end{align*}

Así, $y = x^{k}$ es una solución de la ecuación diferencial homogénea siempre que $k$ sea una solución de la ecuación auxiliar

$$ak(k -1) + bk + c = 0$$

o bien,

$$ak^{2} + (b -a)k + c = 0 \label{5} \tag{5}$$

Hay tres casos distintos a considerar que dependen de si las raíces de esta ecuación auxiliar son reales y distintas, reales e iguales o complejas.

Caso 1: Raíces reales y distintas

Sean $k_{1}$ y $k_{2}$ las raíces reales de (\ref{5}), tales que $k_{1} \neq k_{2}$. Entonces

$$y_{1} = x^{k_{1}} \hspace{1cm} y \hspace{1cm} y_{2} = x^{k_{2}}$$

forman un conjunto fundamental de soluciones. El Wronskiano esta dado como

\begin{align*}
W(x^{k_{1}}, x^{k_{2}}) &= \begin{vmatrix}
x^{k_{1}} & x^{k_{2}} \\
k_{1}x^{k_{1} -1} & k_{2}x^{k_{2} -1} \end{vmatrix} \\
&= k_{2}x^{(k_{1} + k_{2} -1)} -k_{1}x^{(k_{2} + k_{1} -1)}
\end{align*}

Como

$$W(x^{k_{1}}, x^{k_{2}}) = (k_{2} -k_{1}) x^{k_{1} + k_{2} -1} \neq 0$$

$\forall x \in \delta$, entonces la solución general de la ecuación de Cauchy – Euler para $x > 0$, en el caso en el que las raíces son reales y distintas, es

$$y(x) = c_{1}x^{k_{1}} + c_{2}x^{k_{2}} \label{6} \tag{6}$$

Caso 2: Raíces reales repetidas

Si las raíces de (\ref{5}) son repetidas, es decir $k_{1} = k_{2}$, entonces se obtiene sólo una solución particular

$$y = x^{k_{1}} = x^{k_{2}} = x^{k}$$

Cuando las raíces de la ecuación auxiliar (\ref{5}) son iguales, el discriminante necesariamente es cero, es así que de (\ref{5}) se deduce que las raíces deben ser

$$k = -\dfrac{(b -a)}{2a}$$

Cuando estudiamos el método de reducción de orden vimos que conocida una solución no trivial $y_{1}$, una segunda solución $y_{2}$, tal que $y_{1}$ y $y_{2}$ formen un conjunto fundamental de soluciones, puede ser determinada por la expresión

$$y_{2}(x) = y_{1}(x) \int{\dfrac{e^{-\int{P(x) dx}}}{y^{2}_{1}(x)}} \label{7} \tag{7}$$

Para usar este resultado escribamos a la ecuación de Cauchy – Euler en su forma estándar.

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

Identificamos que

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

Vemos que

$$\int{P(x) dx} = \int{\dfrac{b}{ax} dx} = \dfrac{b}{a} \ln(x)$$

Sustituyendo en (\ref{7}) obtenemos lo siguiente

\begin{align*}
y_{2}(x) &= x^{k} \int{\dfrac{e^{-(b/a) \ln(x)}}{x^{2k}} dx} \\
&= x^{k} \int{ \dfrac{x^{-b/a}}{x^{2k}} dx} \\
&= x^{k} \int{\dfrac{x^{-b/a}}{x^{-(b -a)/a}} dx} \\
&= x^{k} \int{ \dfrac{dx}{x}} \\
&= x^{k} \ln (x)
\end{align*}

En el proceso se ha considerado que

$$e^{-(b/a)}\ln (x) = e^{ \ln (x^{-b/a})} = x^{-b/a}$$

y

$$2k = -\dfrac{(b -a)}{a}$$

Entonces, la segunda solución es

$$y_{2}(x) = x^{k} \ln(x)$$

Vemos que

\begin{align*}
W(x^{k} ,x^{k} \ln (x)) &= \begin{vmatrix}
x^{k} & x^{k} \ln (x) \\
kx^{k -1} & kx^{k -1} \ln (x) + x^{k -1} \end{vmatrix} \\
&= kx^{2k -1} \ln (x) + x^{2k -1} -kx^{2k -1} \ln (x) \\
&= x^{2k -1}
\end{align*}

Como

$$W(x^{k} ,x^{k} \ln (x)) = x^{2k -1} \neq 0$$

$\forall x \in \delta$, entonces la solución general de la ecuación de Cauchy – Euler para $x > 0$, en el caso en el que las raíces son iguales, es

$$y(x) = c_{1}x^{k} + c_{2}x^{k} \ln(x) \label{9} \tag{9}$$

Para ecuaciones de orden superior, si $k$ es una raíz de multiplicidad $r$, entonces se puede demostrar que

$$x^{k}, \hspace{0.5cm} x^{k} \ln(x), \hspace{0.5cm} x^{k}(\ln(x))^{2}, \hspace{0.5cm} \cdots, \hspace{0.5cm} x^{k}(\ln(x))^{r -1}$$

son $r$ soluciones linealmente independientes. En correspondencia, la solución general de la ecuación diferencial debe contener una combinación lineal de estas $r$ soluciones.

Caso 3: Raíces complejas conjugadas

Si las raíces de (\ref{5}) son el par conjugado

$$k_{1} = \alpha + i \beta \hspace{1cm} y \hspace{1cm} k_{2} = \alpha -i \beta$$

Donde $\alpha$ y $\beta > 0$ son reales, entonces una solución es

$$y(x) = C_{1}x^{\alpha + i \beta} + C_{2}x^{\alpha -i \beta} \label{10} \tag{10}$$

De tarea moral muestra que

$$W(x^{\alpha + i \beta}, x^{\alpha -i \beta}) = -2i \beta x^{2\alpha -1} \neq 0$$

lo que indica que la solución (\ref{10}) está compuesta por las funciones del conjunto fundamental de soluciones de la ecuación diferencial de Cauchy – Euler.

Tal como lo hicimos en el caso de coeficientes constantes, se desea escribir la solución en términos de funciones reales. Consideremos la identidad

$$x^{i \beta} = (e^{\ln (x^{i \beta})}) = e^{i \beta \ln (x)} \label{11} \tag{11}$$

Usando la fórmula de Euler podemos escribir

$$x^{i \beta} = \cos \left( \beta \ln (x) \right) + i \sin \left( \beta \ln (x) \right) \label{12} \tag{12}$$

De forma similar,

$$x^{-i \beta} = \cos \left( \beta \ln (x) \right) -i \sin \left( \beta \ln (x) \right) \label{13} \tag{13}$$

Si se suman y restan los dos últimos resultados, se obtiene lo siguiente, respectivamente

$$x^{i \beta } + x^{-i \beta } = 2 \cos \left( \beta \ln (x) \right) \hspace{1cm} y \hspace{1cm} x^{i \beta } -x^{i \beta } = 2i \sin \left( \beta \ln (x) \right) \label{14} \tag{14}$$

Debido a que (\ref{10}) es una solución para cualquier valor de las constantes, podemos notar que si elegimos $C_{1} = C_{2} = 1$ y, por otro lado, $C_{1} = 1, C_{2} = -1$, obtenemos las siguientes dos soluciones, respectivamente

$$y_{1}(x) = x^{\alpha}(x^{i \beta} + x^{-i \beta}) \hspace{1cm} y \hspace{1cm} y_{2}(x) = x^{\alpha} (x^{i \beta} -x^{-i \beta}) \label{15} \tag{15}$$

Usando (\ref{14}) podemos escribir

$$y_{1}(x) = 2x^{\alpha} \cos \left( \beta \ln (x) \right) \hspace{1cm} y \hspace{1cm} y_{2}(x) = 2ix^{\alpha} \sin \left( \beta \ln (x) \right) \label{16} \tag{16}$$

De tarea moral muestra que

$$W(x^{\alpha} \cos \left( \beta \ln (x) \right), x^{\alpha } \sin \left(\beta \ln (x) \right) = \beta x^{2\alpha -1} \neq 0$$

Con esto se concluye que

$$y_{1}(x) = x^{\alpha} \cos \left( \beta \ln (x) \right) \hspace{1cm} y \hspace{1cm} y_{2}(x) = x^{\alpha} \sin \left( \beta \ln (x) \right) \label{17} \tag{17}$$

constituyen un conjunto fundamental de soluciones reales de la ecuación diferencial. Así, la solución general de la ecuación de Cauchy – Euler para $x > 0$, en el caso en el que las raíces son complejas conjugadas, es

$$y(x) = x^{\alpha} \left[ c_{1} \cos \left( \beta \ln (x) \right) + c_{2} \sin \left( \beta \ln (x) \right) \right] \label{18} \tag{18}$$

Realicemos algunos ejemplos en los que apliquemos cada caso.

Ejemplo: Resolver la ecuación de Cauchy – Euler

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

Solución: Consideremos la solución $y = x^{k}$, las respectivas derivadas son

$$\dfrac{dy}{dx} = kx^{k -1} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} = k(k -1)x^{k -2}$$

Sustituimos en la ecuación diferencial.

$$x^{2} \left[ k(k -1)x^{k -2} \right] + \dfrac{2}{3} x \left[ kx^{k -1} \right] -\dfrac{2}{9}x^{k} = x^{k} \left[ k(k -1) + \dfrac{2}{3}k -\dfrac{2}{9} \right] = 0$$

Como $x \neq 0$, entonces la ecuación auxiliar es

$$k(k -1) + \dfrac{2}{3}k -\dfrac{2}{9} = 0$$

o bien,

$$k^{2} -\dfrac{1}{3}k -\dfrac{2}{9} = 0$$

Resolviendo para $k$ obtenemos las raíces $k_{1} = \dfrac{2}{3}$ y $k_{2} = -\dfrac{1}{3}$. Como las raíces son reales y distintas, de acuerdo a (\ref{6}), la solución de la ecuación de Cauchy – Euler es

$$y(x) = c_{1}x^{2/3} + c_{2}x^{-1/3}$$

$\square$

Ejemplo: Resolver la ecuación de Cauchy-Euler

$$x^{2} \dfrac{d^{2}y}{dx^{2}} + 3x \dfrac{dy}{dx} + y = 0$$

Solución: Consideremos la solución $y = x^{k}$, las respectivas derivadas son

$$\dfrac{dy}{dx} = kx^{k -1} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} = k(k -1)x^{k -2}$$

Sustituimos en la ecuación diferencial.

$$x^{2} \left[ k(k -1)x^{k -2} \right] + 3x \left[ kx^{k -1} \right] + x^{k} = x^{k} \left[ k(k -1) + 3k + 1 \right] = 0$$

Como $x \neq 0$, entonces la ecuación auxiliar es

$$k(k -1) + 3k + 1 = 0$$

o bien,

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

Resolviendo para $k$ obtenemos las raíces $k_{1} = k_{2} = -1$. Como las raíces son reales repetidas, por (\ref{9}) concluimos que la solución general de la ecuación de Cauchy – Euler es

$$y(x) = c_{1}x^{-1} + c_{2}x^{-1} \ln (x) = \dfrac{1}{x}[c_{1} + c_{2}\ln(x)]$$

$\square$

Ejemplo: Resolver la ecuación de Cauchy-Euler

$$3x^{2} \dfrac{d^{2}y}{dx^{2}} + 6x \dfrac{dy}{dx} + y = 0$$

Solución: Consideremos la solución $y = x^{k}$, las respectivas derivadas son

$$\dfrac{dy}{dx} = kx^{k -1} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} = k(k -1)x^{k -2}$$

Sustituimos en la ecuación diferencial.

$$3x^{2} \left[ k(k -1)x^{k -2} \right] + 6x \left[ kx^{k -1} \right] + x^{k} = x^{k} \left[ 3k(k -1) + 6k + 1 \right] = 0$$

Como $x \neq 0$, entonces la ecuación auxiliar es

$$3k(k -1) + 6k + 1 = 0$$

o bien,

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

Resolviendo para $k$ obtenemos las raíces

$$k_{1} = -\dfrac{1}{2} + i \dfrac{1}{2 \sqrt{3}} \hspace{1cm} y \hspace{1cm} k_{2} = -\dfrac{1}{2} -i \dfrac{1}{2\sqrt{3}}$$

Identificamos que

$$\alpha = -\dfrac{1}{2} \hspace{1cm} y \hspace{1cm} \beta = \dfrac{1}{2\sqrt{3}}$$

Las raíces son complejas conjugadas, de manera que la solución esta dada por (\ref{18}). Así, la solución general de la ecuación de Cauchy – Euler es

$$y(x) = x^{-1/2} \left[ c_{1} \cos \left( \dfrac{1}{2 \sqrt{3}} \ln (x) \right) + c_{2} \sin \left( \dfrac{1}{2 \sqrt{3}} \ln (x) \right) \right]$$

$\square$

Caso no homogéneo

Para resolver la ecuación no homogénea (\ref{3}) podemos aplicar el método de variación de parámetros visto en la entrada anterior, pues basta encontrar el conjunto fundamental de soluciones $\{ y_{1}, y_{2}\}$ de la ecuación homogénea asociada y con ello aplicar la fórmula de la solución particular, esto es

$$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{19} \tag{19}$$

Recordar que la función $g(x)$ se obtiene de la forma estándar de la ecuación diferencial.

Realicemos un ejemplo.

Ejemplo: Usando el método de variación de parámetros, resolver la ecuación de Cauchy – Euler

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

Solución: Debemos hallar el conjunto fundamental de soluciones, así que primero debemos resolver la ecuación homogénea asociada.

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

Consideremos la solución $y = x^{k}$ y sus derivadas

$$\dfrac{dy}{dx} = kx^{k -1} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} = k(k -1)x^{k -2}$$

Sustituimos en la ecuación homogénea asociada.

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

La ecuación auxiliar es

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

De donde $k_{1} = k_{2} = 1$, así la solución complementaria es

$$y_{c}(x) = c_{1} x + c_{2} x \ln (x)$$

Las funciones

$$y_{1}(x) = x \hspace{1cm} y \hspace{1cm} y_{2}(x) = x\ln(x)$$

conforman al conjunto fundamental de soluciones. Para determinar el Wronskiano vamos a considerar la primer derivada de cada solución.

$$\dfrac{dy_{1}}{dx} = 1 \hspace{1cm} y \hspace{1cm} \dfrac{dy_{2}}{dx} = \ln(x) + 1$$

Sustituimos en el Wronskiano

$$W = \begin{vmatrix}
x & x \ln(x) \\
1 & \ln(x) + 1
\end{vmatrix} = x \ln(x) + x -x\ln(x) = x$$

El Wronskiano es

$$W(x) = x$$

Para determinar la función $g$ dividamos entre $x^{2}$ la ecuación diferencial y así escribirla en su forma estándar.

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

Vemos que

$$g(x) = \dfrac{2}{x}$$

Ahora podemos sustituir en la solución particular (\ref{19}).

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

La solución particular es

$$y_{p}(x) = x[\ln (x)]^{2}$$

Por lo tanto, la solución general de la ecuación de Cauchy – Euler será la superposición de ambas soluciones, esto es

$$y(x) = c_{1}x + c_{2}x \ln(x) + x[\ln(x)]^{2}$$

$\square$

Reducción a coeficientes constantes

Las similitudes entre las formas de las soluciones de ecuaciones de Cauchy – Euler y soluciones de ecuaciones con coeficientes constantes no son una coincidencia.

Por ejemplo, cuando las raíces de las ecuaciones auxiliares para

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy =0$$

y

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

son distintas y reales, las soluciones generales respectivas, para $x > 0$, son

$$y(x) = c_{1}e^{k_{1}x} + c_{2}e^{k_{2}x} \hspace{1cm} y \hspace{1cm} y(x) = c_{1}x^{k_{1}} + c_{2}x^{k_{2}} \label{20} \tag{20}$$

Usando la identidad

$$e^{\ln x} = x$$

Para $x> 0$, la segunda solución dada en (\ref{20}) puede expresarse en la misma forma que la primera solución.

$$y(x) = c_{1}e^{k_{1} \ln(x)} + c_{2}e^{k_{2} \ln (x)} = c_{1}e^{k_{1}t} + c_{2}e^{k_{2}t} \label{21} \tag{21}$$

donde $t = \ln (x)$. Este resultado ilustra que cualquier ecuación de Cauchy – Euler se puede escribir como una ecuación con coeficientes constantes haciendo la sustitución $x = e^{t}$ y con esto resolver la nueva ecuación diferencial en términos de la variable $t$, usando los métodos descritos en la entrada correspondiente y una vez obtenida la solución general, sustituir nuevamente $t = \ln (x)$. Este método requiere del uso de la regla de la cadena.

Si se hace la sustitución $x = e^{t}$, (o bien $t = \ln(x)$), aplicando la regla de la cadena obtenemos las siguientes expresiones para las derivadas.

$$\dfrac{dy}{dx} = \dfrac{dy}{dt} \dfrac{dt}{dx} = \dfrac{1}{x} \dfrac{dy}{dt} \label{22} \tag{22}$$

y

$$\dfrac{d^{2}y}{dx^{2}} = \dfrac{d}{dx} \left( \dfrac{1}{x} \dfrac{dy}{dt} \right) = -\dfrac{1}{x^{2}} \dfrac{dy}{dt} + \dfrac{1}{x^{2}} \dfrac{d^{2}y}{dt^{2}} = \dfrac{1}{x^{2}} \left( \dfrac{d^{2}y}{dt^{2}} -\dfrac{dy}{dt} \right) \label{23} \tag{23}$$

Sustituyendo en la ecuación de Cauchy – Euler obtenemos lo siguiente.

\begin{align*}
ax^{2} \dfrac{d^{2}y}{dx^{2}} + bx \dfrac{dy}{dx} + cy &= ax^{2} \left[ \dfrac{1}{x^{2}} \left( \dfrac{d^{2}y}{dt^{2}} -\dfrac{dy}{dt} \right) \right] + bx \left[ \dfrac{1}{x} \dfrac{dy}{dt} \right] + cy \\
&= a \dfrac{d^{2}y}{dt^{2}} + (b -a) \dfrac{dy}{dt} + cy
\end{align*}

Por lo tanto, haciendo la sustitución $x = e^{t}$ reducimos la ecuación de Cauchy – Euler a la ecuación

$$a\dfrac{d^{2}y(t)}{dt^{2}} + (b -a) \dfrac{dy(t)}{dt} + cy(t) = g(t) \label{24} \tag{24}$$

que corresponde a una ecuación diferencial con coeficientes constantes en donde la variable independiente es $t$.

Realicemos un ejemplo.

Ejemplo: Usar el cambio de variable $x= e^{t}$ para convertir la ecuación de Cauchy – Euler

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

en una ecuación de coeficiente constantes y obtener la solución general.

Solución: Consideremos el cambio de variable $x = e^{t}$, usando los resultados (\ref{22}) y (\ref{23}), la ecuación de Cauchy – Euler queda como sigue

$$x^{2}\left[ \dfrac{1}{x^{2}} \left( \dfrac{d^{2}y}{dt^{2}} -\dfrac{dy}{dt} \right) \right] -3x \left[ \dfrac{1}{x} \dfrac{dy}{dt} \right] +13y = 4 + 3e^{t}$$

Esto es,

$$\dfrac{d^{2}y}{dt^{2}} -4 \dfrac{dy}{dt} + 13y = 4 + 3e^{t}$$

Ahora tenemos una ecuación no homogénea con coeficientes constantes. Comencemos por resolver la ecuación homogénea.

$$\dfrac{d^{2}y}{dt^{2}} -4 \dfrac{dy}{dt} + 13y = 0$$

La ecuación auxiliar es

$$k^{2} -4k + 13 = 0$$

Las raíces son

$$k_{1} = 2 + i3 \hspace{1cm} y \hspace{1cm} k_{2} = 2 -i3$$

Identificamos que $\alpha = 2 $ y $\beta = 3$, entonces la solución complementaria, en la variable $t$, es

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

Las funciones correspondientes al conjunto fundamental de soluciones son

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

Las derivadas correspondientes son

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

El Wronskiano esta dado por

$$W = \begin{vmatrix}
e^{2t} \cos(3t) & e^{2t} \sin(3t) \\
2e^{2t} \cos(3t) -3e^{2t} \sin(3t) & 2e^{2t} \sin(3t) + 3e^{2t} \cos(3t)
\end{vmatrix}$$

Calculando el determinante obtendremos

$$W(t) = 3e^{4t}$$

La ecuación diferencial ya se encuentra en su forma estándar, así que la función $g$ es

$$g(t) = 4 + 3e^{t}$$

Ahora podemos sustituir las funciones correspondientes en la solución particular (\ref{19}) para la variable $t$.

\begin{align*}
y_{p}(t) &= -e^{2t} \cos(3t) \int{ \dfrac{e^{2t} \sin(3t) (4 + 3e^{t})} {3e^{4t}} dt} + e^{2t} \sin(3t) \int{ \dfrac{e^{2t} \cos(3t) (4 + 3e^{t})}{3e^{4t}} dt} \\
&= -e^{2t} \cos(3t) \left[ \dfrac{4}{3} \int{ \dfrac{\sin(3t)}{e^{2t}} dt} + \int{ \dfrac{\sin(3t)}{e^{t}} dt }\right] \\
&+ e^{2t} \sin(3t) \left[ \dfrac{4}{3} \int{ \dfrac{\cos(3t)}{e^{2t}} dt} + \int{ \dfrac{\cos(3t)}{e^{t}} dt} \right]
\end{align*}

Las integrales se resuelven con integración por partes. De tarea moral desarrolla el cálculo de cada integral, los resultados correspondientes son

$$\int{ \dfrac{\sin(3t)}{e^{2t}} dt} = -\dfrac{2}{13}e^{-2t} \sin(3t) -\dfrac{3}{13}e^{-2t} \cos(3t)$$

$$\int{ \dfrac{\sin(3t)}{e^{t}} dt} = -\dfrac{3}{10}e^{-t} \cos(3t) -\dfrac{1}{10}e^{-t} \sin(3t)$$

$$\int{ \dfrac{\cos(3t)}{e^{2t}} dt} = -\dfrac{2}{13}e^{-2t} \cos(3t) + \dfrac{3}{13}e^{-2t} \sin(3t)$$

$$ \int{ \dfrac{\cos(3t)}{e^{t}} dt} = \dfrac{3}{10}e^{-t} \sin(3t) -\dfrac{1}{10}e^{-t} \cos(3t)$$

Sustituyendo estos resultados en $y_{p}(t)$ y reduciendo la expresión obtendremos la solución particular

$$y_{p}(t) = \dfrac{4}{13} + \dfrac{3}{10}e^{t}$$

Por tanto, la solución general en términos de la variable $t$ es

$$y(t) = c_{1}e^{2t} \cos(3t) + c_{2}e^{2t} \sin(3t) + \dfrac{4}{13} + \dfrac{3}{10}e^{t}$$

Si regresamos a la variable original $x = e^{t}$ obtenemos finalmente que la solución general de la ecuación de Cauchy – Euler es

$$y(x) = c_{1}x^{2} \cos[3 \ln(x)] + c_{2}x^{2} \sin[3 \ln(x)] + \dfrac{4}{13} + \dfrac{3}{10}x$$

$\square$

Con esto concluimos el estudio de la ecuación de Cauchy – Euler y en general con el estudio de las ecuaciones diferenciales lineales de orden superior con coeficientes constantes y ecuaciones sencillas 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. Resolver las siguientes ecuaciones de Cauchy – Euler.

• $x^{2} \dfrac{d^{2}y}{dx^{2}} -12y = 0$

• $x^{2} \dfrac{d^{2}y}{dx^{2}} + 5x \dfrac{dy}{dx} + 4y = 0$

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

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

• $25x^{2} \dfrac{d^{2}y}{dx^{2}} + 25x \dfrac{dy}{dx} + y = 0$

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

2. Resolver las siguientes ecuaciones de Cauchy – Euler usando el método de variación de parámetros.

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

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

3. Usar el cambio de variable $x = e^{t}$ para convertir las ecuaciones de Cauchy – Euler en ecuaciones diferenciales con coeficientes constantes y resolver la ecuación.

• $x^{2} \dfrac{d^{2}y}{dx^{2}} + 10x \dfrac{dy}{dx} + 8y = x^{2}$

• $x^{2} \dfrac{d^{2}y}{dx^{2}} -4x \dfrac{dy}{dx} + 6y = \ln (x^{2})$

4. Hacer una extensión a orden 3 de la teoría desarrollada en esta entrada y aplicando el método de variación de parámetros obtener la solución general de la siguiente ecuación de Cauchy – Euler.

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

Más adelante…

Con esto concluimos la primera parte de la unidad dos. En la siguiente entrada abordaremos el tema de las oscilaciones mecánicas como ejemplo de aplicación de la teoría que hemos desarrollado hasta este momento.

En entradas posteriores haremos un estudio detallado sobre las ecuaciones diferenciales de orden superior con coeficientes variables y con ello cerraremos la segunda unidad.

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 variación de parámetros
• Siguiente entrada del curso: Oscilaciones mecánicas

• Video relacionado al tema: Ecuación de Euler

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»