1er METODO: COEFICIENTES INDETERMINADOS: MÉTODO DE SUPERPOSICIÓN

INTRODUCCIÓN

Para resolver una EDLNH

𝒂𝒏 𝒚(𝒏) + 𝒂𝒏−𝟏 𝒚(𝒏−𝟏) + ⋯ + 𝒂𝟏 𝒚′ + 𝒂𝟎 𝒚 = 𝒈(𝒙) (1)

Se debe hacer dos cosas:

Encontrar la función o solución complementaria 𝒚𝒄 de la EDLH asociada a (1)
(𝒏) (𝒏−𝟏)
𝒂𝒏 𝒚 + 𝒂𝒏−𝟏 𝒚 + ⋯ + 𝒂𝟏 𝒚′ + 𝒂𝟎 𝒚 = 𝟎
Encontrar alguna solución particular 𝒚𝒑 de la EDLNH (1).
De tal manera que la solución general de la EDL (1) será:
𝒚 = 𝒚𝒄 + 𝒚𝒑
Ya se estudiaron métodos para encontrar 𝒚𝒄 , por lo que los métodos siguientes permitirán
encontrar 𝒚𝒑
MÉTODO DE COEFICIENTES INDETERMINADOS

La idea fundamental de este método, trata de realizar una conjetura acerca de la forma de
𝒚𝒑 que está motivada en la forma que presenta o tiene la función 𝒈(𝒙).
El método se limita a ED lineales como (1) donde
• Los coeficientes 𝒂𝒊 , 𝒊 = 𝟎, 𝟏, … , 𝒏 son cttes y
 𝒈(𝒙) puede ser:

 𝑼𝒏𝒂 𝒇𝒖𝒏𝒄𝒊ó𝒏 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕𝒆 𝑲

 𝑼𝒏𝒂 𝒇𝒖𝒏𝒄𝒊ó𝒏 𝒑𝒐𝒍𝒊𝒏𝒐𝒎𝒊𝒂𝒍
 𝑼𝒏𝒂 𝒇𝒖𝒏𝒄𝒊ó𝒏 𝒆𝒙𝒑𝒐𝒏𝒆𝒏𝒄𝒊𝒂𝒍 𝒆𝜶𝒙
 𝑼𝒏𝒂 𝒇𝒖𝒏𝒄𝒊ó𝒏 𝒔𝒆𝒏𝒐 𝒐 𝒄𝒐𝒔𝒆𝒏𝒐 𝒔𝒆𝒏 𝜷𝒙 𝒐 𝒄𝒐𝒔 𝜷𝒙
 𝑶 𝒔𝒖𝒎𝒂𝒔 𝒚 𝒑𝒓𝒐𝒅𝒖𝒄𝒕𝒐𝒔 𝒇𝒊𝒏𝒊𝒕𝒐𝒔 𝒅𝒆 𝒆𝒔𝒕𝒂𝒔 𝒇𝒖𝒏𝒄𝒊𝒐𝒏𝒆𝒔

A continuación algunos ejemplos de los tipos de 𝒈(𝒙)



Es decir, g(x) es una combinación lineal de funciones de la clase

Donde 𝒏 es un entero
no negativo y
𝜶 𝒚 𝜷 son números
reales

El método no es aplicable a ecuaciones de la forma (1) cuando

Las EDLNH donde 𝑔(𝑥) tenga

una de estas formas, se
estudiarán más adelante

El conjunto de funciones que consiste en 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕𝒆𝒔, 𝒑𝒐𝒍𝒊𝒏𝒐𝒎𝒊𝒐𝒔, 𝒆𝒂𝒙 , 𝒔𝒆𝒏 𝜷𝒙 𝒚 𝒄𝒐𝒔 𝜷𝒙

tienen la notable propiedad de que las derivadas de sus sumas y productos son de nuevo
sumas y productos de 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕𝒆𝒔, 𝒑𝒐𝒍𝒊𝒏𝒐𝒎𝒊𝒐𝒔, 𝒆𝒂𝒙 , 𝒔𝒆𝒏 𝜷𝒙 𝒚 𝒄𝒐𝒔 𝜷𝒙.

Debido a que la combinación lineal de derivadas de la solución particular

(𝒏) (𝒏−𝟏)
𝒂𝒏 𝒚𝒑 + 𝒂𝒏−𝟏 𝒚𝒑 + ⋯ + 𝒂𝟏 𝒚′𝒑 + 𝒂𝟎 𝒚𝒑
 debe ser idéntica a 𝒈(𝒙), parece razonable suponer que 𝒚𝒑 tiene la misma forma que 𝒈(𝒙).
La tabla que sigue muestran algunos ejemplos de 𝒈(𝒙) en (1) junto con la forma
correspondiente de la solución particular.

EJERCICIO – 1
Resolver los ejercicios siguientes:

2do METODO: COEFICIENTES INDETERMINADOS: MÉTODO DEL ANULADOR

La EDLNH
𝒂𝒏 𝒚(𝒏) + 𝒂𝒏−𝟏 𝒚(𝒏−𝟏) + ⋯ + 𝒂𝟏 𝒚′ + 𝒂𝟎 𝒚 = 𝒈(𝒙)

Utilizando el ODL, también se lo puede escribir como sigue

𝒂𝒏 𝑫(𝒏) 𝒚 + 𝒂𝒏−𝟏 𝑫(𝒏−𝟏) 𝒚 + ⋯ + 𝒂𝟏 𝑫𝒚 + 𝒂𝟎 𝒚 = 𝒈(𝒙) (𝟏)

o [𝒂𝒏 𝑫𝒏 + 𝒂𝒏−𝟏 𝑫𝒏−𝟏 + ⋯ + 𝒂𝟏 𝑫 + 𝒂𝟎 ]𝒚 = 𝒈(𝒙) ⟹ 𝑳(𝒚) = 𝒈(𝒙)

𝐿 ⟹ 𝑂𝑝𝑒𝑟𝑎𝑑𝑜𝑟 𝑑𝑖𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎𝑙 𝑙𝑖𝑛𝑒𝑎𝑙 (𝑂𝐷𝐿) 𝑑𝑒 𝑛 − é𝑠𝑖𝑚𝑜 𝑜𝑟𝑑𝑒𝑛

𝒂𝒏 𝑫𝒏 + 𝒂𝒏−𝟏 𝑫𝒏−𝟏 + ⋯ + 𝒂𝟏 𝑫 + 𝒂𝟎 (𝟐)

Utilizando este método del anulador, la 𝒚𝒑 se deduce casi de manera automática una vez
que se encuentra un ODL adecuado que anula a 𝒈(𝒙) en (1).
Pero antes dos conceptos.

1.- FACTORIZACIÓN DE OPERADORES

Si los coeficientes 𝒂𝒊 , 𝒊 = 𝟎, 𝟏, . . . , 𝒏 son cttes R, un ODL (1) se puede factorizar siempre que
el polinomio característico 𝒂𝒏 𝒎𝒏 + 𝒂𝒏−𝟏 𝒎𝒏−𝟏 + ⋯ + 𝒂𝟏 𝒎 + 𝒂𝟎 sea factorizable.
En otras palabras, si 𝒓𝟏 es una raíz de la ecuación auxiliar
𝒂𝒏 𝒎𝒏 + 𝒂𝒏−𝟏 𝒎𝒏−𝟏 + ⋯ + 𝒂𝟏 𝒎 + 𝒂𝟎 = 𝟎
entonces 𝑳 = (𝑫 − 𝒓𝟏 ) 𝑷(𝑫), donde la expresión polinomial 𝑷(𝑫) es un ODL de orden 𝒏 − 𝟏.
Por ejemplo, si se trata a 𝑫 como una cantidad algebraica, entonces el operador diferencial
(𝑫𝟐 + 𝟓𝑫 + 𝟔)𝒚 = (𝑫 + 𝟐)(𝑫 + 𝟑)𝒚 = (𝑫 + 𝟑)(𝑫 + 𝟐)𝒚
Esto muestra una propiedad general:
Los factores de un ODL con coeficientes constantes conmutan.
Una ecuación diferencial tal como
 𝒚′′ + 𝟒𝒚′ + 𝟒𝒚 = 𝟎 𝑝𝑢𝑒𝑑𝑒 𝑠𝑒𝑟 𝑒𝑠𝑐𝑟𝑖𝑡𝑜 𝑐𝑜𝑚𝑜 (𝑫𝟐 + 𝟒𝑫 + 𝟒)𝒚 = 𝟎 𝒐 (𝑫 + 𝟐)(𝑫 + 𝟐)𝒚 = 𝟎 𝒐 (𝑫 + 𝟐)𝟐 𝒚 = 𝟎
2.- OPERADOR ANULADOR.
Si 𝑳 es un ODL con coeficientes cttes y 𝑓(𝑥) es una función suficientemente derivable tal
que
𝑳( 𝒇 (𝒙)) = 𝟎,

entonces se dice que 𝑳 es un anulador de la función 𝑓(𝑥). Por ejemplo:

El operador diferencial 𝑫 𝒂𝒏𝒖𝒍𝒂 𝒂 𝒚 = 𝒌 puesto que 𝑫𝒌 = 𝟎.
El operador diferencial 𝑫𝟐 𝒂𝒏𝒖𝒍𝒂 𝒂 𝒚 = 𝒙 puesto que 𝑫𝟐 𝒙 = 𝟎
El operador diferencial 𝑫𝟑 𝒂𝒏𝒖𝒍𝒂 𝒂 𝒚 = 𝒙𝟐 puesto que 𝑫𝟑 𝒙𝟐 = 𝟎, etcétera.
El operador diferencial 𝑫𝒏 anula a cada una de las funciones

𝟏, 𝒙, 𝒙𝟐 , 𝒙𝟑 , … , 𝒙𝒏−𝟏 (𝟑)

Como una consecuencia inmediata de (3) y el hecho de que la derivación se puede hacer
término a término, un polinomio

𝒄𝟎 + 𝒄𝟏 𝒙 + 𝒄𝟐 𝒙𝟐 + ⋯ + 𝒄𝒏−𝟏 𝒙𝒏−𝟏 (𝟒)

se anula al encontrar un operador que anule la potencia más alta de x.

Las funciones que se anulan por un ODL de n-ésimo orden 𝑳 son simplemente
aquellas funciones que se obtienen de la solución general de la EDH 𝑳(𝒚) = 𝟎.

El operador diferencial (𝑫 − 𝜶)𝒏 anula a cada una de las funciones

𝒆𝜶𝒙 , 𝒙𝒆𝜶𝒙 , 𝒙𝟐 𝒆𝜶𝒙 , … , 𝒙𝒏−𝟏 𝒆𝜶𝒙 (𝟓)

Para ver esto, observe que la ecuación auxiliar de la EDH (𝑫 − 𝜶)𝒏 𝒚 = 𝟎 𝒆𝒔 (𝒎 − 𝜶)𝒏 𝒚 = 𝟎.
Puesto que 𝜶 es una raíz de multiplicidad 𝒏, la solución general es

𝒚 = 𝒄𝟏 𝒆𝜶𝒙 + 𝒄𝟐 𝒙𝒆𝜶𝒙 + 𝒄𝟑 𝒙𝟐 𝒆𝜶𝒙 + ⋯ + 𝒄𝒏−𝟏 𝒙𝒏−𝟏 𝒆𝜶𝒙 (𝟔)

EJEMPLO de operador anulador. Encontrar un operador diferencial que anule la función

dada

Cuando 𝜶 𝒚 𝜷, 𝜷 > 𝟎 son Nros R, la fórmula cuadrática revela que [𝒎𝟐 − 𝟐𝜶𝒎 + (𝜶𝟐 +
𝜷𝟐 )]𝒏 = 𝟎 tiene raíces complejas 𝜶 + 𝒊𝜷 𝒚 𝜶 − 𝒊𝜷, ambas de multiplicidad 𝒏. Del análisis al
final de la sección 4.3, se tiene el siguiente resultado.

𝐸𝑙 𝑜𝑝𝑒𝑟𝑎𝑑𝑜𝑟 𝑑𝑖𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎𝑙 [𝑫𝟐 − 𝟐𝜶𝑫 + (𝜶𝟐 + 𝜷𝟐 )]𝒏 𝑎𝑛𝑢𝑙𝑎 𝑐𝑎𝑑𝑎 𝑢𝑛𝑎𝑑𝑒 𝑙𝑎𝑠 𝑓𝑢𝑛𝑐𝑖𝑜𝑛𝑒𝑠

𝒆𝜶𝒙 𝐜𝐨𝐬 𝜷𝒙 , 𝒙𝒆𝜶𝒙 𝐜𝐨𝐬 𝜷𝒙 , 𝒙𝟐 𝒆𝜶𝒙 𝐜𝐨𝐬 𝜷𝒙 , … , 𝒙𝒏−𝟏 𝒆𝜶𝒙 𝐜𝐨𝐬 𝜷𝒙

𝒆𝜶𝒙 𝐬𝐞𝐧 𝜷𝒙 , 𝒙𝒆𝜶𝒙 𝐬𝐞𝐧 𝜷𝒙 , 𝒙𝟐 𝒆𝜶𝒙 𝐬𝐞𝐧 𝜷𝒙 , … , 𝒙𝒏−𝟏 𝒆𝜶𝒙 𝐬𝐞𝐧 𝜷𝒙 (𝟕)

EJEMPLO de operador anulador. Encontrar un operador diferencial que anule la función

Cuando 𝜶 = 𝟎 𝒚 𝒏 = 𝟏, un caso especial de (7) es

Por ejemplo 𝑫𝟐 + 𝟏𝟔 anulará cualquier combinación lineal de 𝒔𝒆𝒏 𝟒𝒙 𝒚 𝒄𝒐𝒔 𝟒𝒙.

Con frecuencia estamos interesados en anular la suma de dos o más funciones. Como
acabamos de ver en los ejemplos 1 y 2, si 𝑳 es un ODL tal que 𝑳(𝒚𝟏 ) = 𝟎 𝒚 𝑳(𝒚𝟐 ) = 𝟎,
entonces 𝑳 anulará la combinación lineal 𝒄𝟏 𝒚𝟏 (𝒙) + 𝒄𝟐 𝒚𝟐 (𝒙) ⟹ 𝑳[𝒄𝟏 𝒚𝟏 (𝒙) + 𝒄𝟐 𝒚𝟐 (𝒙)] = 𝟎

Esta es una consecuencia directa del principio de superposición. Supongamos ahora que
𝑳𝟏 𝒚 𝑳𝟐 son ODL con coeficientes cttes tales que 𝑳𝟏 anula a 𝒚𝟏 (𝒙) y 𝑳𝟐 anula a 𝒚𝟐 (𝒙), pero
𝑳𝟏 (𝒚𝟐 ) ≠ 𝟎 y 𝑳𝟐 (𝒚𝟏 ) ≠ 𝟎. Entonces el producto de los operadores diferenciales 𝑳𝟏 𝑳𝟐 anula
 la suma 𝒄𝟏 𝒚𝟏 (𝒙) + 𝒄𝟐 𝒚𝟐 (𝒙). Esto se puede demostrar fácilmente, usando la linealidad y el
hecho de que 𝑳𝟏 𝑳𝟐 = 𝑳𝟐 𝑳𝟏 :

𝑳𝟏 𝑳𝟐 (𝒚𝟏 + 𝒚𝟐 ) = 𝑳𝟏 𝑳𝟐 (𝒚𝟏 ) + 𝑳𝟏 𝑳𝟐 (𝒚𝟐 )

= 𝑳𝟐 𝑳𝟏 (𝒚𝟏 ) + 𝑳𝟏 𝑳𝟐 (𝒚𝟐 )
= 𝑳𝟐 [𝑳𝟏 (𝒚𝟏 )] + 𝑳𝟏 [𝑳𝟐 (𝒚𝟐 )] = 𝟎

0 0

Por ejemplo, sabemos de (3) que 𝑫𝟐 𝒂𝒏𝒖𝒍𝒂 𝒂 𝟕 − 𝒙 y de (8) que 𝑫𝟐 − 𝟏𝟔 𝒂𝒏𝒖𝒍𝒂 𝒂 𝒔𝒆𝒏 𝟒𝒙. Por
tanto el producto de operadores 𝑫𝟐 (𝑫𝟐 − 𝟏𝟔) 𝑎𝑛𝑢𝑙𝑎𝑟á 𝑎 𝑙𝑎 𝑐𝑜𝑚𝑏𝑖𝑛𝑎𝑐𝑖ó𝑛 𝑙𝑖𝑛𝑒𝑎𝑙 𝟕 − 𝒙 + 𝟔𝒔𝒆𝒏𝟒𝒙.

NOTA El operador diferencial que anula una función no es único. Vimos en un ejemplo
anterior que 𝑫 + 𝟑 anula a 𝒆−𝟑𝒙 , pero también operadores diferenciales de orden superior
siempre y cuando 𝑫 + 𝟑 sea uno de los factores del operador.

Por ejemplo (𝑫 + 𝟑)(𝑫 + 𝟏), (𝑫 + 𝟑)𝟐 𝒚 𝑫𝟑 (𝑫 + 𝟑) todos anulan a 𝒆−𝟑𝒙 . Como algo natural,
cuando se busca un anulador diferencial para una función 𝒚 = 𝒇(𝒙), se quiere que el
operador de mínimo orden posible haga el trabajo.

COEFICIENTES INDETERMINADOS

Suponga que 𝑳(𝒚) = 𝒈(𝒙) es una EDL con coeficientes cttes y que la función 𝒈(𝒙) consiste
en sumas y productos finitos de las funciones listadas en (3), (5) y (7), es decir, 𝒈(𝒙) es
una combinación lineal de funciones de la forma

𝒌 (𝒄𝒕𝒕𝒆), 𝒙𝒎 , 𝒙𝒎 𝒆𝜶𝒙 , 𝒙𝒎 𝒆𝜶𝒙 𝐜𝐨𝐬 𝜷𝒙 𝒚 𝒙𝒎 𝒆𝜶𝒙 𝐬𝐞𝐧 𝜷𝒙

donde 𝒎 es un entero no negativo y 𝜶 𝒚 𝜷 son Nros R.

Ahora se sabe que una función tal como 𝒈(𝒙) puede ser anulada por un ODL 𝑳𝟏 de menor
orden, que es producto de los operadores 𝑫𝒏 , (𝑫 − 𝜶)𝒏 𝒚 (𝑫𝟐 − 𝟐𝜶𝑫 + 𝜶𝟐 + 𝜷𝟐 )𝒏 .
Al aplicar 𝑳𝟏 a ambos lados de la ecuación 𝑳(𝒚) = 𝒈(𝒙) se obtiene 𝑳𝟏 𝑳(𝒚) = 𝑳𝟏 (𝒈(𝒙)) = 𝟎.
Al resolver la ecuación homogénea de orden superior 𝑳𝟏 𝑳(𝒚) = 0, se descubre la forma de
una solución particular 𝒚𝒑 para la ecuación original no homogénea 𝑳(𝒚) = 𝒈(𝒙).
Entonces sustituimos esta forma supuesta en 𝑳(𝒚) = 𝒈(𝒙) para encontrar una solución
particular explícita. Este procedimiento para determinar 𝒚𝒑 , llamado método de los
coeficientes indeterminados, se ilustra a continuación en varios ejemplos.

Antes de proceder, recuerde que la solución general de una EDLNH 𝑳(𝒚) = 𝒈(𝒙) es 𝒚 = 𝒚𝒄 +
𝒚𝒑 donde 𝒚𝒄 es la función complementaria, es decir, la solución general de la ecuación
homogénea asociada 𝑳(𝒚) = 𝟎. La solución general de cada ecuación 𝑳(𝒚) = 𝒈(𝒙) se define
en el intervalo (−∞, +∞).

RESUMEN DEL MÉTODO - COEFICIENTES INDETERMINADOS: MÉTODO DEL ANULADOR

La ecuación diferencial 𝐿(𝑦) = 𝑔(𝑥) tiene coeficientes constantes y la función 𝒈(𝒙) consiste en sumas
y productos finitos de 𝒄𝒐𝒏𝒔𝒕𝒂𝒏𝒕𝒆𝒔, 𝒑𝒐𝒍𝒊𝒏𝒐𝒎𝒊𝒐𝒔, 𝒇𝒖𝒏𝒄𝒊𝒐𝒏𝒆𝒔 𝒆𝒙𝒑𝒐𝒏𝒆𝒏𝒄𝒊𝒂𝒍𝒆𝒔 𝒆𝜶𝒙 , 𝒔𝒆𝒏𝒐𝒔 𝒚 𝒄𝒐𝒔𝒆𝒏𝒐𝒔

1) Encuentre la función o solución complementaria 𝒚𝒄 para la EDH 𝑳(𝒚) = 𝟎.

2) Opere ambos lados de la EDNH 𝑳(𝒚) = 𝒈(𝒙) con un ODL 𝑳𝟏 que anula la función 𝒈(𝒙).
3) Determine la solución general de la EDH de orden superior 𝑳𝟏 𝑳(𝒚) = 𝟎.
4) Elimine de la solución del paso 3) los términos que se duplican en la solución
complementaria 𝒚𝒄 encontrada en el paso 1). Forme una combinación lineal 𝒚𝒑 de los
términos restantes. Esta es la forma de una solución particular de 𝑳(𝒚) = 𝒈(𝒙).
5) Sustituya 𝒚𝒑 encontrada en el paso 4) en 𝑳(𝒚) = 𝒈(𝒙). Iguale los coeficientes de las distintas
funciones en cada lado de la igualdad y resuelva el sistema resultante de ecuaciones para
determinar los coeficientes desconocidos de 𝒚𝒑 .
6) Con la solución particular encontrada en el paso v), forme la solución general 𝒚 = 𝒚𝒄 + 𝒚𝒑
de la ED dada.

EJERCICIO – 2
VARIACIÓN DE PARÁMETROS

INTRODUCCIÓN
Este es un 3er método para encontrar una solución particular 𝒚𝑷 de una EDL de orden superior. El
procedimiento utilizado es el utilizado para encontrar 𝒚𝑷 de una EDL de 1er estudiado en un
capítulo anterior.
Para adaptar el método de variación de parámetros a una ED de 2do orden

𝒂𝟐 (𝒙)𝒚′′ + 𝒂𝟏 (𝒙)𝒚′ + 𝒂𝟎 (𝒙)𝒚 = 𝒈(𝒙) (𝟏)

Comenzamos por escribir la ecuación en su forma estándar

𝒚′′ + 𝑷(𝒙)𝒚′ + 𝑸(𝒙)𝒚 = 𝒇(𝒙) (𝟐)

Suponemos en (2) que 𝑷(𝒙), 𝑸(𝒙) 𝒚 𝒇(𝒙) son continuas en algún intervalo común I. Por lo estudiado
en secciones anteriores, no hay dificultad para obtener 𝒚𝒄 , de la EDLH asociada de (2), cuando los
coeficientes son constantes.

SUPOSICIONES

Correspondiendo con la suposición 𝒚𝑷 = 𝝁𝟏 (𝒙)𝒚𝟏 (𝒙) utilizado para encontrar una solución particular
𝒚𝑷 de una EDL de 1er orden en el capítulo 2, para la EDL de 2do orden (2) se busca una solución
de la forma

𝒚𝑷 = 𝝁𝟏 (𝒙)𝒚𝟏 (𝒙) + 𝝁𝟐 (𝒙)𝒚𝟐 (𝒙) (𝟑)

donde 𝒚𝟏 y 𝒚𝟐 forman un conjunto fundamental de soluciones en I. Usando la regla del producto

para derivar dos veces a 𝒚𝑷 , se obtiene

𝒚′𝒑 = 𝝁𝟏 𝒚′𝟏 + 𝒚𝟏 𝝁′𝟏 + 𝝁𝟐 𝒚′𝟐 + 𝒚𝟐 𝝁′𝟐

′′ ′ ′ ′′ ′ ′ ′′ ′ ′ ′′ ′ ′
𝒚′′
𝒑 = 𝝁𝟏 𝒚𝟏 + 𝒚𝟏 𝝁𝟏 + 𝒚𝟏 𝝁𝟏 + 𝝁𝟏 𝒚𝟏 + 𝝁𝟐 𝒚𝟐 + 𝒚𝟐 𝝁𝟐 + 𝒚𝟐 𝝁𝟐 + 𝝁𝟐 𝒚𝟐

Sustituyendo (3) y las derivadas anteriores en (2) 𝒚′′ + 𝑷(𝒙)𝒚′ + 𝑸(𝒙)𝒚 = 𝒇(𝒙) y agrupando términos
se obtiene
𝒚′′ ′
𝒑 + 𝑷(𝒙)𝒚𝒑 + 𝑸(𝒙)𝒚𝑷 =
[𝝁𝟏 𝒚′′ ′ ′ ′′ ′ ′ ′′ ′ ′ ′′ ′ ′ ′ ′ ′ ′
𝟏 + 𝒚𝟏 𝝁𝟏 + 𝒚𝟏 𝝁𝟏 + 𝝁𝟏 𝒚𝟏 + 𝝁𝟐 𝒚𝟐 + 𝒚𝟐 𝝁𝟐 + 𝒚𝟐 𝝁𝟐 + 𝝁𝟐 𝒚𝟐 ] + 𝑷(𝒙)[𝝁𝟏 𝒚𝟏 + 𝒚𝟏 𝝁𝟏 + 𝝁𝟐 𝒚𝟐 + 𝒚𝟐 𝝁𝟐 ] +
𝑸(𝒙)[𝝁𝟏 (𝒙)𝒚𝟏 (𝒙) + 𝝁𝟐 (𝒙)𝒚𝟐 (𝒙) ] =
𝝁𝟏 𝒚′′ ′ ′ ′′ ′ ′ ′′ ′ ′ ′′ ′ ′ ′ ′ ′
𝟏 + 𝒚𝟏 𝝁𝟏 + 𝒚𝟏 𝝁𝟏 + 𝝁𝟏 𝒚𝟏 + 𝝁𝟐 𝒚𝟐 + 𝒚𝟐 𝝁𝟐 + 𝒚𝟐 𝝁𝟐 + 𝝁𝟐 𝒚𝟐 + [𝑷(𝒙)𝝁𝟏 𝒚𝟏 + 𝑷(𝒙)𝒚𝟏 𝝁𝟏 + 𝑷(𝒙)𝝁𝟐 𝒚𝟐 + 𝑷(𝒙)𝒚𝟐 𝝁𝟐 ]
′
′′ ′ ′′ ′ ′′
+ [𝑸(𝒙)𝝁𝟏 (𝒙)𝒚𝟏 (𝒙) + 𝑸(𝒙)𝝁𝟐 (𝒙)𝒚𝟐 (𝒙)] = 𝝁𝟏 [𝒚𝟏 + 𝑷𝒚𝟏 + 𝑸𝟏 𝒚𝟏 ] + 𝝁𝟐 [𝒚𝟐 + 𝑷𝒚𝟐 + 𝑸𝒚𝟐 ] + 𝒚𝟏 𝝁𝟏 +
0 0
𝝁′𝟏 𝒚′𝟏 + 𝒚𝟐 𝝁′′
𝟐 + +𝝁′𝟐 𝒚′𝟐+ 𝒚′𝟏 𝝁′𝟏
+ 𝒚′𝟐 𝝁′𝟐
+ 𝑷[𝒚𝟏 𝝁′𝟏 𝒚𝟐 𝝁′𝟐 ]
𝒅 𝒅
= (𝒚 𝝁′ ) + (𝒚 𝝁′ ) + 𝑷[𝒚𝟏 𝝁′𝟏 + 𝒚𝟐 𝝁′𝟐 ] + 𝒚′𝟏 𝝁′𝟏 + 𝒚′𝟐 𝝁′𝟐 =
𝒅𝒙 𝟏 𝟏 𝒅𝒙 𝟐 𝟐
𝒅
[𝒚𝟏 𝝁′𝟏 + 𝒚𝟐 𝝁′𝟐 ] + 𝑷[𝒚𝟏 𝝁′𝟏 + 𝒚𝟐 𝝁′𝟐 ] + 𝒚′𝟏 𝝁′𝟏 + 𝒚′𝟐 𝝁′𝟐 = 𝒇(𝒙) (𝟒)
𝒅𝒙

Lo que se pretende es determinar dos funciones desconocidas 𝝁𝟏 (𝒙) 𝒚 𝝁𝟐 (𝒙), para ello son
necesarias dos ecuaciones. Estas ecuaciones se obtienen con la suposición adicional de que las
funciones 𝝁𝟏 (𝒙) 𝒚 𝝁𝟐 (𝒙), satisfacen 𝒚𝟏 𝝁′𝟏 + 𝒚𝟐 𝝁′𝟐 = 𝟎. Esta suposición no se presenta por sorpresa,
sino que es resultado de los dos primeros términos de (4) puesto que si se requiere que 𝒚𝟏 𝝁′𝟏 +
𝒚𝟐 𝝁′𝟐 = 𝟎, entonces (4) se reduce a 𝒚′𝟏 𝝁′𝟏 + 𝒚′𝟐 𝝁′𝟐 = 𝒇(𝒙). Ahora tenemos nuestras dos ecuaciones
deseadas, a pesar de que sean dos ecuaciones para determinar las derivadas 𝝁′𝟏 y 𝝁′𝟐 . Por la regla
de Cramer, la solución del sistema
𝒚𝟏 𝝁′𝟏 + 𝒚𝟐 𝝁′𝟐 = 𝟎
𝒚′𝟏 𝝁′𝟏 + 𝒚′𝟐 𝝁′𝟐 = 𝒇(𝒙).

puede expresarse en términos de determinantes:

𝑾𝟏 𝒚𝟐 𝒇(𝒙) 𝑾𝟐 𝒚𝟏 𝒇(𝒙)
𝝁′𝟏 = =− 𝒚 𝝁′𝟐 = = (𝟓) 𝒅𝒐𝒏𝒅𝒆
𝑾 𝑾 𝑾 𝑾
𝒚𝟏 𝒚𝟐 𝟎 𝒚𝟐 𝒚𝟏 𝟎
𝑾 = |𝒚′ 𝒚′ |, 𝑾𝟏 = | ′ |, 𝑾𝟐 = | ′ | (𝟔)
𝟏 𝟐 𝒇(𝒙) 𝒚𝟐 𝒚𝟏 𝒇(𝒙)
Las funciones 𝝁′𝟏 y 𝝁′𝟐 . se encuentran integrando los resultados de (5). El determinante W se
reconoce como el Wronskiano de 𝒚𝟏 𝒚 𝒚𝟐 . Por la independencia lineal de 𝒚𝟏 𝒚 𝒚𝟐 en I, se sabe
que 𝑾(𝒚𝟏 (𝒙), 𝒚𝟐 (𝒙)) ≠ 𝟎 ∀ 𝑥 en el intervalo.

RESUMEN DEL MÉTODO

El procedimiento anterior es demasiado largo y complicado para usarse cada vez que se desee
resolver una ED. En este caso resulta más eficaz usar simplemente las fórmulas de (5). Así que
para resolver 𝒂𝟐 𝒚′′ + 𝒂𝟏 𝒚′ + 𝒂𝟎 𝒚 = 𝒈(𝒙), primero se encuentra la solución 𝒚𝒄 = 𝒄𝟏 𝒚𝟏 + 𝒄𝟐 𝒚𝟐 y luego
se calcula el 𝑾(𝒚𝟏 (𝒙), 𝒚𝟐 (𝒙)). Dividiendo entre 𝒂𝟐 , se escribe la ED en forma estándar 𝒚′′ + 𝑷𝒚′ +
𝑾 𝑾
𝑸𝒚 = 𝒇(𝒙) para determinar 𝒇(𝒙). Se encuentra 𝝁𝟏 𝒚 𝝁𝟐 integrando 𝝁′𝟏 = 𝑾𝟏 y 𝝁′𝟏 = 𝑾𝟐, donde 𝑊1 y
𝑊2 se definen como en (6). Una solución particular es 𝒚𝑷 = 𝝁𝟏 𝒚𝟏 + 𝝁𝟐 𝒚𝟐 . Entonces la solución
general de la ecuación es 𝒚 = 𝒚𝒄 + 𝒚𝒑

CONSTANTES DE INTEGRACIÓN

Cuando se calculan las integrales indefinidas de 𝝁𝟏 𝒚 𝝁𝟐, no es necesario introducir algunas

constantes. Esto es porque
𝒚 = 𝒚𝒄 + 𝒚𝒑 = 𝒄𝟏 𝒚𝟏 + 𝒄𝟐 𝒚𝟐 + (𝝁𝟏 + 𝒂𝟏 )𝒚𝟏 + (𝝁𝟐 + 𝒃𝟏 )𝒚𝟐
(𝑪𝟏 + 𝒂𝟏 )𝒚𝟏 + (𝑪𝟐 + 𝒃𝟏 )𝒚𝟐 + 𝝁𝟏 𝒚𝟏 + 𝝁𝟐 𝒚𝟐 = 𝑪𝟏 𝒚𝟏 + 𝑪𝟐 𝒚𝟐 + 𝝁𝟏 𝒚𝟏 + 𝝁𝟐 𝒚𝟐

ECUACIONES DE ORDEN SUPERIOR

El método que se describió para EDNH de 2do orden se puede generalizar a EDL de n-ésimo orden
que se han escrito en forma estándar

𝒚(𝒏) + 𝑷𝒏−𝟏 (𝒙)𝒚(𝒏−𝟏) + ⋯ + 𝑷𝟏 (𝒙)𝒚′ + 𝑷𝟎 (𝒙)𝒚 = 𝒇(𝒙) (9)

Si 𝒚𝒄 = 𝒄𝟏 𝒚𝟏 + 𝒄𝟐 𝒚𝟐 + ⋯ + 𝒄𝒏 𝒚𝒏 es la solución complementaria para (9), entonces una solución

particular es
𝒚𝒑 = 𝝁𝟏 (𝒙)𝒚𝟏 (𝒙) + 𝝁𝟐 (𝒙)𝒚𝟐 (𝒙) + ⋯ + 𝝁𝒏 (𝒙)𝒚𝒏 (𝒙)
Donde los 𝝁′𝒌 , 𝒌 = 𝟏, 𝟐, … , 𝒏 son detreminados por las n ecuaciones
𝒚𝟏 𝝁′𝟏 + 𝒚𝟐 𝝁′𝟐 + ⋯ + 𝒚𝒏 𝝁′𝒏 = 𝟎
𝒚′𝟏 𝝁′𝟏 + 𝒚′𝟐 𝝁′𝟐 + ⋯ + 𝒚′𝒏 𝝁′𝒏 = 𝟎
𝒚′′ ′ ′′ ′ ′′ ′
𝟏 𝝁𝟏 + 𝒚𝟐 𝝁𝟐 + ⋯ + 𝒚𝒏 𝝁𝒏 = 𝟎
∙ ∙ (10)
∙ ∙
∙ ∙
(𝒏−𝟏) ′ (𝒏−𝟏) ′ (𝒏−𝟏) ′
𝒚𝟏 𝝁𝟏 + 𝒚𝟐 𝝁𝟐 + ⋯ + 𝒚𝒏 𝝁𝒏 = 𝒇(𝒙)

Las 1ras 𝒏 − 𝟏 ecuaciones de este sistema, al igual que 𝒚𝟏 𝝁′𝟏 + 𝒚𝟐 𝝁′𝟐 = 𝟎 en (4), son suposiciones
que se hacen para simplificar la ecuación resultante después de que 𝒚𝒑 = 𝝁𝟏 𝒚𝟏 + ⋯ + 𝝁𝒏 𝒚𝒏 se
sustituye en (9). En este caso usando la regla de Cramer se obtiene

𝑾𝒌
𝝁′𝒌 = , 𝒌 = 𝟏, 𝟐, … , 𝒏
𝑾
donde 𝑾 es el Wronskiano de 𝒚𝟏 , 𝒚𝟐 , …. , 𝒚𝒏 𝑦 𝑾𝒌 es el determinante que se obtiene al
remplazar la k-ésima columna del Wronskiano por la columna formada por el lado derecho de
(10), es decir, la columna que consta de (𝟎, 𝟎, . . . , 𝒇(𝒙)). Cuando 𝒏 = 𝟐, se obtiene la ecuación (5).
Cuando 𝒏 = 𝟑, la solución particular 𝒚𝒑 = 𝝁𝟏 𝒚𝟏 + 𝝁𝟐 𝒚𝟐 + 𝝁𝟑 𝒚𝟑 , donde 𝒚𝟏 , 𝒚𝟐 , 𝒚𝟑 constituyen un
conjunto linealmente independiente de soluciones de la EDH asociada y 𝝁𝟏 , 𝝁𝟐 , 𝝁𝟑 se determinan
a partir de
𝑾𝟏 𝑾𝟐 𝑾𝟑
𝝁′𝟏 = , 𝝁′𝟐 = , 𝝁′𝟑 = (𝟏𝟏)
𝑾 𝑾 𝑾

𝟎 𝒚𝟐 𝒚𝟑 𝒚𝟏 𝟎 𝒚𝟑 𝒚𝟏 𝒚𝟐 𝟎 𝒚𝟏 𝒚𝟐 𝒚𝟑
𝑾𝟏 = | 𝟎 𝒚′𝟐 𝒚′𝟑 |, 𝑾𝟐 = | 𝒚′𝟏 𝟎 𝒚′𝟑 |, 𝑾𝟑 = | 𝒚′𝟏 𝒚′𝟐 𝟎 | 𝒚 𝑾 = | 𝒚′𝟏 𝒚′𝟐 𝒚′𝟑 |
𝒇(𝒙) 𝒚′′
𝟐 𝒚′′
𝟑 𝒚′′
𝟏 𝒇(𝒙) 𝒚′′
𝟑 𝒚′′
𝟏 𝒚′′
𝟐 𝒇(𝒙) 𝒚′′
𝟏 𝒚′′
𝟐 𝒚′′
𝟑

Ver problemas 25 y 26 del ejercicio 4.6

EJERCICIO – 3
ECUACIÓN DE CAUCHY-EULER

INTRODUCCIÓN

La relativa facilidad con que pudimos encontrar soluciones explícitas de ecuaciones lineales de
orden superior con coeficientes constantes en las secciones anteriores, en general no se realiza
en ecuaciones lineales con coeficientes variables. En el capítulo 6 veremos que cuando una ED
lineal tiene coeficientes variables, lo mejor que podemos esperar, usualmente, es encontrar una
solución en forma de serie infinita. Sin embargo, el tipo de ecuación diferencial que consideramos
en esta sección es una excepción a esta regla; esta es una ecuación lineal con coeficientes
variables cuya solución general siempre se puede expresar en términos de potencias de x, senos,
cosenos y funciones logarítmicas. Además este método de solució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 Una EDL de la forma

𝒏
𝒅𝒏 𝒚 𝒏−𝟏
𝒅𝒏−𝟏 𝒚 𝒅𝒚
𝒂𝒏 𝒙 + 𝒂𝒏−𝟏 𝒙 + ⋯ + 𝒂𝟏 𝒙 + 𝒂𝟎 𝒚 = 𝒈(𝒙)
𝒅𝒙𝒏 𝒅𝒙𝒏−𝟏 𝒅𝒙
donde los coeficientes 𝒂𝒏 , 𝒂𝒏−𝟏 , … , 𝒂𝟎 son cttes, se conoce como ecuación de Cauchy-Euler
(ECE). La característica observable de este tipo de ecuación es que el grado 𝒌 = 𝒏, 𝒏 − 𝟏, …. , 𝟏, 𝟎
de los coeficientes monomiales 𝒙𝒌 coincide con el orden k de la derivación 𝒅𝒌 𝒚/𝒅𝒙𝒌 :

Mismo Mismo

𝒅𝒏 𝒚 𝒅𝒏−𝟏 𝒚
𝒂𝒏 𝒙𝒏 + 𝒂𝒏−𝟏 𝒙 𝒏−𝟏
+⋯
𝒅𝒙𝒏 𝒅𝒙𝒏−𝟏
Iniciamos el análisis con un examen detallado de las formas de las soluciones generales de la
EDLH de 2do orden
𝒅𝟐 𝒚 𝒅𝒚
𝒂𝒙𝟐 𝟐
+ 𝒃𝒙 + 𝒄𝒚 = 𝟎
𝒅𝒙 𝒅𝒙

La solución de ecuaciones de orden superior se deduce de manera análoga. También, podemos

resolver la ecuación no homogénea ax2y_ _ bxy_ _ cy _ g(x) por variación de parámetros, una
vez que se ha determinado la función complementaria yc.
NOTA El coefi ciente ax2 de y_ es cero en x _ 0. Por lo que, para garantizar que los resultados
fundamentales del teorema 4.1.1 sean aplicables a la ecuación de Cauchy-Euler, centramos
nuestra atención en encontrar soluciones generales defi nidas en el intervalo (0, _). Las soluciones
en el intervalo (__, 0) se obtienen al sustituir t _ _x en la ecuación diferencial. Véanse los
problemas 37 y 38 de los ejercicios 4.7.
MÉTODO DE SOLUCIÓN Se prueba una solución de la forma y _ xm, donde m es un valor que se
debe determinar. Análogo a lo que sucede cuando se sustituye emx en una ecuación lineal con
coefi cientes constantes, cuando se sustituye xm, cada término de una ecuación de Cauchy-Euler
se convierte en un polinomio en m veces xm, puesto que

Por ejemplo, cuando sustituimos y _ xm, la ecuación de segundo orden se transforma en

Así y _ xm es una solución de la ecuación diferencial siempre que m sea una solución de la
ecuación auxiliar

Hay tres casos distintos a considerar que dependen de si las raíces de esta ecuación cuadrática
son reales y distintas, reales e iguales o complejas. En el último caso las raíces aparecen como un
par conjugado.
CASO I: RAÍCES REALES Y DISTINTAS Sean m1 y m2 las raíces reales de (1), tales que m1 _ m2.
Entonces y1 xm1 y y2 xm2 forman un conjunto fundamental de soluciones. Por tanto, la solución
general es
CASO II: RAÍCES REALES REPETIDAS Si las raíces de (l) son repetidas (es decir, m1 _ m2),
entonces se obtiene sólo una solución particular, y _ xm1. Cuando las raíces de la ecuación
cuadrática am2 _ (b _ a)m _ c _ 0 son iguales, el discriminante de los coefi cientes necesariamente
es cero. De la fórmula cuadrática se deduce que las raíces
deben ser m1 _ _(b _ a)_2a.
Ahora se puede construir una segunda solución y2, con la ecuación (5) de la sección 4.2. Primero
se escribe la ecuación de Cauchy-Euler en la forma estándar

y haciendo las identifi caciones P(x) _ b_ax y (b ax) dx (b a) ln x . Así

La solución general es entonces

CASO III: RAÍCES COMPLEJAS CONJUGADAS Si las raíces de (1) son el par conjugado m1 _ a _
ib, m2 _ a _ ib, donde a y b 0 son reales, entonces una solución es

Pero cuando las raíces de la ecuación auxiliar son complejas, como en el caso de las ecuaciones
con coeficientes constantes, se desea escribir la solución sólo en términos de funciones reales.
Observemos la identidad

que, por la fórmula de Euler, es lo mismo que

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

también son soluciones. Como W(xa cos(b ln x), xa sen(b ln x)) _ bx2a_1 _ 0, b 0
en el intervalo (0, _), se concluye que

constituyen un conjunto fundamental de soluciones reales de la ecuación diferencial.

Así la solución general es

REDUCCIÓN A COEFICIENTES CONSTANTES

Las similitudes entre las formas de soluciones de ecuaciones de Cauchy-Euler y soluciones de
ecuaciones lineales con coefi cientes constantes no sólo son una coincidencia. Por ejemplo, cuando
las raíces de las ecuaciones auxiliares para ay_ _ by_ _ cy _ 0 y ax2y_ _ bxy_ _ cy _ 0 son distintas
y reales, las soluciones generales respectivas son
Usando la identidad eln x _ x, x 0, la segunda solución dada en (5) puede expresarse en la
misma forma que la primera solución:

donde t _ ln x. Este último resultado ilustra el hecho de que cualquier ecuación de Cauchy-Euler
siempre se puede escribir de nuevo como una ecuación diferencial lineal con coefi cientes
constantes sustituyendo x _ et. La idea es resolver la nueva ecuación diferencial en términos de
la variable t, usando los métodos de las secciones anteriores y una vez obtenida la solución
general, sustituir nuevamente t _ ln x. Este método, que se ilustró en el último ejemplo, requiere
el uso de la regla de la cadena de la derivación.

EJERCICIO – 4

SOLUCIÓN DE SISTEMAS DE ED LINEALES POR ELIMINACIÓN

INTRODUCCIÓN
Las ecuaciones diferenciales ordinarias simultáneas tienen que ver con dos o más ecuaciones que
contienen derivadas de dos o más variables dependientes (las funciones desconocidas) respecto
a una sola variable independiente. El método de eliminación sistemática para resolver sistemas
de ecuaciones diferenciales con coefi cientes constantes se basa en el principio algebraico de
eliminación de variables. Veremos que la operación análoga de multiplicar una ecuación algebraica
por una constante es operar en una EDO con cierta combinación de derivadas.

ELIMINACIÓN SISTEMÁTICA La eliminación de una incógnita en un sistema de ecuaciones

diferenciales lineales se facilita al rescribir cada ecuación del sistema en notación de operador
diferencial. Recuerde de la sección 4.1 que una sola ecuación lineal

donde las ai, i _ 0, 1, . . . , n son constantes, puede escribirse como

Si el operador diferencial de n-ésimo orden anDn an 1D(n 1) a1D a0 se factoriza en operadores

diferenciales de menor orden, entonces los factores conmutan.
Ahora, por ejemplo, para rescribir el sistema

en términos del operador D, primero se escriben los términos con variables dependientes en un
miembro y se agrupan las mismas variables.

SOLUCIÓN DE UN SISTEMA Una solución de un sistema de ecuaciones diferenciales es un

conjunto de funciones suficientemente derivables x _ f1(t), y _ f2(t), z _ f3(t), etcétera, que
satisface cada ecuación del sistema en algún intervalo común I.
MÉTODO DE SOLUCIÓN Considere el sistema simple de ecuaciones lineales de primer orden

Operando con D la primera ecuación de (1) en tanto que la segunda se multiplica por _ 3 y después
se suma para eliminar y del sistema, se obtiene D2x _ 6x _ 0. Puesto que las raíces de la ecuación
auxiliar de la última ED son m1 16 y m2 16 , se obtiene

Multiplicando la primera ecuación en (1) por 2 mientras que se opera la segunda con D y después
restando, se obtiene la ecuación diferencial para y, D2y _ 6y _ 0.
Inmediatamente se tiene que
Ahora (2) y (3) no satisfacen el sistema (1) para toda elección de c1, c2, c3 y c4 porque el sistema
en sí pone una restricción al número de parámetros en una solución que se puede elegir en forma
arbitraria. Para ver esto, observe que sustituyendo x(t) y y(t) en la primera ecuación del sistema
original (1), después de simplifi car, se obtiene

Puesto que la última expresión es cero para todos los valores de t, debemos tener 16c1 3c3 0 y
16c2 3c4 0. Estas dos ecuaciones nos permiten escribir c3 como un múltiplo de c1 y c4 como un
múltiplo de c2:

Por tanto se concluye que una solución del sistema debe ser

Se recomienda sustituir (2) y (3) en la segunda ecuación de (1) y comprobar que se cumple la
misma relación (4) entre las constantes.

EJERCICIO – 5