Tema 2 (Ecuaciones Diferenciales Ordinarias de Primer Orden)
Tema 2 (Ecuaciones Diferenciales Ordinarias de Primer Orden)
Tema 2 (Ecuaciones Diferenciales Ordinarias de Primer Orden)
Ecuaciones Diferenciales
Ordinarias de Primer Orden
ESQUEMA
- 39 -
40 Ecuaciones Diferenciales
EJERCICIOS
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 41
y′ = 4 y − x
Definición (Forma general o implícita de una EDO de primer orden). Dada una ecuación
diferencial ordinaria de primer orden de incógnita y (x ) , se llama forma general o implícita
de la ecuación a la forma de expresión en la que la variable independiente x, la función y,
y su derivada y' aparecen en un solo miembro de la ecuación, siendo el otro miembro
igual a cero:
F (x, y, y ′) = 0
o en notación diferencial:
dy
F x, y, = 0
dx
y′ + y + 3x2 = 0
42 Ecuaciones Diferenciales
Definición (Forma normal o explícita de una EDO de primer orden). Dada una ecuación
diferencial ordinaria de primer orden de incógnita y (x) , se llama forma normal o explícita
de la ecuación a la forma de expresión en la que la derivada aparece despejada en un
miembro de la ecuación:
y ′ = F (x, y )
o en notación diferencial:
dy
= F (x, y )
dx
donde F es una expresión analítica elemental.
y ′ = 2xy 2 + 2x 2 y + x + y − 2
Definición (Forma diferencial de una EDO de primer orden). Dada una ecuación
diferencial ordinaria de primer orden de incógnita y (x ) , se llama forma diferencial de la
dy
ecuación a la forma de expresión en la que la derivada y ′ = aparece en forma de
dx
suma de diferenciales:
G (x, y )dx + H ( x, y )dy = 0
(y + 3x )dx + 2ydy = 0
2
sen2 (x)
1− t
( )
dt + x 2 + ln(t ) dx = 0
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 43
Observación. Una EDO de primer orden en forma diferencial siempre puede reescribirse
en forma normal, ya que:
dy −G (x, y )
G (x, y )dx + H (x, y )dy = 0 ⇔ = = F (x, y )
dx H (x, y )
Asimismo, en el caso de que la EDO de primer orden sea lineal respecto a la derivada y',
es decir, cuando y' aparezca solamente en forma de potencia de primer grado, también
es posible pasar fácilmente la ecuación a forma normal, tal como se muestra en los
siguientes ejemplos. Precisamente, las EDO de primer orden lineales en y' constituyen
los tipos más básicos a considerar en el estudio de las ecuaciones diferenciales de
primer orden.
x ⋅ y ′ − xy = arcsen( x)
y ′ = F (x, y ) ; y (a) = b
y2
y′ =
x −1
Aplicando el teorema básico de existencia y unicidad, es posible
determinar el dominio D del plano en cada uno de cuyos puntos (x, y ) ∈ D
existe y es única la solución de la EDO. En efecto, se considera la función
y2
F (x, y ) = , con dom(F ) = {(x, y ) ∈ o 2 | x ≠ 1} , así como la derivada
x −1
∂F (x, y ) 2y
parcial = . Se tiene que:
∂y x −1
∂F
• es una función elemental, continua en todo el dominio de
∂y
derivabilidad parcial dom(F y′ ) = {(x, y ) ∈ o 2 | x ≠ 1} .
3
∂F (0,0) F (0,0 + h) − F (0,0) h2 − 0
= lim = lim =
∂y h→ 0 h h→ 0 h
1
3 2 lim+ 3
= 1+ = +∞
h 1 h 0 ∂F (0,0)
= lim = lim 3 = 0 = h→0
1 ⇒ ∃/
h→ 0 h h→0 h 1 ∂y
lim = 1− = −∞
h→0 −
3
h 0
Se tiene que:
• F es una función elemental, continua en todo su dominio
dom(F ) = o 2 .
∂F
• es una función elemental, continua en todo el dominio de
∂y
derivabilidad parcial dom(F y′ ) = o 2 − {(0,0)} .
y ′ = 1 − cos(xy ) ; y (0) = 0
xsen( xy )
∂F (x, y ) si (x, y ) ≠ (0,0)
Por tanto, = 2 1 − cos(xy ) .
∂y
0 si (x, y ) = (0,0)
Continuidad de F en (0,0):
• F (0,0) = 1 − cos(0) = 0
∴ F continua en (0,0)
46 Ecuaciones Diferenciales
∂F
Continuidad de en (0,0):
∂y
∂F (0,0)
• =0
∂y
∂F
∴ continua en (0,0)
∂y
Se trata de uno de los tipos más sencillos de EDO de primer orden, ya que tras una
sencilla manipulación algebraica se convierten en ecuaciones de integración directa.
1
donde G1 (x ) = 3x 2 , G2 ( y ) = y , H1 (x ) = , H2 ( y ) = 2 y .
cos( x )
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 47
x x ′ = et + t + 2
puede ponerse en forma de ecuación de variables separadas:
dx
x = et + t + 2
dt
( )
xdx = et + t + 2 dt
(− e − t − 2)dt + xdx = 0
t
donde G1 (t ) = −e t − t − 2 , G2 (x) = 1 , H1 (t ) = 1 , H2 (x ) = x.
Demostración:
dy
En efecto, llamando y ′ = , resulta:
dx
dy 1 1
y ′ = G (x )H ( y ) ⇔ = G (x)H ( y ) ⇔ dy = G ( x)dx ⇔ −G (x )dx + dy = 0
dx H( y) H( y)
donde F (x) = x + 3 , G ( y ) = tg ( y ) .
H2 ( y ) G (x)
dy = − 1 dx
G2 ( y ) H1 (x)
48 Ecuaciones Diferenciales
−1
y 2dy = dx
5x
2 −1
∫y dy = ∫ 5x dx
1
3
y3 = −1 ln
5
(x )+ C
y 3 = −53 ln( x ) + C
y= 3 −3 ln
5
(x )+ C
Ejemplo: Integrar la EDO de primer orden de variables separadas:
(
− 2xdx + 1 + x 4 e y dy = 0 )
Siguiendo el procedimiento descrito, se tiene:
(1 + x )e dy = 2xdx
4 y
2x
e ydy = dx
1 + x4
y 2x
∫e dy = ∫ 1 + x4 dx
( )
e y = arctg x2 + C
y = ln(arctg(x ) + C ) 2
( )
− e2t dt + 1 + e2t xdx = 0 ; x (0) = 1
(1 + e )xdx = e
2t 2t
dt
e2t
xdx = dt
1 + e2t
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 49
e2t
∫ xdx = ∫ 1 + e2t dt
1 x2
2
= 2
(
1 ln 1 + e2t )+ C
(
x 2 = ln 1 + e 2t + 2C )
Sustituyendo ahora la condición x(0) = 1 en la solución general, se
obtiene una ecuación que permite determinar el valor de la constante de
integración C:
x( 0 ) = 1
(1)2 = ln(1 + 1) + 2C
1 = ln(2) + 2C
C= 1
2
− 12 ln(2)
( )
x2 = ln 1 + e2t + 1 − ln(2)
1
∫ 0'0004 − 0'0002q dq = ∫ dt
−1 ln
0'0002
( 0'0004 − 0'0002q ) = t + C
ln( 0'0004 − 0'0002q ) = (−0'0002)t − 0'0002C
q = 2 − 5000Ce ( −0'0002)t
2 − 5000C = 1
C= 1
5000
q = 2 − e ( −0'0002)t
Se trata de un tipo de EDO de primer orden que se puede reducir a una ecuación de
variables separadas mediante un cambio de variable específico. Previamente a su
definición es necesario recordar el concepto de función homogénea de varias variables.
f (x, y ) = x2 y + xy 2
( )
f (tx, ty ) = (tx)2 (ty ) + (tx)(ty )2 = t 3x2 y + t 3xy 2 = t 3 x2 y + xy 2 = t 3f (x, y )
xy + xz + yz
f (x, y, z ) =
x
t 2 (xy + xz + yz ) xy + xz + yz 1
= = t = t 2 f (x, y, z )
tx x
f (x, y ) = x2 + xy − y
t αf (x, y ) = t α x2 + xy − y = t 2α x2 + t 2α xy − t 2α y
(x 2
)
+ xy dx + 4 y 2dy = 0
( )
G (tx, ty ) = (tx ) 2 + (tx )(ty ) = t 2 x 2 + t 2 xy = t 2 x 2 + xy = t 2G (x, y )
1
G (tx, ty ) = 3 tx − ty = 3 t ( x − y ) = 3 t 3 x − y = t 3 G (x, y )
1
H (tx, ty ) = 3 tx = 3 t 3 x = t 3 H (x, y )
y
F (tx, ty ) = F (x, y ) = F
x
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 53
Demostración:
Considérese en primer lugar que la EDO y ′ = F (x, y ) es homogénea, es decir, puede expresarse en forma
diferencial:
donde G(x,y) y H(x,y) son funciones homogéneas de igual grado α ∈ o . Si se pasa esta ecuación a forma
normal, se tiene que:
−G (x, y )
y′ = = F ( x, y )
H (x, y )
Por tanto,
1
luego F es una función homogénea de grado cero. Además, haciendo t = en la igualdad anterior, se tiene:
x
x y y y
F (x, y ) = F (tx, ty ) = F , = F 1, = F
x x x x
y
Recíprocamente, si la EDO y′ = F (x, y ) verifica que F puede expresarse en la forma F , es decir, F es una
x
función homogénea de grado cero, entonces pasando la ecuación a forma diferencial, resulta:
dy
y ′ = F (x, y ) ⇔ = F (x, y ) ⇔ dy = F (x, y )dx ⇔ −F (x, y )dx + dy = 0
dx
donde G (x, y ) = −F (x, y ) y H (x, y ) = 1 son funciones homogéneas de igual grado α = 0 , luego la EDO y′ = F (x, y )
es homogénea.
3xy
ya que F (x, y ) = 2 es una función homogénea de grado cero:
x − y2
y
Además, F puede expresarse en la forma F :
x
3y 3xy y
3
3xy x x 2 x = F y
F (x, y ) =
2 2
=
2 2
=
2
=
2
x −y x −y y y x
1− 1−
2 2
x x x
54 Ecuaciones Diferenciales
dx t2
=
dt tx
t2
ya que F (t, x ) = es una función homogénea de grado cero:
tx
(rt ) 2 r 2t 2 r2 t2 t2
F (rt, rx) = = = = = F (t, x )
(rt )(rx) r 2tx r 2 tx tx
x
Además, H puede expresarse en la forma F :
t
t2
t2 t2 1 x
F (t, x ) = = = = F
tx tx x t
t2 t
y = tx
G (x, y )dx + H (x, y )dy = 0 ⇒ ⇒ G1 (t )G2 (x)dt + H1 (t )H2 (x)dx = 0
dy = xdt + tdx
Demostración:
y = tx
donde G(x,y) y H(x,y) son funciones homogéneas de igual grado α ∈ o . Si se aplica el cambio ,
dy = xdt + tdx
resulta:
H (t )xdt + (G (t ) + tH (t ))dx = 0
que es una ecuación de variables separadas de la forma G1(t )G2 (x )dt + H1(t )H2 (x )dx = 0 con G1(t ) = H (t ) ,
G2 (x ) = x , H1(t ) = G (t ) + tH (t ) , H2 (x ) = 1 .
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 55
y
4. Se deshace el cambio en la solución sustituyendo t = .
x
(x 2
)
− 3 y 2 dx + 2xydy = 0
y = tx
En primer lugar, se aplica el cambio , obteniéndose:
dy = xdt + tdx
(x 2
)
− 3t 2x2 dx + 2tx2 (xdt + tdx ) = 0
(
2tx 3dt + 1 − t 2 x2dx = 0 )
A continuación, se integra la ecuación de variables separadas obtenida:
(1 − t )x dx = −2tx dt
2 2 3
1 −2t
dx = dt
x 1 − t2
1 −2t
∫ x dx = ∫ 1 − t 2 dt
ln x = ln 1 − t 2 + C
( )
x = eC 1 − t 2
(
x = ±eC 1 − t 2 )
(
x = C 1 − t2 )
Finalmente, se deshace el cambio de variable:
56 Ecuaciones Diferenciales
y2
x = C1 − 2
x
y = tx
En primer lugar, se aplica el cambio , obteniéndose:
dy = xdt + tdx
y y
dy = + cos2 dx
x x
tx tx
xdt + tdx = + cos2 dx
x x
(
xdt + tdx = t + cos2 (t ) dx )
cos2 (t )dx = xdt
A continuación, se integra la ecuación de variables separadas obtenida:
1 1
dx = dt
x cos2 (t )
1 1
∫ x dx = ∫ cos2 (t ) dt
ln( x ) = tg(t ) + C
x = etg(t ) +C
x = ±eCetg(t )
x = Cetg(t )
Finalmente, se deshace el cambio de variable:
y
tg
x= Ce x
xdx = x2dt
A continuación, se integra la ecuación de variables separadas obtenida:
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 57
xdx = x2dt
1
dx = dt
x
1
∫ x dx = ∫ dt
( )
ln x = t + C
1 ⋅ ln(1 ) − C = 1
C = −1
Por tanto, la solución del problema de valor inicial es:
y = x ln( x ) + x
ax + by + c
y ′ = F
a′x + b′y + c′
donde F es una expresión analítica elemental.
ax + by + c
Observación. Nótese que la ecuación y ′ = F no es homogénea, porque la
a′x + b′y + c′
función F no es homogénea de grado cero, debido a la presencia de los términos c y c'
(siempre que al menos uno de ellos no sea nulo). Además, en una ecuación reducible a
homogénea las expresiones ax + by + c , a′x + b′y + c′ son expresiones afines de dos
variables x e y que corresponden a sendas rectas del plano OXY. Precisamente, la
dependencia o independencia lineal entre los vectores directores de estas dos rectas
permite definir los cambios de variable apropiados para su resolución.
ax + by + c
y ′ = F
a′x + b′y + c′
se verifica que:
• Si los vectores directores de las rectas ax + by + c = 0 , a ′x + b ′y + c ′ = 0 son
libres, es decir, si ab′ ≠ a′b , las rectas se cortan en un punto (x0,y0), por lo que
u = x − x0
mediante el cambio de variable la EDO se transforma en una
v = y − y0
ecuación homogénea:
dv au + bv
= F
du a ′u + b ′v
• Si los vectores directores de las rectas ax + by + c = 0 , a ′x + b ′y + c ′ = 0 son
ligados, es decir, si ab′ = a′b , las rectas son paralelas, por lo que mediante el
cambio de variable y = −ba x + u la EDO se transforma en una ecuación de
variables separadas:
du a bu + c
= b
+ F
dx b ′u + c ′
Demostración:
ax + by + c
y ′ = F
a′x + b′y + c′
y considérese en primer lugar el caso ab′ ≠ a′b , tal que las rectas ax + by + c = 0 , a ′x + b ′y + c ′ = 0 se cortan en
u = x − x0
un punto (x0,y0). Si se aplica el cambio , se produce una traslación del origen de coordenadas al
v = y − y0
punto (x0,y0) que elimina los términos c y c' y transforma la EDO en una ecuación homogénea:
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 59
dy ax + by + c u = x − x0 ; du = dx dv a(u + x0 ) + b(v + y0 ) + c
= F ⇒ ⇒ = F ⇒
dx a′x + b′y + c′ v = y − y 0 ; dv = dy du a′(u + x0 ) + b′(v + y0 ) + c′
dv au + bv + ax0 + ay0 + c dv au + bv
⇒ = F ⇒
= F
du a′u + b′v + a′x0 + b′y0 + c′ du a′u + b′v
ya que, por ser (x0,y0) punto de corte de ambas rectas ax + by + c = 0 , a ′x + b ′y + c ′ = 0 , se cumple que
ax0 + by0 + c = 0 , a′x0 + b′y0 + c′ = 0 .
Por otra parte, si ab′ = a′b , donde k = a′ = b' , se aplica el cambio y = −a x + u , con lo que se obtiene:
a b b
ax + by + c ax + by + c y = −a x + u
dy dy b
= F ⇒ = F ⇒ ⇒
′ ′ ′ −a
dx a x + b y + c dx kax + kby + c′ dy = b dx + du
⇒ b
−a dx + du
= F
( b
)
ax + b −a x + u + c
⇒ −a + du = F bu + c ⇒ du = a + F bu + c
dx
b
( )
kax + kb −a x + u + c′
b dx kbu + c′ dx b b′u + c′
dy ax + by + c y − y 0 = t (x − x 0 ) (x − x 0 )dt + tdx a + bt
= F ⇒ ⇒ = F
dx a ′x + b ′y + c ′ dy = ( x − x 0 )dt + tdx dx a ′ + b ′t
y − y0
4. Se deshace el cambio en la solución sustituyendo t = .
x − x0
x − y + 1 = 0
⇒ 2x + 2 = 0 ⇒ x = −1, y = 0
x + y + 1 = 0
y = t (x + 1)
Se aplica el cambio , obteniéndose:
dy = ( x + 1)dt + tdx
x − y +1
dy = dx
x + y +1
x − t (x + 1) + 1
(x + 1)dt + tdx = dx
x + t (x + 1) + 1
(x + 1) − t (x + 1)
(x + 1)dt + tdx = dx
(x + 1) + t (x + 1)
1− t
(x + 1)dt + tdx = dx
1+ t
1 − 2t − t 2
(x + 1)dt = dx
1+ t
A continuación, se integra la ecuación de variables separadas obtenida:
1 1+ t
dx = dt
x +1 1 − 2t − t 2
1 −2 − 2t
∫ x + 1 dx = −1
2 ∫ 1 − 2t − t 2 dt
ln x + 1 = −21 ln 1 − 2t − t 2 + C
( )
eC
x +1 =
1 − 2t − t 2
± eC
x +1 =
1 − 2t − t 2
C
x +1 =
1 − 2t − t 2
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 61
dy ax + by + c y = −a x + u du bu + c
= F ⇒ b ⇒ a + F
− a dx + du
= b
dx a ′x + b ′y + c ′ dy = b dx b ′u + c ′
x − x + u +1
− dx + du = dx
x − x + u −1
u +1
− dx + du = dx
u −1
62 Ecuaciones Diferenciales
2u
du = dx
u −1
A continuación, se integra la ecuación de variables separadas obtenida:
u −1
du = 2dx
u
1
∫ 1 − u du = ∫ 2dx
( )
u − ln u = 2x + C
x + y = e −C e y − x
x + y = ±e −C e y − x
x + y = Ce y − x
Las ecuaciones diferenciales exactas son EDO de primer orden cuya expresión en forma
diferencial coincide con la diferencial (total) de una función de dos variables. Su
resolución puede llevarse a cabo mediante un proceso de integración parcial respecto a
la primera variable, seguido de una derivación parcial respecto a la segunda variable, y
de una nueva integración simple respecto a la segunda variable.
es decir,
∂F (x, y ) ∂F (x, y )
= G ( x, y ) ; = H (x, y )
∂x ∂y
(3x 2
) ( )
+ 4xy 2 dx + 4x 2 y − 3 y 2 dy = 0
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 63
∂F (x, y ) ∂F (x, y )
= 3x 2 + 4xy 2 ; = 4x 2 y − 3 y 2
∂x ∂y
es decir,
( ) (
dF ( x, y ) = 3x 2 + 4xy 2 dx + 4x 2 y − 3 y 2 dy )
Ejemplo: La siguiente EDO de primer orden de incógnita x (t ) es una ecuación
diferencial exacta:
2t 8x 3
dt + dx = 0
1 + t + 2x 4
2
1 + t 2 + 2x 4
(
ya que considerando la función F (t, x ) = ln 1 + t 2 + 2x 4 se cumple que: )
∂F (t, x ) 2t ∂F (t, x) 8x3
= ; =
∂t 1 + t 2 + 2x 4 ∂x 1 + t 2 + 2x 4
es decir,
2t 8x 3
dF (t, x) = dt + dx
1 + t + 2x 4
2
1 + t 2 + 2x 4
∂G (x, y ) ∂H (x, y )
=
∂y ∂x
Demostración:
1. Si la ecuación G (x, y )dx + H(x, y )dy = 0 es una ecuación diferencial exacta en D ⊆ o2 , entonces para todo
∂G (x, y ) ∂H (x, y )
(x, y ) ∈ D se verifica = .
∂y ∂x
En efecto, por definición de ecuación diferencial exacta existe una función F : D ⊆ o2 → o , tal que para todo
(x, y ) ∈ D verifica:
∂F ( x, y ) ∂F (x, y )
= G ( x, y ) ; = H (x, y )
∂x ∂y
64 Ecuaciones Diferenciales
y de aquí:
Ahora bien, como las derivadas parciales de G y H son continuas en (x, y ) ∈ D y coinciden con las derivadas
parciales segundas de F, podemos aplicar el teorema de Schwarz a la función F, de donde resulta para todo
(x, y ) ∈ D :
∂ 2F (x, y ) ∂ 2F (x, y )
=
∂x∂y ∂y∂x
y en consecuencia:
∂G (x, y ) ∂H (x, y )
=
∂y ∂x
∂G (x, y ) ∂H (x, y )
2. Si para todo (x, y ) ∈ D se verifica = , entonces la ecuación G (x, y )dx + H (x, y )dy = 0 es una
∂y ∂x
Se trata de probar que existe una función F : D ⊆ o2 → o , tal que para todo (x, y ) ∈ D verifica
dF (x, y ) = G (x, y )dx + H (x, y )dy . Dado que las funciones G y H son continuas en (x, y ) ∈ D , ambas son
parcialmente continuas en (x, y ) ∈ D respecto a cada variable, luego existe la función F : D ⊆ o2 → o integral
parcial de G respecto a la variable x:
siendo ϕ una función arbitraria que depende sólo de la variable y, equivalente a la constante de integración de
las integrales simples, pues su derivada respecto a x es nula. Derivando respecto a y la función F obtenida,
resulta:
∂F (x, y )
∂y
=
∂
∂y
(∫ G(x, y )dx ) + ddyϕ
Por tanto, para que la función F verifique dF (x, y ) = G (x, y )dx + H (x, y )dy y la EDO sea exacta, debe cumplirse
para (x, y ) ∈ D que:
∂
∂y
(∫ G(x, y)dx) + ddyϕ = H(x, y)
es decir,
dϕ
dy
= H ( x, y ) −
∂
∂y
(∫ G(x, y)dx )
dϕ
y como ϕ es una función arbitraria que depende sólo de la variable y, su derivada también deberá ser
dy
2
∂
H (x, y ) −
∂x
∂
∂y
(∫ G(x, y)dx) = ∂H∂(xx, y ) − ∂∂y∂x (∫ G(x, y )dx) =
∂G (x, y ) ∂H (x, y )
de donde, por el teorema de Schwarz y la hipótesis = , se obtiene:
∂y ∂x
∂2
=
∂H (x, y )
∂x
−
∂x∂y
(∫ G(x, y)dx) = ∂H∂(xx, y ) − ∂∂y ∂∂x (∫ G(x, y )dx ) = ∂H∂(xx, y ) − ∂G∂(xy, y ) = 0
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 65
En conclusión, la función F verifica dF (x, y ) = G (x, y )dx + H (x, y )dy y la EDO es exacta.
∂G (x, y ) ∂H ( x, y )
= cos(x ) ; = cos(x)
∂y ∂x
(2x 2
)
− t dt + 4txdx = 0
∂G (t, x ) ∂H (t, x)
= 4x ; = 4x
∂x ∂t
∂F (x, y )
∂y
=
∂
∂y
(∫ G(x, y)dx ) + ddyϕ = H(x, y)
66 Ecuaciones Diferenciales
dϕ
3. Se despeja en la ecuación formada y se integra respecto a y. La función ϕ
dy
calculada se lleva a la expresión de F(x,y), y con ello se obtiene la solución
general de la EDO exacta.
(x 2
) (
+ 2xy dx + x 2 + 2y 2 dy = 0 )
Se trata efectivamente de una EDO exacta, ya que llamando
G (x, y ) = x 2 + 2xy , H (x, y ) = x 2 + 2y 2 , se cumple que:
∂G (x, y ) ∂H (x, y )
= 2x ; = 2x
∂y ∂x
∫ (x )
2
F (x, y ) = + 2xy dx = 13 x 3 + x2 y + ϕ( y )
∂F (x, y )
∂y
=
∂ 1 3
∂y 3
(
x + x2y +
dϕ
dy
= x2 +
dϕ
dy
)
dϕ
x2 + = x 2 + 2y 2
dy
dϕ
Por último, se despeja en la ecuación formada y se integra respecto a
dy
la variable y:
dϕ
= 2y 2
dy
ϕ′( y ) = 2y 2
2
∫ ϕ′( y)dy = ∫ 2y dy
ϕ( y ) = 2
3
y3 + C
(e x
)
+ 3 y dx + (3x + cos( y ))dy = 0
∂G (x, y ) ∂H (x, y )
= 3; =3
∂y ∂x
∂F (x, y ) ∂ dϕ dϕ
= (3xy + sen( y )) + = 3y +
∂x ∂x dx dx
dϕ
3y + = e x + 3y
dx
dϕ
Por último, se despeja en la ecuación formada y se integra respecto a
dx
la variable x:
dϕ
= ex
dx
ϕ′(x) = e x
x
∫ ϕ′(x)dx = ∫ e dx
ϕ(x) = e x + C
3xy + sen( y ) + e x + C = 0
∂H (x, y )
= cos(xy ) − 3xysen(xy ) − x2 y 2 cos(xy )
∂x
68 Ecuaciones Diferenciales
F (x, y ) = ∫ (y cos(xy ) − xy
2
)
sen(xy ) dx =
∴ F (x, y ) = xy cos(xy ) + ϕ( y )
∂F (x, y ) ∂ dϕ dϕ
= (xy cos(xy )) + = x cos(xy ) − x2 ysen(xy ) +
∂y ∂y dy dy
dϕ
x cos(xy ) − x2 ysen( xy ) + = x cos(xy ) − x2 ysen( xy )
dy
dϕ
Por último, se despeja en la ecuación formada y se integra respecto a
dy
la variable y:
dϕ
=0
dy
ϕ′( y ) = 0
∫ ϕ′( y )dy = ∫ 0 ⋅ dy
ϕ( y ) = C
()
Sustituyendo ahora la condición y 3π = 1 en la solución general, se obtiene
una ecuación que permite determinar el valor de la constante de
integración C:
()
y 3π = 1
C = −6π
xy cos(xy ) − π =0
6
∂G (x, y ) ∂H ( x, y )
≠
∂y ∂x
∂ ∂
(µ(x, y )G (x, y )) = (µ(x, y )H (x, y ))
∂y ∂x
esto es, EDO exactas que pueden resolverse por el método visto de integración-
derivación-integración. La utilidad de esta técnica está en que la solución general de la
ecuación original y de la transformada en exacta es esencialmente la misma, ya que al
multiplicar la EDO por el factor integrante sólo se añadirán algunas soluciones a su
solución general, concretamente el conjunto de soluciones particulares que hacen
µ( x, y ) = 0 (siempre que existan, y no coincidan con las de la EDO original), y
recíprocamente sólo se eliminarán algunas soluciones de la EDO original, aquellas para
las cuales el factor integrante no esté definido (si existen).
∂G (x, y ) ∂H ( x, y )
≠
∂y ∂x
∂ ∂
(µ(x, y )G (x, y )) = (µ(x, y )H (x, y ))
∂y ∂x
(2x 2
) (
y + y 2 dx + x 3 + 3
2
)
xy dy = 0
∂G ( x, y ) ∂H (x, y )
= 2x2 + 2 y ≠ = 3x2 + 32 y
∂y ∂x
( ) (
xy 2x 2 y + y 2 dx + xy x 3 + 3
2
)
xy dy = 0
(2x 3 2
) ( )
y + xy 3 dx + x 4 y + 32 x 2 y 2 dy = 0
tsen(tx ) 2
(sen(tx ) − t cos(x ))dt + + t sen(x ) dx = 0
x
tsen(tx )
ya que, llamando G (t, x ) = sen(tx ) − t cos( x) , H (t, x) = + t 2sen(x) ,
x
se cumple que:
∂G (t, x ) ∂H (t, x ) sen(tx )
= t cos(tx ) + tsen(x ) ≠ = + t cos(tx ) + 2tsen( x )
∂x ∂t x
1
La función µ(t, x ) = es un factor integrante de la ecuación diferencial, ya
t
que al multiplicarla por él, se obtiene la EDO exacta:
1 1 tsen(tx ) 2
(sen(tx ) − t cos(x ))dt + + t sen( x ) dx = 0
t t x
sen(tx ) sen(tx )
− cos( x) dt + + tsen(x) dx = 0
t x
puesto que se cumple:
∂ ∂
(µ(t, x )G (t, x)) = (µ(t, x)H (t, x)) = cos(tx) + sen(x)
∂x ∂t
Sea una EDO de primer orden no exacta G (x, y )dx + H (x, y )dy = 0 , es decir, una ecuación
diferencial tal que:
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 71
∂G (x, y ) ∂H ( x, y )
≠
∂y ∂x
∂G ( x, y ) ∂H ( x, y ) ∂µ(x, y ) ∂µ(x, y )
µ(x, y ) − µ( x, y ) = H ( x, y ) − G ( x, y )
∂y ∂x ∂x ∂y
∂G ( x, y ) ∂H ( x, y ) ∂µ(x, y ) ∂µ( x, y )
µ( x, y ) − = H ( x, y ) − G ( x, y )
∂ y ∂ x ∂x ∂y
que es la ecuación diferencial general de cálculo del factor integrante. Esta ecuación es
una ecuación en derivadas parciales, normalmente mucho más difícil de resolver que la
EDO original. De ahí que resulte inviable considerar el cálculo general de factores
integrantes para ecuaciones diferenciales ordinarias. Sin embargo, la EDD de cálculo del
factor integrante es resoluble en algunos casos particulares, por ejemplo, si se posee
alguna información adicional respecto al factor integrante µ, tal como que µ dependa
sólo de la variable x, o sólo de la variable y, o de una expresión concreta de dos variables
como la suma x+y o el producto xy. Precisamente, estos casos son los que se van a
considerar en el cálculo práctico de factores integrantes.
µ(x) = Ce ∫ f ( x )dx
Demostración:
En efecto, sea µ un factor integrante dependiente sólo de la variable x que convierte a la ecuación en una
ecuación diferencial exacta de la forma µ(x )G (x, y )dx + µ(x )H (x, y )dy = 0 . Por definición de ecuación diferencial
exacta se debe cumplir:
∂
(µ(x)G (x, y )) = ∂ (µ(x)H(x, y ))
∂y ∂x
y de aquí:
72 Ecuaciones Diferenciales
∂G (x, y ) dµ ∂H (x, y )
µ( x ) = H (x, y ) + µ( x )
∂y dx ∂x
∂G (x, y ) ∂H ( x, y ) dµ
µ( x ) − µ( x ) = H (x, y )
∂y ∂x dx
∂G (x, y ) ∂H (x, y ) dµ
µ(x ) − = H (x, y )
∂y ∂x dx
∂G (x, y ) ∂H (x, y )
−
1 dµ ∂y ∂x
⋅ =
µ(x ) dx H (x, y )
1 dµ
y como ⋅ es una función que depende sólo de la variable x, de la igualdad anterior resulta la condición
µ( x ) dx
referida:
∂G (x, y ) ∂H (x, y )
−
∂y ∂x
= f (x )
H (x, y )
1 dµ
⋅ = f (x )
µ(x ) dx
µ′(x )
= f (x )
µ( x )
µ′(x )
∫ µ(x) dx = ∫ f (x)dx
( )
ln µ(x ) = ∫ f (x )dx + C
µ(x ) = Ce ∫ f ( x )dx
2 y
x
2
xy − dx + 2x y + 1 dy = 0 ( )
En primer lugar, se comprueba que la EDO admite un factor integrante de
y
la forma µ(x ) . Llamando G (x, y ) = xy 2 − , H (x, y ) = 2x2 y + 1 , se cumple
x
que:
∂G (x, y ) ∂H (x, y ) 1 1
− 2xy − − 4 xy − − 2xy
∂y ∂x x x − 1 − 2x2 y −1
= = = = = f (x )
H ( x, y ) 2
2x y + 1 2
2x y + 1 2
x 2x y + 1 x ( )
El factor integrante es:
−1
∫ dx −ln( x ) C 1
µ( x ) = Ce ∫ f ( x )dx = Ce x = Ce = = {C = 1} =
x x
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 73
µ( y ) = Ce ∫ f ( y )dy
Demostración:
En efecto, sea µ un factor integrante dependiente sólo de la variable y que convierte a la ecuación en una
ecuación diferencial exacta de la forma µ( y )G ( x, y )dx + µ( y )H (x, y )dy = 0 . Por definición de ecuación diferencial
exacta se debe cumplir:
∂
(µ( y )G (x, y )) = ∂ (µ( y )H(x, y ))
∂y ∂x
y de aquí:
dµ ∂G (x, y ) ∂H ( x, y )
G (x, y ) + µ( y ) = µ( y )
dy ∂y ∂x
dµ ∂H (x, y ) ∂G (x, y )
G (x, y ) = µ( y ) − µ( y )
dy ∂x ∂y
dµ ∂H (x, y ) ∂G (x, y )
G (x, y ) = µ( y ) −
dy ∂x ∂y
∂H (x, y ) ∂G (x, y )
−
1 dµ ∂x ∂y
⋅ =
µ( y ) dy G (x, y )
1 dµ
y como ⋅ es una función que depende sólo de la variable y, de la igualdad anterior resulta la condición
µ( y ) dy
referida:
∂H (x, y ) ∂G (x, y )
−
∂x ∂y
= f ( y)
G (x, y )
1 dµ
⋅ = f (y)
µ( y ) dy
µ′( y )
= f ( y)
µ( y )
µ′( y )
∫ µ( y ) dy = ∫ f ( y )dy
( )
ln µ( y ) = ∫ f ( y )dy + C
µ( y ) = ece ∫ f ( y )dy
74 Ecuaciones Diferenciales
µ( y ) = ±ece ∫ f ( y )dy
µ( y ) = Ce ∫ f ( y )dy
(xy 2
) (
− y 3 dx + 1 − xy 2 dy = 0 )
En primer lugar, se comprueba que la EDO admite un factor integrante de
la forma µ(y ) . Llamando G (x, y ) = xy 2 − y 3 , H (x, y ) = 1 − xy 2 , se cumple
que:
∂H (x, y ) ∂G (x, y )
−
∂x ∂y − y 2 − 2xy + 3 y 2 − 2xy + 2 y 2
= = =
G (x, y ) xy 2 − y 3 xy 2 − y 3
2y ( y − x) −2
= 2
= = f ( y)
− y ( y − x) y
µ(w ) = Ce ∫ f (w )dw
Demostración:
En efecto, sea µ un factor integrante dependiente de una expresión w(x,y) que convierte a la ecuación en una
ecuación diferencial exacta de la forma µ(w (x, y ))G (x, y )dx + µ(w (x, y ))H (x, y )dy = 0 . Por definición de ecuación
diferencial exacta se debe cumplir:
∂
(µ(w (x, y ))G (x, y )) = ∂ (µ(w (x, y ))H(x, y ))
∂y ∂x
y de aquí:
∂G (x, y ) ∂H (x, y )
−
1 dµ ∂y ∂x
⋅ =
µ(w (x, y )) dw H (x, y ) ∂w (x, y ) − G (x, y ) ∂w (x, y )
∂x ∂y
1 dµ
y como ⋅ es una función que depende de la expresión w(x,y), de la igualdad anterior resulta la
µ(w (x, y )) dw
condición referida:
∂G ( x, y ) ∂H (x, y )
−
∂y ∂x
= f (w (x, y ))
∂w (x, y ) ∂w (x, y )
H (x, y ) − G (x, y )
∂x ∂y
1 dµ
⋅ = f (w )
µ(w ) dw
µ′(w )
= f (w )
µ(w )
µ′(w )
∫ µ(w ) dw = ∫ f (w )dw
( )
ln µ(w ) = ∫ f (w )dw + C
µ(w ) = Ce ∫ f (w )dw
(y + xy )dx + (x − x y )dy = 0
2 2
4xy 4xy −2 −2
= 2 2 2 2
= 2 2
= = = f (w (x, y ))
xy − x y − xy − x y − 2x y xy w
−2
∫ dw −2ln(w ) C C 1
µ(w ) = Ce ∫ f (w )dw = Ce w = Ce = 2
= 2 2
= {C = 1} = 2 2
w x y x y
76 Ecuaciones Diferenciales
∂G (x, y ) ∂H (x, y )
1. Se calculan las derivadas parciales , y se estudia la existencia
∂y ∂x
del factor integrante considerando tres casos:
∂G (x, y ) ∂H (x, y )
−
∂y ∂x
• Si = f (x ) , entonces el factor integrante es de la forma
H (x, y )
µ(x) = Ce ∫ f ( x )dx .
∂H (x, y ) ∂G (x, y )
−
∂x ∂y
• Si = f ( y ) , entonces el factor integrante es de la forma
G (x, y )
µ( y ) = Ce ∫ f ( y )dy .
• Si es posible encontrar una expresión sencilla w(x,y), tal como x+y, x−y,
∂G ( x, y ) ∂H (x, y )
−
∂y ∂x
xy, x/y, para la cual se verifica = f (w (x, y )) ,
∂w (x, y ) ∂w (x, y )
H (x, y ) − G (x, y )
∂x ∂y
entonces el factor integrante es de la forma µ(w ) = Ce ∫ f (w )dw .
2. Se multiplica la EDO original G (x, y )dx + H (x, y )dy = 0 por el factor integrante
µ( x, y ) , obteniéndose la EDO exacta:
∂F (x, y )
•
∂y
=
∂
∂y
(∫ µ(x, y)G(x, y)dx ) + ddyϕ = µ(x, y)H(x, y)
• ϕ( y ) = ∫ µ(x, y )H (x, y ) −
∂
∂y
(∫ µ(x, y)G(x, y)dx )dy
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 77
∂G (x, y ) 4x ∂H (x, y )
= 1 + 2xy + ≠ = 3 + 4xy
∂y y3 ∂x
4x
1 + 2xy +
− 3 − 4xy
y3
= =
2x ∂w (x, y )
( ∂x
)
2 ∂w ( x, y )
3x + 2x y 2
− y + xy −
y 2 ∂y
y 3 + 2xy 4 + 4 x − 3 y 3 − 4 xy 4
= =
2x ∂w (x, y )
3
( 2 ∂w ( x, y )
y 3x + 2x y
∂x
)
3
− y y + xy −
2
y 2 ∂y
− 2 y 3 − 2xy 4 + 4 x
=
(3xy 3
+ 2x 2 y 4 ) ∂w (x, y )
∂x
(
− y 4 + xy 5 − 2xy
∂w (x, y )
∂y
)
∂w ∂w
• Caso factor integrante µ(x ) ( w = x , = 1, = 0 ):
∂x ∂y
− 2 y 3 − 2xy 4 + 4 x
=
(3xy 3
) (
+ 2x 2 y 4 (1) − y 4 + xy 5 − 2xy (0) )
− 2y 3 − 2xy 4 + 4x
= ≠ f (x)
3xy 3 + 2x 2 y 4
∂w ∂w
• Caso factor integrante µ(y ) ( w = y , = 0, = 1 ):
∂x ∂y
− 2 y 3 − 2xy 4 + 4 x
=
(3xy 3
) (
+ 2x 2 y 4 (0) − y 4 + xy 5 − 2xy (1) )
2y 3 + 2xy 4 − 4x 2
= = = f ( y)
(y 3
+ xy 4
− 2x y ) y
2
∫ y dy 2 ln( y )
µ( y ) = Ce ∫ f ( y )dy = Ce = Ce = Cy 2 = {C = 1} = y 2
(y 3
) (
+ xy 4 − 2x dx + 3xy 2 + 2x 2 y 3 dy = 0 )
Se integra µ( y )G (x, y ) = y 3 + xy 4 − 2x parcialmente respecto a x:
∫ (y )
3
F (x, y ) = + xy 4 − 2x dx = xy 3 + 1
2
x 2 y 4 − x 2 + ϕ( y )
∂F (x, y )
∂y
=
∂
∂y
xy 3 + ( 1
2
x2y 4 − x2 + ) dϕ
dy
= 3xy 2 + 2x 2 y 3 +
dϕ
dy
dϕ
3xy 2 + 2x 2 y 3 + = 3xy 2 + 2x 2 y 3
dy
dϕ
Por último, se despeja en la ecuación formada y se integra respecto a
dy
la variable y:
dϕ
=0
dy
ϕ′( y ) = 0
xy 3 + 12 x 2 y 4 − x 2 + C = 0
2
(x )
(2xsen(x) cos(x) + xysen(x) )dx + xsen − 2x cos( x) dy = 0
y
Se trata efectivamente de una EDO no exacta, ya que llamando
xsen2 (x )
G (x, y ) = 2xsen(x ) cos( x) + xysen(x) , H (x, y ) = − 2x cos( x ) , se
y
cumple que:
∂G (x, y )
= xsen(x) ≠
∂y
∂G (x, y ) ∂H (x, y )
−
∂y ∂x
=
∂w ( x, y ) ∂w (x, y )
H (x, y ) − G (x, y )
∂x ∂y
∂w ∂w
• Caso factor integrante µ(x ) ( w = x , = 1, = 0 ):
∂x ∂y
y y ∂w − y ∂w 1
• Caso factor integrante µ ( w = , = , = ):
x x ∂x x 2 ∂y x
y
Se multiplica la EDO original por el factor integrante µ(x, y ) = ,
x
obteniéndose la EDO exacta:
∫ (2ysen(x) cos(x) + y )
2
F (x, y ) = sen(x ) dx = ysen2 (x ) − y 2 cos( x) + ϕ( y )
∂F (x, y )
∂y
=
∂
∂y
(
ysen2 (x) − y 2 cos( x) +
dϕ
dy
)
= sen2 (x) − 2y cos( x) +
dϕ
dy
dϕ
sen2 (x ) − 2y cos( x) + = sen2 (x) − 2y cos(x )
dy
dϕ
=0
dy
dϕ
Por último, se despeja en la ecuación formada y se integra respecto a
dy
la variable y:
dϕ
=0
dy
ϕ′( y ) = 0
ysen2 (x ) − y 2 cos( x ) + C = 0
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 81
2 + 1 dx + x − 1 dy = 0 ; y (1) = 1
xy 2 y y3
Se trata efectivamente de una EDO no exacta, ya que llamando
1 x 1
G (x, y ) = 2 + 2 , H (x, y ) = − 3 , se cumple que:
xy y y
∂G (x, y ) − 2 ∂H (x, y ) 1
= 3 ≠ =
∂y xy ∂x y
−2 1 −2 1 − 2 − xy 2
− −
xy 3 y xy 3 xy 3 y
1 1
= = = 2
= = = f (w (x, y ))
1 1 2 − 2 − xy xy w
x − 2 − 2x − 2 − 2 −x
y y y y2
1 x
2xy + dx + x2 − dy = 0
y y 2
1
Se integra µ(x, y )G (x, y ) = 2xy + parcialmente respecto a x:
y
1 x
F (x, y ) = ∫ 2xy + dx = x2 y + + ϕ( y )
y y
∂F (x, y ) ∂ 2 x dϕ x dϕ
= x y + + = x2 − 2 +
∂y ∂y y dy y dy
x dϕ x
x2 − 2
+ = x2 − 2
y dy y
82 Ecuaciones Diferenciales
dϕ
Por último, se despeja en la ecuación formada y se integra respecto a
dy
la variable y:
dϕ
=0
dy
ϕ′( y ) = 0
1
(1)2 (1) + +C = 0
1
2+C = 0
C = −2
Por tanto, la solución del problema de valor inicial es:
x
x2 y + −2 = 0
y
Teorema (Multiplicidad del factor integrante). Dada una EDO de primer orden no exacta
G (x, y )dx + H (x, y )dy = 0 , si µ1(x, y ) y µ2 ( x, y ) son dos factores integrantes distintos de la
misma, entonces la solución general de la ecuación es:
µ1(x, y )
=C
µ2 ( x, y )
Demostración:
∂ µ1(x, y ) ∂ µ1(x, y )
dx + dy = 0
∂x µ2 (x, y ) ∂y µ2 (x, y )
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 83
Por tanto, se trata de probar que la ecuación diferencial anterior es equivalente a la original
G (x, y )dx + H (x, y )dy = 0 , es decir, poniendo ambas ecuaciones en forma normal debe verificarse que:
∂µ ( x, y ) ∂µ (x, y )
µ1(x, y ) 2 − µ2 (x, y ) 1
dy − G (x, y ) ∂x ∂x
= =
dx H (x, y ) ∂µ ( x, y ) ∂µ (x, y )
µ2 (x, y ) 1 − µ1(x, y ) 2
∂y ∂y
En efecto, como µ1(x, y ) y µ2 (x, y ) son dos factores integrantes de la ecuación original, al multiplicarla por
cada uno de ellos se obtienen sendas ecuaciones diferenciales exactas:
∂ ∂
(µ1(x, y )G (x, y )) = (µ1(x, y )H (x, y ))
∂y ∂x
∂ ∂
(µ2 (x, y )G (x, y )) = (µ2 (x, y )H (x, y ))
∂y ∂x
de donde resulta:
multiplicando la primera de estas igualdades por µ2 (x, y ) , la segunda por µ1(x, y ) , y restándolas miembro a
miembro, se obtiene:
−G
y despejando resulta:
H
∂µ (x, y ) ∂µ (x, y )
µ1(x, y ) 2 − µ2 (x, y ) 1
− G (x, y ) ∂x ∂x
=
H (x, y ) ∂µ (x, y ) ∂µ (x, y )
µ2 ( x, y ) 1 − µ1(x, y ) 2
∂y ∂y
y y
Ejemplo: Dada la EDO y cos − x dx + − x cos dy = 0 , deducir la solución
x x
general a partir de un factor integrante dependiente sólo de x, y un factor
y
integrante dependiente de .
x
Se trata efectivamente de una EDO no exacta, ya que llamando
y y
G (x, y ) = y cos − x , H (x, y ) = − x cos , se cumple que:
x x
y
2 cos
= x = − 2 = f (x )
y x
− x cos
x
y
2 cos
x y
= = 2 cos = 2 cos(w ) = f (w ( x, y ))
y y y y x
cos − cos + 1
x x x x
Efectivamente, la EDO admite un factor integrante dependiente del
cociente y/x, de la forma:
y
2sen
= Ce ∫
2 cos(w )dw
µ2 (w ) = Ce ∫ f (w )dw = Ce2sen(w ) = {C = 1} = e2sen(w ) = e x
1
x2 =C
y
2sen
e x
y
−2sen
1 x
e =C
x2
a0 (x) b( x )
donde G (x ) = y H (x ) = . Si b(x ) = 0 (respectivamente H ( x ) = 0 ), la ecuación
a1(x) a1(x )
lineal se llama homogénea, y en caso contrario, b( x) ≠ 0 (respectivamente H ( x) ≠ 0 ), la
ecuación lineal se llama completa. Además, si a0 y a1 son constantes (respectivamente G
es constante) la ecuación lineal se llama ecuación lineal de coeficientes constantes, y en
caso contrario, si a0 o a1 dependen de x (respectivamente G depende de x), la ecuación
lineal se llama ecuación lineal de coeficientes variables.
y′ + 5y = 0
dy
tg (x ) + xy = 0
dx
es una ecuación lineal homogénea de coeficientes variables.
dx
+ x = 2t
dt
es una ecuación lineal completa de coeficientes constantes.
86 Ecuaciones Diferenciales
y
y′ + = ln(1 + x )
x
es una ecuación lineal completa de coeficientes variables.
Procedimiento (Resolución de una EDO lineal homogénea). Para resolver una ecuación
lineal homogénea y ′ + G (x) ⋅ y = 0 (EDO de variables separadas) se procede de la
siguiente forma:
1. Se separan los diferenciales en la igualdad, poniendo dy a la izquierda y dx a la
derecha:
dy
+ G (x) ⋅ y = 0
dx
dy
= −G (x ) ⋅ y
dx
dy = −G (x) ⋅ ydx
1
dy = −G (x )dx
y
y = eCe ∫ −G ( x )dx
y = Ce ∫ −G ( x )dx
2xy
dy = 2
dx
x + 4x + 5
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 87
1 2x
dy = 2 dx
y x + 4x + 5
1 2x
∫ y dy = ∫ x2 + 4x + 5 dx
2x + 4 − 4 2x + 4 4
ln( y ) = ∫ x2 + 4x + 5 dx = ∫ x2 + 4x + 5 dx − ∫ x2 + 4x + 5 dx =
(
= ln x2 + 4x + 5 − 4 ∫) 1
(x + 2)2 + 1
(
dx = ln x2 + 4x + 5 − 4arctg(x + 2) + C )
( )
∴ ln( y ) = ln x2 + 4x + 5 − 4arctg(x + 2) + C
( )
y = eC x2 + 4x + 5 e−4arctg( x +2)
( )
y = C x2 + 4 x + 5 e−4arctg( x +2)
y ′ − x2e2x y = 0
dy = x2e2x ydx
1
dy = x2e2xdx
y
1 2 2x
∫ y dy = ∫ x e dx
u = x2 dv = e2xdx
ln( y ) = ∫ x2e2xdx = 1 2x
= 1 x2e2x
2
− ∫ xe2xdx =
du = 2x v = 2 e
u = x dv = e2xdx
= = 1 x2e2x − 1 xe2x + ∫ 12 e2xdx =
1 2x 2 2
du = 1 v = 2 e
= 1 x2e2x
2
− 1 xe2x
2
+ 1 e2x
4
+C = (
1 x2
2
− 1x
2
+ 1
4
)e 2x
+C
∴ ln( y ) = ( 1 x2
2
− 1x
2
+ 1
4
)e
2x
+C
y = eCe 2
(1 x2 − 12 x + 14 )e2x
y = Ce 2
(1 x2 − 12 x + 14 )e2x
Ejemplo: Resolver el problema de valor inicial:
dx
+
tx
dt sen2 (t )
= 0; x (2π ) = 2
88 Ecuaciones Diferenciales
dx −tx
=
dt sen2 (t )
−tx
dx = dt
sen2 (t )
1 −t
dx = dt
x sen2 (t )
1 −t
∫ x dx = ∫ sen2 (t ) dt
−1
−t u=t dv = dt
ln( x ) = ∫ dt = sen 2
(t ) = tcotg(t ) − ∫ cotg(t )dt =
sen2 (t ) du = dt v = cotg(t )
cos(t )
= tcotg(t ) − ∫ dt = tcotg(t ) − ln(sen(t )) + C
sen(t )
eCetcotg(t )
x =
sen(t )
Cetcotg(t )
x=
sen(t )
x (2π ) = 2
Ce 2
π cotg π
2
()
2=
sen 2π ()
Ce0
2=
1
C=2
Por tanto, la solución del problema de valor inicial es:
2etcotg(t )
x=
sen(t )
Fr
50a = 490 − 2v
donde a es la aceleración de la caída del cuerpo. De aquí, como a = v ′ ,
resulta la ecuación lineal completa de coeficientes constantes y término
independiente constante que determina la velocidad de descenso del
cuerpo:
50v ′ = 490 − 2v
además, dado que el cuerpo parte del reposo se trata de resolver el
problema de valor inicial:
50v ′ = 490 − 2v ; v (0) = 0
La solución general de la ecuación es:
dv
25 = 245 − v
dt
−1 −1 dt
dv = 25
245 − v
90 Ecuaciones Diferenciales
−1
∫ 245 − v dv = ∫ 25
−1 dt
(
ln 245 − v = ) −1 t
25
+C
−1 t
245 − v = eCe 25
−1 t
245 − v = Ce 25
−1 t
v = 245 − Ce 25
Sustituyendo ahora la condición v (0) = 0 en la solución general, se
obtiene el valor de la constante de integración C:
v (0) = 0
−1 ( 0 )
0 = 245 − Ce 25
0 = 245 − C
C = 245
Por tanto, la velocidad de caída del cuerpo en cada instante t es:
−1 t
v = 2451 − e 25
−1 t
x′ = 2451 − e 25
y considerando que la distancia recorrida en el instante inicial es nula se
trata de resolver el problema de valor inicial:
−1 t
x′ = 2451 − e 25 ; x (0) = 0
La solución general de la ecuación es:
dx −1 t
= 2451 − e 25
dt
−1 t
dx = 2451 − e 25 dt
−1 t
∫ dx = ∫
2451 − e 25 dt
−1 t
x = 245 t + 25e 25 + C
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 91
−1 ( 0 )
0 = 245 0 + 25e 25 + C
0 = 6125 + C
C = −6125
Por tanto, la distancia recorrida por el cuerpo en cada instante t es:
−1 t
x = 245 t + 25e 25 − 6125
Altitud inicial: Dado que el cuerpo tarda en llegar al nivel del mar 25
segundos, la altitud inicial coincide con la distancia total recorrida por el
cuerpo en 25 s, es decir:
( )
x = 245 25 + 25e−1 − 6125 ≅ 2253'26 m
El siguiente teorema afirma que la solución general de una ecuación lineal completa
y ′ + G ( x) ⋅ y = H (x ) puede obtenerse sumando a la solución general de la correspondiente
ecuación lineal homogénea y ′ + G (x) ⋅ y = 0 una solución particular de la ecuación lineal
completa.
Teorema (Solución general de la ecuación lineal completa). Dada una ecuación lineal
completa y ′ + G (x) ⋅ y = H (x ) , si y = f (x, C) es la solución general de la correspondiente
ecuación lineal homogénea y ′ + G (x) ⋅ y = 0 , e y = g(x ) es una solución particular no
trivial de la ecuación lineal completa, entonces la solución general de la ecuación lineal
completa es:
y = f (x, C) + g(x)
Demostración:
Sea y = f (x, C) la solución general de la ecuación lineal homogénea, e y = g(x) una solución particular de la
ecuación lineal completa. Se trata de probar que la solución general de la ecuación lineal completa es
y = f (x, C) + g (x) . Derivando esta expresión, se obtiene:
y ′ = f ′(x, C) + g ′(x)
Si se llevan las expresiones de la solución general de la ecuación lineal completa y de su derivada a la ecuación
lineal completa y ′ + G (x ) ⋅ y = H (x ) , resulta:
0 + H (x ) = H (x)
H (x) = H (x )
es decir, la familia uniparamétrica de funciones y = f (x, C) + g (x) verifica la ecuación lineal completa, luego
constituye su solución general, como se quería probar.
Ejemplo: Integrar la siguiente ecuación lineal completa, sabiendo que una solución
particular es y = x 2 :
( )
xy ′ + x 2 − 1 y = x 4 + x 2
(
xy ′ + x 2 − 1 y = 0 )
x
dy
dx
( )
+ x2 − 1 y = 0
xdy = (1 − x )ydx 2
1 1 − x2
dy = dx
y x
1 1 − x2
∫ y dy = ∫ x
dx
2
1
( ) ∫ 1 −xx
ln y = dx = ∫ − x dx = ln x − ( ) 1
2
x2 + C
x
( )
∴ ln y = ln x − ( ) 1
2
x2 + C
− 1 x2
y = eC x e 2
− 1 x2
y = ±eC xe 2
− 1 x2
y = Cxe 2
Las ecuaciones lineales completas son, en general, ecuaciones no exactas, pero pueden
convertirse en exactas e integrarse mediante el método del factor integrante. Además, en
este caso el cálculo del factor integrante es sencillo, ya que siempre es un factor que
depende sólo de la variable independiente de la ecuación, lo que se prueba en la
siguiente proposición.
µ(x) = Ce ∫ G ( x )dx
Demostración:
dy
+ G (x ) ⋅ y = H (x )
dx
(G(x) ⋅ y − H(x))dx + dy = 0
luego, de acuerdo con la condición necesaria y suficiente para que una ecuación diferencial ordinaria admita
un factor integrante dependiente sólo de x, se tiene que:
∂
(G(x) ⋅ y − H(x)) − ∂ (1)
∂y ∂x G (x ) − 0
= = G ( x ) = f (x )
1 1
es decir, toda ecuación lineal completa admite un factor integrante de la forma µ(x ) . Para hallar el factor
integrante, basta aplicar la fórmula:
Observación. Nótese que la expresión del factor integrante de una ecuación lineal
completa, µ(x) = Ce ∫ G ( x )dx , coincide con la inversa de la solución general de la ecuación
lineal homogénea correspondiente, y = Ce ∫ −G ( x )dx .
Procedimiento (Resolución de una EDO lineal completa mediante el método del factor
integrante). Para resolver una ecuación lineal completa y ′ + G (x) ⋅ y = H (x ) mediante un
factor integrante se procede de la siguiente forma:
1. Se calcula el factor integrante dependiente sólo de x asociado a la ecuación lineal
mediante la fórmula:
µ(x) = Ce ∫ G ( x )dx
∂F (x, y ) dϕ
• = µ ′(x) ⋅ y + = µ(x)(G (x) ⋅ y − H (x, y ) )
∂x dx
Ejemplo: Integrar la siguiente EDO lineal completa por el método del factor
integrante:
y
y′ + + x3 = 0
x
Expresando la ecuación en forma estándar, queda:
1
y′ + y = −x 3
x
1
donde G (x ) = , H (x) = −x 3 . Esta EDO lineal completa admite un factor
x
integrante dependiente de la variable x, de la forma:
1
∫ x dx ( )
ln x
µ( x ) = Ce ∫ G ( x )dx = Ce = Ce = C x = {C = 1} = x
xy ′ + y = − x 4
dy
x + y = −x 4
dx
(x 4
)
+ y dx + xdy = 0
∂F (x, y )
=
∂
(xy ) + dϕ = y + dϕ
∂x ∂x dx dx
dϕ
y+ = x4 + y
dx
dϕ
Por último, se despeja en la ecuación formada y se integra respecto a
dx
la variable x:
dϕ
= x4
dx
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 95
ϕ′(x ) = x 4
4
∫ ϕ′(x)dx = ∫ x dx
ϕ(x) = 1
5
x5 + C
xy + 1
5
x5 + C = 0
C
y = − 15 x 4 −
x
Ejemplo: Resolver el siguiente problema de valor inicial asociado a una EDO lineal
completa por el método del factor integrante:
(x 2
)
+ 1 y ′ + 4 xy = x ; y (0 ) = −1
(x 2
)
2
( ) (
+ 1 y ′ + 4x x 2 + 1 y = x x 2 + 1 )
(x + 2x + 1) dy
4
dx
+ (4x y + 4xy ) = x + x
2 3 3
∫ (x ) ( )
4
F (x, y ) = + 2x 2 + 1 dy = x 4 + 2x 2 + 1 y + ϕ(x )
∂F (x, y )
∂x
=
∂
∂x
((
x 4 + 2x 2 + 1 y +
dϕ
dx
= 4x 3 + 4x y +
dϕ
dx
)) ( )
dϕ
4 x 3 y + 4 xy + = 4x 3 y + 4 xy − x 3 − x
dx
dϕ
Por último, se despeja en la ecuación formada y se integra respecto a
dx
la variable x:
96 Ecuaciones Diferenciales
dϕ
= −x 3 − x
dx
ϕ′(x ) = − x 3 − x
∫ ϕ′(x)dx = ∫ (− x )
3
− x dx
ϕ(x) = −1 x 4 − 1 x2 + C
4 2
(x 4
)
+ 2x 2 + 1 y − 1
4
x4 − 1
2
x2 + C = 0
1 +C
1 4
y= −
4
(x 2
+1 )
2
1 − 1 − C = −1
4 4
C =1
Por tanto, la solución del problema de valor inicial es:
5
1 4
y= −
4
(x 2
+1 )
2
Este método, de nombre cuando menos curioso (variación de constantes), fue ideado por
Joseph Louis Lagrange partiendo de una hipótesis en cierta forma similar a la planteada
en el teorema de la solución general de la ecuación lineal completa. La hipótesis consiste
en suponer que la ecuación lineal completa y ′ + G ( x) ⋅ y = H (x ) tiene una solución general
y = f ( x,C) de la misma forma que la solución general de la correspondiente ecuación
lineal homogénea y ′ + G (x ) ⋅ y = 0 , es decir, y = f (x,C) = Ce ∫ −G ( x )dx , pero tomando la
constante de integración C como una función que depende de la variable x.
Desarrollando esta idea, la proposición siguiente establece la expresión de la solución
general de la ecuación lineal completa.
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 97
Demostración:
solución general y = Ce ∫ −G ( x )dx de la correspondiente ecuación lineal homogénea, pero tomando la constante
de integración C como una función que depende de la variable x, es decir, una expresión de la forma
y = C (x )e ∫ −G ( x )dx . Derivando esta expresión, se obtiene:
Si se llevan las expresiones de la solución general de la ecuación lineal completa y de su derivada a la ecuación
lineal completa y ′ + G (x ) ⋅ y = H (x ) , resulta:
C′(x )e ∫ −G ( x )dx = H (x )
H (x)
C′(x ) =
e ∫ −G ( x )dx
C′(x ) = H (x )e ∫ G ( x )dx
C (x ) = ∫ H (x )e ∫ G ( x )dxdx
y = Ce ∫ −G ( x )dx
y = C(x )e ∫ −G ( x )dx
(C′(x)e ∫ −G ( x )dx ) (
− C (x )G ( x)e ∫ −G ( x )dx + G (x ) C (x )e ∫ −G ( x )dx = H (x ) )
C ′(x ) = H (x )e ∫ G ( x )dx
4. Se integra C'(x), y se lleva el factor resultante C(x) a la solución general de la
ecuación lineal:
dy
+ 2xy = 0
dx
dy = −2xydx
1
dy = −2xdx
y
1
∫ y dy = ∫ − 2xdx
ln( y ) = − x2 + C
2
y = eCe − x
2
y = Ce− x
2
C (x ) = 2e x + C
2 2 2
y = 2e x + C e− x = 2 + Ce − x
E
−+
L R
10 ⋅ i + 1 ⋅ i′ = 12sen(500t )
50
di
+ 500 ⋅ i = 600sen(500t )
dt
En primer lugar, se resuelve la ecuación lineal homogénea asociada:
di
+ 500 ⋅ i = 0
dt
1
di = −500dt
i
1
∫ i di = ∫ − 500dt
ln( i ) = −500t + C
i = eCe−500t
100 Ecuaciones Diferenciales
i = Ce−500t
Se deriva esta solución general tomando la constante de integración C
como una función dependiente de t:
i = C (t )e −500t
di
= C′(t )e −500t − 500C (t )e−500t
dt
Se llevan ambas expresiones a la ecuación lineal completa y se despeja
C'(x):
u = sen(500t ) dv = 600e500t dt
C(t ) = ∫ 600sen(500t )e500t dt = =
du = 500 cos ( 500t) v = 65 e500t
u = cos(500t ) dv = 600e500t dt
= =
du = −500sen( 500t) v = 65 e500t
∴ C (t ) = ∫ 600sen(500t )e500t dt = 3
5
e500t (sen(500t ) − cos(500t )) + C
y ′ + G (x ) ⋅ y = H (x) ⋅ y p (p ∈ o)
y ′ + y = xy 4
y ′ + sen(x ) y = (x + 1) y 3
dx x
+ = tx 2
dt t
1
es una ecuación de Bernoulli con G (t ) = , H (t ) = t , p = 2 .
t
y ′ + G (x ) ⋅ y = H (x) ⋅ y 0 ⇔ y ′ + G (x ) ⋅ y = H (x)
Demostración:
y′ + G (x) ⋅ y = H (x) ⋅ y p
dy
+ G (x) ⋅ y = H (x ) ⋅ y p
dx
dy = H (x ) ⋅ y p − G (x) ⋅ y dx
(1 − p) y − pdy = (1 − p) y − p H (x ) ⋅ y p − G (x ) ⋅ y dx
u = y1− p
Si a esta ecuación se le aplica el cambio , resulta:
du = (1 − p) y − pdy
du
+ (1 − p)G (x ) ⋅ u = (1 − p)H (x)
dx
102 Ecuaciones Diferenciales
u′ + (1 − p)G (x ) ⋅ u = (1 − p)H (x )
y − y2
y′ − =
x x3
−1 −1
Se trata de una EDO de Bernoulli con G (x ) = , H (x ) = , p = 2.
x x3
u = y −1
Aplicando el cambio , resulta la ecuación lineal completa:
du = − y − 2dy
u 1
u′ + =
x x3
−1
ϕ′(x ) =
x2
−1
∫ ϕ′(x)dx = ∫ x 2 dx
1
ϕ( x) = +C
x
La solución general de la ecuación lineal completa es:
1
xu + +C = 0
x
x
y=
−1
−C
x
y ′ + y = xy 3
u′ − 2u = −2x
En este caso, se aplica el método de variación de constantes. Primero, se
resuelve la ecuación lineal homogénea asociada:
u′ − 2u = 0
du
− 2u = 0
dx
104 Ecuaciones Diferenciales
du = 2udx
1
du = 2dx
u
1
∫ u du = ∫ 2dx
( )
ln u = 2x + C
u = eC e 2x
u = Ce 2x
Tomando la constante de integración C como una función dependiente de
x, se deriva esta solución general:
u = C(x )e 2x
u′ = C ′(x )e 2x + 2C(x )e 2x
u = x dv = −2e −2x dv
C(x ) = ∫ − 2xe − 2x dx = =
du = dx e − 2x
(
u = xe −2x − 12 e −2x + C e 2x )
u = x − 12 + Ce 2x
y ′ = a( x ) y 2 + b(x ) y + c(x )
y′ = y2 + y + 3
y ′ = xy 2 + x 2 y + 3x 3
dx
= t ⋅ x 2 + (t − 3)x + t 3
dt
1
y = g(x) +
u
y ′ = a(x ) y 2 + b(x ) y + c(x ) ⇒ ⇒ u′ + (2a(x )g(x ) + b(x ) )u = −a(x )
1
dy = g ′(x )dx − du
u 2
106 Ecuaciones Diferenciales
Demostración:
y′ = a( x) y 2 + b( x) y + c(x )
1
y = g(x ) +
u
Si a esta ecuación se le aplica el cambio 1 , donde y = g(x) es una solución particular de
dy = g ′(x)dx − du
u 2
la misma, resulta:
2
u′ 1 1
g ′(x ) − = a(x ) g(x ) + + b(x) g(x ) + + c(x )
2 u u
u
u′ 1 1 1
g ′(x ) − = a(x ) g 2 (x ) + 2g(x ) + + b(x ) g(x ) + + c(x )
2 u 2 u
u u
u′ 1 1 1
g ′(x ) − = a(x)g 2 (x ) + b( x)g( x) + c(x ) + a(x ) 2g(x ) + + b(x )
u2 u u2 u
y como y = g(x) es solución particular de ecuación se cumple que g ′(x ) = a(x)g 2 (x ) + b(x)g (x) + c(x ) , luego:
u′ 1 1 1
− = a(x ) 2g(x ) + + b(x )
u2 u u2 u
u′ = a(x)(−2g (x )u − 1) − b(x)u
1
y = g(x ) +
u
⇒ u′ + (2a( x )g( x ) + b(x ) )u = −a(x )
2
y ′ = a(x ) y + b(x ) y + c(x ) ⇒
1
dy = g ′(x )dx − du
u2
y ′ = y 2 − 2xy + 1 + x 2
y′ = m
m = m2 x 2 − 2mx 2 + 1 + x 2
( )
m2 x 2 + m − 2x 2 − 1 + 1 + x 2 = 0
m=
2x 2 + 1 ± (− 2x − 1)
2 2
(
− 4x 2 1 + x 2 )=
2
2x
1
2x 2 + 1 ± 4 x 4 + 4x 2 + 1 − 4x 2 − 4x 4 2x 2 + 1 ± 1 1 +
= = = x2
2x 2 2x 2 1
u ′ = −1
Se trata de una EDO de variables separadas, de integración directa:
du
= −1
dx
du = −dx
∫ du = ∫ − dx
u = −x + C
Finalmente, se deshace el cambio en la solución obtenida, sustituyendo
1
u= :
y−x
1
= −x + C
y−x
1
y=x+
− x +C
(1 + x )y ′ + 2xy
3 2
+ x2y + 1 = 0
− 2x x2 1
y′ = y2 − y−
3 3
1+ x 1+ x 1 + x3
−2x − x2 −1
que es una EDO de Riccati con a(x ) = , b(x ) =
. , c(x) =
3 3
1+ x 1+ x 1 + x3
En primer lugar, se debe obtener la solución particular, para lo cual se
deriva y = mx , y se llevan y e y' a la ecuación para así despejar m:
y′ = m
− 2x 3 x3 1
m= m2 − m−
3 3
1+ x 1+ x 1 + x3
2x 3 2x 3 + 1 1
m2 + m+ =0
3 3
1+ x 1+ x 1+ x3
(
2x 3 m2 + 2x 3 + 1 m + 1 = 0 )
m=
− 2x 3 − 1 ± (2x 3
+1 )2
− 8x 3
=
− 2x 3 − 1 ± 4x 6 − 4x 3 + 1
=
4x 3 4x 3
=
− 2x 3 − 1 ± (2x 3
−1 )2
=
− 2x 3 − 1 ± 2x 3 − 1 ( ) = 2−x1 3
4x 3 4x 3 − 1
4x 2 x2
u′ + − u = 2x
1 + x3 1 + x3 1 + x3
3x 2 2x
u′ + u=
3
1+ x 1 + x3
En este caso, se aplica el método de variación de constantes. Primero, se
resuelve la ecuación lineal homogénea asociada:
3x 2
u′ + u=0
1 + x3
du 3x 2
+ u=0
dx 1 + x 3
1 − 3x 2
du = dx
u 1 + x3
1 − 3x 2
∫u du = ∫ 1 + x 3 dx
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 109
ln u = − ln 1 + x 3 + C
( )
eC
u =
1 + x3
C
u=
1 + x3
Tomando la constante de integración C como una función dependiente de
x, se deriva esta solución general:
C( x )
u=
1 + x3
C ′(x) 3x 2C(x )
u′ = −
1 + x3 (1 + x ) 3 2
C ′(x) = 2x
∫ C ′(x)dx = ∫ 2xdx
C( x ) = x 2 + C
x2 + C
u=
1 + x3
1
Finalmente, se deshace el cambio de variable sustituyendo u = , y se
y+x
obtiene la solución general de la ecuación de Riccati:
1 x2 + C
=
y + x 1 + x3
1 + x3
y= −x
x2 + C
110 Ecuaciones Diferenciales
y = xtg(2 y ′) + cos( y ′ + π)
x′
x = tx ′ 2 +
x′ + 1
x′
es una ecuación de Lagrange con f (x ′) = x ′2 , g (x′) = .
x′ + 1
y = xy ′2 + xy ′ + 3 y ′
Demostración:
y = xf (t ) + g (t )
dt dt
y ′ = t = f (t ) + xf ′(t ) + g ′(t )
dx dx
dt
t = f (t ) + (xf ′(t ) + g ′(t ) )
dx
dt
de donde despejando , se obtiene:
dx
dt t − f (t )
=
dx xf ′(t ) + g ′(t )
e invirtiendo esta ecuación, lo que equivale a tomar x como variable dependiente y t como variable
independiente, se obtiene:
dx xf ′(t ) + g ′(t )
=
dt t − f (t )
dx f ′(t ) g ′(t )
− x=
dt t − f (t ) t − f (t )
( ) (
y = xf h−1( x,C) + g h−1(x,C) )
Ejemplo: Hallar la solución general de la ecuación de Lagrange:
y = xy ′2 + y ′ + 5
112 Ecuaciones Diferenciales
y = xt 2 + t + 5
dt dt
y ′ = t = t 2 + 2xt +
dx dx
dt
t = t 2 + (2xt + 1)
dx
t − t2 dt
=
2xt + 1 dx
dx 2xt + 1
=
dt t − t2
dx 2 1
− x=
dt 1 − t t − t2
En este caso, se aplica el método de variación de constantes. Primero, se
resuelve la ecuación lineal homogénea asociada:
dx 2
− x=0
dt 1 − t
1 2
dx = dt
x 1− t
1 2
∫ x dx = ∫ 1 − t dt
ln( x ) = −2ln(1 − t ) + C
eC
x =
(1 − t )2
C
x=
(1 − t )2
dx C′(t ) 2C(t )
= 2
+
dt (1 − t ) (1 − t )3
1− t
C′(t ) =
t
Se integra C'(t), y se lleva el resultado a la solución general de la ecuación
lineal:
1− t
∫ C′(t )dt = ∫ t
dt
C(t ) = ln( t ) − t + C
x=
( )
ln t − t + C
(1 − t )2
(1 − t )2
y = x (1 + y ′) + y ′2
y = x (1 + t ) + t 2
dt dt
y′ = t = 1 + t + x + 2t
dx dx
dt
0 = 1 + (x + 2t )
dx
−1 dt
=
x + 2t dx
dx
= −x − 2t
dt
dx
+ x = −2t
dt
En este caso, se aplica el método de variación de constantes. Primero, se
resuelve la ecuación lineal homogénea asociada:
dx
+x=0
dt
1
dx = −dt
x
114 Ecuaciones Diferenciales
1
∫ x dx = ∫ − dt
( )
ln x = −t + C
x = eC e −t
x = Ce −t
Tomando la constante de integración C como una función dependiente de
t, se deriva esta solución general:
x = C(t )e −t
dx
= C ′(t )e − t − C(t )e − t
dt
Se llevan ambas expresiones a la ecuación lineal completa y se despeja
C'(t):
C ′(t ) = −2te t
(
x = (−2t + 2)e t + C e −t )
x = −2t + 2 + Ce −t
y = xy ′ + arctg( y ′)
y = xy ′ − 2 y ′
1
x = tx ′ +
x′
1
es una ecuación de Clairaut con g (x′) = .
x′
Demostración:
y = xt + g (t )
dt dt
y′ = t = t + x + g ′(t )
dx dx
dt dt
0=x + g ′(t )
dx dx
y de aquí, resulta:
(x + g′(t )) dt =0
dx
y = xy ′ + y ′ − y ′2
y = xt + t − t 2
dt dt dt
y′ = t = t + x + − 2t
dx dx dx
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 117
dt
(x + 1 − 2t ) =0
dx
y de aquí, pueden considerarse dos casos:
dt
• = 0 ⇒ t = C ⇒ y = Cx + C − C 2 (Solución general).
dx
x = 2t − 1 x = 2t − 1
• x + 1 − 2t = 0 ⇒ x = 2t − 1 ⇒ 2 ⇒ 2
y = (2t − 1)t + t − t y = t
(Solución singular).
En el estudio diferencial de las funciones reales de una variable se define el ángulo entre
dos curvas funcionales derivables y = f (x ) , y = g(x) en un punto de corte (x, y ) entre
ambas, como el ángulo θ que forman entre sí las rectas tangentes a f y g en dicho punto,
el cual viene determinado por la expresión:
θ = α − β = arctg(f ′(x ) ) − arctg(g ′(x ) )
y = g(x)
θ
y y = f (x)
α
β
X
O x
(
tg α + π
2
) = F (−x1, y )
y = g (x )
π
2
y y = f ( x,C)
α+ π
α 2
O x X
−1
y′ =
F (x, y )
x2 + y 2 = C2
2x + 2 yy ′ = 0
−x
y′ =
y
Según esto, la recta tangente a la circunferencia de la familia que pasa
−x
por un punto arbitrario (x,y) tiene como pendiente y ′ = en ese punto.
y
Por tanto, la EDO de primer orden que define a la familia uniparamétrica
de trayectorias y = g(x, K ) que verifican la condición de ortogonalidad
planteada es:
−1 y
y′ = =
F (x, y ) x
Se trata de una EDO de primer orden de variables separadas:
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 119
y
y′ =
x
dy y
=
dx x
dy dx
=
y x
1 1
∫ y dy = ∫ x dx
( )
ln y = ln x + K( )
y =e
( )
ln x + K
y = ±e K x
y = Kx
y = Kx
2 2 2
x + y =C π
2
y
α+ π
2
α
O x X
Al igual que en el caso anterior, supóngase que y = f (x,C ) es una familia uniparamétrica
de curvas planas, y que y ′ = F (x, y ) es la EDO de primer orden asociada a dicha familia.
De acuerdo con el significado geométrico de la derivada, es evidente que F(x,y)
determina la pendiente de la recta tangente a la curva de la familia que pasa por el
punto (x,y). Esta recta tangente forma un ángulo α respecto al eje OX dado por:
tg (α) = y ′ = F (x, y )
Por tanto, la recta tangente a una trayectoria y = g (x ) que corte a la curva y = f (x,C ) en
el punto (x,y) con un ángulo θ ≠ π (trayectoria oblicua) tendrá una pendiente igual a la
2
tangente del ángulo suma α + θ , es decir:
120 Ecuaciones Diferenciales
y = g (x )
θ
y y = f ( x,C)
α+θ
α
O x X
Figura 2.7: Trayectoria y = g (x ) oblicua a cada curva y = f (x,C ) con ángulo constante
θ≠ π .
2
F (x, y ) + tg(θ)
Por definición de derivada se cumple que g ′(x ) = tg(α + θ) = , de donde se
1 − F (x, y )tg(θ)
obtiene la EDO de primer orden en forma normal que define a la familia uniparamétrica
de trayectorias y = g(x, K ) que cortan a cada curva y = f (x,C ) con ángulo constante
θ≠ π :
2
F (x, y ) + tg(θ)
y′ =
1 − F (x, y )tg(θ)
y′ = C
y = Cx y
C = y
⇒ x ⇒ y′ =
′
y = C C = y ′ x
y
+ tg(θ)
F (x, y ) + tg(θ)
y′ = = x
1 − F (x, y )tg (θ) y
1 − tg(θ)
x
Se trata de una EDO de primer orden homogénea:
y
+ tg(θ)
y′ = x
y
1 − tg(θ)
x
y = tx
Para resolverla, se aplica el cambio , obteniéndose:
dy = xdt + tdx
y
+ tg(θ)
dy
= x
dx 1 − y tg(θ)
x
y
+ tg(θ)
dy = x dx
y
1 − tg(θ)
x
tx
+ tg(θ)
xdt + tdx = x dx
tx
1 − tg(θ)
x
t + tg(θ)
xdt + tdx = dx
1 − t ⋅ tg(θ)
( )
(1 − t ⋅ tg(θ))xdt + t − t 2 ⋅ tg(θ) dx = (t + tg(θ))dx
(1 + t )⋅ tg(θ)dx = (1 − t ⋅ tg(θ))xdt
2
( )
1 arctg (t ) − 1 ln 1+t 2 +K
x = e tg ( θ) 2
1 arctg (t )
± eK e tg (θ)
x=
1 + t2
122 Ecuaciones Diferenciales
1 arctg (t )
Ke tg (θ)
x=
1+ t2
Finalmente, se deshace el cambio de variable:
1 arctg y
tg ( θ) x
Ke
x=
y2
1+
x2
Y
1 arctg y
tg ( θ) x
Ke
x=
y2 y = Cx
1+
x2 θ
y
α+θ
α
O x X
1 arctg y
tg ( θ) x
Ke
Figura 2.8: Trayectoria x = oblicua a cada recta y = Cx con
y2
1+
x2
ángulo constante θ ≠ 2π .
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 123
EJERCICIOS
a. y′ = x4 + y 4 . Solución: D = o2
b. y′ = 3 x + y .
Solución: D = {(x, y ) ∈ o 2 | y ≠ x }
c. y′ = x y .
Solución: D = {(x, y ) ∈ o 2 | y > 0}
dx 3
d. = tx .
dt
Solución: D = {(x, y ) ∈ o 2 | x ≠ 0}
a. ( )
y ′ = ln 1 + x2 y 2 ; y (0) = 0 .
b. y′ = 3
xy 2 ; y (0) = 0 .
c. y′ = y 3 ; y (0) = 0 .
d. y ′ = 3 x 3 + y 3 ; y (0) = 0 .
a. ( )
(x − 1) y 2 − 1 dx + x2dy = 0 . Solución:
1
x
y = th ln( x ) + + C
b. y ′ = x ln( x ) cos2 ( y ) .
Solución: (
y = arctg 12 x 2 ln( x ) − 1
4
x2 + C )
1 − cos (x)
c. y ′sen(x) = y . Solución: y=C
1 + cos(x )
1
d. x 2 ydy = 1 − y 2 dx . Solución: 1 − y2 = +C
x
1 x − 1 sen(2x )
b. y ′ = sen2 (x )(1 + y ) ; y (0) = 0 . Solución: y = e2 4 −1
x2
a. (y 2
)
− xy dx +
x+y
dy = 0 .
1 5
b. dx + dy = 0 .
4 3
x +x y x + y2
2
x − 3y
c. y ′ = cos .
2x + 2y
3 2
dy 2x − xy
d. = .
dx x + 2y
y
4 xy − y 2 x
b. y′ = . Solución: x =C
2x 2 y
2−
x
y 2 − 4x 2 C
c. y′ = . Solución: x=
2xy 2
y
+4
x2
C
d. 2xyy ′ = x 2 − y 2 . Solución: x=
3
3y 2
−1
x2
y
a. xtg + y dx − xdy = 0 ; y (1) = 6π . Solución: y = x ⋅ arcsen 12 x ( )
x
b. ( x+y + )
x − y dx + ( x−y − )
x + y dy = 0 ; y (1) = 1 .
1
Solución: x=
y2
1+ 1−
x2
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 125
x
e
−x t
2 t x
dx (t + x) e + x2 − e 1+
c. = ; x (1) = 1 . Solución: t= e2e t
dt tx
x2 + y 2
− 5
y
xy 2e 2 xe
d. y′ = ; y (1) = 2 . Solución: x=
x + y x2 + y 2
2 y
b. (x − 2y + 4)dy = (x − y + 1)dx .
C
Solución: x −2 =
2( y − 3) 2( y − 3)2
1− +
x −2 (x − 2) 2
2x − 2y + 4
c. y′ = .
2x − 2 y − 1
5
2( y − x) − 4
−5 − 8
Solución: x= ln 1 − +
16 2( y − x) + 1 2
2( y − x ) − 4
1 −
2( y − x) + 1
5
5 ln 1 + 2( y − x ) − 4
+ 16 − 8 +C
2( y − x ) + 1 2( y − x) − 4
1+
2( y − x ) + 1
dx t +x+2
d. = ; x(1) = 1 .
dt 2t + 2x − 1
Solución: t = 23 (t + x) − 5 ln
9
(3t + 3x + 1 ) + 59 ln(7) − 13
9. Verificar si las siguientes ecuaciones diferenciales son exactas:
1 ln(x )
a. 2xe y + dx + x 2e y − dy = 0 .
xy y 2
2 y 1 y
a. y − 2 dx + 2xy + dy = 0 . Solución: xy 2 + +C = 0
x x x
1 1 y x
b. + − dx + y cos( y ) + dy = 0 .
x x2 x2 + y 2 x2 + y 2
1 x
Solución: ln( x ) − − arctg + ysen( y ) + cos( y ) + C = 0
x y
c. (x 3
) (
− yexy dx + 2y − xexy dy = 0 . ) Solución: 1 x4
4
− e xy + y 2 + C = 0
1 dρ θ
d. − θ dθ = ρ − e .
2
− 3 + 4ρ − ρ
Solución: eθ − θρ + arcsen(ρ − 2) + C = 0
− ex − y
a. y′ = ; y (0) = 1 .
2 + x + ye y
Solución: ex + xy + 2 y + ( y − 1)e y − 3 = 0
b.
1− y
2
+ cos( x ) − 2xy
2
(
dy = y + ysen(x ) dx ; y (0) = 1 . )
1− y
Solución: − xy 2 + y cos(x ) + arcsen( y ) + 1 − y 2 − 1 − π
2
=0
c. (x + y )2dx + x2 + 2xy + 2 y − 21 dy = 0 ; y (1) = 1 .
1+ y
Solución: 1 (x
3
+ y )3 − 1
3
( )
y 3 + ln 1 + y 2 − arctg( y ) + π
4
− ln(2) − 7
3
=0
d.
( )
1 − x2 + y 3 dx + 3xy 2 + 1 dy = 0 ; y (0) = 0 .
Solución: x3 y + y + 1 arcsen( x ) + 1 x
2 2
1 − x2 = 0
a. (x 2
) (
− y 2 − x − y dx + y 2 − x2 dy = 0 , ) sabiendo que admite un factor
integrante dependiente de x+y.
Solución: 1 x2 − xy − x + 1 y2 + C = 0
2 2
b. (1 − x y )dx + x (y − x )dy = 0 ,
2 2
sabiendo que admite un factor integrante
dependiente sólo de x.
1
Solución: 1
2
y 2 − xy − +C = 0
x
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 127
y2
c. y′ = , sabiendo que admite un factor integrante dependiente de
x2 + xy
x2y.
y
Solución:
x
( )
+ ln y + C = 0
d. (3y 2
) ( )
− x dx + 2 y 3 − 6xy dy = 0 , sabiendo que admite un factor integrante
dependiente de x+y2.
− 2y 2 1
Solución: + +C = 0
(x + y ) 2 2 x + y2
a. (
6xydx + 4 y + 9x2 dy = 0 . ) Solución: 3x 2 y 3 + y 4 + C = 0
b. (x 4
) (
ln( x ) − 2xy 3 dx + 3x2 y 2 dy = 0 . )
Solución: y= 3
− x 3 ln(x) + x 3 − Cx 2
1
−3
2
y2 + −1
d. xy ′ = x .
y
Solución: 1
2
x3y2 + 1
3
x3 − 1
2
x2 + C = 0
14. Resolver los siguientes problemas de valor inicial mediante factores integrantes:
a. (2x 3
) (
y + 3x 2 dx + 4x 2 y 2 + x 4 dy = 0 ; y (0) = 1 .)
2
Solución: x y + 3x + 4
3
y3 − 4
3
=0
3xy 2
b. (x 2
− xy + y 2 dx + ) x+y
dy = 0 ; y (1) = 1 .
5
Solución: y= 3− 1
4
x3 +
4x
y
c. xych(x ) + sh( y ) dx − ych( y )dy = 0 ; y (1) = 0 .
x
Solución: y = argsh(xsh(x ) − xsh(1) )
x π
d. 2arctg( y )dx + 2
dy = 0 ; y (1) = 1 . Solución: y = tg
1+ y 4x 2
128 Ecuaciones Diferenciales
−4 y 2x + 2y
15. Dada la EDO dx + dy = 0 , deducir la solución general a partir de
2
(x − y ) (x − y ) 2
un factor integrante dependiente sólo de y, y un factor integrante dependiente de
x−y.
x−y
Solución: =C
y
dy
b. = y 1 + x2 .
dx
1 x 1+ x2 + 1 argsh( x )
Solución: y = Ce 2 2
1 + 2x − x2
c. y′ + y = 0 .
x
1+ 2x − x2 − arcsen 1 x − 1
Solución: y = Ce 2 2
1
1−
−e x
d. x 2 y ′ = (1 − x ) y ; y (1) = −1 . Solución: y=
x
1 + 3x 3 −1
c. y ′ − 3x 2 y = , sabiendo que una solución particular es y = .
2 x
x
3 1
Solución: y = Ce x −
x
d. y ′ − 2th(x ) y = −sh(x ) , sabiendo que una solución particular es y = ch(x ) .
Solución: y = C ⋅ ch2 (x ) + ch(x )
2xy
b. y′ − = x + 1.
1 + x2
Solución: (
y = 1 + x2 1
2
)( ( )
ln 1 + x 2 + arctg(x ) − C )
dy 3x2 − 1
c. − 3 y = x+3.
dx x −x
Solución: ( )( ( ) (
y = x 3 − x − 3 ln x + ln x + 1 + 2 ln x − 1 − C ) ( ) )
dx tg2 (t ) tg(t ) − t − C
d. +x = . Solución: x=
dt et et
19. Resolver los siguientes problemas de valor inicial asociados a ecuaciones lineales
completas mediante el método del factor integrante:
y
a. y′ − = x2 ; y (1) = 0 . Solución: y = 12 x 3 − 12 x
x
1
b. y ′ + ytg(x ) = ; y (0) = 1 . Solución: y = sen(x ) + cos( x )
cos(x )
dx
d. − x = sen(2t ) ; x (0) = 0 .
dt
Solución: x = −51 (sen(2t ) + 2 cos(2t ) ) + 2
5
et
y
b. 1 − x2 ⋅ y ′ − =1.
arcsen(x )
Solución: ( ( )
y = ln arcsen(x ) + C arcsen(x ) )
dy 3
c. + 3x 2 y = x 2 . Solución: y = 13 + Ce − x
dx
d.
dy
+
y
dx 2 x
−4 = 0. Solución: ( )
y = 8 x − 1 + Ce − x
21. Resolver los siguientes problemas de valor inicial asociados a ecuaciones lineales
completas mediante el método de variación de constantes:
y
a. y′ − = xarctg(x ) ; y (1) = 1 .
x
Solución: ( ) (
y = x 2arctg (x ) − 12 x ln 1 + x 2 + 1 + 12 ln(2) − 4π x )
(2x + 1) y
b. y′ + = e −2x ; y (1) = 1
2
.
x
−2x
Solución: (
y = 12 xe − 2x + 12 e 2 − 1 )e x
130 Ecuaciones Diferenciales
dy 2
c. + 2xy = x ; y (0) = 2 . Solución: y = 12 + 3
2
e−x
dx
dy y 1
d. + = ; y (1) = 1 . Solución: y =1
dx x2 x2
dx 4
b. t + x = t 3x 3 ; x (1) = 4 . Solución: x=
dt t (t − 2) 2
−t −C
x =
a. y = 2xy ′ + ln( y ′) . Solución: t2
y = − 2t − 2C + ln(t )
t
−1 C
3 x = 2 +
b. y = − xy ′ − . Solución: t t
y′ y = − 2 − C t
t
c. y = xy ′ + 1 + y ′ .
−1
x =
Solución: y = Cx + 1 + C (general); 2 1+ t (singular)
−t
y = + 1+ t
2 1+ t
d. y = xy ′ − 2 − y ′ .
x = 1
Solución: y = Cx − 2 − C (general); (singular)
y = −2
24. Obtener las trayectorias ortogonales a las siguientes familias de curvas planas:
1 x2 + 1 y 2
a. Familia de elipses C 2 x 2 + y 2 = 1 . Solución: y = Ke 2 2
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 131
K
b. Familia de hipérbolas x 2 − y 2 = C . Solución: y=
x
25. Obtener las trayectorias oblicuas que cortan a las siguientes familias de curvas
planas con un ángulo constante θ = 4π :
a. Familia de parábolas y = Cx 2 .
3 arctg 4 y + 1
Ke 17 17 x 17
Solución: x=
2y 2 y
+ +1
2 x
x
b. Familia de circunferencias x 2 + y 2 = C .
K
Solución: x=
y
arctg
y2 x
+1 ⋅e
x2