Cap 1 Ecuaciones Diferenciales 1er Orden

Capítulo 1 Ecuaciones Diferenciales de Primer Orden M.Sc. Ing.
Juan Carlos Quispe Apaza
CAPITULO 1
ECUACIONES DIFERENCIALES DE PRIMER ORDEN
𝒅𝒚 𝒅𝟑 𝒚 𝝏𝒇 𝝏𝟑 𝒇 𝝏𝒇
1.1 Ecuación Diferencial 𝒙𝟐 + 𝟐𝒙 − 𝒅𝒙 = 𝟖 + 𝒅𝒙𝟑 ; 𝒙𝟐 + 𝟐𝒙 − 𝝏𝒙 = 𝟖 + 𝝏𝒚𝟑 − 𝝏𝒛
Ecuación matemática que presenta a una función como incógnita (variable dependiente) y a
sus derivadas.
1.1.1 La ecuación diferencial ordinaria
Una ecuación diferencial ordinaria, es aquella ecuación diferencial en la que la función incógnita
(variable dependiente) solo depende de una variable independiente, tiene la forma:
𝐹(𝑥, 𝑦, 𝑦 ′ , 𝑦 ′′ , 𝑦 ′′′ , … … , 𝑦 (𝑛) ) = 0 … (1.1)
En una ecuación diferencial la incógnita a calcularse viene a ser la función 𝑦 = 𝑓 (𝑥 ) ?
𝑑3𝑦 𝑑2𝑦
+ 4𝑥 𝑙𝑛 𝑥 + 2 = cos 𝑥
𝑑𝑥 3 𝑑𝑥
𝑦′′′ + 4𝑥 𝑙𝑛 𝑥 + 𝑦′′ = cos 𝑥
Ejemplo 1: Dé una solución de la ecuación diferencial
𝑦 ′ − 𝑦 = 0 … (𝛼)
Tanteando al solución de (𝛼)
𝑦 = 𝑒𝑥 → 𝑦′ = 𝑒 𝑥
Entonces: 𝑦 ′ − 𝑦 = 0 → 𝑒 𝑥 − 𝑒 𝑥 = 0 → 0 = 0 verifica.
Por tanto 𝑦 = 𝑒 𝑥 es solución de (𝛼)
También verifica 𝑦 = 0 porque 𝑦 ′ = 0, entonces 𝑦 ′ − 𝑦 = 0 → 0 − 0 = 0 cumple.
Rpta.: 𝑦 = 𝑒𝑥
Prueba: 𝑦 ′ = 𝑒 𝑥 reemplazando en la ecuación diferencial (𝛼)
𝑦′ − 𝑦 = 0 → 𝑒𝑥 − 𝑒𝑥 = 0 → 0=0
Ejemplo 2: Dé una solución de la ecuación diferencial
𝑦′′ + 𝑦 = 0 … (𝛽)
Rpta.: 𝑦 = cos 𝑥
Capítulo 1 Ecuaciones Diferenciales de Primer Orden M.Sc. Ing. Juan Carlos Quispe Apaza
Prueba: 𝑦 ′ = − sin 𝑥 , 𝑦 ′′ = − cos 𝑥 𝑒𝑛 (𝛽)
𝑦 ′′ + 𝑦 = 0 → − cos 𝑥 + cos 𝑥 = 0 → 0=0
Sin embargo, también existe otra solución:
𝑦 = sin 𝑥 → 𝑦′ = cos 𝑥 𝑦 ′′ = − sin 𝑥 𝑒𝑛 (𝛽)

𝑦 ′′ + 𝑦 = 0 → − sin 𝑥 + sin 𝑥 = 0 → 0=0
En realidad, en el capítulo 3 se verá que si dos soluciones linealmente independientes son

solución de una ecuación diferencial de segundo orden también los será una combinación lineal
de las mismas, para el problema la solución general viene dada por:
𝑦 = 𝛼 sin 𝑥 + 𝛽 cos 𝑥 → 𝑦 ′ = 𝛼 cos 𝑥 − 𝛽 sin 𝑥 ; 𝑦 ′′ = −𝛼 sin 𝑥 − 𝛽 cos 𝑥 en (𝛽)
𝑦 ′′ + 𝑦 = 0 → −𝛼 sin 𝑥 − 𝛽 cos 𝑥 + 𝛼 sin 𝑥 + 𝛽 cos 𝑥 = 0
1.1.2 Ecuación diferencial en derivadas parciales
Ecuación en la que la función incógnita (variable dependiente) depende de más de una variable
independiente.
𝜕 2 𝑓 𝜕 2 𝑓 𝜕 2 𝑓 1 𝜕𝑓
+ + =
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2 𝑘 𝜕𝑡
1.2 Clasificación según el orden y el grado de una ecuación diferencial
El orden de una ecuación diferencial ordinaria es la derivada de mayor orden en la ecuación

diferencial.
El grado de una ecuación diferencial es el grado algebraico de la derivada de mayor orden en la

ecuación diferencial.
a) 𝑦 ′′′ + 3𝑥(𝑦 ′ )6 = 𝑒 𝑥 De orden 3 y grado 1

1
b) (𝑦 ′′ )3 + 𝑥 (𝑦 ′ )4 = (𝑦 ′′′ )2 De orden 3 y grado 2
4 5
𝑑5 𝑦 𝑑2 𝑦 𝑑𝑦
c) ( ) + sin 𝑡 ( ) = De orden 5 y grado 4
𝑑𝑡 5 𝑑𝑡 2 𝑑𝑡
d) (𝑦 ′ )2 + cos 𝑡 𝑦 = tan 𝑡 De orden 1 y grado 2
Problema 1 Determine el orden y el grado de la ecuación diferencial:
(|𝑥 − 2| + |𝑥 + 2|)𝑦 ′′′ + (|𝑥 − 2| − |𝑥 + 2|)𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥
Solución:
x  2 x 20 x  2 x 2

x 2  → x 2 
(x  2) x 20 2  x x 2
x  2 x 20 x  2 x  2

x 2  → x 2 
(x  2) x20 x  2 x  2
I1 I2 I3
-2 2
i) Para: −∞ < 𝒙 < −𝟐
(|𝑥 − 2| + |𝑥 + 2|)𝑦 ′′′ + (|𝑥 − 2| − |𝑥 + 2|)𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥

(2 − 𝑥 + (−𝑥 − 2))𝑦 ′′′ + (2 − 𝑥 − (−𝑥 − 2))𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥
−2𝑥𝑦 ′′′ + 4𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥 De tercer orden y grado 1
ii) Para: −𝟐 ≤ 𝒙 < 𝟐
(|𝑥 − 2| + |𝑥 + 2|)𝑦 ′′′ + (|𝑥 − 2| − |𝑥 + 2|)𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥

(2 − 𝑥 + (𝑥 + 2))𝑦 ′′′ + (2 − 𝑥 − (𝑥 + 2))𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥
4𝑦 ′′′ − 2𝑥𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥 De tercer orden y grado 1
iii) Para: 𝒙 ≥ 𝟐
(𝑥 − 2 + 𝑥 + 2)𝑦 ′′′ + (𝑥 − 2 − (𝑥 + 2))𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥
2𝑥𝑦 ′′′ − 4𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥 De tercer orden y grado 1
1.3 Solución de una ecuación diferencial
La función 𝑦(𝑥) es solución de la ecuación diferencial
𝐹(𝑥, 𝑦, 𝑦 ′ , 𝑦 ′′ , 𝑦 ′′′ , … … , 𝑦 (𝑛) ) = 0 … (1.2)
𝑦(𝑥) es solución si satisface la ecuación (1) en un intervalo “I” donde la misma 𝑦(𝑥) y sus
derivadas son continuas y diferenciables, ahora es vital que se distinga entre el dominio de la
función 𝑦(𝑥) y el intervalo I donde se satisface la ecuación (I), generalmente 𝐼 ⊂ 𝐷𝑜𝑚𝑦(𝑥) .
Ejemplo: Verifique si la función 𝑦 = 𝐶1 𝑒 2𝑥 + 𝐶2 𝑥𝑒 2𝑥 es solución de la ecuación diferencial:
𝑦 ′′ − 4𝑦 ′ + 4𝑦 = 0
Prueba:
𝑦 = 𝐶1 𝑒 2𝑥 + 𝐶2 𝑥𝑒 2𝑥 → 𝑦 ′ = 2𝐶1 𝑒 2𝑥 + 𝐶2 (1 ⋅ 𝑒 2𝑥 + 𝑥2𝑒 2𝑥 ) = (2𝐶1 + 𝐶2 )𝑒 2𝑥 + 2𝐶2 𝑥𝑒 2𝑥

𝑦 ′′ = 2(2𝐶1 + 𝐶2 )𝑒 2𝑥 + 2𝐶2 (1 ⋅ 𝑒 2𝑥 + 𝑥2𝑒 2𝑥 ) = (4𝐶1 + 4𝐶2 )𝑒 2𝑥 + 4𝐶2 𝑥𝑒 2𝑥
Reemplazando en la ecuación diferencial
𝑦 ′′ − 4𝑦 ′ + 4𝑦 = 0
(4𝐶1 + 4𝐶2 )𝑒 2𝑥 + 4𝐶2 𝑥𝑒 2𝑥 − 4[(2𝐶1 + 𝐶2 )𝑒 2𝑥 + 2𝐶2 𝑥𝑒 2𝑥 ] + 4(𝐶1 𝑒 2𝑥 + 𝐶2 𝑥𝑒 2𝑥 ) = 0

(4𝐶1 + 4𝐶2 − 8𝐶1 − 4𝐶2 + 4𝐶1 )𝑒 2𝑥 + (4𝐶2 − 8𝐶2 + 4𝐶2 )𝑥𝑒 2𝑥 = 0
(0)𝑒 2𝑥 + (0)𝑥𝑒 2𝑥 = 0 → 0 = 0 𝑣𝑒𝑟𝑖𝑓𝑖𝑐𝑎
1.3.1 Regla de Leibniz
Derivación bajo el signo integral

𝑣(𝑥)
𝐹(𝑥 ) = ∫ 𝑔(𝑡) 𝑑𝑡
𝑢(𝑥)
𝐹′ (𝑥 ) = 𝑔[𝑣(𝑥)] ⋅ 𝑣 ′ (𝑥 ) − 𝑔[𝑢(𝑥 )] ⋅ 𝑢′ (𝑥 ) … (1.3)
Problema 2 Compruebe si la función:

𝑥
cos 𝑡
𝑦 = √𝑥 ∫ 𝑑𝑡
4 √𝑡
Es solución de la ecuación diferencial: 2𝑥𝑦 ′ − 𝑦 = 2𝑥 cos 𝑥
Solución: Calculando la derivada de la función y reemplazando la misma en la ecuación

diferencial.
𝑥 𝑥
′
1 cos 𝑡 cos 𝑥 cos 4 1 cos 𝑡
𝑦 = ∫ 𝑑𝑡 + √𝑥 ⋅ [ ⋅1− ⋅ 0] → 𝑦′ = ∫ 𝑑𝑡 + cos 𝑥
2√𝑥 4 √𝑡 √𝑥 √4 2√𝑥 4 √𝑡
En la ecuación diferencial:
𝑥 𝑥
′
1 cos 𝑡 cos 𝑡
2𝑥𝑦 − 𝑦 = 2𝑥 cos 𝑥 → 2𝑥 [ ∫ 𝑑𝑡 + cos 𝑥] − √𝑥 ∫ 𝑑𝑡 = 2𝑥 cos 𝑥
2√𝑥 4 √𝑡 4 √𝑡
𝑥 𝑥
cos 𝑡 cos 𝑡
√𝑥 ∫ 𝑑𝑡 + 2𝑥 cos 𝑥 − √𝑥 ∫ 𝑑𝑡 = 2𝑥 cos 𝑥
4 √𝑡 4 √𝑡
𝟐𝒙 𝐜𝐨𝐬 𝒙 = 𝟐𝒙 𝐜𝐨𝐬 𝒙 𝑣𝑒𝑟𝑖𝑓𝑖𝑐𝑎
Muchas veces se conoce a la función solución de una ecuación diferencial, la pregunta es que:
¿A partir de la solución se pueda llegar a la ecuación diferencial?
Ecuación Diferencial Solución
Problema 3 Halle la ecuación diferencial cuya solución es la función:
𝑦 = 𝐴√1 + 𝑥 2 + 𝐵𝑥
Solución: Para hallar la ecuación diferencial cuya solución general es la función presentada,
debemos eliminar sus constantes generadas por la integración, el número de constantes nos indica
que el orden de la ecuación diferencial, como aquí se ven dos contantes la ecuación diferencial
que originó esta solución debe ser una ecuación diferencial de segundo orden. Una forma de
eliminar estas constantes es aislar la misma para posteriormente eliminarla con derivación.
1
2𝑥 ⋅ 𝑥 − √1 + 𝑥 2 ⋅ 1
𝑦 𝐴√1 + 𝑥 2 𝑑 𝑥𝑦 ′ − 𝑦 ⋅ 1 2√1 + 𝑥 2
= +𝐵 / → = 𝐴[ ]+0
𝑥 𝑥 𝑑𝑥 𝑥2 𝑥2
𝑥2 − 1 − 𝑥2 𝐴 𝑑
𝑥𝑦 ′ − 𝑦 = 𝐴 ⋅ → 𝑥𝑦 ′ − 𝑦 = − → (𝑥𝑦 ′ − 𝑦)√1 + 𝑥 2 = −𝐴 /
√1 + 𝑥 2 √1 + 𝑥 2 𝑑𝑥
1
(1 ⋅ 𝑦 ′ + 𝑥𝑦 ′′ − 𝑦 ′ )√1 + 𝑥 2 + (𝑥𝑦 ′ − 𝑦) ⋅ 2𝑥 = 0
2√1 + 𝑥 2
(𝑥𝑦 ′ − 𝑦)𝑥
𝑥𝑦 ′′ √1 + 𝑥 2 + =0 / ⋅ √1 + 𝑥 2 → 𝑥𝑦 ′′ (1 + 𝑥 2 ) + 𝑥 2 𝑦 ′ − 𝑥𝑦 = 0 / ÷ 𝑥
√1 + 𝑥2
(𝟏 + 𝒙𝟐 )𝒚′′ + 𝒙𝒚′ − 𝒚 = 𝟎
Problema 4 Halle la ecuación diferencial que tiene por solución a la función:
𝑦 = 𝐴𝑥 4 + 𝐵𝑥 −4 + 3𝑥 5
Solución: A diferencia de la forma del anterior ejercicio esta vez trataremos de eliminar las
constantes empleando derivación y desarrollo algebraicos de sumas y restas.
𝑑
𝑦 = 𝐴𝑥 4 + 𝐵𝑥 −4 + 3𝑥 5 /⋅ 𝑥 4 → 𝑥 4 𝑦 = 𝐴𝑥 8 + 𝐵 + 3𝑥 9 / 𝑑𝑥
4𝑥 3 𝑦 + 𝑥 4 𝑦 ′ = 8𝐴𝑥 7 + 27𝑥 8 / ⋅ 𝑥 −3
4𝑦 + 𝑥𝑦 ′ = 8𝐴𝑥 4 + 27𝑥 5 … (1)
Derivando nuevamente (1) respecto de "𝑥"
12𝑥 2 𝑦 + 4𝑥 3 𝑦 ′ + 4𝑥 3 𝑦 ′ + 𝑥 4 𝑦 ′′ = 56𝐴𝑥 6 + 216𝑥 7 →

12𝑦 + 8𝑥𝑦 ′ + 𝑥 2 𝑦 ′′ = 56𝐴𝑥 4 + 216𝑥 5 … (2)
Trabajando con (1) y (2)
4y xy' 8Ax 4 27x 5 / 7 28 y 7x y ' 56Ax 4 189x 5

12y 8x y ' x 2y '' 56Ax 4 216x 5 / 1 12y 8x y ' x 2y '' 56Ax 4 216x 5
16𝑦 − 𝑥𝑦 ′ − 𝑥 2 𝑦 ′′ = −27𝑥 5
𝑹𝒑𝒕𝒂.: 𝒙𝟐 𝒚′′ + 𝒙𝒚′ − 𝟏𝟔𝒚 = 𝟐𝟕𝒙𝟓
1.4 Existencia y Unicidad
Sea la ecuación diferencial 𝑦 ′ = 𝑓(𝑥, 𝑦) sujeta a la condición de valor inicial 𝑦(𝑥0 ) = 𝑦0

y
R
y0
x0 x
𝜕𝑓
Si 𝑓(𝑥, 𝑦) y 𝜕𝑦
son continuas en la región R, entonces existe una única solución, que pasa por
(𝑥0 , 𝑦0 ).
Si 𝑓(𝑥, 𝑦) es continua en R, entonces existe al menos una solución en su interior

𝜕𝑓
Si es continua, entonces existe una sola solución que pasa por el punto (𝑥0 , 𝑦0 ).
𝜕𝑦
1.5 Tipos de soluciones
Solución General: Una ecuación diferencial de primer orden se resuelve integrando, lo cual genera
una constante obteniendo la denominada curva solución (curva integral) que es una solución
monoparamétrica de soluciones, es decir que contiene una constante.
Es decir si se tiene la ecuación diferencial: 𝐹 (𝑥, 𝑦, 𝑦 ′ ) = 0 … (1.4)
Su solución será de la forma: 𝐺 (𝑥, 𝑦, 𝐶 ) = 0 … (1.5)
C1
C2
C3
x
C4
Solución Particular: No posee parámetro (constante) se halla para algún valor inicial de la
solución general.
Solución Singular: Es una solución que verifica la ecuación diferencial pero no proviene de ningún
valor de la constante de la solución general, se obtienen por otros procedimientos, se presenta en
raras ocasiones, una de estas generalmente es la solución trivial 𝒚 = 𝟎
1.6 Problemas de Valor Inicial y Problemas de valor de frontera

Problemas de valor inicial vs problema de valor de frontera
𝐹(𝑥, 𝑦, 𝑦 ′ , 𝑦 ′′ , 𝑦 ′′′ , … . . , 𝑦 (𝑛) ) = 0 𝐹(𝑥, 𝑦, 𝑦 ′, 𝑦 ′′ , 𝑦 ′′′ , … . . , 𝑦 (𝑛) ) = 0 Con

las condiciones son en el mismo punto las condiciones son en distintos valores
𝑦(𝑥0 ) = 𝑦0 𝑦(𝑥0 ) = 𝑦0
𝑦 ′ (𝑥0 ) = 𝑦0′ 𝑦(𝑥1 ) = 𝑦1
𝑦 ′′ (𝑥0 ) = 𝑦0′′ 𝑦(𝑥2 ) = 𝑦2
… ….
… ….
(𝑛−1)
𝑦 (𝑛−1) (𝑥0 ) = 𝑦0 𝑦( 𝑥𝑛−1 ) = 𝑦𝑛−1
Ejemplo: Movimiento rectilíneo acelerado Ejemplo: Viga
Fx mg sin
x (t ) x2
x1 y2
W mg
y1
𝑥(0) = 0 (cae del reposo) 𝑦(𝑥1 ) = 𝑦1
𝑥 ′ (0) = 0 (parte del reposo) 𝑦(𝑥2 ) = 𝑦2
𝑥 ′′ (0) = 𝑔 sin 𝛼 (parte del reposo sin rozamiento) 𝑦(𝑥3 ) = 𝑦3
En el plano inclinado se aplican las leyes de Newton en la dirección del movimiento:
∑ 𝐹𝑚𝑜𝑣 = 𝑚𝑎 → 𝐹𝑥 = 𝑚𝑎 → 𝑚𝑔𝑠𝑖𝑛 𝛼 = 𝑚𝑎 → 𝒂 = 𝑔 sin 𝛼
La diferencia esencial ente los problemas de valor inicial y los de frontera es que en los problemas
de condiciones iniciales todos las condiciones son en un punto en particular 𝑥0 en el ejemplo en
𝑥0 = 0.
1.7 Métodos de Resolución de Ecuaciones diferenciales Ordinarias de Primer Orden
𝐹(𝑥, 𝑦, 𝑦 ′ ) = 0
Existen dos formas:
i) La forma diferencial ii) Forma Normal

𝑑𝑦
𝑃 (𝑥, 𝑦)𝑑𝑥 + 𝑄 (𝑥, 𝑦)𝑑𝑦 = 0 𝑑𝑥
= 𝑓(𝑥, 𝑦)
Empezaremos desarrollando técnicas de solución en la forma diferencial.
1.7.1 Ecuación diferencial de variable separable
𝑃 (𝑥, 𝑦)𝑑𝑥 + 𝑄 (𝑥, 𝑦)𝑑𝑦 = 0 …. (1.6)
Por este método es posible factorizar 𝑃(𝑥, 𝑦) = 𝑝1 (𝑥 )𝑝2 (𝑦) ; 𝑄 (𝑥, 𝑦) = 𝑞1 (𝑥 )𝑞2 (𝑦), es decir se
puede escribir ambas expresiones 𝑃 y 𝑄 como productos de factores que dependan exclusivamente
de una sola variable.
En (𝛼 ) 𝑝1 (𝑥 )𝑝2 (𝑦)𝑑𝑥 + 𝑞1 (𝑥 )𝑞2 (𝑦)𝑑𝑦 = 0 / ÷ 𝑝2 (𝑦)𝑞1 (𝑥 )
𝑝1 (𝑥 ) 𝑞2 (𝑦)
𝑑𝑥 + 𝑑𝑦 = 0 / ∫
𝑞1 (𝑥 ) 𝑝2 (𝑦)
𝑝1 (𝑥 ) 𝑞2 (𝑦)
∫ 𝑑𝑥 + ∫ 𝑑𝑦 = ∫ 0
𝑞1 (𝑥 ) 𝑝2 (𝑦)
𝑭(𝒙) + 𝑮(𝒚) = 𝑪 … (1.7)
Problema 5 Resolver:
1 + 𝑦2
𝑥𝑦 𝑑𝑦 − 𝑑𝑥 = 0
1 + 𝑥2
Solución:
1 + 𝑦2
𝑑𝑥 − 𝑥𝑦
⏟ 𝑑𝑦 = 0 / ÷ (1 + 𝑦 2 )𝑥
1⏟+ 𝑥 2
𝑄
𝑃
1 𝑦
∫ 2
𝑑𝑥 − ∫ 𝑑𝑦 = ∫ 0
𝑥(1 + 𝑥 ) 1 + 𝑦2
2𝑥 1 2𝑦
∫ 𝑑𝑥 − 2 ∫ 𝑑𝑦 = ∫ 0
2𝑥 2 (1 2
+𝑥 ) 1 + 𝑦2
Para la primera integral: 𝑢 = 𝑥 2 → 𝑑𝑢 = 2𝑥𝑑𝑥
1 1 1 1 1 −1 1
∫ 𝑑𝑢 − ln(1 + 𝑦 2 ) = 𝐶 → ∫( + ) 𝑑𝑢 − ln(1 + 𝑦 2 ) = 𝐶 / ⋅ 2
2 𝑢(1 + 𝑢) 2 2 𝑢 1+𝑢 2
𝑢
ln|𝑢| − ln|1 + 𝑢| − ln|1 + 𝑦 2 | = 2𝐶 → ln | |=𝐾
(1 + 𝑢)(1 + 𝑦 2 )
Retornando: A
𝑥2 𝑥2
ln ( )=𝐾 → = 𝑒𝐾
(1 + 𝑥 2 )(1 + 𝑦 2 ) (1 + 𝑥 2 )(1 + 𝑦 2 )
𝒙𝟐 = 𝑨(𝟏 + 𝒙𝟐 )(𝟏 + 𝒚𝟐 )
Graficando la solución
Llevando la ecuación diferencial a la forma normal.
1 + 𝑦2 1 + 𝑦2 𝑑𝑦 1 + 𝑦2 ′
1 + 𝑦2
𝑥𝑦 𝑑𝑦 − 𝑑𝑥 = 0 → 𝑥𝑦 𝑑𝑦 = 𝑑𝑥 → = → 𝑦 =
1 + 𝑥2 1 + 𝑥2 𝑑𝑥 𝑥𝑦(1 + 𝑥 2 ) 𝑥𝑦(1 + 𝑥 2 )
campo direccional de la ecc. dif .
curva solución
x 2  A 1  x 2  1  y 2 
1.7.2 Ecuación diferencial homogénea
Extendiendo el concepto de homogeneidad de un polinomio al concepto de función
𝑃(𝑥, 𝑦) = 6𝑥 2 𝑦 5 + 𝑥 5 𝑦 2 + 𝑦 7
Extendiendo este concepto para funciones, la función 𝑓(𝑥, 𝑦) será homogénea de grado "𝑛" si
cumple:
𝑓(𝑡𝑥, 𝑡𝑦) = 𝑡 𝑛 𝑓(𝑥, 𝑦) → 𝑓(𝑥, 𝑦) es homogénea de grado "𝑛"
Para el polinomio del ejemplo:
𝑃(𝑡𝑥, 𝑡𝑦) = 6(𝑡𝑥)2 (𝑡𝑦)5 + (𝑡𝑥)5 (𝑡𝑦)2 + (𝑡𝑦)7

= 6𝑡 7 𝑥 2 𝑦 5 + 𝑡 7 𝑥 5 𝑦 2 + 𝑡 7 𝑦 7 = 𝑡 7 (6𝑥 2 𝑦 5 + 𝑥 5 𝑦 2 + 𝑦 7 )
𝑃(𝑡𝑥, 𝑡𝑦) = 𝑡 7 𝑃(𝑥, 𝑦)
Más ejemplos:
𝑦
𝑓(𝑥, 𝑦) = 3𝑥 2 𝑦 − 𝑥 3 cos ( ) − 𝑦 3
𝑥
𝑡𝑦 𝑦
𝑓(𝑡𝑥, 𝑡𝑦) = 3(𝑡𝑥 )2 (𝑡𝑦) − (𝑡𝑥)3 cos ( ) − (𝑡𝑦)3 = 𝑡 3 (3𝑥 2 𝑦 − 𝑥 3 cos ( ) − 𝑦 3 ) = 𝑡 3 𝑓(𝑥, 𝑦)
𝑡𝑥 𝑥
Es homogénea de grado 3
Procedimiento de Solución de ecuaciones diferenciales homogéneas
𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄 (𝑥, 𝑦)𝑑𝑦 = 0 …. (𝛼)

Para la resolución se tiene cumplir que 𝑃(𝑥, 𝑦) y 𝑄(𝑥, 𝑦) sean homogéneas del mismo grado
1 𝑦 1 𝑛 𝑦
𝑃(𝑡𝑥, 𝑡𝑦) = 𝑡 𝑛 𝑃(𝑥, 𝑦) C.V. 𝑡 = 𝑥 → 𝑃 (1, 𝑥 ) = (𝑥) 𝑃(𝑥, 𝑦) → 𝑃(𝑥, 𝑦) = 𝑥 𝑛 𝑃(1, 𝑥 )
1 𝑦 1 𝑛 𝑦
𝑄 (𝑡𝑥, 𝑡𝑦) = 𝑡 𝑛 𝑄(𝑥, 𝑦) C.V. 𝑡 = 𝑥 → 𝑄 (1, 𝑥 ) = (𝑥) 𝑄(𝑥, 𝑦) → 𝑄 (𝑥, 𝑦) = 𝑥 𝑛 𝑄(1, 𝑥 )
En (𝛼)
𝑦 𝑦
𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄(𝑥, 𝑦)𝑑𝑦 = 0 → 𝑥 𝑛 𝑃 (1, ) 𝑑𝑥 + 𝑥 𝑛 𝑄 (1, ) 𝑑𝑦 = 0
𝑥 𝑥
𝑦 𝑦
𝑥 𝑛 [𝑃 (1, ) 𝑑𝑥 + 𝑄 (1, ) 𝑑𝑦] = 0 𝑟𝑒𝑐𝑜𝑟𝑑𝑎𝑛𝑑𝑜 𝑎⋅𝑏 =0 ↔𝑎 = 0∨𝑏 = 0
𝑥 𝑥
𝑦 𝑦
𝑃 (1, ) 𝑑𝑥 + 𝑄 (1, ) 𝑑𝑦 = 0
𝑥 𝑥
𝑦
𝐶. 𝑉. 𝑢 = 𝑥 → 𝑦 = 𝑢𝑥 → 𝑑𝑦 = 𝑥𝑑𝑢 + 𝑢𝑑𝑥
𝑃(1, 𝑢)𝑑𝑥 + 𝑄 (1, 𝑢)(𝑥𝑑𝑢 + 𝑢𝑑𝑥) = 0

[⏟𝑃 (1, 𝑢) + 𝑢𝑄(1, 𝑢)] 𝑑𝑥 + 𝑥𝑄
⏟ (1, 𝑢) 𝑑𝑢 = 0
𝑀 𝑁
𝑑𝑥 𝑄(1, 𝑢)
∫ +∫ 𝑑𝑢 = ∫ 0
𝑥 𝑃(1, 𝑢) + 𝑢𝑄(1, 𝑢)
𝑸(𝟏, 𝒖)
𝐥𝐧|𝒙| + ∫ 𝒅𝒖 = 𝑪 … (1.8)
𝑷(𝟏, 𝒖) + 𝒖𝑸(𝟏, 𝒖)
Ya es de variable separable:
En resumen, si la ecuación diferencial es homogénea, se pueden hacer uno de esto dos

cambios: 𝑦 = 𝑢𝑥 → 𝑑𝑦 = 𝑥𝑑𝑢 + 𝑢𝑑𝑥 ó 𝑥 = 𝑢𝑦 → 𝑑𝑥 = 𝑦𝑑𝑢 + 𝑢𝑑𝑦
La elección por lo general brindará un ahorro en operaciones del tipo algebraico e incluso en la
integración.
Problema 6 Resolver
𝑥𝑦 ′ = 𝑦 + √𝑦 2 − 𝑥 2
Solución: Llevándola a la forma diferencial
𝑑𝑦
𝑥 = 𝑦 + √𝑦 2 − 𝑥 2 → 2 2
⏟ + √𝑦 − 𝑥 ) 𝑑𝑥 − ⏟
(𝑦 𝑥 𝑑𝑦 = 0
𝑑𝑥 𝑄
𝑃
(𝑦 ⏟ 2 − 𝑥 2 ) 𝑑𝑥 − ⏟
⏟ + √𝑦 𝑥 𝑑𝑦 = 0 es una ecuación diferencial homogénea:
⏟1° 1° ⏟
1°
𝑃 𝑄
CV. 𝑦 = 𝑢𝑥 → 𝑑𝑦 = 𝑥𝑑𝑢 + 𝑢𝑑𝑥

(𝑢𝑥 + √(𝑢𝑥)2 − 𝑥 2 ) 𝑑𝑥 − 𝑥(𝑥𝑑𝑢 + 𝑢𝑑𝑥 ) = 0 → (𝑢𝑥 + √𝑥 2 (𝑢2 − 1)) 𝑑𝑥 − 𝑥(𝑥𝑑𝑢 + 𝑢𝑑𝑥) = 0
𝑥 [(𝑢 + √𝑢2 − 1) 𝑑𝑥 − (𝑥𝑑𝑢 + 𝑢𝑑𝑥)] = 0 → (𝑢 + √𝑢2 − 1) 𝑑𝑥 − (𝑥𝑑𝑢 + 𝑢𝑑𝑥 ) = 0
(𝑢 + √𝑢2 − 1 − 𝑢) 𝑑𝑥 − 𝑥𝑑𝑢 = 0
√𝑢2 − 1 𝑑𝑥 − 𝑥𝑑𝑢 = 0
𝑑𝑥 1
∫ −∫ 𝑑𝑢 = ∫ 0
𝑥 √𝑢2 − 1
ln|𝑥| − ln |𝑢 + √𝑢2 − 1| = 𝐾
ln 𝐶
𝑥
ln | | = ln 𝐶
𝑢 + √𝑢2 − 1
𝑦 𝑦
𝑥 = 𝐶(𝑢 + √𝑢2 − 1) retornando a la variable original: 𝑥 = 𝐶 ( + √( )2 − 1)
𝑥 𝑥
𝒙𝟐 = 𝑪 (𝒚 + √𝒚𝟐 − 𝒙𝟐 )
campo direccional de la ecc. dif .
curva solución

x 2  C y  y2  x 2 
Problema 7 Resolver: (II-2018)
(𝑦 4 − 𝑥 + √𝑥𝑦 2 )𝑑𝑥 − 4𝑦 5 √𝑥𝑑𝑦 = 0
Solución: no es homogénea, pero intentaremos un cambio de variable
4 2
(𝑦 ⏟ 5 √𝑥 𝑑𝑦 = 0
⏟ − 𝑥 + √𝑥𝑦 ) 𝑑𝑥 − 4𝑦
𝑃(𝑥,𝑦) 𝑄(𝑥,𝑦)
2
((𝑦 2 )2 − (√𝑥) + √𝑥𝑦 2 ) 𝑑𝑥 − 4√𝑥 (𝑦 2 )2 𝑦𝑑𝑦 = 0 / ÷ 2√𝑥
2 𝑑𝑥
((𝑦 2 )2 − (√𝑥) + √𝑥𝑦 2 ) − (𝑦 2 )2 2𝑦𝑑𝑦 = 0
2√𝑥
1
C.V. 𝑡 = 𝑦 2 → 𝑑𝑡 = 2𝑦𝑑𝑦 𝑧 = √𝑥 → 𝑑𝑧 = 2 𝑑𝑥
√𝑥
𝑡 2 𝑑𝑡 = 0
(𝑡 2 − 𝑧 2 + 𝑧𝑡) 𝑑𝑧 − ⏟
⏟
𝑀 𝑁
Es homogénea de grado 2: C.V. 𝑡 = 𝑢𝑧 → 𝑑𝑡 = 𝑧𝑑𝑢 + 𝑢𝑑𝑧
(𝑢2 𝑧 2 − 𝑧 2 + 𝑧𝑢𝑧)𝑑𝑧 − 𝑢2 𝑧 2 (𝑧𝑑𝑢 + 𝑢𝑑𝑧) = 0

𝑧 2 [(𝑢2 − 1 + 𝑢)𝑑𝑧 − 𝑢2 (𝑧𝑑𝑢 + 𝑢𝑑𝑧)] = 0 → (𝑢2 − 1 + 𝑢 − 𝑢3 )𝑑𝑧 − 𝑢2 𝑧𝑑𝑢 = 0
𝑑𝑧 𝑢2
∫ −∫ 2 𝑑𝑢 = ∫ 0
𝑧 𝑢 + 𝑢 − 1 − 𝑢3
𝑢2 + 𝑢 − 1 − 𝑢3 = 𝑢(𝑢 + 1) − (1 + 𝑢3 ) = 𝑢(𝑢 + 1) − (1 + 𝑢)(1 − 𝑢 + 𝑢2 )
= (1 + 𝑢)(𝑢 − 1 + 𝑢 − 𝑢2 )
𝑢2
ln|𝑧| − ∫ =𝐾
(1 + 𝑢)(2𝑢 − 1 − 𝑢2 )
𝑢2 𝑢2
ln|𝑧| + ∫ =𝐾 → ln|𝑧| + ∫ =𝐾
(1 + 𝑢)(𝑢2 − 2𝑢 + 1) (1 + 𝑢)(𝑢 − 1)2
ln|𝑧| + ∫ ( + + )
𝑢 + 1 𝑢 − 1 (𝑢 − 1)2
1 1
𝑢2 4 2 𝑢2 𝐴 𝐵 𝐶
= + + → = + +
(𝑢 + 1)(𝑢 − 1)2 𝑢 + 1 𝑢 − 1 (𝑢 − 1)2 (𝑢 + 1)(𝑢 − 1)2 𝑢 + 1 𝑢 − 1 (𝑢 − 1)2
𝑑
𝑢2 = 𝐴(𝑢 − 1)2 + 𝐵 (𝑢 + 1)(𝑢 − 1) + 𝐶 (𝑢 + 1) / 𝑑𝑢
2𝑢 = 2𝐴(𝑢 − 1) + 𝐵[(𝑢 − 1) + (𝑢 + 1)] + 𝐶 / |𝑢=1
1 1 3
2 ⋅ 1 = 2 ⋅ 4 ⋅ 0 + 𝐵 (0 + 2) + → 𝐵=
2 4
1 3 1
4 4 2
ln|𝑧| + ∫ ( + + ) 𝑑𝑢 = 𝐾
𝑢 + 1 𝑢 − 1 (𝑢 − 1)2
1 3 1 1
ln|𝑧| + ln|𝑢 + 1| + ln|𝑢 − 1| − ⋅ =𝐾 /⋅4
4 4 2 𝑢−1
2 2
ln|𝑧 4 (𝑢 + 1)(𝑢 − 1)3 | = 4𝐾 + → 𝑧 4 (𝑢 + 1)(𝑢 − 1)3 = 𝑒 4𝐾+𝑢−1
𝑢−1
2
𝑧 4 (𝑢 + 1)(𝑢 − 1)3 = 𝑒 4𝐾 𝑒 𝑢−1 retornando a la variable original
2
3 𝑦2
𝑡 𝑦2 4 𝑦2 𝑦2 −1
𝑢= = → (√𝑥) ( + 1) ( − 1) = 𝐶𝑒 √𝑥
𝑧 √𝑥 √𝑥 √𝑥
2√𝑥 3
𝑦 2 + √𝑥 𝑦 2 − √𝑥
2
𝑥 ( )( ) = 𝐶𝑒 𝑦2−√𝑥
√𝑥 √𝑥
𝟐√𝒙
𝟑
(𝒚𝟐 + √𝒙)(𝒚𝟐 − √𝒙) = 𝑪𝒆𝒚𝟐 −√𝒙
Problema 8 Resolver la ecuación diferencial (I-2021)
𝑦+𝑥
𝑥𝑦 ′ − 𝑦 = (𝑦 + 𝑥 )𝑙𝑛 ( )
𝑥
Solución: Es homogénea de grado 1
𝑦 = 𝑡𝑥 → 𝑑𝑦 = 𝑥𝑑𝑡 + 𝑡𝑑𝑥 o 𝑦 ′ = 𝑥𝑡 ′ + 𝑡
𝑡𝑥 + 𝑥 𝑥(𝑡 + 1)
𝑥(𝑥𝑡 ′ + 𝑡) − 𝑡𝑥 = (𝑡𝑥 + 𝑥 )𝑙𝑛 ( ) ∶ 𝑥 2 𝑡 ′ + 𝑥𝑡 − 𝑡𝑥 = 𝑥 (𝑡 + 1)𝑙𝑛 ( )
𝑥 𝑥
𝑥(𝑡 + 1)
𝑥𝑡 ′ = (𝑡 + 1)𝑙𝑛 ( ) ∶ 𝑥𝑡 ′ = (𝑡 + 1)𝑙𝑛(𝑡 + 1)
𝑥
𝑑𝑡 𝑑𝑥 𝑑𝑡 𝑑𝑥
= → ∫ =∫
(𝑡 + 1)𝑙𝑛(𝑡 + 1) 𝑥 (𝑡 + 1)𝑙𝑛(𝑡 + 1) 𝑥
ln[ln(𝑡 + 1)] = 𝑙𝑛𝑥 + 𝑙𝑛𝐶 → ln(𝑡 + 1) = 𝐶𝑥 → 𝑡 + 1 = 𝑒 𝐶𝑥
𝑦
Retornando a la variable original: = 𝑒 𝐶𝑥 − 1 → 𝑦 = 𝑥(𝑒 𝐶𝑥 − 1)
𝑥
𝒚 = 𝒙(𝒆𝑪𝒙 − 𝟏)
1.7.3 Ecuación diferencial reducible a homogénea
1.7.3.1 Ecuación diferencial Isobárica
𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄 (𝑥, 𝑦)𝑑𝑦 = 0 … (𝛼)
Recordando que la ecuación (𝛼) es homogénea si 𝑃(𝑡𝑥, 𝑡𝑦) = 𝑡 𝑛 𝑃(𝑥, 𝑦) y 𝑄 (𝑡𝑥, 𝑡𝑦) = 𝑡 𝑛 𝑄(𝑥, 𝑦),
sin embargo, las ecuaciones isóbaras generalizan a las ecuaciones homogéneas por cuanto los
coeficientes 𝑃(𝑥, 𝑦) y 𝑄(𝑥, 𝑦) no son funciones homogéneas del mismo grado. Se busca una
transformación que convierta la ecuación (𝛼) en homogénea si es que no lo es.
Una ecuación Diferencial es Isóbara si cumple:
𝑃(𝑡𝑥, 𝑡 𝑚 𝑦) = 𝑡 𝑛 𝑃(𝑥, 𝑦)
… (1.9)
𝑄(𝑡𝑥, 𝑡 𝑚 𝑦) = 𝑡 𝑛−𝑚+1 𝑄(𝑥, 𝑦)
1
Si empelamos el cambio: 𝑡 = 𝑥
𝑦 1 𝒚
𝑃 (1, 𝑚
) = 𝑛 𝑃(𝑥, 𝑦) → 𝑷(𝒙, 𝒚) = 𝒙𝒏 𝑷 (𝟏, )
𝑥 𝑥 𝒙𝒎
𝑦 1 𝒚
𝑄 (1, ) = 𝑄(𝑥, 𝑦) → 𝑸(𝒙, 𝒚) = 𝒙𝒏−𝒎+𝟏 𝑸 (𝟏, )
𝑥𝑚 𝑥 𝑛−𝑚+1 𝒙𝒎
Reemplazando estas últimas expresiones en (𝛼): 𝑃 (𝑥, 𝑦)𝑑𝑥 + 𝑄 (𝑥, 𝑦)𝑑𝑦 = 0
𝒚 𝒚
𝒙𝒏 𝑷 (𝟏, ) 𝒅𝒙 + 𝒙𝒏−𝒎+𝟏
𝑸 (𝟏, ) 𝒅𝒚 = 𝟎 … (1.10)
𝒙𝒎 𝒙𝒎
El cambio de variable será:
𝒚
𝒖= … (1.11)
𝒙𝒎
𝑦 = 𝑥𝑚𝑢 → 𝑑𝑦 = 𝑚𝑥 𝑚−1 𝑢𝑑𝑥 + 𝑥 𝑚 𝑑𝑢
𝑥 𝑛 [𝑃(1, 𝑢)𝑑𝑥 + 𝑥 1−𝑚 𝑄 (1, 𝑢)(𝑚𝑥 𝑚−1 𝑢𝑑𝑥 + 𝑥 𝑚 𝑑𝑢)] = 0
[𝑃(1, 𝑢) + 𝑚𝑄(1, 𝑢)𝑢]𝑑𝑥 + 𝑄 (1, 𝑢)𝑥𝑑𝑢 = 0
Ya es de variable separable:
𝒅𝒙 𝑸(𝟏, 𝒖)
∫ +∫ 𝒅𝒖 = 𝑪 … (1.12)
𝒙 𝑷(𝟏, 𝒖) + 𝒎𝑸(𝟏, 𝒖)𝒖
Nota: Las ecuaciones Homogéneas tienen una invariancia de escala frente a la transformación:
(𝑥, 𝑦) → (𝑡𝑥, 𝑡𝑦)
En cambio, una ecuación isobárica tiene una invariancia de escala frente al cambio:
(𝑥, 𝑦) → (𝑡𝑥, 𝑡 𝑚 𝑦)
𝒚
El cambio como bien se vio es: 𝒖 = 𝒙𝒎, que convierte a la ecuación diferencial en una ecuación
de variable separable, pero primero se debe identificar 𝑚, en consecuencia el camino más
recomendable es primero hallar "𝑚" y con ello transformar primero la ecuación diferencial en
homogénea. Entonces el cambio sugerido primero es:
𝟏 1
𝑧 = 𝒙𝒎 → 𝒙 = 𝒛𝒎 o mejor aún 𝑥 = 𝑧 𝛼 donde 𝛼 =
𝑚
En conclusión, se sugiere primero transformar a la ecuación en homogénea con cualquiera de

los cambios siguientes:
𝑥 = 𝑧 𝛼 → 𝑑𝑥 = 𝛼𝑧 𝛼−1 𝑑𝑧 … (1.13)
𝑦 = 𝑧 𝛼 → 𝑑𝑦 = 𝛼𝑧 𝛼−1 𝑑𝑧 … (1.14)
Donde "𝛼" es una constante que debe calcularse para generar la ecuación homogénea
4𝑥𝑦 2 𝑑𝑥 + (3𝑥 2 𝑦 − 1)𝑑𝑦 = 0 ; 𝑦(1) = 3
Solución:
4𝑥𝑦 2 𝑑𝑥 + ⏟
⏟ (3𝑥 2 𝑦 − 1) 𝑑𝑦 = 0
𝑃 𝑄
Tratando por la forma isobárica: C.V. 𝑥 = 𝑧 𝑚 → 𝑑𝑥 = 𝑚𝑧 𝑚−1 𝑑𝑧

4𝑧 𝑚 𝑦 2 𝑚𝑧 𝑚−1 𝑑𝑧 + (3𝑧 2𝑚 𝑦 − 1)𝑑𝑦 = 0 → 4𝑚𝑦 2 𝑧 2𝑚−1 𝑑𝑧 + (3𝑧

⏟ ⏟ 2𝑚 𝑦 − ⏟
1 ) 𝑑𝑦 = 0 … (∗)
2𝑚+1 2𝑚+1 0
Sabemos para que esta ecuación sea homogénea debe cumplirse que:
1
2𝑚 + 1 = [2𝑚 + 1 = 0] → 𝑚=−
2
𝟏 1 3
Por lo tanto el cambio debe ser: 𝒙 = 𝒛−𝟐 → 𝑑𝑥 = − 𝑧 − 2𝑑𝑧 en (∗)
2
1 1 1
4 (− ) 𝑦 2 𝑧 2(−2)−1 𝑑𝑧 + (3𝑧 2(−2)𝑦 − 1) 𝑑𝑦 = 0
2
−2𝑦 2 𝑧 −2 𝑑𝑧 + (3𝑧 −1 𝑦 − 1)𝑑𝑦 = 0
𝑧 −2 (−2𝑦 2 𝑑𝑧 + (3𝑦𝑧 − 𝑧 2 )𝑑𝑦) = 0
(3𝑦𝑧 − 𝑧 2 )𝑑𝑦 − 2𝑦 2 𝑑𝑧 = 0 ya es homogénea
C.V. 𝑧 = 𝑢𝑦 → 𝑑𝑧 = 𝑦𝑑𝑢 + 𝑢𝑑𝑦
(3𝑦𝑢𝑦 − 𝑢2 𝑦 2 )𝑑𝑦 − 2𝑦 2 (𝑦𝑑𝑢 + 𝑢𝑑𝑦) = 0 → 𝑦 2 [(3𝑢 − 𝑢2 )𝑑𝑦 − 2(𝑦𝑑𝑢 + 𝑢𝑑𝑦)] = 0

(3𝑢 − 𝑢2 − 2𝑢)𝑑𝑦 − 2𝑦𝑑𝑢 = 0 → (𝑢 − 𝑢2 )𝑑𝑦 − 2𝑦𝑑𝑢 = 0
𝑑𝑦 2
∫ −∫ 𝑑𝑢 = ∫ 0
𝑦 𝑢(1 − 𝑢)
2 −2 2
ln|𝑦| + ∫ 𝑑𝑢 = 𝐾 → ln|𝑦| + ∫ ( + ) 𝑑𝑢 = 𝐾
𝑢(𝑢 − 1) 𝑢 𝑢−1
𝑦(𝑢 − 1)2
ln|𝑦| − 2 ln|𝑢| + 2 ln|𝑢 − 1| = ln 𝐶 → ln | | = ln 𝐶
𝑢2
𝑧 𝑥 −2 1
Retornando 𝑢 = = =
𝑦 𝑦 𝑥2𝑦
2
1
𝑦( 2 − 1)
𝑥 𝑦
=𝐶 → 𝑦(1 − 𝑥 2 𝑦)2 = 𝐶
1 2
( 2 )
𝑥 𝑦
Para la constante usamos la condición inicial: 𝑦(1) = 3 → 𝑥=1 ⟹𝑦=3
3(1 − 12 ⋅ 3)2 = 𝐶 → 𝐶 = 12
𝒚(𝟏 − 𝒙𝟐 𝒚)𝟐 = 𝟏𝟐
Problema 10 Resolver
(5𝑥 3 𝑦 + 2𝑦 2 )𝑑𝑥 + (4𝑥𝑦 + 3𝑥 4 )𝑑𝑦 = 0 ; 𝑦(1) = 2
Solución: Verificando si la ecuación diferencial es una ecuación isóbara (puede reduciré a una
ecuación homogénea) bajo el cambio 𝑦 = 𝑧 𝛼
⏟ 3 𝑦 + 2𝑦
(5𝑥 ⏟2 ) 𝑑𝑥 + (4𝑥𝑦 ⏟4 ) 𝑑𝑦 = 0
⏟ + 3𝑥
4° 2° 2° 4°
Cambio de variable 1: 𝑦 = 𝑧 𝛼 → 𝑑𝑦 = 𝛼𝑧 𝛼−1 𝑑𝑧
(5𝑥 3 𝑧 𝛼 + 2𝑧 2𝛼 )𝑑𝑥 + (4𝑥𝑧 𝛼 + 3𝑥 4 )𝛼𝑧 𝛼−1 𝑑𝑧 = 0
⏟ 3 𝑧 𝛼 + 2𝑧
(5𝑥 ⏟ 2𝛼
⏟ 2𝛼−1 + ⏟
) 𝑑𝑥 + (4𝛼𝑥𝑧 3𝛼𝑥 4 𝑧 𝛼−1 ) 𝑑𝑧 = 0 … (∗)
3+𝛼 2𝛼 2𝛼 3+𝛼
[3 + 𝛼 = 2𝛼 ] = [2𝛼 = 3 + 𝛼 ]
𝛼=3
Por lo tanto el cambio de variable 1 es: 𝒚 = 𝒛𝟑
Sustituyendo en (∗)
(5𝑥 3 𝑧 3 + 2𝑧 6 )𝑑𝑥 + (12𝑥𝑧 5 + 9𝑥 4 𝑧 2 )𝑑𝑧 = 0
Es homogénea, cambio de variable 2: 𝑧 = 𝑢𝑥 → 𝑑𝑧 = 𝑥𝑑𝑢 + 𝑢𝑑𝑥
(5𝑥 3 𝑢3 𝑥 3 + 2𝑥 6 𝑢6 )𝑑𝑥 + (12𝑥𝑢5 𝑥 5 + 9𝑥 4 𝑢2 𝑥 2 )(𝑥𝑑𝑢 + 𝑢𝑑𝑥 ) = 0

𝑥 6 [(5𝑢3 + 2𝑢6 )𝑑𝑥 + (12𝑢5 + 9𝑢2 )(𝑥𝑑𝑢 + 𝑢𝑑𝑥 )] = 0
(14𝑢6 + 14𝑢3 )𝑑𝑥 + (12𝑢5 + 9𝑢2 )𝑥𝑑𝑢 = 0
𝑑𝑥 12𝑢3 + 9 𝑑𝑥 4𝑢3 + 3
∫ +∫ 𝑑𝑢 = ∫ 0 → ∫ +∫ 3𝑢2 𝑑𝑢 = ∫ 0
𝑥 14𝑢(𝑢3 + 1) 𝑥 14𝑢3 (𝑢3 + 1)
Integrando con el cambio de variable: 𝑡 = 𝑢3 → 𝑑𝑡 = 3𝑢2 𝑑𝑢
𝑑𝑥 4𝑡 + 3 𝑑𝑥 1 3 1
∫ +∫ 𝑑𝑡 = ∫ 0 → ∫ + 14 ∫ [ + ] 𝑑𝑡 = ∫ 0
𝑥 14𝑡(𝑡 + 1) 𝑥 𝑡 𝑡+1
1
ln 𝑥 + 14[3 ln 𝑡 + ln(𝑡 + 1)] = ln 𝑘 → ln[𝑥 14 𝑡 3 (𝑡 + 1)] = ln 𝑘 14 → 𝑥 14 𝑡 3 (𝑡 + 1) = 𝐶
1
𝑧 𝑦3
Retornando de todos los cambios de variable: 𝑡 = 𝑢3 , 𝑢=𝑥= 𝑥
1 9 1 3
𝑦3 𝑦3
𝑥 14 𝑡 3 (𝑡 + 1) = 𝐶 → 𝑥 14 (𝑢3 )3 (𝑢3 + 1) = 𝐶 → 𝑥 14 ( ) [( ) + 1] = 𝐶 → 𝑥 14
𝑥 𝑥
𝑥 2 𝑦 3 (𝑦 + 𝑥 3 ) = 𝐶 evaluando en la condición inicial : 12 ⋅ 23 (2 + 13 ) = 𝐶 → 𝐶 = 24
𝒙𝟐 𝒚𝟑 (𝒚 + 𝒙𝟑 ) = 𝟐𝟒
1.7.3.2 Ecuación diferencial de Jacobi Tiene la forma:
𝒂𝒙 + 𝒃𝒚 + 𝒄
𝒚′ = 𝒇 ( ) … (1.15)
𝒅𝒙 + 𝒆𝒚 + 𝒇
Esta puede llevarse a la forma diferencial y tanto 𝑃(𝑥, 𝑦) y 𝑄(𝑥, 𝑦) son rectas
𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄 (𝑥, 𝑦)𝑑𝑦 = 0 … (𝛼)

(𝐴𝑥 + 𝐵𝑦 + 𝐶 ) 𝑑𝑥 + ⏟
⏟ (𝐷𝑥 + 𝐸𝑦 + 𝐹) 𝑑𝑦 = 0
𝑷 𝑸
Hay dos casos:
i) Rectas paralelas: Se transforma en variable separable con un cambio de variable.
(𝐴𝑥 + 𝐵𝑦 + 𝐶 )𝑑𝑥 + (𝐷𝑥 + 𝐸𝑦 + 𝐹)𝑑𝑦 = 0

(𝐴𝑥 + 𝐵𝑦 + 𝐶 )𝑑𝑥 + [𝛼(𝐴𝑥 + 𝐵𝑦) + 𝐹]𝑑𝑦 = 0
El cambio deberá ser: 𝒖 = 𝑨𝒙 + 𝑩𝒚
ii) Rectas no paralelas: Las rectas deben intersectarse, luego desfasarse a un nuevo sistema
de ejes, cuyo punto de intersección sea el origen, el resultado será una ecuación homogéneas
de grado 1, pues ambas rectas pasan por el origen y carecen de termino independiente.
y Ax  By  C  0
y0
Dx  Ey  F  0
x
x0
ℎ = 𝑥 − 𝑥0 𝑑ℎ = 𝑑𝑥
C.V. →
𝑘 = 𝑦 − 𝑦0 𝑑𝑘 = 𝑑𝑦
Problema 11 Resolver: (I-2019)
(3𝑦 2 − 𝑥 + 1)𝑑𝑥 + (2𝑦 3 − 6𝑥𝑦 − 2𝑦)𝑑𝑦 = 0
Solución:
(3𝑦 2 − 𝑥 + 1)𝑑𝑥 + (𝑦 2 − 3𝑥 − 1)2𝑦𝑑𝑦 = 0
C.V. 𝑧 = 𝑦 2 → 𝑑𝑧 = 2𝑦𝑑𝑦
(⏟3𝑧 − 𝑥 + 1) 𝑑𝑥 + (⏟𝑧 − 3𝑥 − 1) 𝑑𝑧 = 0
𝑃 𝑄
Son dos rectas, resolviendo el sistema de ecuaciones:
3z  x  1  0 1 1
z  3x  1  0 →
⏟ 𝑥 = −2 , 𝑧 = −2
 𝑟𝑒𝑠𝑜𝑙𝑣𝑖𝑒𝑛𝑑𝑜
1
ℎ = 𝑥 − (− 2) 𝑑ℎ = 𝑑𝑥
1
→
𝑘 = 𝑧 − (− ) 𝑑𝑘 = 𝑑𝑧
2
1 1 1 1
[3 (𝑘 − ) − (ℎ − ) + 1] 𝑑ℎ + [(𝑘 − ) − 3 (ℎ − ) − 1] 𝑑𝑘 = 0
2 2 2 2
(3𝑘 − ℎ)𝑑ℎ + (𝑘 − 3ℎ)𝑑𝑘 = 0
Ya es homogénea de grado 1
C.V. ℎ = 𝑢𝑘 → 𝑑ℎ = 𝑘𝑑𝑢 + 𝑢𝑑𝑘
(3𝑘 − 𝑢𝑘 )(𝑘𝑑𝑢 + 𝑢𝑑𝑘 ) + (𝑘 − 3𝑢𝑘 )𝑑𝑘 = 0

𝑘[(3 − 𝑢)(𝑘𝑑𝑢 + 𝑢𝑑𝑘 ) + (1 − 3𝑢)𝑑𝑘 ] = 0 → (3𝑢 − 𝑢2 + 1 − 3𝑢)𝑑𝑘 + (3 − 𝑢)𝑘𝑑𝑢 = 0
𝑑𝑘 𝑢−3
∫ +∫ 𝑑𝑢 = ∫ 0
𝑘 (𝑢 − 1)(𝑢 + 1)
−1 2
ln|𝑘| + ∫ ( + ) 𝑑𝑢 = 𝐾
𝑢−1 𝑢+1
ln|𝑘| − ln|𝑢 − 1| + 2 ln|𝑢 + 1| = ln 𝐶
1 1
𝑘 (𝑢 + 1)2 ℎ 𝑥+2 𝑥+2
ln | | = ln 𝐶 → 𝑟𝑒𝑡𝑜𝑟𝑛𝑎𝑛𝑑𝑜 𝑢= = =
𝑢−1 𝑘 𝑧 + 1 𝑦2 + 1
2 2
2
1
1 𝑥 +
(𝑦 2 + 2) ( 2 + 1)
2 1
𝑦 +2
=𝐶 → (𝑥 + 𝑦 2 + 1)2 = 𝐶 (𝑥 − 𝑦 2 )
1
𝑥+2
2 1−1
𝑦 +2
(𝒙 + 𝒚𝟐 + 𝟏)𝟐 = 𝑪(𝒙 − 𝒚𝟐 )
1.7.4 Ecuación Diferencial Exacta
Sea la ecuación diferencial en su forma diferencial:
𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄(𝑥, 𝑦)𝑑𝑦 = 0 … (𝛼)
Si (I) proviene de la diferencial total de primer orden de una función escalar "𝑓" tal que:
𝑓(𝑥, 𝑦) = 𝐶 … (1.16)
Aplicando el diferencial de primer orden sobre (1.16)
𝑑 𝑓(𝑥, 𝑦) = 𝑑 𝐶
𝜕𝑓 𝜕𝑓
𝑑𝑥 + 𝑑𝑦 = 0 … (1.17)
𝜕𝑥 𝜕𝑦
En consecuencia 𝑓 es un campo escalar conservativo, 𝑓 también recibe el nombre de función

potencial.
𝜕𝑓 𝜕𝑓
𝑃(𝑥, 𝑦) = ; 𝑄(𝑥, 𝑦) =
𝜕𝑥 𝜕𝑦
𝜕 𝜕 𝜕𝑓 𝜕2 𝑓 𝜕 𝜕 𝜕𝑓 𝜕2 𝑓
(𝑃) = ( )= 𝑦 (𝑄 ) = ( )=
𝜕𝑦 𝜕𝑦 𝜕𝑥 𝜕𝑦𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑥𝜕𝑦
Pero sabemos que:
𝜕2 𝑓 𝜕2 𝑓 𝜕𝑃 𝜕𝑄
= → =
𝜕𝑦𝜕𝑥 𝜕𝑥𝜕𝑦 𝜕𝑦 𝜕𝑥
Condición que se cumple cuando "𝑓" y sus derivadas parciales son continuas en una región 𝑎 <
𝑥 < 𝑏 , 𝑐 < 𝑦 < 𝑑. En nuestra hipótesis de cobertura asumimos que todas las funciones que
analizamos son suficientemente continuas y diferenciables como para garantizar la validez de las
operaciones que se realizan en ellas.
"𝑓" se construye mediante integración parcial.
𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄(𝑥, 𝑦)𝑑𝑦 = 0 … (𝛼)
Será exacta si se cumple la condición de la diferencial exacta (Condición de Euler):
𝜕𝑃 𝜕𝑃
=
𝜕𝑦 𝜕𝑥
Resolución:
⏟(𝑥, 𝑦)𝑑𝑥 + 𝑄
𝑃 ⏟(𝑥, 𝑦)𝑑𝑦 = 0
𝜕𝑓 𝜕𝑓
𝑑𝑥 𝑑𝑦
𝜕𝑥 𝜕𝑦
𝜕𝑓
= 𝑃(𝑥, 𝑦) → 𝜕𝑓 = 𝑃(𝑥, 𝑦)𝜕𝑥 → 𝑑𝑓 = 𝑃(𝑥, 𝑦)𝑑𝑥 / ∫
𝜕𝑥
∫ 𝑑𝑓 = ∫ 𝑃(𝑥, 𝑦)𝑑𝑥 + Φ(𝑦)
𝑓 = ∫ 𝑃(𝑥, 𝑦)𝑑𝑥 + Φ(𝑦) … (𝛽 )
Derivando (𝛽) respecto de "𝑦"
𝜕𝑓 𝜕 𝜕
= [∫ 𝑃(𝑥, 𝑦)𝑑𝑥 + Φ(𝑦)] → 𝑄= [∫ 𝑃(𝑥, 𝑦)𝑑𝑥] + Φ′ (𝑦)
𝜕𝑦 𝜕𝑦 𝜕𝑦
𝑑Φ(𝑦) 𝜕 𝜕
= 𝑄 − [∫ 𝑃(𝑥, 𝑦)𝑑𝑥] → ∫ 𝑑Φ(𝑦) = ∫ [𝑄 − (∫ 𝑃(𝑥, 𝑦)𝑑𝑥)] 𝑑𝑦
𝑑𝑦 𝜕𝑦 𝜕𝑦
𝜕
Φ(𝑦) = ∫ [𝑄 − (∫ 𝑃(𝑥, 𝑦)𝑑𝑥)] 𝑑𝑦 + 𝐶 … (𝛾)
𝜕𝑦
(𝛾) en (𝛽 )
𝝏
𝒇 = ∫ 𝑷(𝒙, 𝒚)𝒅𝒙 + ∫ [𝑸 − (∫ 𝑷(𝒙, 𝒚)𝒅𝒙)] 𝒅𝒚 + 𝑪 … (1.18)
𝝏𝒚
1.7.4.1 Técnica de diferenciales
El objetivo es el que mediante transformaciones se lleguen a diferenciales totales de primer orden,

para efectuar la cancelación de la integral con el diferencial (técnica de la búsqueda y aplicación
de la antiderivada)
1. 𝒅(𝑢 + 𝑣 ) = 𝒅𝑢 + 𝒅𝑣
2. 𝒅(𝑢 ⋅ 𝑣 ) = 𝑣𝒅𝑢 + 𝑢𝒅𝑣
𝑢 𝑣 𝒅𝑢−𝑢 𝒅𝑣
3. 𝒅 ( ) =
𝑣 𝑣2
4. 𝒅(𝑐 ) = 0
Problema 12 Resolver (𝑥 + 𝑦)𝑑𝑥 + (2𝑦 + 𝑥 )𝑑𝑦 = 0
Solución:
(⏟𝑥 + 𝑦) 𝑑𝑥 + (⏟2𝑦 + 𝑥 ) 𝑑𝑦 = 0
𝑃 𝑄
𝜕𝑃 𝜕𝑄
=1 ; =1
𝜕𝑦 𝜕𝑥
(⏟𝑥 + 𝑦) 𝑑𝑥 + (⏟2𝑦 + 𝑥 ) 𝑑𝑦 = 0 … (𝛼)
𝜕𝑓 𝜕𝑓
𝜕𝑥 𝜕𝑦
𝜕𝑓
=𝑥+𝑦 → 𝜕𝑓 = (𝑥 + 𝑦)𝜕𝑥 → 𝑑𝑓 = (𝑥 + 𝑦)𝑑𝑥 / ∫
𝜕𝑥
𝑥2
∫ 𝑑𝑓 = ∫(𝑥 + 𝑦) 𝑑𝑥 → 𝑓= + 𝑥𝑦 + Φ(𝑦) … (𝛽 )
2
Derivando (𝛽 ) respecto de "𝑦"

𝜕𝑓
= 𝑥 + Φ′ (𝑦) … (𝛾) igualando (𝛼 ) = (𝛾)
𝜕𝑦
2𝑦 + 𝑥 = 𝑥 + Φ′ (𝑦)
𝑑Φ(𝑦)
Φ′ (𝑦) = 2𝑦 → = 2𝑦 → ∫ 𝑑Φ(𝑦) = ∫ 2𝑦 𝑑𝑦 → Φ(𝑦) = 𝑦 2 − 𝐶 𝑒𝑛 (𝛼)
𝑑𝑦
𝑥2 𝑥2
𝑓= + 𝑥𝑦 + 𝑦 2 − 𝐶 = 0 → + 𝑥𝑦 + 𝑦 2 = 𝐶
2 2
𝒙𝟐
+ 𝒙𝒚 + 𝒚𝟐 = 𝑪
𝟐
Empleando el método de diferenciales

𝑥𝑑𝑥 + 𝑦𝑑𝑥
⏟ + 2𝑦𝑑𝑦 + 𝑥𝑑𝑦
⏟ =0 → (𝑦𝑑𝑥 + 𝑥𝑑𝑦) + 𝑥𝑑𝑥 + 2𝑦𝑑𝑦 = 0 / ∫
// //
𝒙𝟐 𝟐𝒚𝟐 𝒙𝟐
1 ∫ 𝒅(𝒙𝒚) + ∫ 𝒙 𝒅𝒙 + ∫ 𝟐𝒚 𝒅𝒚 = ∫ 𝟎 → 𝒙𝒚 + 𝟐
+ 𝟐
=𝑪 → 𝟐
+ 𝒙𝒚 + 𝒚𝟐 = 𝑪
Problema 13 Resolver: (invierno 2019)
[𝑦 + 𝑥 2 (𝑥𝑦 + 2)]𝑑𝑥 + [𝑥 + 4𝑦 3 (𝑥𝑦 + 2)]𝑑𝑦 = 0
Solución: Distribuyendo para buscar diferenciales.
𝑦𝑑𝑥 + 𝑥 2 (𝑥𝑦 + 2)𝑑𝑥 + 𝑥𝑑𝑦 + 4𝑦 3 (𝑥𝑦 + 2)𝑑𝑦 = 0

𝑦𝑑𝑥 + 𝑥𝑑𝑦 + (𝑥𝑦 + 2)(𝑥 2 𝑑𝑥 + 4𝑦 3 𝑑𝑦) = 0 / ÷ (𝑥𝑦 + 2)
𝑦𝑑𝑥 + 𝑥𝑑𝑦
+ 𝑥 2 𝑑𝑥 + 4𝑦 3 𝑑𝑦 = 0 / ∫
𝑥𝑦 + 2
1 1 𝑦𝑑𝑥+𝑥𝑑𝑦
Sabemos que 𝑑 [ln(𝑥𝑦 + 2)] = ⋅ 𝑑(𝑥𝑦 + 2) = ⋅ (𝑑(𝑥𝑦) + 𝑑(2)) =
𝑥𝑦+2 𝑥𝑦+2 𝑥𝑦+2
∫ 𝑑[ln(𝑥𝑦 + 2)] + ∫ 𝑥 2 𝑑𝑥 + ∫ 4𝑦 3 𝑑𝑦 = ∫ 0
𝒙𝟑
𝐥𝐧(𝒙𝒚 + 𝟐) + + 𝒚𝟒 = 𝑪
𝟑
1.7.4.2 Independencia de la Trayectoria: En la integral
∫ 𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄 (𝑥, 𝑦)𝑑𝑦

𝐶
𝜕𝑃 𝜕𝑄
Al cumplirse la condición de una diferencial exacta: = es independiente de la trayectoria,
𝜕𝑦 𝜕𝑥
en consecuencia, no depende de la trayectoria el resultado solo depende de las condiciones inicial

y final, por comodidad se toman rectas horizontales y verticales.
y P (x , y )
C2
C1
y0
x0 x x
∫ 𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄(𝑥, 𝑦)𝑑𝑦 = ∫ 𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄(𝑥, 𝑦)𝑑𝑦 + ∫ 𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄 (𝑥, 𝑦)𝑑𝑦
𝐶 𝐶1 𝐶2
Para la curva 𝐶1 : 𝑦 = 𝑦0 recta horizontal 𝑑𝑦 = 0
Para la curva 𝐶2 : 𝑥 = 𝑥 vertical 𝑑𝑥 = 0 (𝑃(𝑥, 𝑦) es un punto general cualquiera)

∫ 𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄 (𝑥, 𝑦)𝑑𝑦 = ∫ 𝑃(𝑥, 𝑦0 )𝑑𝑥 + 𝑄(𝑥, 𝑦0 ) ⋅ 0 + ∫ 𝑃(𝑥, 𝑦)0 + 𝑄(𝑥, 𝑦)𝑑𝑦
𝐶 𝐶1 𝐶2
𝑥 𝑦
∫ 𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄(𝑥, 𝑦)𝑑𝑦 = ∫ 𝑃 (𝑥, 𝑦0 )𝑑𝑥 + ∫ 𝑄(𝑥, 𝑦)𝑑𝑦
𝐶 𝑥0 𝑦0
El punto 𝑃0 (𝑥0 , 𝑦0 ) es arbitrario, puede elegirse a conveniencia, aprovechamos este resultado para
resolver la ecuación diferencial exacta:
𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄 (𝑥, 𝑦)𝑑𝑦 = 0
Integramos ambos miembros, llegando a la fórmula:

𝒙 𝒚
∫ 𝑷(𝒙, 𝒚𝟎 )𝒅𝒙 + ∫ 𝑸(𝒙, 𝒚)𝒅𝒚 = 𝑪 … (1.19)
𝒙𝟎 𝒚𝟎
Si se emplea el otro camino se llega también a la fórmula:

𝒙 𝒚
∫ 𝑷(𝒙, 𝒚)𝒅𝒙 + ∫ 𝑸(𝒙𝟎 , 𝒚)𝒅𝒚 = 𝑪 … (1.20)
𝒙𝟎 𝒚𝟎
Problema 14 Hallar el valor de "𝐾" de modo que la ecuación diferencial sea exacta y luego
resuelva dicha ecuación.
(𝑦 3 + 𝑘𝑥𝑦 4 − 2𝑥 )𝑑𝑥 + (3𝑥𝑦 2 + 20𝑥 2 𝑦 3 )𝑑𝑦 = 0
Solución: Aplicando la condición de diferencial exacta
𝜕𝑃 𝜕𝑄
(𝑦 3 + 𝑘𝑥𝑦 4 − 2𝑥 )𝑑𝑥 + ⏟
⏟ (3𝑥𝑦 2 + 20𝑥 2 𝑦 3 ) 𝑑𝑦 = 0 → =
𝜕𝑦 𝜕𝑥
𝑃 𝑄
3𝑦 2 + 4𝑘𝑥𝑦 3 = 3𝑦 2 + 40𝑥𝑦 3 → 𝑘 = 10 en la ecuación diferencial
(𝑦 3 + 10𝑥𝑦 4 − 2𝑥 )𝑑𝑥 + (3𝑥𝑦 2 + 20𝑥 2 𝑦 3 )𝑑𝑦 = 0
Empleando independencia de la trayectoria, tomando como 𝑃0 (𝑥0 , 𝑦0 ) = 𝑃0 (0,0)

𝒙 𝒚
∫ 𝑷(𝒙, 𝒚𝟎 )𝒅𝒙 + ∫ 𝑸(𝒙, 𝒚)𝒅𝒚 = 𝑪
𝒙𝟎 𝒚𝟎
𝒙 𝑦
∫ (0 + 10𝑥 ⋅ 0 − 2𝑥 )𝑑𝑥 + ∫ (3𝑥𝑦 2 + 20𝑥 2 𝑦 3 )𝑑𝑦 = 𝐶
3 4
𝟎 0
𝑦
−𝑥 2 |0𝑥 + (𝑥𝑦 3 + 5𝑥 2 𝑦 4 )|0 = 𝐶
𝟓𝒙𝟐 𝒚𝟒 + 𝒙𝒚𝟑 − 𝒙𝟐 = 𝑪
1.7.5 Factor Integrante
Es notorio que las ecuaciones diferenciales exactas son raras, la exactitud depende de su
equilibrio preciso mismo que se puede destruir con facilidad. Si la Ecuación diferencial:
𝑃(𝑥, 𝑦)𝑑𝑥 + 𝑄 (𝑥, 𝑦)𝑑𝑦 = 0
No es exacta, puede serlo si multiplicamos esta por una función 𝑢 = 𝑢(𝑥, 𝑦) con la propiedad
de que ahora esta nueva expresión ya es una diferencial exacta:
𝑢𝑃(𝑥, 𝑦) 𝑑𝑥 + ⏟
⏟ 𝑢𝑄 (𝑥, 𝑦) 𝑑𝑦 = 0
𝑀 𝑁
Sabemos que debe cumplirse:
𝜕𝑀 𝜕𝑁
= … (1)
𝜕𝑦 𝜕𝑥
 Caso1: Si asumimos que 𝑢 = 𝑢(𝑥)
𝜕 𝜕𝑃
(𝑢𝑃) = 𝑢 …. (2)
𝜕𝑦 𝜕𝑦
𝜕 𝜕𝑢 𝜕𝑄
(𝑢𝑄 ) = 𝑄+𝑢 …. (3)
𝜕𝑥 𝜕𝑥 𝜕𝑥
Reemplazando (2) y (3) en (1)
𝜕𝑃 𝜕𝑄
𝜕𝑃 𝜕𝑢 𝜕𝑄 𝜕𝑢 𝜕𝑃 𝜕𝑄 𝑑𝑢 ( 𝜕𝑦 − 𝜕𝑥 )
𝑢 = 𝑄+𝑢 → 𝑄 = 𝑢( − ) → = 𝑑𝑥 / ∫
𝜕𝑦 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑥 𝑢 𝑄
𝜕𝑃 𝜕𝑄 𝜕𝑃 𝜕𝑄
𝑑𝑢 ( − ) ( − )
∫ = ∫ 𝜕𝑦𝑄 𝜕𝑥 𝑑𝑥 → ln|𝑢| = ∫ 𝜕𝑦 𝜕𝑥
𝑑𝑥
𝑢 𝑄
𝝏𝑷 𝝏𝑸
(𝝏𝒚 − 𝝏𝒙 )
∫⏟ 𝑸 𝒅𝒙
𝒖(𝒙) = 𝒆 𝒇(𝒙) … (1.21)
 Caso 2: Si asumimos que 𝑢 = 𝑢(𝑦), trabajando de forma similar, el factor integrante

que depende de "𝑦" será.
𝝏𝑸 𝝏𝑷
( 𝝏𝒙 − 𝝏𝒚 )
∫⏟ 𝑷 𝒅𝒚
𝒖(𝒚) = 𝒆 𝒇(𝒚) … (1.22)
 Caso 3: Si asumimos que 𝑢 = 𝑢(𝑧) ∧ 𝑧 = 𝑓(𝑥, 𝑦)
x
u z
y
𝜕 𝜕
(𝑢𝑃) = 𝑢𝑧 𝑧𝑦 𝑃 + 𝑢𝑃𝑦 ; (𝑢𝑄 ) = 𝑢𝑧 𝑧𝑥 𝑄 + 𝑢𝑄𝑥
𝜕𝑦 𝜕𝑥
Igualando tenemos: 𝑢𝑧 𝑧𝑦 𝑃 + 𝑢𝑃𝑦 = 𝑢𝑧 𝑧𝑥 𝑄 + 𝑢𝑄𝑥 → 𝑢𝑧 (𝑃𝑧𝑦 − 𝑄𝑧𝑥 ) = 𝑢(𝑄𝑥 − 𝑃𝑦 )
𝑑𝑢 (𝑄𝑥 − 𝑃𝑦 ) (𝑄𝑥 − 𝑃𝑦 )
= 𝑑𝑧 / ∫ → ln|𝑢| = ∫ 𝑑𝑧
𝑢 (𝑃𝑧𝑦 − 𝑄𝑧𝑥 ) (𝑃𝑧𝑦 − 𝑄𝑧𝑥 )
𝝏𝑸 𝝏𝑷
( − )
𝝏𝒙 𝝏𝒚
∫ 𝝏𝒛 𝝏𝒛𝒅𝒛
𝑷 −𝑸
⏟𝝏𝒚 𝝏𝒙
𝒖(𝒛) = 𝒆 𝒇(𝒛) … (1.23)
Probar con los siguientes cambios: 𝑧 = 𝑥 + 𝑦 , 𝑧 = 𝑥𝑦 , 𝑧 = 𝑥 2 + 𝑦 2
4𝑥𝑦 2 𝑑𝑥 + (3𝑥 2 𝑦 − 1)𝑑𝑦 = 0 ; 𝑦(1) = 3
Solución: Verificando si se trata de una ecuación diferencial exacta
4𝑥𝑦 2 𝑑𝑥 + ⏟
⏟ (3𝑥 2 𝑦 − 1) 𝑑𝑦 = 0
𝑃 𝑄
𝜕𝑃 𝜕𝑄
= 8𝑥𝑦 ; = 6𝑥𝑦
𝜕𝑦 𝜕𝑥
Buscamos un factor integrante:

𝜕𝑃 𝜕𝑄 𝜕𝑄 𝜕𝑃
( − ) ( − )
𝜕𝑦 𝜕𝑥 𝜕𝑥 𝜕𝑦
∫⏟ 𝑄 𝑑𝑥 ∫⏟ 𝑃 𝑑𝑦
𝑢(𝑥 ) = 𝑒 𝑓(𝑥) , 𝑢(𝑦) = 𝑒 𝑓(𝑦)
Elegimos el factor integrante que depende de “𝑦”

(6𝑥𝑦−8𝑥𝑦) 1 1 1 1
∫ 𝑑𝑦 1 −
= 𝑦 −2
4𝑥𝑦 2 − ∫ 𝑑𝑦 − ⋅ln 𝑦 ln 𝑦 2
𝑢(𝑦) = 𝑒 =𝑒 2 𝑦 =𝑒 2 =𝑒
Multiplicamos a la ecuación diferencial original por nuestro factor integrante

1
4𝑥𝑦 2 𝑑𝑥 + (3𝑥 2 𝑦 − 1)𝑑𝑦 = 0 / ⋅ 𝑦 −2
3 1 1
4𝑥𝑦 2 𝑑𝑥 + (3𝑥 2 𝑦 2 − 𝑦 −2 ) 𝑑𝑦 = 0
Verificando la condición de una diferencial exacta

3 1 1
4𝑥𝑦 2 𝑑𝑥 + (3𝑥 2 𝑦 2 − 𝑦 −2 ) 𝑑𝑦 = 0
⏟ ⏟
𝑀
𝑁
𝜕𝑀 𝜕𝑁
= → 6𝑥𝑦1/2 = 6𝑥𝑦1/2 𝑣𝑒𝑟𝑖𝑓𝑖𝑐𝑎
𝜕𝑦 𝜕𝑥
Calculando la función de la que proviene

3 1 1
4𝑥𝑦
⏟ 2 𝑑𝑥 + (3𝑥 𝑦 − 𝑦 −2 ) 𝑑𝑦 = 0 … (𝛼)
2 2
⏟
𝜕𝑓
𝜕𝑥 𝜕𝑓
𝜕𝑦
𝜕𝑓 1 1 1 1 1 1
= 3𝑥 2 𝑦 2 − 𝑦 −2 → 𝑑𝑓 = (3𝑥 2 𝑦 2 − 𝑦 −2 ) 𝑑𝑦 / ∫ → ∫ 𝑑𝑓 = ∫ (3𝑥 2 𝑦 2 − 𝑦 −2 ) 𝑑𝑦
𝜕𝑦
3 1
𝑦2 𝑦2 3 1
𝑓 = 3𝑥 2 − + Φ(𝑥) → 𝑓 = 2𝑥 2 𝑦 2 − 2𝑦 2 + Φ(𝑥) … (𝛽)
3 1
2 2
𝜕𝑓 3
Derivando (𝛽) respecto de la variable "𝑥" 𝜕𝑥
= 4𝑥𝑦 2 + Φ′ (𝑥 ) … (𝛾)
Igualando (𝛼 ) y (𝛾)
3 3 𝑑Φ(𝑥)
4𝑥𝑦 2 = 4𝑥𝑦 2 + Φ′ (𝑥 ) → Φ′ (𝑥 ) = 0 → =0 → ∫ 𝑑Φ(𝑥) = ∫ 0
𝑑𝑥
Φ(𝑥) = −𝐾 en (𝛼)
3 1 3 1 1 𝐾
𝑓 = 2𝑥 2 𝑦 2 − 2𝑦 2 − 𝐾 = 0 → 2𝑥 2 𝑦 2 − 2𝑦 2 = 𝐾 → 𝑦 2 (𝑥 2 𝑦 − 1) = / ( )2
2
𝒚(𝒙𝟐 𝒚 − 𝟏)𝟐 = 𝑪
Problema 16 Resuelva la ecuación diferencial
(𝑥𝑦 + 2𝑥𝑦 ln2 𝑥 + 𝑥 ln 𝑥 )𝑑𝑦 + (2𝑦 2 ln 𝑥 + 𝑦)𝑑𝑥 = 0
Solución:
Tratando de darle la forma de una diferencial exacta
(2𝑦 2 ln 𝑥 + 𝑦) 𝑑𝑥 + ⏟
⏟ (𝑥𝑦 + 2𝑥𝑦 ln2 𝑥 + 𝑥 ln 𝑥 ) 𝑑𝑦 = 0
𝑃 𝑄
𝜕𝑃
= 4𝑦 ln 𝑥 + 1
𝜕𝑦
𝜕𝑄 1 𝜕𝑄
= 𝑦 + 2𝑦 (1 ⋅ ln2 𝑥 + 𝑥2 ln 𝑥 ⋅ ) + ln 𝑥 + 1 → = 𝑦 + 2𝑦 ln2 𝑥 + 4𝑦 ln 𝑥 + ln 𝑥 + 1
𝜕𝑥 𝑥 𝜕𝑥
𝜕𝑃 𝜕𝑄
− 𝑦+2𝑦 ln2 𝑥+ln 𝑥 −(𝑦+2𝑦 ln2 𝑥+ln 𝑥)
𝜕𝑦 𝜕𝑥
∫ 𝑄 𝑑𝑥 ∫ 𝑑𝑥 ∫ 𝑑𝑥 1 −1 1
𝑢(𝑥 ) = 𝑒 = 𝑒 𝑥𝑦+2𝑥𝑦 ln2 𝑥+𝑥 ln 𝑥 = 𝑒 𝑥(𝑦+2𝑦 ln2 𝑥+ln 𝑥) = 𝑒 − ∫𝑥𝑑𝑥 = 𝑒 ln 𝑥 =
𝑥
(2𝑦 2 ln 𝑥 + 𝑦)𝑑𝑥 + (𝑥𝑦 + 2𝑥𝑦 ln2 𝑥 + 𝑥 ln 𝑥 )𝑑𝑦 = 0 / ⋅ 1
𝑥
2𝑦 2 ln 𝑥 𝑦
( (𝑦 + 2𝑦 ln2 𝑥 + ln 𝑥 ) 𝑑𝑦 = 0 … (1)
+ ) 𝑑𝑥 + ⏟
⏟ 𝑥 𝑥 𝜕𝑓
𝜕𝑓 𝜕𝑦
𝜕𝑥
𝜕𝑓
= 𝑦 + 2𝑦 ln2 𝑥 + ln 𝑥 → ∫ 𝑑𝑓 = ∫(𝑦 + 2𝑦 ln2 𝑥 + ln 𝑥 ) 𝑑𝑦
𝜕𝑦
𝑦2
𝑓= + 𝑦 2 ln2 𝑥 + 𝑦 ln 𝑥 + Φ(𝑥) … (2)
2
Derivando (2) respecto de "𝑥" y luego comparando con (1)
𝜕𝑓 2𝑦 2 ln 𝑥 𝑦 2𝑦 2 ln 𝑥 𝑦 2𝑦 2 ln 𝑥 𝑦
= + + Φ′ (𝑥 ) → ′( )
+ +Φ 𝑥 = +
𝜕𝑥 𝑥 𝑥 𝑥 𝑥 𝑥 𝑥
𝑑Φ
Φ′ (𝑥 ) = 0 → =0 → ∫ 𝑑Φ = ∫ 0 𝑑𝑥 → Φ = −𝐶
𝑑𝑥
𝑦2
Sustituyendo en (2) 𝑓= 2
+ 𝑦 2 ln2 𝑥 + 𝑦 ln 𝑥 − 𝐶 = 0
𝒚𝟐
+ 𝒚𝟐 𝐥𝐧𝟐 𝒙 + 𝒚 𝐥𝐧 𝒙 = 𝑪
𝟐
5 1 52 1 7 3
Problema 17 Resolver la ecuación diferencial: (2 x 2 y  y )dx  ( x 2  x y 2 )dy  0
2 2
Si se conoce que admite un factor integrante de la forma    ( xm y n )
Solución:
Primera Forma: Tomando directamente un factor integrante de la forma 𝑢 = 𝑥 𝑚 𝑦 𝑛
1 1
(2𝑥 5/2 𝑦 + 𝑦 5/2 ) 𝑑𝑥 − ( 𝑥 7/2 + 𝑥𝑦 3/2 ) 𝑑𝑦 = 0 / ⋅ 𝑥𝑚𝑦𝑛
2 2
5 1 5 1 7 3
(2𝑥 2+𝑚 𝑦1+𝑛 + 𝑥 𝑚 𝑦 2+𝑛 ) 𝑑𝑥 + (− 𝑥 2+𝑚 𝑦 𝑛 − 𝑥 1+𝑚 𝑦 2+𝑛 ) 𝑑𝑦 = 0
⏟ 2 ⏟ 2
𝑃 𝑄
𝜕𝑃 5 1 5 3 𝜕𝑄 1 7 5 3
= 2(1 + 𝑛)𝑥 2+𝑚 𝑦 𝑛 + 2 ( + 𝑛) 𝑥 𝑚 𝑦 2+𝑛 ; = −2 ( + 𝑚) 𝑥 2+𝑚 𝑦 𝑛 − (1 + 𝑚)𝑥 𝑚 𝑦 2+𝑛
𝜕𝑦 2 𝜕𝑥 2
𝜕𝑃 𝜕𝑄
Comparando coeficientes para que se cumpla: =
𝜕𝑦 𝜕𝑥
7
2(1 + 𝑛) = −12 (2 + 𝑚) … (1) 𝑛 9
{1 5 de (2) 𝑚 = − 2 − 4 … (3)
( + 𝑛) = −(1 + 𝑚) … (2)
2 2
7 𝑛 9 1
(3) en (1) 2(1 + 𝑛) = −12 ( + − − ) → 2(1 + 𝑛) = − (14 − 2𝑛 − 9)
2 2 4 8
𝟑 𝟑
16 + 16𝑛 = −5 + 2𝑛 → 14𝑛 = −21 → 𝒏 = −𝟐 𝑒𝑛 (3) 𝒎 = −𝟐
El factor integrante es: 𝑢 = 𝑥 −3/2 𝑦 −3/2
1 1
(2𝑥 5/2 𝑦 + 𝑦 5/2 ) 𝑑𝑥 − ( 𝑥 7/2 + 𝑥𝑦 3/2 ) 𝑑𝑦 = 0 / ⋅ 𝑥 −3/2 𝑦 −3/2
2 2
1 1 3
(2𝑥𝑦 −1/2 + 𝑥 −3/2 𝑦) 𝑑𝑥 + (− 𝑥 2 𝑦 −2 − 𝑥 −1/2 ) 𝑑𝑦 = 0
⏟ 2 ⏟ 2
𝑃 𝑄
Empleando la fórmula de la independencia de la trayectoria, tomando en forma arbitraria y

cuidando que haya división entre cero tomamos (𝑥0 , 𝑦0 ) = (1,1)
𝑥 𝑦
∫ 𝑃(𝑥, 𝑦0 ) 𝑑𝑥 + ∫ 𝑄(𝑥, 𝑦) 𝑑𝑦 = 𝐶
𝑥0 𝑦0
𝑥 𝑦
1 1 3
∫ (2𝑥 ⋅ 1−1/2 + 𝑥 −3/2 ⋅ 1) 𝑑𝑥 + ∫ (− 𝑥 2 𝑦 −2 − 𝑥 −1/2 ) 𝑑𝑦 = 𝐶
1 2 1 2
1 𝑥 1 𝑦
− − 1 𝑥 𝑦
1 𝑥 2 1 𝑦 2
2
(𝑥 + 2 ⋅ 1 )| + (− 2 𝑥 2 1 −𝑥 −1/2
𝑦)| = 𝐶 → (𝑥 2 − 𝑥 −2 )| + (𝑥 2 𝑦 −1/2 − 𝑥 −1/2 𝑦)|1 = 𝐶
−2 −2 1
1 1
1 1 1 1 1
(𝑥 2 − 𝑥 −2 ) − (12 − 1−2 ) + (𝑥 2 𝑦 −2 − 𝑥 −2 𝑦) − (𝑥 2 − 𝑥 −2 ) = 𝐶
𝒙𝟐 𝒚
− =𝑪
√𝒚 √𝒙
Problema 18 Resolver (𝑥 2 𝑦 3 + 2𝑥𝑦 2 + 𝑦)𝑑𝑥 + (𝑥 3 𝑦 2 − 2𝑥 2 𝑦 + 𝑥 )𝑑𝑦
Solución: Procedimiento propuesto por Univ. Jhovani Brian Casas Huanca, empleo de
diferenciales, agrupando y transformando convenientemente la expresión:
𝑥 2 𝑦 3 𝑑𝑥 + 2𝑥𝑦 2 𝑑𝑥 + 𝑦𝑑𝑥 + 𝑥 3 𝑦 2 𝑑𝑦 − 2𝑥 2 𝑦𝑑𝑦 + 𝑥𝑑𝑦 = 0

𝑥 2 𝑦 2 (𝑦𝑑𝑥 + 𝑥𝑑𝑦) + (𝑦𝑑𝑥 + 𝑥𝑑𝑦) + 2𝑥𝑦(𝑦𝑑𝑥 − 𝑥𝑑𝑦) = 0
[(𝑥𝑦)2 + 1](𝑦𝑑𝑥 + 𝑥𝑑𝑦) + 2𝑥𝑦(𝑦𝑑𝑥 − 𝑥𝑑𝑦) = 0 / ÷ (𝑥𝑦)2
(𝑥𝑦)2 + 1 𝑥𝑦 𝑦𝑑𝑥 − 𝑥𝑑𝑦
[ ] (𝑦𝑑𝑥 + 𝑥𝑑𝑦) + 2 2 ( ) = 0 / ÷ (𝑥𝑦)2
(𝑥𝑦) 2 𝑥 𝑦2
1 1 𝑦𝑑𝑥 − 𝑥𝑑𝑦
[1 + 2
] (𝑦𝑑𝑥 + 𝑥𝑑𝑦) + 2 𝑥 ( )=0
(𝑥𝑦) 𝑦2
𝑦
1 1 𝑥
∫ [1 + ] 𝑑 (𝑥𝑦 ) + 2 ∫ 𝑥 𝑑 ( ) = ∫0
(𝑥𝑦)2 𝑦
𝑦
𝟏 𝒙
𝒙𝒚 − + 𝟐 𝐥𝐧 ( ) = 𝑪
𝒙𝒚 𝒚
1.8 Métodos de Resolución de Ecuaciones diferenciales de primer orden en la forma normal

𝑑𝑦
= 𝑓 (𝑥, 𝑦) ó 𝑦 ′ = 𝑓(𝑥, 𝑦)
𝑑𝑥
1.8.1 Ecuación Lineal
𝒅𝒚
+ 𝑃(𝑥)𝒚 = 𝑄(𝑥 ) … (1)
𝒅𝒙
Existen los siguientes métodos de resolución:
 Factor integrante
 Variación de parámetros
 Método de Bellman
i) Factor integrante: Buscamos a un factor integrante "𝑢" tal que:
𝑑𝑦 𝑑𝑦
+ 𝑃(𝑥 )𝑦 = 𝑄 (𝑥 ) /⋅ 𝑢 → 𝑢⋅ + 𝑢 ⋅ 𝑃(𝑥 )𝑦 = 𝑢 ⋅ 𝑄 (𝑥 )
𝑑𝑥 𝑑𝑥
𝑢 ⋅ 𝑦′ + ⏟
𝑢 ⋅ 𝑃(𝑥 ) 𝑦 = 𝑢 ⋅ 𝑄(𝑥 ) … (2)
//
Sabemos que: (𝑢𝑦)′ = 𝑢′ 𝑦 + 𝑢𝑦 ′ → (𝑢𝑦)′ = 𝑢𝑦 ′ + 𝑢⏟′ 𝑦 … (3)

//
Comparando (2) con (3)
𝑑𝑢 𝑑𝑢
𝑢′ = 𝑢𝑃(𝑥 ) → = 𝑢𝑃(𝑥 ) / ∫ = ∫ 𝑃(𝑥) 𝑑𝑥 → ln 𝑢 = ∫ 𝑃(𝑥) 𝑑𝑥
𝑑𝑥 𝑢
𝒖(𝒙) = 𝒆∫ 𝑷(𝒙)𝒅𝒙 … (1.566)
𝑢 ⋅ 𝑦 ′ + 𝑢 ⋅ 𝑃(𝑥 )𝑦 = 𝑢 ⋅ 𝑄 (𝑥 )
⏟
𝑑
(𝑢⋅𝑦)
𝑑𝑥
𝑑(𝑢𝑦) = 𝑢𝑄 (𝑥 )𝑑𝑥 → ∫ 𝑑(𝑢𝑦) = ∫ 𝑢𝑄(𝑥 )𝑑𝑥
𝒖𝒚 = ∫ 𝒖𝑸(𝒙)𝒅𝒙 + 𝑪 … (1.24)
ii) Por Variación de Parámetros:
𝑑𝑦
+ 𝑃(𝑥)𝑦 = 𝑄(𝑥 ) … (𝛼)
𝑑𝑥
Buscamos una solución para la ecuación homogénea de (𝛼), la misma resulta haciendo 𝑄 (𝑥 ) = 0
𝑑𝑦 𝑑𝑦 𝑑𝑦
+ 𝑃(𝑥 )𝑦 = 0 → = −𝑃(𝑥 )𝑑𝑥 → ∫ = − ∫ 𝑃(𝑥 ) 𝑑𝑥 → ln 𝑦 = − ∫ 𝑃(𝑥 ) 𝑑𝑥 + 𝑘
𝑑𝑥 𝑦 𝑦
𝑦 = 𝑒 − ∫ 𝑃(𝑥)𝑑𝑥+𝑘 = 𝑒 − ∫ 𝑃(𝑥)𝑑𝑥 𝑒 𝑘 → 𝑦 = 𝐾𝑒 − ∫ 𝑃(𝑥)𝑑𝑥 … (𝛽)
Ahora asumimos que la constante es función de la variable independiente, es decir:
𝑦 = 𝐾(𝑥 )𝑒 − ∫ 𝑃(𝑥)𝑑𝑥 … (𝛾)
Derivando (3)
𝑦 ′ = 𝐾 ′(𝑥 )𝑒 − ∫ 𝑃(𝑥)𝑑𝑥 + 𝐾(𝑥)𝑒 − ∫ 𝑃(𝑥)𝑑𝑥 [−𝑃(𝑥)] reemplazando en (𝛼)
𝐾 ′ (𝑥 )𝑒 − ∫ 𝑃(𝑥)𝑑𝑥 + 𝐾(𝑥 )𝑒 − ∫ 𝑃(𝑥)𝑑𝑥 [−𝑃(𝑥 )] + 𝑃(𝑥 )𝐾(𝑥 )𝑒 − ∫ 𝑃(𝑥)𝑑𝑥 = 𝑄(𝑥 )

𝑑𝐾(𝑥 )
𝐾 ′ (𝑥 )𝑒 − ∫ 𝑃(𝑥)𝑑𝑥 = 𝑄(𝑥 ) → = 𝑒 ∫ 𝑃(𝑥)𝑑𝑥 𝑄(𝑥 ) → 𝑑𝐾(𝑥 ) = 𝑒 ∫ 𝑃(𝑥)𝑑𝑥 𝑄(𝑥 )𝑑𝑥 →
𝑑𝑥
∫ 𝑑𝐾(𝑥 ) = ∫ 𝑒 ∫ 𝑃(𝑥)𝑑𝑥 𝑄(𝑥 )𝑑𝑥 → 𝐾(𝑥 ) = ∫ 𝑒 ∫ 𝑃(𝑥)𝑑𝑥 𝑄(𝑥 )𝑑𝑥 + 𝐶 𝑒𝑛 (𝛾)
𝑦 = [∫ 𝑒 ∫ 𝑃(𝑥)𝑑𝑥 𝑄(𝑥 )𝑑𝑥 + 𝐶] 𝑒 − ∫ 𝑃(𝑥)𝑑𝑥
𝒚 = 𝒆− ∫ 𝑷(𝒙)𝒅𝒙 [∫ 𝒆∫ 𝑷(𝒙)𝒅𝒙 𝑸(𝒙)𝒅𝒙 + 𝑪] … (1.25)
iii) Método de Bellman: En realidad el método es el de sustitución de Lagrange

𝑑𝑦
+ 𝑃(𝑥)𝑦 = 𝑄(𝑥 ) … (𝛼)
𝑑𝑥
Se considera que "y" es el producto de dos funciones de "x".
𝑦 = 𝑢(𝑥 )𝑣(𝑥 ) … (1.26)

𝑦 = 𝑢𝑣 → 𝑦 ′ = 𝑢′ 𝑣 + 𝑢𝑣´
Reemplazando esta última expresión en (𝛼)
𝑢′ 𝑣 + 𝑢𝑣´ + 𝑃(𝑥 )𝑢𝑣 = 𝑄 (𝑥 )

𝑢𝑣 ′ + [𝑢′ + 𝑃(𝑥)𝑢]𝑣 = 𝑄 (𝑥 ) … (𝛽)
𝑑𝑦
Asumiendo que "𝑢" es solución de la ecuación diferencial homogénea: 𝑑𝑥
+ 𝑃(𝑥 )𝑦 = 0
𝑑𝑢 𝑑𝑢
𝑢′ + 𝑃(𝑥 )𝑢 = 0 → = −𝑃(𝑥 )𝑢 → = −𝑃(𝑥 )𝑑𝑥 → ln 𝑢 = − ∫ 𝑃(𝑥 ) 𝑑𝑥
𝑑𝑥 𝑢
𝑢(𝑥 ) = 𝑒 − ∫ 𝑃(𝑥)𝑑𝑥 … (1.27)
𝑑𝑣 𝑄(𝑥)
Retomando (1.56) 𝑢𝑣 ′ + [⏟𝑢′ + 𝑃(𝑥)𝑢] 𝑣 = 𝑄 (𝑥 ) → 𝑢𝑣 ′ = 𝑄 (𝑥 ) → 𝑑𝑥
= 𝑢
0
𝑑𝑣 𝑄 (𝑥 ) 𝑑𝑣 𝑄 (𝑥 ) 𝑑𝑣
= → = − ∫ 𝑃(𝑥)𝑑𝑥 → = 𝑒 ∫ 𝑃(𝑥)𝑑𝑥 𝑄(𝑥 ) → ∫ 𝑑𝑣 = ∫ 𝑒 ∫ 𝑃(𝑥)𝑑𝑥 𝑄(𝑥 ) 𝑑𝑥
𝑑𝑥 𝑢 𝑑𝑥 𝑒 𝑑𝑥
𝑣(𝑥 ) = ∫ 𝑒 ∫ 𝑃(𝑥)𝑑𝑥 𝑄(𝑥 ) 𝑑𝑥 + 𝐶 … (1.28)
Reemplazando (1.27) y (1.28) en (1.26)
𝒚 = 𝒆− ∫ 𝑷(𝒙)𝒅𝒙 [∫ 𝒆∫ 𝑷(𝒙)𝒅𝒙 𝑸(𝒙) 𝒅𝒙 + 𝑪] … (1.28)
𝒅𝒚 2
+ 𝒚 = 𝑥2
𝒅𝒙 𝑥
Solución: Empleando el método de factor integrante
2
𝑢 = 𝑒 ∫𝑥𝑑𝑥 = 𝑒 2 ln 𝑥 = 𝑥 2
𝒅𝒚 2 𝑑 2
+ 𝒚 = 𝑥2 / ⋅ 𝑥2 → 𝑥 2 𝒚′ + 2𝑥𝒚 = 𝑥 4 → (𝑥 𝑦) = 𝑥 4 → 𝑑(𝑥 2 𝑦) = 𝑥 4 𝑑𝑥
𝒅𝒙 𝑥 𝑑𝑥
2 4 2
𝑥5
∫ 𝑑(𝑥 𝑦) = ∫ 𝑥 𝑑𝑥 → 𝑥 𝑦= +𝐶 →
5
𝟏 𝟑 𝑪
𝒚= 𝒙 + 𝟐
𝟓 𝒙
Problema 20 Resolver:
1
𝑦′ =
𝑥 sin 𝑦 + 2 sin 2𝑦
Solución: Dándole la forma de ecuación lineal
𝑑𝑦 1 𝑑𝑥 𝑑𝑥
= → = 𝑥 sin 𝑦 + 2 sin 2𝑦 → − sin 𝑦 𝑥 = 2 sin 2𝑦
𝑑𝑥 𝑥 sin 𝑦 + 2 sin 2𝑦 𝑑𝑦 𝑑𝑦
Empleando la fórmula de variación de parámetros
𝑑𝑥
+ 𝑃(𝑦)𝑥 = 𝑄(𝑥 ) → 𝑥 = 𝑒 − ∫ 𝑃(𝑦)𝑑𝑦 [∫ 𝑒 ∫ 𝑃(𝑦)𝑑𝑦 𝑄(𝑦) 𝑑𝑦 + 𝐶]
𝑑𝑦
𝑥 = 𝑒 − ∫ − sin 𝑦𝑑𝑦 [∫ 𝑒 ∫ − sin 𝑦𝑑𝑦 2 sin 2𝑦 𝑑𝑦 + 𝐶] = 𝑒 − cos 𝑦 [4 ∫ 𝑒 cos 𝑦 sin 𝑦 cos 𝑦 𝑑𝑦 + 𝐶]
Cambio de variable: 𝑡 = cos 𝑦 → 𝑑𝑡 = − sin 𝑦 𝑑𝑦
𝑥 = 𝑒 − cos 𝑦 [−4 ∫ 𝑒 𝑡 𝑡 + 𝐶] = 𝑒 − cos 𝑦 [−4𝑒 𝑡 (𝑡 − 1) + 𝐶 ] = 𝑒 − cos 𝑦 [−4𝑒 cos 𝑦 (cos 𝑦 − 1) + 𝐶 ]
𝒙 = −𝟒(𝐜𝐨𝐬 𝒚 − 𝟏) + 𝑪𝒆− 𝐜𝐨𝐬 𝒚
Problema 21 Resuelva la ecuación diferencial

(1 seny )dx 2y cosy x (secy tgy ) dy
Solución: Llevando a la forma normal
dx dx secy tgy 2y cosy

(1 seny ) x (secy tgy ) 2y cosy x ecuación lineal
dy dy 1 seny 1 seny
dx 1 seny 2y cosy dx 2y cosy

x x sec y
dy cos y(1 seny ) 1 seny dy 1 seny
sec y dy tgy )
Factor integrante: u e e ln(secy sec y tg y
2y cosy
d (sec y tg y ) x (sec y tg y )dy (sec y tg y ) x 2ydy
1 seny
(sec y tg y ) x y2 C
Segunda forma: Factor diferencial exacta, factor integrante
P Q
(1 seny )dx x (secy tgy ) 2y cosy dy 0 cos y ; secy tgy
y x
secy tgy coy 1 seny cos2 y (1 seny ) seny

dy dy dy
1 seny cosy(1 seny ) cosy(1 seny )
u e e e e ln(secy) sec y
(1 seny ) sec ydx x (secy tgy ) sec y 2y dy 0

1 seny
f (1 seny ) sec ydx f x (1 seny ) sec y (y ) x( ) (y ) ..(1)
cos y
f
x cos y sec y (1 seny ) sec ytgy '(y )
y
f seny cos2 y sen 2y seny
x 1 (1 seny ) '(y ) x '(y )
y 2 2
cos y cos y
f 1 seny f 1 seny
x '(y ) igualando x 2y '(y) 2y '(y ) y2
y 2 y 2
cos y cos y
en (1)
1 seny
f x( ) y2 C
cos y
1.8.2 Ecuación diferencial de Bernoulli
Nicolaus Senior
1623-1708
Jacob I Nicolaus I Johannes I

1654-1705 1662-1716 1667-1748
Nicolaus II Nicolaus III Daniel Johannes II

1687-1759 1695-1726 1700-1782 1710-1790
Johannes III Jacob II

1746-1807 1759-1789
La familia Bernoulli produjo ocho matemáticos, muchos de ellos sobresalientes, Los hermanos
Jacob I (hermano mayor) y Johannes I (hermano menor) asistieron a una conferencia de
Gottfried Wilhelm Leibniz y quedaron fascinados por su nuevo e innovador método propuesto,
aún no bautizado como cálculo diferencial el integral, desde entonces se convirtieron en sus
seguidores y sentían antipatía por Newton, los tres potenciaron esta nueva rama de la
matemática más adelante conocida como análisis matemático.
Jacob I Bernoulli descubrió al número "𝒆" en un curioso e hipotético caso de estudio del interés
compuesto en un banco que paga un tasa del 100% anual, este ejercicio le oriento a
identificar que era mejor sacar el dinero a los seis meses porque como el banco paga al 100%
anual la ganancia sería mayor si se asume que se colocan 100 Bs al banco al finalizar el año se
tendrá 100 Bs de interés en total 200 Bs, pero si se saca a medio año el dinero se tendrá 50 Bs
de ganancia de interés en total se tendría los 100 Bs iniciales más los 50 Bs total 150 Bs pero
que sin embargo aún hay medio año (6 meses) para seguir ganando intereses de los 150 Bs que
ahora se colocarán en el banco, cuyo interés será la mitad de los 150 Bs por estar colocados solo
los restantes 6 meses entonces se ganarán 75 Bs más dando un total de 225 Bs, como se puede
apreciar las ganancias son mayores, si se procede colocar el dinero al en la mañana y retirarlo al
finalizar el día se generarían más ganancias.
1 año 100(1 + 1) = 200
Medio año (dos colocaciones) 1 2

100 (1 + ) = 225
2
Trimestre (cuatro colocaciones) 1 4
100 (1 + ) = 244,14
4
Al día (365 colocaciones) 1 365
100 (1 + ) = 271,45
365
Sin embargo notó que por más que ingresará y retirará el dinero al minuto, hay un límite de
usura que precisamente es el número 𝑒.
1 n
𝑒 = lim (1 + )
n→∞ n
Johannes I Bernoulli, matemático tal vez más prolífico que su hermano mayor Jacob I, siempre
sintió que vivió a su sombra sin embargo no era la realidad, Johannes de un gran talento fue
profesor del joven Márquez Guillaume de L’Hopital, con el firmó el convenio de enviarle sus
descubrimientos en Matemática y que él podía hacer lo que quisiera con ellos a cambio de un
salario regular, uno de los resultados de acuerdo produjo el siguiente descubrimiento por parte
de Johannes:
𝑓(𝑥) 𝑓′(𝑥)
𝑙𝑖𝑚 = 𝑙𝑖𝑚
𝑥→𝑎 𝑔(𝑥) 𝑥→𝑎 𝑔′(𝑥)
Que sin desmerecer la gran capacidad y talento matemático de L´Hopital quien entrego esta
valiosa regla a las matemáticas fue Johannes Bernoulli, pero que sin embargo al ser publicada
por el Márquez de L´Hopital se quedó con ese nombre como: Regla de L´Hopital. Posteriormente
L´Hopital aclaró que la regla no era suya sin embargo ya todo el mundo la conocía por su
nombre, situación que produjo enorme ira en Johannes Bernoulli que siempre tuvo un carácter
receloso con todo el mundo incluso con su hermano y su hijo Daniel al sentir que ellos iban
siendo reconocidos y premiados por la sociedad sin importar que a él también se lo reconocía,
pero él lo quería todo para él, demando y acuso a L´Hopital de plagió, pero las cosas ya estaban
hechas.
Daniel Bernoulli, un brillante Matemático dela tercera generación gran amigo de Euler, sus
trabajos abarcan calculo, ecuaciones diferenciales, probabilidades, pero es más recordado por su
trabajo de la conservación de la energía en hidrodinámica y mecánica de fluidos, reconocido como
fundador de la física matemática y el último de los más famosos Bernoulli.
El tópico que veremos concerniente a las ecuaciones diferenciales fue propuesto por Jacob I
Bernoulli referida a la ecuación 𝒚′ + 𝑷(𝒙)𝒚 = 𝑸(𝒙)𝒚𝒏 , cuya solución fue propuesta por su
hermano Johannes con la sustitución 𝒛 = 𝒚𝟏−𝒏 que reduce la ecuación a una ecuación lineal.
Procedimiento
Posee la forma:
𝒅𝒚
+ 𝑃(𝑥 )𝒚 = 𝑄 (𝑥 ) ⋅ 𝑦 𝑛 … (1.29)
𝒅𝒙
Se debe eliminar 𝑦 𝑛 multiplicando ambos miembro por 𝑦 −𝑛
𝒅𝒚 𝒅𝒚
+ 𝑃(𝑥 )𝒚 = 𝑄 (𝑥 ) ⋅ 𝑦 𝑛 / ⋅ 𝑦 −𝑛 → 𝑦 −𝑛 𝟏−𝒏
+ 𝑃(𝑥 ) 𝒚⏟ = 𝑄(𝑥 ) … (𝛼)
𝒅𝒙 𝒅𝒙
𝒖
Luego el cambio de variables es en:
𝑑
𝑢 = 𝑦1−𝑛 /
𝑑𝑥
𝑑𝑢 𝑑𝑦 𝑑𝑦 1 𝑑𝑢
= (1 − 𝑛)𝑦 −𝑛 → 𝑦 −𝑛 = ⋅ 𝑒𝑛 (𝛼)
𝑑𝑥 𝑑𝑥 𝑑𝑥 1 − 𝑛 𝑑𝑥
1 𝑑𝑢
⋅ + 𝑃(𝑥 )𝑢 = 𝑄(𝑥 ) / ⋅ (1 − 𝑛)
1 − 𝑛 𝑑𝑥
𝒅𝒖
+ (1 − 𝑛)𝑃(𝑥 )𝒖 = 𝑄(𝑥 )
⏟
𝒅𝒙
𝒆𝒄𝒖𝒂𝒄𝒊ó𝒏 𝒍𝒊𝒏𝒆𝒂𝒍
4𝑥𝑦 2 𝑑𝑥 + (3𝑥 2 𝑦 − 1)𝑑𝑦 = 0 ; 𝑦(1) = 3
Solución: generando la forma de Bernoulli, pero tomando como variable dependiente a "𝑥"
𝒅𝒙
+ 𝑃(𝑦)𝒙 = 𝑄(𝑦) ⋅ 𝑥 𝑛
𝒅𝒚
𝑑𝑥 𝒅𝒙 3 1 𝑑𝑥 3 2 1
4𝑥𝑦 2 + 3𝑥 2 𝑦 − 1 = 0 → + 𝒙 = 2 ⋅ 𝑥 −1 / ⋅ 𝑥 → 𝑥 + ⏟𝑥 = 2
𝑑𝑦 𝒅𝒚 4𝑦 4𝑦 𝑑𝑦 4𝑦 𝑢 4𝑦
Cambio de variable: 𝑢 = 𝑥 2
𝑑 𝑑𝑢 𝑑𝑥
𝑢 = 𝑥2 / → = 2𝑥 en la ecuación diferencial
𝑑𝑦 𝑑𝑦 𝑑𝑦
1 𝑑𝑢 3 1 𝑑𝑢 3 1 3 3 1
− ∫2𝑦 𝑑𝑦 ∫2𝑦 𝑑𝑦
+ 𝑢= 2 → + 𝑢= 2 → 𝑢=𝑒 [∫ 𝑒 𝑑𝑦 + 𝐶]
2 𝑑𝑦 4𝑦 4𝑦 𝑑𝑦 2𝑦 2𝑦 2𝑦 2
1
3 3 1 −2 1
3 3 3 1 𝑦2
𝑢= 𝑒 −2 ln 𝑦 [∫ 𝑒 2 ln 𝑦 𝑑𝑦 + 𝐶] = 𝑦 [ ∫ 𝑦 2 ⋅ 𝑦 −2 𝑑𝑦 + 𝐶] = 𝑦 −2 [ ⋅ + 𝐶]
2𝑦 2 2 2 1
2
Retornando a la variable original:
3 1 3 1 1
𝑥 2 = 𝑦 −2 [𝑦 2 + 𝐶] → 𝑥 2𝑦2 = 𝑦2 + 𝐶 → 𝑦 2 (𝑥 2 𝑦 − 1) = 𝐶 / ( )2
𝑦(𝑥 2 𝑦 − 1)2 = 𝐾
Evaluando: 3(12 ⋅ 3 − 1)2 = 𝐾 → 𝐾 = 12
𝒚(𝒙𝟐 𝒚 − 𝟏)𝟐 = 𝟏𝟐
𝑑𝑦 𝑥
sec 2 𝑦 + tan 𝑦 sec 𝑦 + = −𝑥 sec 2 𝑦
𝑑𝑥 cos 𝑦
Solución: Transformando la ecuación diferencial
1 𝑑𝑦 sin 𝑦 𝑥 𝑥
2
+ 2
+ =− / ⋅ cos 2 𝑦
cos 𝑦 𝑑𝑥 cos 𝑦 cos 𝑦 cos 2 𝑦
𝑑𝑦 𝑑𝑦
+ sin 𝑦 + 𝑥 cos 𝑦 + 𝑥 = 0 → + sin 𝑦 + 𝑥(cos 𝑦 + 1) = 0
𝑑𝑥 𝑑𝑥
Empleando la identidad de ángulo medio: cos 𝑦 = cos 2(𝑦2) = cos 2 𝑦2 − sin2 𝑦2
𝑑𝑦 𝑦
+ 2 sin 𝑦2 cos2 + 2𝑥 cos2 𝑦2 = 0 / ÷ cos2 𝑦2 → sec 2 𝑦2 𝒚′ + 2 tan 𝑦2 + 2𝑥 = 0
𝑑𝑥
𝑑 1
Cambio de variable: 𝑢 = tan 𝑦2 / 𝑑𝑥 → 𝑢′ = sec 2 𝑦2 ⋅ 2 𝑦′
2𝑢′ + 2𝑢 = −2𝑥 → 𝒖′ + 𝑢 = −𝑥 ecuación lineal, por variación de parámetros:
𝑢 = 𝑒 − ∫ 𝑑𝑥 [∫ 𝑒 ∫ 𝑑𝑥 (−𝑥 )𝑑𝑥 + 𝐶] → 𝑢 = 𝑒 −𝑥 [− ∫ 𝑒 𝑥 𝑥 𝑑𝑥 + 𝐶 ] → 𝑢 = 𝑒 −𝑥 [−𝑒 𝑥 (𝑥 − 1) + 𝐶 ]
𝑢 = (1 − 𝑥 ) + 𝐶𝑒 −𝑥 retornando a la variable original

𝒚
𝐭𝐚𝐧 ( ) = (𝟏 − 𝒙) + 𝑪𝒆−𝒙
𝟐
1.8.3 Ecuación diferencial de Riccati
Dentro el conjunto de matemáticos italianos notables resalta Jacobo Francesco conde de

Riccati, natural de Venecia (1676-1754) quien dio a conocer los trabajos de Newton en Italia,
interesado por la acústica, Riccati quiso resolver ciertas ecuaciones diferenciales de segundo
orden, considerando curvas cuyos radios de curvatura dependen sólo de las ordenadas, llego a la
ecuación:
𝑑 2𝑥 𝑑2𝑦
𝑚
𝑑𝑦 2
𝑥 = +( )
𝑑𝑝2 𝑑𝑝2 𝑑𝑝
Que con un cambio de variable apropiado la redujo a una ecuación de primer orden:
𝑚
𝑑𝑞 𝑑𝑥 𝑢2
𝑥 = +
𝑑𝑥 𝑑𝑥 𝑞
Llamando a "𝑞" a una potencia de "𝑥", 𝑥 𝑛 por ejemplo se obtuvo la forma:
𝑑𝑢 𝑢2
+ = 𝑛𝑥 𝑚+𝑛−1
𝑑𝑥 𝑥 𝑛
Operando se llegó a:
𝑑𝑢
= 𝐴(𝑥 ) + 𝐵(𝑥 )𝑢 + 𝐶 (𝑥 )𝑢2
𝑑𝑥
La ecuación de Riccati fue importante porque permitió reducir las ecuaciones de segundo orden
a ecuaciones diferenciales de primer orden.
Procedimiento:
𝒅𝒚
+ 𝑃(𝑥 )𝒚𝟐 + 𝑄 (𝑥 )𝒚 = 𝑅(𝑥 ) … (1.30)
𝒅𝒙
Para resolver esta ecuación diferencial es necesario conocer una solución particular de la misma
que deberá obtenerse por tanteo, pero de acuerdo a los coeficientes variables estas generalmente
son:
 Al conocer una solución particular 𝒚𝟏 de (1.30) se efectúa el cambio:
1
𝑦 = 𝑦1 + … (1.31)
𝑧
𝑧′
Derivando y reemplazando en (I) 𝑦 ′ = 𝑦1′ − 𝑧 2
𝑧′ 1 2 1
𝑦1′ − 2 + 𝑃(𝑥 ) (𝑦1 + ) + 𝑄 (𝑥 ) (𝑦1 + ) = 𝑅 (𝑥 )
𝑧 𝑧 𝑧
𝑧′ 2𝑦1 1 1
𝑦1′ − 2
+ 𝑃(𝑥 ) (𝑦12 + + 2 ) + 𝑄 (𝑥 ) (𝑦1 + ) = 𝑅(𝑥 )
𝑧 𝑧 𝑧 𝑧
𝑧′ 2𝑦1 1 1
(𝑦1′ + 𝑃(𝑥 )𝑦12 + 𝑄 (𝑥 )𝑦1 ) − 2
+ 𝑃(𝑥 ) ( + 2 ) + 𝑄(𝑥 ) = 𝑅(𝑥 )
𝑧 𝑧 𝑧 𝑧
𝑧′ 2𝑦1 𝑧 + 1 𝑄 (𝑥 )
𝑅(𝑥 ) − + 𝑃 (𝑥 ) [ ] + = 𝑅(𝑥 ) / ⋅ (−𝑧 2 )
𝑧2 𝑧2 𝑧
𝑧 ′ − [2𝑦1 𝑃(𝑥 ) + 𝑄(𝑥 )]𝑧 = 𝑃(𝑥)
⏟
𝑒𝑐𝑢𝑎𝑐𝑖ó𝑛 𝑙𝑖𝑛𝑒𝑎𝑙
Sugerencia, para los tanteos de la solución particular, según los coeficientes de la variable 𝑦, 𝑦′
se presenta la siguiente tabla:
Coeficientes Solución Particular
𝑥 , 𝑥2 , 𝑥3 , 𝑥4 , … 𝑦𝑝 = 𝑎0 + 𝑎1 𝑥
𝑦𝑝 = 𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2
𝑦𝑝 = 𝑎𝑥 𝑛
1 1 1 𝑎1
, , 𝑦𝑝 = 𝑎0 +
𝑥 𝑥2 𝑥4 𝑥
𝑦𝑝 = 𝐴 sin 𝑥
sin 𝑥 , cos 𝑥 𝑦𝑝 = 𝐴 cos 𝑥
𝑦𝑝 = 𝐴 sin 𝑥 + 𝐵 cos 𝑥
tan 𝑥 , sec 𝑥 𝑦𝑝 = 𝐴 tan 𝑥
𝑦𝑝 = 𝐴 sec 𝑥
𝑦𝑝 = 𝐴 tan 𝑥 + 𝐵 sec 𝑥
𝑒 𝑥 , 𝑒 2𝑥 , 𝑒 3𝑥 , … 𝑦𝑝 = 𝐴𝑒 𝛼𝑥
1.7.4 Formulas asociadas si se conocen más de una solución particular
i) 2 soluciones particulares: Si se conocen 𝑦1 , 𝑦2 como soluciones particulares
𝑦 ′ + 𝑃(𝑥 )𝑦 2 + 𝑄 (𝑥 )𝑦 = 𝑅(𝑥 ) … (𝐼)
Ambas soluciones cumplen la ecuación diferencial
𝑦1′ + 𝑃(𝑥 )𝑦12 + 𝑄 (𝑥 )𝑦1 = 𝑅 (𝑥 ) … (1′)

𝑦2′ + 𝑃(𝑥 )𝑦22 + 𝑄 (𝑥 )𝑦2 = 𝑅(𝑥 ) … (2′)
Restando en forma separada (I) con (1’) y (2’)
′
𝑦 ′ − 𝑦1′
𝑦 − 𝑦1′ + 𝑃(𝑥 )[𝑦 2 − 𝑦12 ] + 𝑄(𝑥 )[𝑦 − 𝑦1 ] = 0 → + 𝑃(𝑥 )[𝑦 + 𝑦1 ] + 𝑄(𝑥 ) = 0 … (3′)
𝑦 − 𝑦1
𝑦 ′ − 𝑦2′
𝑦 ′ − 𝑦2′ + 𝑃(𝑥 )[𝑦 2 − 𝑦22 ] + 𝑄 (𝑥 )[𝑦 − 𝑦2 ] = 0 → + 𝑃(𝑥 )[𝑦 + 𝑦2 ] + 𝑄 (𝑥 ) = 0 … (4′)
𝑦 − 𝑦2
Restando (3’) y (4’)
𝑦 ′ − 𝑦1′ 𝑦 ′ − 𝑦2′ 𝑑(𝑦 − 𝑦1 ) 𝑑(𝑦 − 𝑦2 )

− = −𝑃(𝑥 )[𝑦1 − 𝑦2 ] / → ∫ −∫ = − ∫ 𝑃(𝑥 )[𝑦1 − 𝑦2 ] 𝑑𝑥
𝑦 − 𝑦1 𝑦 − 𝑦2 𝑦 − 𝑦1 𝑦 − 𝑦2
𝑦 − 𝑦1
ln|𝑦 − 𝑦1 | − ln|𝑦 − 𝑦2 | = − ∫ 𝑃(𝑥 )[𝑦1 − 𝑦2 ] 𝑑𝑥 + 𝑘 → ln | | = − ∫ 𝑃(𝑥 )[𝑦1 − 𝑦2 ] 𝑑𝑥 + 𝑘
𝑦 − 𝑦2
𝒚 − 𝒚𝟏
= 𝑪𝒆− ∫ 𝑷(𝒙)[𝒚𝟏−𝒚𝟐]𝒅𝒙 … (1.32)
𝒚 − 𝒚𝟐
ii) 3 soluciones particulares: Si se conocen 𝑦1 , 𝑦2 , 𝑦3 como soluciones particulares
Si se conocen 3 soluciones particulares podemos usar la fórmula (1.32)

𝑦 − 𝑦1
= 𝐶1 𝑒 − ∫ 𝑃(𝑥)[𝑦1−𝑦2]𝑑𝑥 … (5′)
𝑦 − 𝑦2
Evaluando 𝑦3 en la fórmula (5’)

𝑦3 − 𝑦1
= 𝐶2 𝑒 − ∫ 𝑃(𝑥)[𝑦1−𝑦2]𝑑𝑥 … (6′)
𝑦3 − 𝑦2
Dividiendo (5’) entre (6’)

𝒚 − 𝒚𝟏 𝒚𝟑 − 𝒚𝟏
= 𝑪( ) … (1.33)
𝒚 − 𝒚𝟐 𝒚𝟑 − 𝒚𝟐
Problema 24: Resuelva la ecuación diferencial: (I-2005)
𝑦 ′ = 𝑒 2𝑥 𝑦 2 − 2𝑦 − 9𝑒 −2𝑥 ; 𝑦(0) = 4
Solución: Tanteando una solución particular: 𝑦𝑝 = 𝐴𝑒 𝛼𝑥 → 𝑦𝑝′ = 𝐴𝛼𝑒 𝛼𝑥 en la ecuación:
𝐴𝛼𝑒 𝛼𝑥 = 𝑒 2𝑥 (𝐴𝑒 𝛼𝑥 )2 − 2𝐴𝑒 𝛼𝑥 − 9𝑒 −2𝑥 → 𝐴𝛼𝑒 𝛼𝑥 = 𝐴2 𝑒 2𝑥+2𝛼𝑥 − 2𝐴𝑒 𝛼𝑥 − 9𝑒 −2𝑥 →
𝐴2 𝑒 2𝑥+2𝛼𝑥 − (2𝐴 + 𝐴𝛼)𝑒 𝛼𝑥 − 9𝑒 −2𝑥 = 0 por comparación
2𝐴 + 𝐴𝛼 = 0 → 𝐴(2 + 𝛼 ) = 0 → 𝛼 = −2
𝐴2 𝑒 2𝑥+2𝛼𝑥 − 9𝑒 −2𝑥 = 0 → 𝐴2 𝑒 2𝑥+2𝛼𝑥 = 9𝑒 −2𝑥 → 𝛼 = −2 ∧ 𝐴 = ±3
Nota: 𝐴 no puede ser cero porque eliminaría a la solución particular, además de no cumplir con
la ecuación.
′
En consecuencia, se tienen dos soluciones particulares. 𝑦𝑝1 = 3𝑒 −2𝑥 ; ′
𝑦𝑝2 = −3𝑒 −2𝑥
Tomando la primera solución particular el cambio de variable será:
1 𝑧′
𝑦 = 3𝑒 −2𝑥 + … (1) → 𝑦 ′ = −6𝑒 −2𝑥 −
𝑧 𝑧2
𝑧′ 1 2 1
−6𝑒 −2𝑥 − 2
= 𝑒 2𝑥
(3𝑒 −2𝑥
+ ) − 2 (3𝑒 −2𝑥 + ) − 9𝑒 −2𝑥 →
𝑧 𝑧 𝑧
−2𝑥
𝑧′ 2𝑥 −4𝑥
6𝑒 −2𝑥 1 1
−6𝑒 − 2 = 𝑒 (9𝑒 + + 2 ) − 2 (3𝑒 −2𝑥 + ) − 9𝑒 −2𝑥 →
𝑧 𝑧 𝑧 𝑧
−2𝑥
𝑧′ −2𝑥
6 𝑒 2𝑥 2
−6𝑒 − 2 = 9𝑒 + + 2 − 6𝑒 −2𝑥 − − 9𝑒 −2𝑥 →
𝑧 𝑧 𝑧 𝑧
𝑧 ′ 4 𝑒 2𝑥
− 2 = + 2 /⋅ (−𝑧 2 ) → 𝑧 ′ + 4𝑧 = −𝑒 2𝑥 /⋅ 𝑢 = 𝑒 ∫ 4𝑑𝑥 = 𝑒 4𝑥
𝑧 𝑧 𝑧
𝑒 4𝑥 𝑧 ′ + 4𝑒 4𝑥 𝑧 = −𝑒 6𝑥 → 𝑑(𝑒 4𝑥 𝑧) = −𝑒 6𝑥 / ∫ → ∫ 𝑑(𝑒 4𝑥 𝑧) = − ∫ 𝑒 6𝑥 𝑑𝑥
𝑒 6𝑥 1
𝑒 4𝑥 𝑧 = − +𝐶 → 𝑧 = − 𝑒 2𝑥 + 𝐶𝑒 −4𝑥
6 6
Reemplazando en (1)
1
𝑦 = 3𝑒 −2𝑥 +
1
− 𝑒 2𝑥 + 𝐶𝑒 −4𝑥
6
1 1 7
Evaluando en la condición inicial 𝑦(0) = 4 → 4 = 3 + 1 → − +𝐶 =1 → 𝐶 =
−6+𝐶 6 6
𝟔
𝒚 = 𝟑𝒆−𝟐𝒙 +
𝟕𝒆−𝟒𝒙 − 𝒆𝟐𝒙
Veamos cómo nos hubiera ido si usábamos la fórmula para la solución de una ecuación diferencial
de Riccati al tener dos soluciones particulares.
𝑦 − 𝑦1
𝑦 ′ + 𝑃(𝑥 )𝑦 2 + 𝑄(𝑥 )𝑦 = 𝑅 (𝑥 ) → 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛: = 𝐶𝑒 − ∫ 𝑃(𝑥)[𝑦1 −𝑦2]𝑑𝑥
𝑦 − 𝑦2
𝑦 ′ = 𝑒 2𝑥 𝑦 2 − 2𝑦 − 9𝑒 −2𝑥 ; 𝑦(0) = 4 → 𝑦 ′ − 𝑒 2𝑥 𝑦 2 + 2𝑦 = −9𝑒 −2𝑥 ; 𝑦(0) = 4
′
Las soluciones particulares: 𝑦𝑝1 = 3𝑒 −2𝑥 ; ′
𝑦𝑝2 = −3𝑒 −2𝑥
𝑦 − 3𝑒 −2𝑥 2𝑥 −2𝑥 −2𝑥 𝑦 − 3𝑒 −2𝑥

−2𝑥
= 𝐶𝑒 − ∫ −𝑒 [3𝑒 −(−3𝑒 )]𝑑𝑥 → = 𝐶𝑒 ∫ 6𝑑𝑥
𝑦 + 3𝑒 𝑦 + 3𝑒 −2𝑥
𝑦 − 3𝑒 −2𝑥 4−3 1
−2𝑥
= 𝐶𝑒 6𝑥 𝑒𝑣𝑎𝑙𝑢𝑎𝑛𝑑𝑜 =𝐶 → 𝐶=
𝑦 + 3𝑒 4+3 7
𝒚 − 𝟑𝒆−𝟐𝒙 𝟏 𝟔𝒙
= 𝒆
𝒚 + 𝟑𝒆−𝟐𝒙 𝟕
Sin embargo, en la mayoría de los casos mediante transformaciones (operaciones algebraicas) es

posible llegar a la forma de la solución del primer método elegido.
−2𝑥 6𝑥 4𝑥
21𝑒 −2𝑥 + 3𝑒 4𝑥 21𝑒 −2𝑥 + 3𝑒 4𝑥 𝑒 −4𝑥
7𝑦 − 21𝑒 = 𝑦𝑒 + 3𝑒 → 𝑦= → 𝑦= ⋅ −4𝑥
7 − 𝑒 6𝑥 7 − 𝑒 6𝑥 𝑒
21𝑒 −6𝑥 − 3 + 3 + 3 3𝑒 −2𝑥 (7𝑒 −4𝑥 − 𝑒 2𝑥 ) + 6
𝑦= =
7𝑒 −4𝑥 − 𝑒 2𝑥 7𝑒 −4𝑥 − 𝑒 2𝑥
𝟔
𝒚 = 𝟑𝒆−𝟐𝒙 +
𝟕𝒆−𝟒𝒙 − 𝒆𝟐𝒙
Problema 25 Resuelva la ecuación diferencial: (I-2019)
dy
(2x 4 2x 3 ) (x 6 x 5 )y 2 (2x 4 8x 2 )y 4x 8 0
dx
3
Si se conoce que tiene una solución de la forma: y x
4
Solución: y ' 3 x
(2x 4 2x 3 )( 3 x 4
) (x 6 x 5 )( 2x 6
) (2x 4 8x 2 )( x 3
) 4x 8 0
1 2 2 1 1 2 1
6 6 x x 2 x 8 x 4x 8 0 ( 6 8) (2 )x (2 4)x 0
3 1 4 z'
Verifica para 2 Ecuación de Riccati y 2x y' 6x
z z2
3
z' x 1 1
(2x 4 2x 3 )( 6x 4
) (x 6 x 5 )(4x 6
4 ) (2x 4 8x 2 )(2x 3
) 4x 8 0
2 z 2 z
z z
1 z' 4x 3 4x 2 x6 x5 2x 4 8x 2
12 12x (2x 4 2x 3 ) 4 4x 1
4x 16x 1
4x 8 0
z2 z z2 z
2x 3 (x 1)z ' (2x 4 4x 3 4x 2 )z x 5 (x 1) 0
x2 2x 2 x 2 (x 1)
z' z Ecuación. Lineal
x2 x 2(x 1)
x2 2x 2 x2 x x 2 x 2 2 1
frac. par . 1 1
2 2 x (x 1) x x 1
x x x x
2 1 2 1
(1 )dx (1 )dx x 2 (x 1) x2 x 1 x 2 (x 1)
z e x x 1 e x x 1 dx C ex e x
dx C
2(x 1) x 1 x2 2(x 1)
x2 x 1 x2 x
z ex e x
C ex e x
C
x 1 2 x 1 2
x3 x2 x 3 2Cx 2e x en el cambio
z Ce x ( )
2(x 1) x 1 2(x 1)
3 2(x 1)
y 2x
x 3
2Cx 2e x
Problema 26 Resolver la ecuación diferencial (I-2021):
𝑥 2 𝑦 ′ − 𝑥 2 𝑦 2 + (2𝑥𝑒 𝑥 − 𝑥 2 )𝑦 + 𝑒 𝑥 − 𝑒 2𝑥 = 0
𝑒 𝑎𝑥
Sabiendo que una solución primicial es de la forma: 𝑦1 =
𝑥
Solución: Reemplazamos la solución particular propuesta para hallar el valor de "𝑎":
𝑒 𝑎𝑥 𝑎𝑒 𝑎𝑥 𝑥 − 𝑒 𝑎𝑥
𝑦1 = → 𝑦1′ =
𝑥 𝑥2
𝑎𝑒 𝑎𝑥 𝑥 − 𝑒 𝑎𝑥 𝑒 𝑎𝑥 2 𝑒 𝑎𝑥
𝑥2 − 𝑥 2
( ) + (2𝑥𝑒 𝑥
− 𝑥 2)
+ 𝑒 𝑥 − 𝑒 2𝑥 = 0
𝑥2 𝑥 𝑥
𝑎𝑒 𝑎𝑥 𝑥 − 𝑒 𝑎𝑥 − 𝑒 2𝑎𝑥 + (2𝑒 𝑥 − 𝑥 )𝑒 𝑎𝑥 + 𝑒 𝑥 − 𝑒 2𝑥 = 0
Igualando términos semejantes:
(𝑎 − 1)𝑥𝑒 𝑎𝑥 + (𝑒 𝑥 − 𝑒 𝑎𝑥 ) + [2𝑒 (𝑎+1)𝑥 − 𝑒 2𝑎𝑥 − 𝑒 2𝑥 ] = 0
Para que la ecuación verifique el valor de la incógnita debe ser: 𝒂 = 𝟏

𝑒𝑥
Entonces la solución particular es: 𝑦1 = 𝑥
Efectuando el cambio de variable y reemplazando en la ecuación diferencial:
𝒆𝒙 𝟏 𝑒 𝑥𝑥 − 𝑒 𝑥 𝑧′
𝒚= + … (𝛼) → 𝑦′ = − 2
𝒙 𝒛 𝑥2 𝑧
𝑒 𝑥𝑥 − 𝑒 𝑥 𝑧′ 𝑒𝑥 1 2 𝑒𝑥 1
𝑥 2( − ) − 𝑥 2
( + ) + (2𝑥𝑒 𝑥
− 𝑥 2)
( + ) + 𝑒 𝑥 − 𝑒 2𝑥 = 0
𝑥2 𝑧2 𝑥 𝑧 𝑥 𝑧
Simplificando términos semejantes:
𝑧′ 𝑒 2𝑥 𝑒𝑥 1 2𝑥𝑒 𝑥 𝑥 2
𝑒 𝑥𝑥 − 𝑒 𝑥 − 𝑥2 − 𝑥 2
( + 2 + ) + 2𝑒 2𝑥
− 𝑥𝑒 𝑥
+ − + 𝑒 𝑥 − 𝑒 2𝑥 = 0
𝑧2 𝑥2 𝑥𝑧 𝑧 2 𝑧 𝑧
𝑧′ 𝑥2 𝑥2 𝑧2
−𝑥 2 2
− 2− = 0 / ⋅ (−𝑥 2)
𝑧 𝑧 𝑧
𝑧′ + 𝑧 = −1 Ecuación lineal o de variables separables
𝑑𝑧
+ (𝑧 + 1) = 0
𝑑𝑥
𝑑𝑧
∫ + ∫ 𝑑𝑥 = ∫ 0 → 𝑙𝑛(𝑧 + 1) + 𝑥 = 𝐶 → 𝑧 = 𝐶𝑒 −𝑥 − 1
𝑧+1
𝑒𝑥 1
Reemplazando en (𝛼) 𝑦= + 𝐶𝑒 −𝑥−1
𝑥
𝒆𝒙 𝒆𝒙
𝒚= +
𝒙 𝑪 − 𝒆𝒙
1.8 Envolventes y Soluciones Singulares
Antes de estudiar las ecuaciones no resueltas en la derivada, haremos una breve digresión
geométrica. Consideremos el haz de curvas de la figura, cuya ecuación es 𝜑(𝑥, 𝑦, 𝐶 ) = 0. Como
puede observarse, existe una curva, E, que no pertenece a la familia (y, por tanto, no satisface
dicha ecuación uniparamétrica), pero que es tangente en cada punto a una curva de la familia.
Una curva con esta propiedad se llama envolvente del haz de curvas. Como sabemos, al derivar
la ecuación del haz y eliminar el parámetro C obtenemos la ecuación diferencial de la
familia F (x, y, y′) = 0 que, a diferencia de lo que ocurre con la ecuación finita, sí será satisfecha
por la envolvente, ya que en cada punto (x, y) su pendiente y′ es la misma que la de la curva
del haz a la que allí es tangente. Será por tanto una solución singular tanto en el sentido dado
hasta ahora a ese concepto, puesto que no pertenece a la solución general como en el más
restrictivo usado a veces de ser una solución con la propiedad de que en todos sus puntos se
infringe la propiedad de unicidad por pasar por ellos más de una solución con la misma
tangente.
Como veremos, ésta es una forma típica, aunque no exclusiva, en que aparecen las soluciones
singulares de ecuaciones no resueltas en la derivada. Para que podamos comprobar este extremo,
veamos cómo puede calcularse la ecuación de la envolvente. El punto P de la figura está en dos
curvas del haz correspondientes a los valores C y C + ∆C del parámetro; por tanto, satisface
las ecuaciones de ambas curvas 𝜑(𝑥, 𝑦, 𝐶 ) = 0 y 𝜑(𝑥, 𝑦, 𝐶 + Δ𝐶 ) = 0, o combinaciones lineales
independientes de estas dos ecuaciones como
𝜑(𝑥, 𝑦, 𝐶 ) = 0 𝜑(𝑥, 𝑦, 𝐶 ) = 0 … (1′ )

{𝜑(𝑥, 𝑦, 𝐶 + 𝛥𝐶 ) − 𝜑(𝑥, 𝑦, 𝐶 ) → {𝑑
=0 𝜑(𝑥, 𝑦, 𝐶 ) = 0 … (2′ )
𝛥𝐶 𝑑𝐶
En el límite ∆C → 0 las dos curvas tenderán a confundirse y el punto P irá a la envolvente que,
por tanto, satisfará las ecuaciones que se obtienen tomando dicho límite.
La ecuación finita de la envolvente se obtiene eliminando el parámetro C entre las dos ecuaciones
anteriores, la envolvente no satisface la ecuación de la familia para ningún valor fijo de C,
𝜕2 𝜑
suponiendo además de la regularidad que ≠ 0 puede en principio usarse la ecuación (2’) para
𝜕𝐶 2
obtener el valor C(x, y) que corresponde a la curva de la familia tangente a la envolvente en
(𝑥, 𝑦). Sustituyendo ese valor variable de la constante en (1’) se obtiene la ecuación implícita de
la envolvente.
1.9 Ecuaciones no resueltas con respecto a la Primera Derivada
Si se quiere resolver una ecuación diferencial F(x, y, y′) = 0, puede intentarse despejar p de la
ecuación finita F(x, y, p) = 0 y resolver las ecuaciones y′ = p(x, y) con 𝑭(𝒙, 𝒚, 𝒑(𝒙, 𝒚)) = 𝟎
∀ (𝑥, 𝑦) que aparezcan, al tiempo que se tiene cuidado con las soluciones singulares que puedan
perderse o ganarse en el proceso.
Ejemplo: Resuelva:
(𝑦 ′)2 − (𝑥 + 𝑦)𝑦 ′ + 𝑥𝑦 = 0
Efectuando el cambio de variable: 𝑦 ′ = 𝑝, llegamos a una ecuación algebraica.
𝒑𝟐 − (𝑥 + 𝑦)𝒑 + 𝑥𝑦 = 0
Resolviéndola y retornando a la variable original: (𝒑 − 𝑥 )(𝒑 − 𝑦) = 0
′
𝒙𝟐
𝒑 − 𝑥 = 0 → 𝒑 = 𝑥 → 𝒚 = 𝑥 → 𝑑𝑦 = 𝑥𝑑𝑥 → ∫ 𝑑𝑦 = ∫ 𝑥 𝑑𝑥 → 𝒚= +𝑪
𝟐
𝑑𝑦
𝒑 − 𝑦 = 0 → 𝒑 = 𝑦 → 𝒚′ = 𝑦 → ∫ = ∫ 𝑑𝑥 → ln 𝑦 = 𝑥 + 𝐶1 → 𝒚 = 𝑪𝒆𝒙
𝑦
Ambas soluciones verifican la ecuación diferencial.
1.9.1 Ecuación de Clairaut
La ecuación de Clairaut tiene la forma: 𝒚 = 𝒙𝒚′ + 𝒇(𝒚′ ) … (1.34)
Para resolverla se emplea el cambio: 𝑦 ′ = 𝑝 … (𝛼)
𝑦 = 𝑥𝑝 + 𝑓(𝑝) … (𝛽)
Derivando respecto de "𝑥" obtenemos:
𝑦⏟′ = 𝑝 + 𝑥𝑝′ + 𝑓 ′ (𝑝)𝑝′ → 𝑝 = 𝑝 + 𝑥𝑝′ + 𝑓 ′ (𝑝)𝑝′ → 0 = [𝑥 + 𝑓 ′ (𝑝)]𝑝′ = 0

𝑝
𝑑𝑝
⏟+ 𝑓 ′(𝑝) = 0 ∨ ⏟
Se obtienen las dos posibilidades: 𝑥
𝑑𝑥
=0
1° 2°
Trabajando inicialmente con 2°:
𝑑𝑝
= 0 → 𝑑𝑝 = 0 → ∫ 𝑑𝑝 = ∫ 0 → 𝑝 = 𝐶 … (𝛾)
𝑑𝑥
Reemplazando (𝛾) en (𝛽) hallamos la solución general, que es una familia de rectas:
𝒚 = 𝐶𝒙 + 𝑓 (𝐶 ) … (1.35)
Trabajando con 1°: Se tendrá directamente una solución paramétrica
𝑥 + 𝑓 ′ (𝑝) = 0 → 𝑥 = −𝑓 ′ (𝑝)
Que con (9999) forman la solución paramétrica de la ecuación diferencial:
𝑥 = −𝑓 ′ (𝑝)
{ … (1.36)
𝑦 = 𝑥𝑝 + 𝑓(𝑝)
Eliminando de (1.36) precisamente la solución singular y es precisamente la ecuación de la

envolvente de la solución general (1.36)
𝑦 = −𝑝𝑓 ′ (𝑝) + 𝑓 (𝑝) … (1.37)

1.9.2 Ecuación de Lagrange
La ecuación de Lagrange tiene la forma: 𝒚 = 𝒙𝒇(𝒚′) + 𝒈(𝒚′ ) … (1.38)
Para resolverla se emplea el cambio: 𝑦 ′ = 𝑝 … (𝛼)
𝑦 = 𝑥𝑓(𝑝) + 𝑔(𝑝) … (𝛽)
Derivando respecto de "𝑥" obtenemos una ecuación lineal para "𝑥" respecto de "𝑝"
𝑦 ′ = 𝑓(𝑝) + 𝑥𝑓 ′ (𝑝)𝑝′ + 𝑔′ (𝑝)𝑝′

𝑑𝑝
𝑝 − 𝑓 (𝑝) = [𝑥𝑓 ′ (𝑝) + 𝑔′ (𝑝)]𝑝′ → 𝑝 − 𝑓(𝑝) = [𝑥𝑓 ′ (𝑝) + 𝑔′ (𝑝)]
𝑑𝑥
𝒅𝒙 𝑓 ′ (𝑝) 𝑔′ (𝑝)
+ 𝒙=
𝒅𝒑 𝑓(𝑝) − 𝑝 𝑝 − 𝑓(𝑝)
Resolviendo obtendremos: 𝑥 = Φ(𝑝, 𝐶) sustituyendo en (𝛽) obtenemos finalmente la solución

paramétrica, donde el parámetro viene a ser "𝑝"
𝒙 = Φ(𝑝, 𝐶 )
{ … (1.39)
𝒚 = Φ(𝑝, 𝐶 )𝑓(𝑝) + 𝑔(𝑝)
En el caso de Lagrange puede ocurrir que existan soluciones singulares, sin embargo, no es una
regla como el caso de Clairaut.
Alexis Claude Clairaut (1713-1765), nacido en Paris Francia, entre los muchos problemas
que estudió destaca el de los tres cuerpos. En su estudio del movimiento de la Luna, usó métodos
de cálculo de soluciones singulares. Calculó la vuelta del cometa Halley en 1759 y participó en
la expedición dirigida por Maupertuis para comprobar la predicción teórica de Newton de que la
Tierra es un esferoide oblato.
Joseph Louis Lagrange (1736-1813). Destacó en estudios sobre análisis, teoría de números,
mecánica celeste y mecánica analítica, de la que puede ser considerado como principal fundador.
Recogió todo lo que sobre mecánica se había hecho desde Newton en su obra maestra: Mécanique
Analytique (1788). Se le recuerda por sus ecuaciones del movimiento, las coordenadas
generalizadas, el concepto de energía potencial y los multiplicadores que llevan su nombre.
′
𝑎𝑦 ′
𝑦 = 𝑥𝑦 +
√1 + (𝑦 ′ )2
Solución: Podemos apreciar que se acomoda al caso de Clairaut
Efectuamos el cambio de variable: 𝑦 ′ = 𝑝 posteriormente derivando la ecuación respecto de "𝑥".

𝑎𝑝
𝑦 = 𝑥𝑝 + … (𝛼)
√1 + (𝑝)2
2𝑝𝑝′
𝑝′ √1 + 𝑝2 − 𝑝
2√1 + 𝑝2 1 + 𝑝2 − 𝑝2 ′
𝑦 ′ = 𝑝 + 𝑥𝑝′ + 𝑎 → 𝑝 = 𝑝 + 𝑥𝑝′ + 𝑎 [ 3 ]𝑝 →
1 + 𝑝2 (1 + 𝑝2 )2
[ ]
𝑎 ′
𝑎
0 = [𝑥 + 3] 𝑝 𝑥+ ⏟′ = 0
3 =0 ∨ 𝑝
(1 + 𝑝2 )2 2
⏟ (1 + 𝑝 )2 2°
1°
𝑑𝑝
Trabajando con (2°) = → ∫ 𝑑𝑝 = 0 → 𝑝 = 𝐶 en (𝛼) obtenemos la solución general
𝑑𝑥
𝑎𝐶
𝑺𝒐𝒍𝒖𝒄𝒊ó𝒏 𝑮𝒆𝒏𝒆𝒓𝒂𝒍: 𝒚 = 𝐶𝒙 + … (𝛽)
√1 + 𝐶 2
Trabajando con (1°) Hallamos la solución paramétrica conjuntamente con (𝛼)
𝑎 𝑎
𝑥=− 3
𝑥=− 3 … (1′)
(1 +𝑝2 )2 (1 + 𝑝2 )2
𝑎
𝑎𝑝 →
𝑎𝑝3
𝑦=− 3 𝑝 + 𝑦= 3 … (2′)
{ (1 + 𝑝2 )2 √1 + 𝑝2 { (1+𝑝 22 )
Tratando de eliminar el parámetro "𝑝" para hallar al envolvente.
Dividiendo (2′) entre (1′)
𝑦 3 𝑦 𝑎 𝑎
𝑝3 = − → 𝑝 = −√ 𝑒𝑛 (1′ ) → 𝑥=− 3 =− 3 →
𝑥 𝑥 2 2 2 2 2
3 𝑦 𝑥3 + 𝑦3
[1 + (√𝑥 ) ] ( 2 )
𝑥3
3
𝑎𝑥 2 2 2 2
𝑥=− 3 → (𝑥 3 + 𝑦3) = −𝑎 / ( )3
2 2 2
(𝑥 3 + 𝑦3)
Se obtiene la envolvente de la solución general, que en este caso es la Astroide:
𝟐 𝟐 𝟐
𝒙𝟑 + 𝒚𝟑 = 𝒂𝟑
aC
y  Cx 
1 C2
x
2 2 2
x3  y3  a3
𝑦 = 2𝑥𝑦 ′ + cos(𝑦 ′ )
Solución: La ecuación diferencial es una ecuación de diferencial de Lagrange:
Efectuamos el cambio de variable: 𝑦 ′ = 𝑝 posteriormente derivando la ecuación respecto de "𝑥".

𝑑
𝑦 = 2𝑥𝑝 + cos(𝑝) … (𝛼 ) / → 𝑦 ′ = 2𝑝 + 2𝑥𝑝′ − sin 𝑝 𝑝′ → 𝑝 = 2𝑝 + 2𝑥𝑝′ − sin 𝑝 𝑝′
𝑑𝑥
𝑑𝑝 𝑑𝑥 𝑑𝑥 2 sin 𝑝
−𝑝 = (2𝑥 − sin 𝑝)𝑝′ → −𝑝 = (2𝑥 − sin 𝑝) 𝑑𝑥 → −𝑝 𝑑𝑝 = 2𝑥 − sin 𝑝 → +𝑝𝑥 =
𝑑𝑝 𝑝
Resolviendo la ecuación lineal:

2
∫𝑝𝑑𝑝 sin 𝑝
2 sin 𝑝
− ∫𝑝𝑑𝑝 2
𝑥=𝑒 [∫ 𝑒 ( ) 𝑑𝑝 + 𝐶] → 𝑥 = 𝑒 −2 ln 𝑝 [∫ 𝑒 ln 𝑝 𝑑𝑝 + 𝐶]
𝑝 𝑝
1 1 sin 𝑝 cos 𝑝 𝐶
𝑥= (∫ 𝑝 sin 𝑝 𝑑𝑝 + 𝐶) = (−𝑝 cos 𝑝 + sin 𝑝 + 𝐶 ) → 𝑥= − + 2
𝑝2 𝑝2 𝑝2 𝑝 𝑝
Reemplazando en (𝛼)
sin 𝑝 2𝐶
𝑦 = 2𝑥𝑝 + cos 𝑝 → 𝑦 = 2 − 2 cos 𝑝 + + cos 𝑝
𝑝 𝑝
La solución es paramétrica:
sin 𝑝 cos 𝑝 𝐶
𝑥= − + 2
𝑝2 𝑝 𝑝
sin 𝑝 2𝐶
𝑦=2 − cos 𝑝 +
{ 𝑝 𝑝
Bibliografía
1. Rey Pastor Julio Análisis Matemático Volumen III

2. Aguirregabiria Juan Ecuaciones diferenciales para estudiantes de Física
3. Collete Jean Paul Historia de las Matemáticas Vol. II
4. Boyer Carl. Historia de la Matemática
5. Kreider - Kuller Ecuaciones Diferenciales
6. Zill Dennis Ecuaciones Diferenciales
7. Espinoza Eduardo Ecuaciones Diferenciales
8. Cruz Llanos Hernán (†) Apuntes de Clases de Ecuaciones Diferenciales
9. Espinoza Ríos Guillermo Apuntes de Clases de Ecuaciones Diferenciales
10. Lliulli Edgar Problemas Resueltos de Ecuaciones Diferenciales
11. Plantel Docente de Ecc. dif Apuntes de Clases de Mat-207
Facultad de Ingeniería UMSA

Cap 1 Ecuaciones Diferenciales 1er Orden

Cargado por

Copyright:

Formatos disponibles

Cap 1 Ecuaciones Diferenciales 1er Orden

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Cap 1 Ecuaciones Diferenciales 1er Orden

Cargado por

Copyright:

Formatos disponibles

Capítulo 1 Ecuaciones Diferenciales de Primer Orden M.Sc. Ing.

Juan Carlos Quispe Apaza

ECUACIONES DIFERENCIALES DE PRIMER ORDEN

1.1.1 La ecuación diferencial ordinaria

𝐹(𝑥, 𝑦, 𝑦 ′ , 𝑦 ′′ , 𝑦 ′′′ , … … , 𝑦 (𝑛) ) = 0 … (1.1)

En una ecuación diferencial la incógnita a calcularse viene a ser la función 𝑦 = 𝑓 (𝑥 ) ?

Ejemplo 1: Dé una solución de la ecuación diferencial

Tanteando al solución de (𝛼)

Por tanto 𝑦 = 𝑒 𝑥 es solución de (𝛼)

También verifica 𝑦 = 0 porque 𝑦 ′ = 0, entonces 𝑦 ′ − 𝑦 = 0 → 0 − 0 = 0 cumple.

Prueba: 𝑦 ′ = 𝑒 𝑥 reemplazando en la ecuación diferencial (𝛼)

Ejemplo 2: Dé una solución de la ecuación diferencial

Prueba: 𝑦 ′ = − sin 𝑥 , 𝑦 ′′ = − cos 𝑥 𝑒𝑛 (𝛽)

𝑦 ′′ + 𝑦 = 0 → − cos 𝑥 + cos 𝑥 = 0 → 0=0

Sin embargo, también existe otra solución:

𝑦 = sin 𝑥 → 𝑦′ = cos 𝑥 𝑦 ′′ = − sin 𝑥 𝑒𝑛 (𝛽)

En realidad, en el capítulo 3 se verá que si dos soluciones linealmente independientes son

𝑦 = 𝛼 sin 𝑥 + 𝛽 cos 𝑥 → 𝑦 ′ = 𝛼 cos 𝑥 − 𝛽 sin 𝑥 ; 𝑦 ′′ = −𝛼 sin 𝑥 − 𝛽 cos 𝑥 en (𝛽)

𝑦 ′′ + 𝑦 = 0 → −𝛼 sin 𝑥 − 𝛽 cos 𝑥 + 𝛼 sin 𝑥 + 𝛽 cos 𝑥 = 0

1.1.2 Ecuación diferencial en derivadas parciales

1.2 Clasificación según el orden y el grado de una ecuación diferencial

El orden de una ecuación diferencial ordinaria es la derivada de mayor orden en la ecuación

El grado de una ecuación diferencial es el grado algebraico de la derivada de mayor orden en la

a) 𝑦 ′′′ + 3𝑥(𝑦 ′ )6 = 𝑒 𝑥 De orden 3 y grado 1

Problema 1 Determine el orden y el grado de la ecuación diferencial:

(|𝑥 − 2| + |𝑥 + 2|)𝑦 ′′′ + (|𝑥 − 2| − |𝑥 + 2|)𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥

x  2 x 20 x  2 x 2

x  2 x 20 x  2 x  2

i) Para: −∞ < 𝒙 < −𝟐

(|𝑥 − 2| + |𝑥 + 2|)𝑦 ′′′ + (|𝑥 − 2| − |𝑥 + 2|)𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥

−2𝑥𝑦 ′′′ + 4𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥 De tercer orden y grado 1

ii) Para: −𝟐 ≤ 𝒙 < 𝟐

(|𝑥 − 2| + |𝑥 + 2|)𝑦 ′′′ + (|𝑥 − 2| − |𝑥 + 2|)𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥

4𝑦 ′′′ − 2𝑥𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥 De tercer orden y grado 1

(𝑥 − 2 + 𝑥 + 2)𝑦 ′′′ + (𝑥 − 2 − (𝑥 + 2))𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥

2𝑥𝑦 ′′′ − 4𝑦 ′′ + 3𝑥𝑒 𝑥 𝑦 = tan 𝑥 De tercer orden y grado 1

1.3 Solución de una ecuación diferencial

La función 𝑦(𝑥) es solución de la ecuación diferencial

𝐹(𝑥, 𝑦, 𝑦 ′ , 𝑦 ′′ , 𝑦 ′′′ , … … , 𝑦 (𝑛) ) = 0 … (1.2)

Ejemplo: Verifique si la función 𝑦 = 𝐶1 𝑒 2𝑥 + 𝐶2 𝑥𝑒 2𝑥 es solución de la ecuación diferencial:

𝑦 = 𝐶1 𝑒 2𝑥 + 𝐶2 𝑥𝑒 2𝑥 → 𝑦 ′ = 2𝐶1 𝑒 2𝑥 + 𝐶2 (1 ⋅ 𝑒 2𝑥 + 𝑥2𝑒 2𝑥 ) = (2𝐶1 + 𝐶2 )𝑒 2𝑥 + 2𝐶2 𝑥𝑒 2𝑥

Reemplazando en la ecuación diferencial

(4𝐶1 + 4𝐶2 )𝑒 2𝑥 + 4𝐶2 𝑥𝑒 2𝑥 − 4[(2𝐶1 + 𝐶2 )𝑒 2𝑥 + 2𝐶2 𝑥𝑒 2𝑥 ] + 4(𝐶1 𝑒 2𝑥 + 𝐶2 𝑥𝑒 2𝑥 ) = 0

1.3.1 Regla de Leibniz

Derivación bajo el signo integral

𝐹′ (𝑥 ) = 𝑔[𝑣(𝑥)] ⋅ 𝑣 ′ (𝑥 ) − 𝑔[𝑢(𝑥 )] ⋅ 𝑢′ (𝑥 ) … (1.3)

Problema 2 Compruebe si la función:

Es solución de la ecuación diferencial: 2𝑥𝑦 ′ − 𝑦 = 2𝑥 cos 𝑥

Solución: Calculando la derivada de la función y reemplazando la misma en la ecuación

Ecuación Diferencial Solución

Problema 3 Halle la ecuación diferencial cuya solución es la función:

Problema 4 Halle la ecuación diferencial que tiene por solución a la función:

Derivando nuevamente (1) respecto de "𝑥"

12𝑥 2 𝑦 + 4𝑥 3 𝑦 ′ + 4𝑥 3 𝑦 ′ + 𝑥 4 𝑦 ′′ = 56𝐴𝑥 6 + 216𝑥 7 →

Trabajando con (1) y (2)

4y xy' 8Ax 4 27x 5 / 7 28 y 7x y ' 56Ax 4 189x 5

𝑹𝒑𝒕𝒂.: 𝒙𝟐 𝒚′′ + 𝒙𝒚′ − 𝟏𝟔𝒚 = 𝟐𝟕𝒙𝟓

1.4 Existencia y Unicidad

Sea la ecuación diferencial 𝑦 ′ = 𝑓(𝑥, 𝑦) sujeta a la condición de valor inicial 𝑦(𝑥0 ) = 𝑦0