UNIVERSIDAD TECNOLOGICA

DE LOS ANDES

ESCUELA PROFESIONAL DE INGENIERIA CIVIL

FILIAL CUSCO

ECUACIONES DIFERENCIALES

Alfredo Julio Escalante Auccapure

Cusco Perú
IC ECUACIONES DIFERENCIALES UTEA

4.6.3 Generalización del método de variación de parámetros .

Dada la ecuación diferencial lineal no homogénea, escrita en su forma canónica

𝑦 (𝑛) + 𝑎𝑛−1 (𝑥)𝑦 (𝑛−1) + ⋯ + 𝑎2 (𝑥)𝑦 ′′ + 𝑎1 (𝑥)𝑦 ′ + 𝑎0 (𝑥)𝑦 = 𝑓(𝑥) … (4.6.3.1)

Generalizando el proceso de variación de parámetros que se hizo para una ecuación

diferencial lineal no homogénea de segundo orden, se propone los siguientes pasos:

a) Determinar la solución homogénea 𝑦𝑐 .

Hallamos 𝑦1 , 𝑦2 , … , 𝑦𝑛 soluciones linealmente independientes de la ecuación homogénea
asociada, es decir:

𝑦 (𝑛) + 𝑎𝑛−1 (𝑥)𝑦 (𝑛−1) + ⋯ + 𝑎2 (𝑥)𝑦 ′′ + 𝑎1 (𝑥)𝑦 ′ + 𝑎0 (𝑥)𝑦 = 0

Entonces la solución homogénea o complementaria es

𝑦𝑐 = 𝑐1 𝑦1 + 𝑐2 𝑦2 + ⋯ + 𝑐𝑛 𝑦𝑛

b) Hallamos el wronskiano.
𝑦1 𝑦2 ⋯ 𝑦𝑖 ⋯ 𝑦𝑛
𝑦1′ 𝑦2′ ⋯ 𝑦𝑖′ ⋯ 𝑦𝑛′
𝑊=| ⋮ ⋮ ⋮ ⋮ |
(𝑛−1) (𝑛−1) (𝑛−1) (𝑛−1)
𝑦1 𝑦2 ⋯ 𝑦𝑖 ⋯ 𝑦𝑛

Si reemplazamos la 𝑖-ésima columna del wronskiano 𝑊 con la columna

0
0
⋮
0
(1)

obtendremos las sub-determinantes 𝑊𝑖 .

c) Calculamos los valores de 𝑢𝑖 con la fórmula:

𝑊𝑖
𝑢𝑖 = ∫ 𝑓(𝑥) 𝑑𝑥, ∀𝑖 = 1,2, … , 𝑛
𝑊

d) Hallamos la solución particular 𝑦𝑝 .

Por el método de variación de parámetros la solución particular será:

𝑦𝑝 = 𝑢1 𝑦1 + 𝑢2 𝑦2 + ⋯ + 𝑢𝑛 𝑦𝑛

e) Finalmente escribimos la solución general

AJEA 2
IC ECUACIONES DIFERENCIALES UTEA

𝑦 = 𝑐1 𝑦1 + 𝑐2 𝑦2 + ⋯ + 𝑐𝑛 𝑦𝑛 + 𝑢1 𝑦1 + 𝑢2 𝑦2 + ⋯ + 𝑢𝑛 𝑦𝑛

Aplicación:
16
Resolver: (𝐷 + 1)3 𝑦 = (2𝑥+3)𝑒 𝑥
4.7 Método de operadores inversos.
Operadores.
En esta parte nos referimos específicamente al operador diferencial

𝐿(𝐷) = 𝑎𝑛 𝐷𝑛 + 𝑎𝑛−1 𝐷(𝑛−1) + ⋯ + 𝑎1 𝐷 + 𝑎0

y a las propiedades de algunos de sus casos particulares.

Propiedades.
1°) Si la función 𝑓(𝑥) es de clase 𝑛 en todo el intervalo 𝐼, es decir 𝑓(𝑥) ∈ 𝐶 𝑛 (𝐼) y 𝑎 ∈ ℝ (ó
𝑎 ∈ ℂ) entonces:

i) 𝐷𝑛 (𝑒 𝑎𝑥 𝑓(𝑥)) = 𝑒 𝑎𝑥 (𝐷 + 𝑎)𝑛 𝑓(𝑥)

ii) 𝐷𝑛 (𝑒 𝑎𝑥 ) = 𝑎𝑛 𝑒 𝑎𝑥

2°) El siguiente teorema, llamado teorema básico de los operadores, nos permite sacar una
exponencial que está dentro de un operador.

Teorema 1.
Si la función 𝑓(𝑥) es de clase 𝑛 en todo el intervalo 𝐼, es decir 𝑓(𝑥) ∈ 𝐶 𝑛 (𝐼), 𝐿(𝐷) es un
operador diferencial lineal y 𝑎 ∈ ℝ (ó 𝑎 ∈ ℂ), entonces:
i) 𝐿(𝐷)𝑒 𝑎𝑥 = 𝐿(𝑎)𝑒 𝑎𝑥
ii) 𝐿(𝐷)(𝑒 𝑎𝑥 𝑓(𝑥)) = 𝑒 𝑎𝑥 𝐿(𝐷 + 𝑎)𝑓(𝑥)

Nota: Para introducir una función exponencial dentro de un operador, utilizamos la siguiente
formula

𝑒 𝑎𝑥 𝐿(𝐷)𝑓(𝑥) = 𝐿(𝐷 − 𝑎)(𝑒 𝑎𝑥 𝑓(𝑥))

Aplicación.
1. Demostrar que: (𝐷2 + 4)(𝐷 + 5)6 (𝑥 5 𝑒 −5𝑥 ) = 0

2. Demostrar que: (𝐷 + 5)𝑛 𝑥 𝑛 𝑒 −5𝑥 = 𝑛! 𝑒 −5𝑥

∴ (𝐷 + 5)𝑛 𝑥 𝑛 𝑒 −5𝑥 = 𝑛! 𝑒 −5𝑥

Operadores inversos.
En la ecuación diferencial lineal de coeficientes constantes
AJEA 3
IC ECUACIONES DIFERENCIALES UTEA

𝐿(𝐷)𝑦 = 𝑓(𝑥)

donde 𝐿(𝐷) = 𝑎𝑛 𝐷𝑛 + 𝑎𝑛−1 𝐷𝑛−1 + ⋯ + 𝑎1 𝐷 + 𝑎0 se denomina operador diferencial lineal

de 𝑛-ésimo orden de coeficientes constantes.

Definición.- Dada la ecuación diferencial lineal de coeficientes constantes 𝐿(𝐷)𝑦 = 𝑓(𝑥),

1
definimos el operador inverso de 𝐿(𝐷), como el operador 𝐿−1 (𝐷) = 𝐿(𝐷) tal que 𝐿−1 (𝐷)𝑓(𝑥)
es una solución particular de la ecuación diferencial dada, es decir:

1
𝑦𝑝 = 𝐿−1 (𝐷)𝑓(𝑥) = [𝑓(𝑥)]
𝐿(𝐷)

Teorema 2.
Si 𝐿(𝐷) = (𝐷 − 𝑎)𝑛 y 𝑓(𝑥) es una función continua, entonces una solución particular de la
ecuación diferencial (𝐷 − 𝑎)𝑛 𝑦 = 𝑓(𝑥) es

𝑦𝑝 = 𝑒 𝑎𝑥 ∫ (𝑒 𝑎𝑥 ∫ ⋯ (𝑒 𝑎𝑥 ∫ 𝑒 −𝑎𝑥 𝑓(𝑥) 𝑑𝑥) ⋯ 𝑑𝑥) 𝑑𝑥

⏟
𝑛 𝑣𝑒𝑐𝑒𝑠

Observaciones:
1
De la formula 𝑦𝑝 = 𝐿(𝐷) [𝑓(𝑥)] y el teorema último podemos afirmar lo siguiente:

i) Si 𝐿(𝐷) = 𝐷 − 𝑎 entonces:
1
𝑦𝑝 = [𝑓(𝑥)] = 𝑒 𝑎𝑥 ∫ 𝑒 −𝑎𝑥 𝑓(𝑥)𝑑𝑥
𝐷−𝑎
ii) Si 𝐿(𝐷) = (𝐷 − 𝑎)𝑛 entonces:

1
𝑦𝑝 = [𝑓(𝑥)] = 𝑒 𝑎𝑥 ∫ (𝑒 𝑎𝑥 ∫ ⋯ (𝑒 𝑎𝑥 ∫ 𝑒 −𝑎𝑥 𝑓(𝑥) 𝑑𝑥) ⋯ 𝑑𝑥) 𝑑𝑥
(𝐷 − 𝑎)𝑛

iii) Si 𝐿(𝐷) = (𝐷 − 𝑎1 )(𝐷 − 𝑎2 ) ⋯ (𝐷 − 𝑎𝑛 ) entonces:

𝑦𝑝 = 𝑒 𝑎𝑛 𝑥 ∫ 𝑒 −𝑎𝑛𝑥 ⋯ (𝑒 𝑎2 𝑥 ∫ 𝑒 −𝑎2 𝑥 (𝑒 𝑎1 𝑥 ∫ 𝑒 −𝑎1 𝑥 𝑓(𝑥)𝑑𝑥) 𝑑𝑥) ⋯ 𝑑𝑥

1 𝐵 2 𝐵 𝑛 𝐵
iv) Si 𝐿(𝐷) = 𝐷−𝑎 + 𝐷−𝑎 + ⋯ + 𝐷−𝑎 entonces:
1 2 𝑛

1 1 1
𝑦𝑝 = 𝐵1 [𝑓(𝑥)] + 𝐵2 [𝑓(𝑥)] + ⋯ + 𝐵𝑛 [𝑓(𝑥)]
𝐷 − 𝑎1 𝐷 − 𝑎2 𝐷 − 𝑎𝑛

= 𝐵1 𝑒 𝑎1 𝑥 ∫ 𝑒 −𝑎1 𝑥 𝑓(𝑥)𝑑𝑥 + 𝐵2 𝑒 𝑎2 𝑥 ∫ 𝑒 −𝑎2 𝑥 𝑓(𝑥)𝑑𝑥 + ⋯ + 𝐵𝑛 𝑒 𝑎𝑛 𝑥 ∫ 𝑒 −𝑎𝑛𝑥 𝑓(𝑥)𝑑𝑥

AJEA 4
IC ECUACIONES DIFERENCIALES UTEA

Aplicación.
3. Determinar la solución particular de: 𝐷(𝐷 − 2)(𝐷 + 3)𝑦 = 𝑒 −3𝑥

No sorprenderse con las dos soluciones particulares diferentes que se encontraron para la
misma ecuación diferencial ya que tiene muchas soluciones particulares.

Teorema 3.
Si la función 𝑓(𝑥) es de clase 𝑛 en todo el intervalo 𝐼, es decir 𝑓(𝑥) ∈ 𝐶 𝑛 (𝐼) y 𝑎 ∈ ℝ (ó 𝑎 ∈
ℂ) entonces una solución particular de la ecuación diferencial:
𝑒 𝑎𝑥
a) 𝐿(𝐷)𝑦 = 𝑒 𝑎𝑥 será 𝑦𝑝 = 𝐿(𝑎) , si 𝐿(𝑎) ≠ 0

1
b) 𝐿(𝐷)𝑦 = 𝑒 𝑎𝑥 𝑓(𝑥) será 𝑦𝑝 = 𝑒 𝑎𝑥 𝐿(𝐷+𝑎) [𝑓(𝑥)]
Observación.
Estas se refieren a la parte a) del Teorema 3, y consisten en:

1°) Si 𝐿(𝑎) = 0, por el teorema del resto en algebra se sabe que 𝐿(𝐷) = (𝐷 − 𝑎)𝐿1 (𝐷)
entonces la solución particular será
1 𝑒 𝑎𝑥
𝑦𝑝 = ( )
𝐷 − 𝑎 𝐿1 (𝑎)
siempre que 𝐿1 (𝑎) ≠ 0.

2°) Si 𝐿(𝐷) = (𝐷 − 𝑎)𝑘 , entonces es fácil demostrar que la solución particular será

1 𝑎𝑥
𝑥 𝑘 𝑎𝑥
𝑦𝑝 = 𝑒 = 𝑒
(𝐷 − 𝑎)𝑘 𝑘!
Aplicación.
4. Determinar la solución general de: (𝐷 + 7)3 𝑦 = (𝑥 2 + 3𝑥)𝑒 −7𝑥

Teorema 4.
Si 𝐿(𝐷) = 𝐷 + 𝑎, 𝑎 ≠ 0 y 𝑓(𝑥) = 𝑃𝑛 (𝑥) es una función polinomial de grado 𝑛, entonces

1 1 𝐷 𝐷2 𝐷3 𝐷𝑛
[𝑓(𝑥)] = [1 − + 2 − 3 + ⋯ + (−1)𝑛 𝑛 ] 𝑓(𝑥)
𝐷+𝑎 𝑎 𝑎 𝑎 𝑎 𝑎

Aplicación.
5. Determinar la solución particular de: (𝐷 − 1)𝑦 = (𝑥 4 + 3)𝑒 5𝑥

Teorema 5.
Si 𝐿(𝐷) = 𝑎𝑛 𝐷𝑛 + 𝑎𝑛−1 𝐷(𝑛−1) + ⋯ + 𝑎1 𝐷 + 𝑎0 es un operador diferencial lineal de
coeficientes constantes de orden 𝑛 y 𝑓(𝑥) es un polinomio de grado 𝑘, entonces una solución
particular de la ecuación diferencial 𝐿(𝐷)𝑦 = 𝑓(𝑥) es de la forma:

AJEA 5
IC ECUACIONES DIFERENCIALES UTEA

1
𝑦𝑝 = [𝑓(𝑥)] = (𝑏0 + 𝑏1 𝐷 + 𝑏2 𝐷2 + 𝑏3 𝐷3 + ⋯ + 𝑏𝑟 𝐷𝑟 )𝑓(𝑥); 𝑟 ≤ 𝑘
𝐿(𝐷)
donde
1 1
= = 𝑏0 + 𝑏1 𝐷 + 𝑏2 𝐷2 + 𝑏3 𝐷3 + ⋯ + 𝑏𝑟 𝐷𝑟 + ⋯
𝐿(𝐷) 𝑎0 + 𝑎1 𝐷 + 𝑎2 𝐷2 + ⋯ + 𝑎𝑛 𝐷𝑛

Aplicación.
𝑑4 𝑦 𝑑2 𝑦
6. Determinar la solución general de: 𝑑𝑥 4 + 5 𝑑𝑥 2 + 4𝑦 = 𝑥 5 − 5𝑥 + 2

Lema.
Si 𝑎 y 𝐿(−𝑎2 ) son números reales, entonces:
i) 𝐷2𝑛 𝑠𝑒𝑛(𝑎𝑥 + 𝑏) = (−𝑎2 )𝑛 𝑠𝑒𝑛(𝑎𝑥 + 𝑏)
ii) 𝐷2𝑛 𝑐𝑜𝑠(𝑎𝑥 + 𝑏) = (−𝑎2 )𝑛 𝑐𝑜𝑠(𝑎𝑥 + 𝑏)

Teorema 6.
1
Si 𝑎, 𝐿(−𝑎2 ) y 𝐿(−𝑎2 ) son números reales, entonces:

a) 𝐿(𝐷2 )𝑠𝑒𝑛(𝑎𝑥 + 𝑏) = 𝐿(−𝑎2 )𝑠𝑒𝑛(𝑎𝑥 + 𝑏)

b) 𝐿(𝐷2 )𝑐𝑜𝑠(𝑎𝑥 + 𝑏) = 𝐿(−𝑎2 )𝑐𝑜𝑠(𝑎𝑥 + 𝑏)

1 𝑠𝑒𝑛(𝑎𝑥+𝑏)
c) 𝐿(𝐷2 ) 𝑠𝑒𝑛(𝑎𝑥 + 𝑏) = , siempre que 𝐿(−𝑎2 ) ≠ 0
𝐿(−𝑎2 )

1 𝑐𝑜𝑠(𝑎𝑥+𝑏)
d) 𝐿(𝐷2 ) 𝑐𝑜𝑠(𝑎𝑥 + 𝑏) = , siempre que 𝐿(−𝑎2 ) ≠ 0.
𝐿(−𝑎2 )

Aplicación.
𝑑4 𝑦 𝑑2 𝑦
7. Determinar la solución particular de: 𝑑𝑥 4 + 5 𝑑𝑥 2 + 7𝑦 = 𝑠𝑒𝑛𝑥 𝑐𝑜𝑠𝑥

Observación.
El operador diferencial lineal de orden 𝑛 con coeficientes constantes se puede escribir de la
forma 𝐿(𝐷) = 𝐿1 (𝐷2 ) + 𝐷𝐿2 (𝐷2 )

Teorema 7.
Si 𝑎, 𝐿1 (𝐷2 ) y 𝐿2 (𝐷2 ) son números reales, entonces:

a) 𝐿(𝐷)𝑠𝑒𝑛(𝑎𝑥 + 𝑏) = 𝐿1 (−𝑎2 )𝑠𝑒𝑛(𝑎𝑥 + 𝑏) + 𝐷𝐿2 (−𝑎2 )𝑠𝑒𝑛(𝑎𝑥 + 𝑏)

b) 𝐿(𝐷)𝑐𝑜𝑠(𝑎𝑥 + 𝑏) = 𝐿1 (−𝑎2 )𝑐𝑜𝑠(𝑎𝑥 + 𝑏) + 𝐷𝐿2 (−𝑎2 )𝑐𝑜𝑠(𝑎𝑥 + 𝑏)

1 1 1
c) 𝐿(𝐷) 𝑠𝑒𝑛(𝑎𝑥 + 𝑏) = 𝐿 2 )+𝐷𝐿 2) 𝑠𝑒𝑛(𝑎𝑥 + 𝑏) = 𝐿 2 )+𝐷𝐿 2) 𝑠𝑒𝑛(𝑎𝑥 + 𝑏)
1 (𝐷 2 (𝐷 1 (−𝑎 2 (−𝑎

AJEA 6
IC ECUACIONES DIFERENCIALES UTEA

1 1 1
d) 𝐿(𝐷) 𝑐𝑜𝑠(𝑎𝑥 + 𝑏) = 𝐿 (𝐷 2 )+𝐷𝐿 (𝐷 2 )
𝑐𝑜𝑠(𝑎𝑥 + 𝑏) = 𝐿 (−𝑎 2 )+𝐷𝐿 (−𝑎2 )
𝑐𝑜𝑠(𝑎𝑥 + 𝑏)
1 2 1 2

Aplicación.
8. Determinar la solución particular de: (𝐷3 + 5𝐷2 + 𝐷 + 5)𝑦 = 𝑒 −𝑥 𝑠𝑒𝑛3𝑥

Aplicación.
9. Hallar la solución particular de la ecuación diferencial
(𝐷2 + 𝑎2 )𝑛 𝑦 = 𝑐𝑜𝑠(𝑎𝑥) + 𝑖𝑠𝑒𝑛(𝑎𝑥)

1 𝑥𝑛 𝑛𝜋
[𝑐𝑜𝑠𝑎𝑥] = 𝑐𝑜𝑠 (𝑎𝑥 − 2 )
(𝐷2 +𝑎 )
2 𝑛 (2𝑎)𝑛 𝑛!
1 𝑥𝑛 𝑛𝜋
[𝑠𝑒𝑛𝑎𝑥] = 𝑠𝑒𝑛 (𝑎𝑥 − 2 )
(𝐷2 + 𝑎2 )𝑛 (2𝑎)𝑛 𝑛!

Observación.
Otra fórmula que nos abrevia el cálculo de la solución particular es:

1 𝐿′ (𝐷) 1
𝑦𝑝 = [𝑥𝑓(𝑥)] = [𝑥 − ]( 𝑓(𝑥))
𝐿(𝐷) 𝐿(𝐷) 𝐿(𝐷)

donde 𝑓(𝑥) es una función.

10. Determinar la solución particular de: (𝐷2 + 8𝐷 + 32)4 𝑦 = 𝑒 −4𝑥 𝑠𝑒𝑛4𝑥

Nota:
El método de los operadores inversos también sirve para calcular integrales.

11. Calcular la integral ∫ 𝑥 5 𝑒 3𝑥 𝑑𝑥

12. Calcular la integral ∫ 𝑥 3 𝑐𝑜𝑠2𝑥𝑑𝑥

13. Calcular la integral ∫ 𝑒 3𝑥 𝑠𝑒𝑛4𝑥𝑑𝑥

14. Resolver la ecuación diferencial: (𝐷2 − 4𝐷 + 3)𝑦 = 2𝑥𝑒 3𝑥 + 3𝑒 𝑥 𝑐𝑜𝑠2𝑥

15. Resolver la ecuación diferencial: (𝐷2 + 𝐷)(𝐷 − 2)3 𝑦 = 𝑒 2𝑥 + 𝑒 −𝑥 𝑠𝑒𝑛3𝑥

16. Resolver: 𝐷2 (𝐷4 − 4𝐷3 + 6𝐷2 − 4𝐷 + 1)𝑦 = (𝑥 2 + 1)(1 − 𝑒 𝑥 )

17. Resolver la ecuación diferencial: 𝑦 𝑣𝑖𝑖𝑖 − 2𝑦 𝑖𝑣 + 𝑦 = 𝑐𝑜𝑠4𝑥 + 𝑠𝑒𝑛6𝑥

4.8 Ecuación diferencial ordinaria de Euler-Cauchy.

Definición.- La ecuación diferencial lineal

AJEA 7
IC ECUACIONES DIFERENCIALES UTEA

𝑑𝑛 𝑦 𝑑𝑛−1 𝑦 𝑑2𝑦 𝑑𝑦
𝑎𝑛 𝑥 𝑛 𝑛
+ 𝑎 𝑛−1 𝑥 𝑛−1
𝑛−1
+ ⋯ + 𝑎 2 𝑥 2
2
+ 𝑎1 𝑥 + 𝑎0 𝑦 = 𝑓(𝑥) … (𝜑)
𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑑𝑥
donde 𝑎𝑛 ≠ 0, … , 𝑎1 , 𝑎0 son constantes y 𝑓(𝑥) es continua, se llama ecuación diferencial
lineal de Euler-Cauchy.

Con la sustitución 𝑥 = 𝑒 𝑤 ó 𝑤 = 𝑙𝑛𝑥 se convierte la ecuación diferencial de Euler-Cauchy

(que por cierto es de coeficientes variables) en una ecuación diferencial lineal de coeficientes
constantes.

Se prueba que

𝑑𝑛 𝑦
𝑥𝑛 = 𝐷𝑤 (𝐷𝑤 − 1)(𝐷𝑤 − 2) ⋯ (𝐷𝑤 − (𝑛 − 1))𝑦
𝑑𝑥 𝑛

Aplicando la formula demostrada a la ecuación (𝜑), esta queda como

[𝑎𝑛 𝐷𝑤 (𝐷𝑤 − 1)(𝐷𝑤 − 2) ⋯ (𝐷𝑤 − (𝑛 − 1)) + ⋯ + 𝑎2 𝐷𝑤 (𝐷𝑤 − 1) + 𝑎1 𝐷𝑤 + 𝑎0 ]𝑦

= 𝑓(𝑒 𝑤 )

la ecuación diferencial lineal de coeficientes constantes no homogénea.

La forma más general de la ecuación diferencial de Euler-Cauchy es la ecuación diferencial

Euler-Legendre.

Definición.- La ecuación diferencial lineal

𝑑𝑛 𝑦 𝑑 𝑛−1 𝑦 𝑑2 𝑦 𝑑𝑦
𝑎𝑛 (𝑝𝑥 + 𝑞)𝑛 𝑛
+ 𝑎 𝑛−1 (𝑝𝑥 + 𝑞) 𝑛−1
𝑛−1
+ ⋯ + 𝑎 2 (𝑝𝑥 + 𝑞) 2
2
+ 𝑎1 (𝑝𝑥 + 𝑞)
𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑑𝑥
+ 𝑎0 𝑦 = 𝑓(𝑥)

donde 𝑎𝑛 ≠ 0, … , 𝑎1 , 𝑎0 son constantes y 𝑓(𝑥) es continua, se llama ecuación diferencial

lineal de Euler-Legendre.

La sustitución 𝑝𝑥 + 𝑞 = 𝑒 𝑤 ó 𝑤 = 𝑙𝑛(𝑝𝑥 + 𝑞) convierte la ecuación diferencial de Euler-

Legendre en una ecuación diferencial lineal de coeficientes constantes. Esta sustitución nos
permite obtener la formula
𝑑𝑛 𝑦
(𝑝𝑥 + 𝑞)𝑛 𝑛 = 𝑝𝑛 𝐷𝑤 (𝐷𝑤 − 1)(𝐷𝑤 − 2) ⋯ (𝐷𝑤 − (𝑛 − 1))𝑦, ∀𝑛 = 1,2, …
𝑑𝑥

AJEA 8