Ecuaciones diferenciales de primer orden

2
Forma general: F (x, y, y 0 ) = 0, ó en forma más explicita y 0 = f (x, y)

Resolver una ecuación diferencial de primer orden es hallar el conjunto de todas sus

soluciones.

Es fácil vericar que en muchos casos una familia de curvas monoparamétricas con-

tiene todas las soluciones de una ecuación diferencial de primer orden. En tal caso se

dice que es la solución general. Al especicar el valor de la constante, se obtiene una

solución particular, y si existen otras soluciones que no se pueden obtener a partir de la

ecuación monoparamétrica, éstas se llaman soluciones singulares.

No existe una técnica sistemátuicca para resolver todas las ecuaciones de primer orden,

pero si existen métodos para resolver ciertas clases de ellas.

2.1. Resolución de ecuaciones diferenciales por integra-

ción directa
Las más sencillas de las ecuaciones diferenciales de primer orden son las que tienen la

forma:

y 0 = f (x)

10
2.1. RESOLUCIÓN DE ECUACIONES DIFERENCIALES POR INTEGRACIÓN DIRECTA

donde f sólo depende de x. Estas se resuelven por integración directa.

dy
= f (x) =⇒ dy = f (x) dx
dx
Z Z
=⇒ dy = f (x) dx

=⇒ y = F (x) + C

Ejemplos 2.1.1 Resolver la ecuación diferencial dada:

1. y 0 = 8x − 7
Solución:
dy
= 8x − 7 =⇒ dy = (8x − 7) dx
dx
Z Z
dy = (8x − 7) dx

8x2
y= − 7x + C
2
y = 4x2 − 7x + C


dy
2. = e−3x + e6x
dx
Solución:

dy = e−3x + e6x dx

Z Z
e−3x + e6x dx


dy =

1 1
y = − e−3x + e6x + C
3 6


11
2.2. ECUACIÓN DIFERENCIAL DE VARIABLES SEPARABLES

2.2. Ecuación diferencial de variables separables

dy
Una ecuación diferencial de primer orden de la forma: = f (x, y) es de variables
dx
separables si existen dos funciones g (x) y h (y), tales que

f (x, y) = g (x) · h (y)

Es decir si f (x, y) puede descomponerse en dos factores, uno en x y otro en y. En tal

caso, la ecuación que resulta se resuelve del modo siguiente:

dy
= g (x) · h (y)
dx
dy
=⇒ = g (x) dx
h (y)
Z Z
dy
=⇒ = g (x) dx
h (y)
=⇒ H(y) = G(x) + C

Restricción: Al separar las variables, se debe tener presente que los diferenciales no
pueden quedar en el denominador.

Ejemplos 2.2.1 Resolver la ecuación diferencial dada:

dy
1. = 5xy + 3y
dx
Solución:
dy
= y (5x + 3)
dx
Z Z
dy
= (5x + 3) dx
y
5x2
ln y = + 3x + C1
2
5 2
y = e( 2 x +3x)+C1
5 2
y = e( 2 x +3x) .eC1

Pero eC1 = C (constante). Luego:

12
2.2. ECUACIÓN DIFERENCIAL DE VARIABLES SEPARABLES

5 2
y = Ce( 2 x +3x)


dy p
2. = x 1 − y2
dx
Solución:
dy
p = xdx
1 − y2

Z Z
dy
p = xdx
1 − y2
x2
sen−1 y = +C
2
 2 
x
y = sen +C
2


dy xy + 3x − y − 3
3. =
dx xy − 2x + 4y − 8
Solución:
dy (y + 3) (x − 1)
=
dx (y − 2) (x + 4)
y−2 x−1
Z Z
dy = dx
y+3 x+4
Z   Z  
5 5
1− dy = 1− dx
y+3 x+4

y − 5 ln |y + 3| = x − 5 ln |x + 4| + C

13
2.2. ECUACIÓN DIFERENCIAL DE VARIABLES SEPARABLES

2.2.1. Problemas con valores iniciales

Una ecuación diferencial de primer orden de la forma y 0 = f (x, y) ; y(x0 ) = y0 se

denomina un problema con valores iniciales. En este caso, la solución es una solución par-

ticular, aquella que satisface la condición inicial.

Ejemplos 2.2.2 Resolver los siguientes problemas con valores iniciales:

y
1. y0 = ; y(−3) = 4
x
Solución:
dy y dy dx
= =⇒ =
dx x y x

Z Z
dy dx
=
y x

ln |y| = ln |x| + ln |c|

eln|y| = eln|cx|

Solución general:

y = cx

Aplicando la condición inicial y(−3) = 4

Sustituyendo:

4
4 = c · (−3) =⇒ c = −
3
Entonces, la solución particular es:

4
y=− x
3

1
2. y0 = √ ; y(0) = 0
1 − x2

14
2.2. ECUACIÓN DIFERENCIAL DE VARIABLES SEPARABLES

Solución:
dy 1
=√
dx 1 − x2
Z Z
1
dy = √ dx
1 − x2
∴ y = sen−1 x + c Solución general

Además: x = 0; y = 0

Sustituyendo:

0 = sen−1 (0) + c

c=0

Entonces:

y = sen−1 x Solución particular 

3. y 0 = (3x − 1)5 ; y(2) = 4

Solución:
dy
= (3x − 1)5
dx
Z Z
dy = (3x − 1)5 dx

1
∴ y= (3x − 1)6 + c Solución general
18
Además: x = 2; y = 4

Sustituyendo:

1 1
4= (3(2) − 1)6 + c =⇒ 4 = (5)6 + c
18 18
15, 625
∴ c=−
18

15
2.2. ECUACIÓN DIFERENCIAL DE VARIABLES SEPARABLES

Entonces:

1 15, 625
y= (3x − 1)6 − Solución particular
18 18


4. y 0 = 3x2 (y 2 + 1) ; y(0) = 1

Solución:
dy dy
= 3x2 y 2 + 1 ⇒ 2 = 3x2 dx


dx y +1
Z Z
dy
= 3x2 dx
y2 + 1

y = tan x3 + c


Solución general

y(0) = 1

Sustituyendo:

1 = tan(c)

π
∴ c=
4
 π
y = tan x3 + Solución particular
4


2.2.2. Teorema de existencia y unicidad

Sea y 0 = f (x, y) una ED de primer orden. El teorema de existencia y unicidad establece
∂f (x, y)
que si f (x, y) y son funciones continuas sobre un intervalo I = [a, b] entonces
∂y
existe una y solo una solución de la ecuación que satisface la condición inicial y(x0 ) = y0 .

Es decir, existe una sola solución que pasa por el punto p(x0 , y0 ).

16
2.3. FUNCIONES HOMOGÉNEAS

2.2.3. Interpretación geométrica

Desde el punto de vista geométrico, la ecuación diferencial de primer orden y 0 = f (x, y),
expresa la pendiente de las grácas para las soluciones, en cada punto del plano donde

ellas están denidas. Es decir, la ecuación diferencial de primer orden representa un campo

de direcciones (ver sección 2.1.1, pág.35, Dennis Zill).

2.3. Funciones homogéneas

Si una función f (x, y) tiene la propiedad

f (tx, ty) = tα f (x, y)

para algún número real α, entonces se dice que f (x, y) es una función homogénea de grado
α.
Ejemplos 2.3.1 Verique si cada función es o no homogénea:
1. f (x, y) = x4 + y 4

Solución:

f (tx, ty) = (tx)4 + (ty)4

= t4 x4 + t4 y 4

= t4 x4 + y 4

= t4 f (x, y)

Por tanto, f (x, y) es homogénea de 4to grado. 

2. g(x, y) = x4 + y 4 − 3

Solución:

g(tx, ty) = (tx)4 + (ty)4 − 3

= t4 x4 + t4 y 4 − 3

Por tanto, g(x, y) no es homogénea. 

17
2.4. ECUACIÓN HOMOGÉNEA DE PRIMER ORDEN

2.4. Ecuación homogénea de primer orden

La ecuación diferencial de primer orden y 0 = f (x, y) es homogénea si f (x, y) es una

función homogénea de grado cero, es decir si

f (tx, ty) = f (x, y)

dy 3x2 + 2y 2
Ejemplos 2.4.1 Verique si la ecuación diferencial = es o no homogé-
dx 5xy
nea
Solución:
3x2 + 2y 2
f (x, y) =
5xy
3(tx)2 + 2(ty)2 3t2 x2 + 2t2 y 2
f (tx, ty) = =
5(tx)(ty) 5t2 xy

3x2 + 2y 2
= = f (x, y)
5xy

f (x, y) es una función homogénea de grado cero, por tanto, la ecuación diferencial es

homogénea. 

Notas:
1. Si la función f (x, y) es homogénea de grado cero, entonces se cumple que: f (x, y) =
f (y/x)

Así, para la función f (x, y) del ejemplo anterior:

3x2 2y 2
2
3x + 2y 2
2
+ 2
f (x, y) = = x x
5xy 5xy
x2
3 + 2 (y/x)2
= = f (y/x)
5y/x


2. Una ED de primer orden escrita en la forma diferencial M (x, y)dx + N (x, y)dy = 0
se dice que es homogénea si ambas funciones M y N son funciones homogéneas del mismo
grado.

18
2.5. SOLUCIÓN ED HOMOGÉNEA DE PRIMER ORDEN

2.5. Solución ED homogénea de primer orden

Supongamos que la ecuación diferencial y 0 = f (x, y) es homogénea. En tal caso, la

podemos expresar por:

dy
= f (y/x)
dx
Efectuando la sustitución:  y
 v =


 x
 y = vx


dy = vdx + xdv


vdx + xdv
= f (v)
dx
vdx + xdv = f (v)dx

xdv = (f (v) − v) dx
dv dx
= Ecuación de variables separables
f (v) − v x

Ejemplos 2.5.1 Determine si cada ecuación diferencial es o no homogénea, si lo es

resolverla:

1. (x2 + y 2 ) dx + (x2 − xy) dy = 0

Solución:
dy x2 + y 2
=− 2
dx x − xy

Las funciones coecientes M (x, y) = x2 + y 2 y N (x, y) = x2 − xy son homogéneas

de grado dos, por tanto, la ecuación diferencial es homogénea.

 y
 v =
=⇒ y = vx
x
dy = vdx + xdv


19
2.5. SOLUCIÓN ED HOMOGÉNEA DE PRIMER ORDEN

Sustituyendo:

 2
x + (vx)2

vdx + xdv
=− 2
dx x − x(vx)

1 + v2
 
vdx + xdv = dx
v−1
1 + v2
 
xdv = − v dx
v−1
 
v+1
xdv = dx
v−1
v−1
Z Z
dx
dv =
v+1 x
Z   Z
2 dx
1− dv =
v+1 x

v − 2 ln |v + 1| = ln |x| + ln |c|

v = ln (v + 1)2 cx
 

ev = (v + 1)2 cx

Sustituyendo:

y 2
ey/x = +1 cx Solución general
x


2. (x + y) y 0 = (x − y)

Solución:
dy x−y
=
dx x+y

por tanto, es homogénea de grado cero.

 y
 v = =⇒ y = vx
x
dy = vdx + xdv


20
2.5. SOLUCIÓN ED HOMOGÉNEA DE PRIMER ORDEN

Sustituyendo:

 
vdx + xdv x − vx
=
dx x + vx
   
vdx + xdv 1−v 1−v
= =⇒ vdx + xdv = dx
dx 1+v 1+v
 
1−v
xdv = − v dx
1+v
 2 
v + 2v − 1
xdv = − dx
1+v
Z Z
1+v dx
2
dv = −
v + 2v − 1 x
1
ln |v 2 + 2v − 1| = − ln |x| + ln |c1 |
2

c
v 2 + 2v − 1 =
x2
 y 2 y c
+2 −1= 2
x x x
y 2 2y c
+ − 1 =
x2 x x2

Multiplicando por x2 :

y 2 + 2xy − x2 = c Solución general

Observación 2.5.1 En general, una ED de primer orden y primer grado puede expre-
sarse como:
dy a1 x + a2 y + a3
= (1)
dx b1 x + b2 y + b3
a a a
a) Cuando 1 = 2 = 3 = K , la ecuación (1) se resolverá por transposición de
b1 b2 b3
términos.
a a
b) Cuando 1 = 2 = K (Siendo a3 6= 0 o b3 6= 0 ), entonces la ecuación (1) se
b1 b2
resolverá mediante sustitución de suma de variables.
a1 a2 a a
c) Si sucede que a3 = b3 = 0 y además 6= (Los cocientes 1 y 2 no son
b1 b2 b1 b2

21
2.5. SOLUCIÓN ED HOMOGÉNEA DE PRIMER ORDEN

proporcionales), entonces la ecuación (1) es homogénea.

a a a
d) Si sucede que los cocientes 1 , 2 y 3 no son proporcionales (pudiendo ser a3 = 0
b1 b2 b3
o b3 = 0, pero no simultáneamente), entonces la ecuación (1) puede ser reducida a
homogénea.

2.5.1. Procedimiento para convertir una ED no homogénea a ho-

mogénea
Supongamos que la ecuación diferencial (1) no es homogénea. El procedimiento para

convertirla en homogénea consiste en realizar una traslación del origen de coordenadas,

para ello realizaremos la siguiente sustitución:

y
t x = u + h; dx = du

y = t + k; dy = dt

Donde (h, k) son constantes a determinar. Sustituyendo en la ecuación (1):

dt a1 (u + h) + a2 (t + k) + a3
=
du b1 (u + h) + b2 (t + k) + b3

dt a1 u + a2 t + (a1 h + a2 k + a3 )
=
du b1 u + b2 t + (b1 h + b2 k + b3 )

Para determinar h y k, resolver el sistema de ecuaciones lineales:


a h + a k + a = 0
1 2 3
b 1 h + b 2 k + b 3 = 0

Luego,

dt a1 u + a2 t
=
du b 1 u + b2 t

22
2.5. SOLUCIÓN ED HOMOGÉNEA DE PRIMER ORDEN

La ecuación resultante es homogénea. Para resolverla efectuamos la sustitución:

 t
 v =
u



 t = vu



dt = vdu + udv

Ejemplo 2.5.1.1 Resolver la ED no homogénea (x + y − 1) dx + (x − y − 3) dy = 0

Solución:

x = u + h =⇒ dx = du

y = t + k =⇒ dy = dt

Sustituyendo en el ED:

[(u + h) + (t + k) − 1] du + [(u + h) − (t + k) − 3] dt = 0

[u + t + (h + k − 1)] du + [u − t + (h − k − 3)] dt = 0

[u − t + (h − k − 3)] dt = − [u + t + (h + k − 1)] du

dt u + t + (h + k − 1)
=−
du u − t + (h − k − 3)

h−k−3 = 0
h+k−1 = 0
2h − 4 = 0

h=2
k = −h + 1 = −2 + 1 = −1

Sustituyendo en la ED:

dt u+t dt u+t
=− =⇒ = Ecuación homogénea
du u−t du t−u

23
2.5. SOLUCIÓN ED HOMOGÉNEA DE PRIMER ORDEN

Haciendo:

t
= v =⇒ t = uv =⇒ dt = udv + vdu
u
 
u+t
⇒ dt = du
t−u
 
u + uv
udv + vdu = du
uv − u
 
1+v
udv + vdu = du
v−1
1 + 2v − v 2
 
udv = du
v−1

Separando variables:
 
v−1 du
dv =
1 + 2v − v 2 u

Integrando:

v−1
Z Z
du
dv =
1 + 2v − v 2 u
1
− ln |1 + 2v − v 2 | = ln |u| + ln |c|
2
u2 + 2tu − t2 = c

u=x−h⇒u=x−2

t=y−k ⇒t=y+1

(x − 2)2 + 2 (y + 1) (x − 2) − (y − 1)2 = c Solución general

24
2.6. DIFERENCIAL TOTAL

2.6. Diferencial total

1. Función de una variable

y = f (x)

⇒ dy = f 0 (x)dx

2. Función de dos variable

φ = f (x, y)
∂f ∂f
⇒ df (x, y) = dx + dy
∂x ∂y

3. Función de tres variable

φ = f (x, y, z)
∂f ∂f ∂f
⇒ df (x, y, z) = dx + dy + dz
∂x ∂y ∂z

2.7. Ecuación diferencial exacta del primer orden

La ecuación diferencial M (x, y)dx + N (x, y)dy = 0 es exacta, si existe una función

f (x, y) tal que:

df (x, y) = M dx + N dy

Esto es si:
∂f ∂f
M= ; N=
∂x ∂y

25
2.7. ECUACIÓN DIFERENCIAL EXACTA DEL PRIMER ORDEN

2.7.1. Criterio de exactitud

Teorema 2.7.1.1 La ecuación diferencial M (x, y)dx + N (x, y)dy = 0 es exacta, si:
∂f ∂f
M= ; N=
∂x ∂y

Demostración: Si M dx + N dy = 0 es exacta, entonces:

∂f ∂f
M= ; N=
∂x ∂y

Ejemplo 2.7.1.1 Resolver las siguientes ED

 
x y
1) (cos x + ln y)dx = + e dy
y
Solución: Primer método:
∂M 1
M = cos x + ln y ⇒ =
∂y y
x ∂N 1
N= + ey ⇒ =
y ∂x y
∂M ∂N
= . La ED es exacta. Luego,
∂y ∂x

Z Z
f (x, y) = M dx + h(y) = (cos x + ln y)dx + h(y)

= sen x + x ln y + h(y) (∗)

Para determinar la función h(y), partimos de la condición:

∂f x x
= N ⇒ + h0 (y) = + ey ⇒ h0 (y) = ey
∂y y y

h(y) = ey .

Sustituyendo en (∗):

f (x, y) = sen x + x ln y + ey = c Solución general

26
2.7. ECUACIÓN DIFERENCIAL EXACTA DEL PRIMER ORDEN

Segundo método:

Z Z  
x y
f (x, y) = N dy + h(x) = + e dy + h(x)
y

= x ln y + ey + h(x) (∗)

Para determinar la función h(x), partimos de la condición:

∂f
= M ⇒ ln y + h0 (x) = cos x + ln y ⇒ h0 (y) = cos x
∂x
h(y) = sen x.

Sustituyendo en (∗):

f (x, y) = sen x + x ln y + ey = c Solución general


ydx − xdy 1
2) + dy = 0
(x + y)2 y
Solución:
 
y 1 x
dx + − dy = 0
(x + y)2 y (x + y)2

y ∂M (x + y)2 (1) − 2y(x + y) x−y

M= 2 ⇒ = 4
=
(x + y) ∂y (x + y) (x + y)3
1 x ∂N (x + y)2 (1) − 2x(x + y) x−y
N= − 2
⇒ = − 4
=
y (x + y) ∂x (x + y) (x + y)3
∂M ∂N
= . La ED es exacta. Así,
∂y ∂x

Z Z
y
f (x, y) = M dx + h(y) = dx + h(y)
(x + y)2
y
=− + h(y) (∗)
x+y

27
2.7. ECUACIÓN DIFERENCIAL EXACTA DEL PRIMER ORDEN

∂f (x + y)(1) − y(1) 0 1 x
=N ⇒− + h (y) = −
∂y (x + y)2 y (x + y)2
x 1 x 1
⇒− 2
+ h0 (y) = − 2
⇒ h0 (y) =
(x + y) y (x + y) y

h(y) = ln y

Sustituyendo en (∗):
y
f (x, y) = ln y − =c Solución general
x+y

3) [y cos(xy) − sen x] + x cos(xy)dy = 0

Solución:
∂M
M = y cos(xy) − sen x ⇒ = cos(xy) − xy sen(xy)
∂y
∂N
N = x cos(xy) ⇒ = cos(xy) − xy sen(xy)
∂x
∂M ∂N
= . La ED es exacta. Así,
∂y ∂x

Z Z
f (x, y) = N dy + h(x) = x cos(xy) dx + h(x)

= sen(xy) + h(x) (∗)

∂f
= N ⇒ y cos(xy) + h0 (x) = y cos(xy) − sen x ⇒ h0 (x) = sen x
∂y

h(x) = cos x

Sustituyendo en (∗):

f (x, y) = sen(xy) + cos x = c Solución general

28
2.7. ECUACIÓN DIFERENCIAL EXACTA DEL PRIMER ORDEN

2.7.2. Factor integrante

Supongamos que la ecuación diferencial M dx + N dy = 0 no es exacta. una función
µ(x, y) es un factor integrante de M dx + N dy = 0 si µM dx + µN dy = 0 es una ecuación
diferencial exacta.

2.7.3. Búsqueda de un factor integrante

En algunas ocasiones podemos obtener un factor integrante ajustando la ecuación a

fórmulas de integración conocidas, pero en general, no existe un método para la obtención

de un factor integrante, excepto en el caso de que este factor dependa de una sola de las

variables.

2.7.4. Factor integrante que depende de una de las variables

Supongamos que

M dx + N dy = 0 (∗)

es una ED no exacta y µ(x, y) un factor integrante de (∗). Entonces la ecuación:

µM dx + µN dy = 0 (∗∗)

es exacta, esto es:

∂ ∂
(µM ) = (µN )
∂y ∂x

Por lo que, µ debe satisfacer la EDP

∂M ∂µ ∂N ∂µ
µ +M = µ +N
∂y ∂y ∂x ∂x
!
∂M ∂N ∂µ ∂µ
µ − = N −M
∂y ∂x ∂x ∂y

Supongamos que µ depende exclusivamente de la variable x (No en absoluto de la

29
2.7. ECUACIÓN DIFERENCIAL EXACTA DEL PRIMER ORDEN

∂µ ∂µ dµ
variable y ), entonces =0 y se convierte en . Así,
∂y ∂x dx
!
∂M ∂N dµ
µ − = N
∂y ∂x dx

∂M ∂N
−
1 dµ ∂y ∂x
=
µ dx N

∂M ∂N
−
∂y ∂x d ln µ
Si f (x) = ⇒ = f (x). Por tanto,
N dx
Z 
µ(x) = exp f (x) dx factor integrante

De igual manera, si µ depende solo de y, entonces el factor integrante es

 Z 
µ(y) = exp − g(y) dx

∂M ∂N
−
∂y ∂x
donde g(y) =
M

Ejemplo 2.7.4.1 Resolver las siguientes ED:

1) (x2 + y 2 + 1)dx + x(x − 2y)dy = 0

Solución:
∂M
M = x2 + y 2 + 1 ⇒ = 2y
∂y
∂N
N = x2 − 2xy ⇒ = 2x − 2y
∂x

∂M ∂N
6= . La ED no es exacta. Determinar un factor integrante.
∂y ∂x

∂M ∂N
−
∂y ∂x −2x + 4y 2
= =−
N 2x − 2y x

30
2.7. ECUACIÓN DIFERENCIAL EXACTA DEL PRIMER ORDEN

2
f (x) = −
x
factor integrante:

Z 
2
µ(x) = exp − dx = x−2
x

Multiplicando la ED por el factor integrante µ(x) = x−2 :

(x−2 + x−2 y 2 + 1)dx + (1 − 2x−1 y)dy = 0

∂M
M = x−2 + x−2 y 2 + 1 ⇒ = 2x−2 y
∂y
∂N
N = 1 − 2x−1 y ⇒ = 2x−2 y
∂x
∂M ∂N
= . La ED es exacta.
∂y ∂x

Z
f (x, y) = M dx + h(y)

= x − x−1 y 2 − x−1 + h(y) (∗)

∂f
= N ⇒ −2x−1 y + h0 (y) = 1 − 2x−1 y ⇒ h(y) = y
∂y

Sustituyendo en (∗) y multiplicando por x:

x2 + y 2 + xy − 1 = cx Solución general

2) y(2x − y + 1)dx + x(3x − 4y + 3)dy = 0

31
2.7. ECUACIÓN DIFERENCIAL EXACTA DEL PRIMER ORDEN

Solución:
∂M
M = 2xy − y 2 + y ⇒ = 2x − 2y + 1
∂y
∂N
N = 3x2 − 4xy + 3x ⇒ = 6x − 4y + 3
∂x
∂M ∂N
6= . La ED no es exacta. Determinar un factor integrante.
∂y ∂x

∂M ∂N
−
∂y ∂x 4x − 2y + 2 2
=− 2
= − = g(y)
M 2xy − y + y y

factor integrante:

 Z  Z 
2 2
µ(y) = exp − g(y) dy = = exp dy = y
y

Multiplicando la ED por el factor integrante µ(y) = y 2 :

y 3 (2x − y + 1)dx + xy 2 (3x − 4y + 3)dy = 0

∂M
M = 2xy 3 − y 4 + y 3 ⇒ = 6xy 2 − 4y 3 + 3y 2
∂y
∂N
N = 3x2 y 2 − 4xy 3 + 3xy 2 ⇒ = 6xy 2 − 4y 3 + 3y 2
∂x
∂M ∂N
= . La ED es exacta.
∂y ∂x

Z
f (x, y) = M dx + h(y)

= x2 y 3 − xy 4 + xy 3 + h(y) (∗)

∂f
= N ⇒ 3x2 y 2 − 4xy 3 + 3xy 2 + h0 (y) = 3x2 y 2 − 4xy 3 + 3xy 2 ⇒ h(y) = c
∂y

32
2.8. ECUACIÓN DIFERENCIAL LINEAL DE PRIMER ORDEN

Sustituyendo en (∗):

xy 3 (x − y + 1) = c Solución general

2.8. Ecuación diferencial lineal de primer orden

Una ED que se pueda escribir en la forma

dy
a(x) + b(x)y = r(x)
dx
donde a(x), b(x) y r(x) son funciones denidas en el intervalo (α, β), es llamada una

ecuación diferencial de primer orden.

Supongamos que a(x) 6= 0 en (α, β), entonces obtenemos la forma estandar:

dy
+ p(x)y = Q(x) (∗)
dx

b(x) r(x)
donde: p(x) = , Q(x) = .
a(x) a(x)

p(x) y Q(x) son funciones continuas en el intervalo (α, β).

Por el momento, supongamos que existe un factor integrante positivo µ > 0, que solo

depende de x, para (∗)

 
dy
µ(x) + p(x)y = µ(x)Q(x) (∗∗)
dx

debe ser una ecuación diferencial exacta.

dy
µ(x) + p(x)µ(x)y = µ(x)Q(x)
dx
Esta ecuación se puede escribir en la forma:

M dx + N dy = 0

33
2.8. ECUACIÓN DIFERENCIAL LINEAL DE PRIMER ORDEN

donde: M = µpy − µQ; N =µ

p, Q, y µ son funciones solo de x. Si (∗) es exacta, entonces:

∂M ∂N
=
∂y ∂x

Por lo que, µ satisface la ecuación:

dµ
µp =
dx
Z 
µ(x) = exp p(x) dx factor integrante

Multiplicando (∗), por µ(x)

p(x) dx dy
R R R
p(x) dx p(x) dx
e +e p(x)y = e Q(x)
dx
d h R p(x) dx i R
e y = e p(x) dx Q(x)
dx
Integrando:

R
Z R
p(x) dx p(x) dx
e y= e Q(x) dx

R
Z R
− p(x) dx p(x) dx
y=e e Q(x) dx Ecuación de Leibniz

Ejemplos 2.8.1 Resolver las siguientes ED:

2x
1) y0 + y = 9x2
x ln x
Solución:
2x
p(x) = ; Q(x) = 9x2
x ln x

Factor integrante:

Z  Z 
2x
µ(x) = exp p(x) dx = exp dx = ln x
x ln x

34
2.8. ECUACIÓN DIFERENCIAL LINEAL DE PRIMER ORDEN

Multiplicando la ED por el factor integrante µ(x) = ln x:

2x
ln xy 0 + y = 9x2 ln x
x

d
[y ln x] = 9x2 ln x
dx
Z
1
y= 9x2 ln x dx
ln x
x3 (3 ln x − 1) + c
y= Solución general
ln x

dy 2x 4
2) + y=
dx 1 + x 2
(1 + x2 )2
Solución:
2x 4
p(x) = ; Q(x) =
1+x 2
(1 + x2 )2

Factor integrante:

Z 
2x
µ(x) = exp dx = 1 + x2
1 + x2

Multiplicando la ED por el factor integrante µ(x) = 1 + x2 :

dy 4
1 + x2 + 2xy =
dx 1 + x2

d  4
y 1 + x2 =


dx 1 + x2
Z
1 4
y= dx
(1 + x2 ) 1 + x2
4 tan−1 x 2 −1


y= + c 1 + x Solución general
1 + x2


35
2.9. ECUACIONES REDUCIBLES A LA FORMA LINEAL

ECUACIÓN DE BERNOULLI

2.9. Ecuaciones reducibles a la forma lineal

Ecuación de Bernoulli
Una ED que se pueda escribir en la forma

dy
+ p(x)y = Q(x)y n (∗)
dx
donde n 6= 0, n 6= 1, p(x) y Q(x) son funciones continuas [a, b], es llamada ecuación de

Bernoulli.

La ecuación (∗) se puede reducir a una ecuación diferencial lineal como sigue. Multi-

plicando (∗) por y −n , tenemos:

dy
y −n + p(x)y 1−n = Q(x) (∗∗)
dx
Haciendo el cambio de variable

u = y 1−n

Implica:

du dy
= (1 − n) y −n
dx dx
dy 1 du
y −n =
dx 1 − n dx

Sustituyendo en (∗∗):

1 du
+ p(x)u = Q(x)
1 − n dx
du
+ (1 − n)p(x)u = (1 − n)Q(x) Ecuación lineal en u
dx
Ejemplos 2.9.1 Resolver las siguientes ED:
dy 2
1) + y = 6x4 y 2
dx x
Solución: Multiplicando la ED por y−2
dy 2 −1
y −2 + y = 6x4
dx x

36
2.9. ECUACIONES REDUCIBLES A LA FORMA LINEAL

ECUACIÓN DE BERNOULLI

Haciendo el cambio de variable u = y −1

du 2
− u = −6x4 Ecuación lineal en u
dx x

2
p(x) = − ; Q(x) = −6x4
x

Factor integrante:

Z 
2
µ(x) = exp − dx = x−2
x

Multiplicando la ED por el factor integrante µ(x) = x−2 :

du
x−2 − 2x−3 u = −6x2
dx

d  −2 
x u = −6x2
dx
Z
−2
x u = −6x2 dx = −2x3 + c

u = −2x5 + cx2

Sustituyendo u = y −1

1
y= Solución general
x2 (c − 2x3 )


dy 3
2) − y = 6x2 y 1/3 ln x
dx 2x
2
Solución: Multiplicando la ED por y−1/3
3
2 −1/3 dy 1
y − y 2/3 = 4x2 ln x
3 dx x

Haciendo el cambio de variable u = y 2/3

du 1
− u = 4x2 ln x Ecuación lineal en u
dx x

37
2.10. DIFERENCIALES NOTABLES EXACTAS: COMBINACIONES INTEGRABLES

1
p(x) = − ; Q(x) = 4x2 ln x
x

Factor integrante:

Z 
1
µ(x) = exp − dx = x−1
x

Multiplicando la ED por el factor integrante µ(x) = x−1 :

du
x−1 − x−2 u = 4x ln x
dx

d  −1 
x u = 4x ln x
dx
Z
−1
x u = 4x ln x dx

u = x3 (2 ln x − 1) + cx

Sustituyendo u = y 2/3

y 2/3 = x3 (2 ln x − 1) + cx Solución general

2.10. Diferenciales notables exactas: Combinaciones in-

tegrables
1) d (xy) = xdy + ydx
 
x ydx − xdy
2) d =
y y2
 y  xdy − ydx
3) d =
x x2
 2
x y(2x)dx − x2 dy
4) d =
y y2
5) d (xp y q ) = xp d (y q ) + y q d (xp ) = xp−1 y q−1 (qxdy + pydx)

38
2.10. DIFERENCIALES NOTABLES EXACTAS: COMBINACIONES INTEGRABLES

y
h  y i d xdy − ydx
6) d tan−1 = xy 2 =
x x2 + y 2
1+
x
Ejemplos 2.10.1 Halle la solución de la ED usando diferenciales notables por lo me-
nos en una parte de la ecuación:

1) 3x4 dx + ydx = −xdy

Solución:

3x4 dx = −d (xy)

Integrando

3 5
x = −xy + c Solución general
5


2) y (xy + 1) dx = xdy

Solución:
ydx − xdy
xdx = −
y2
 
x
xdx = −d
y

Integrando

1 2 x
x =− +c Solución general
2 y

3) ydx + [y (x2 + y 2 ) − x] dy = 0

Solución:
ydx − xdy
+ ydy = 0
x2 + y 2
    
−1 x 1 2
d tan +d y =0
y 2

39
2.11. SOLUCIONES DE ED POR SUSTITUCIÓN

Integrando

 
−1 x 1
tan + y2 = c Solución general
y 2

2.11. Soluciones de ED por sustitución

Normalmente resolvemos una ecuación diferencial reconociéndola dentro de una cierta

clase de ecuaciones (Separables, lineales, exactas, etc.) y después aplicamos un procedi-

miento, que consiste en pasos matemáticos especícos para el tipo de ecuación que nos

conducen a la solución de la misma.

Cuando tenemos una ecuación diferencial que no pertenece a una de las clases citadas,

para resolverla puede ser de utilidad transformarla en otra ecuación diferencial mediante

una sustitución apropiada:

y
a) Homogénea: v=
x
b) Bernoulli: u = y 1−n
dy
c) Reducción a separación de variables: Una ED de la forma = f (Ax+By+C), B 6=
dx
0, se puede reducir siempre a una ecuación con variables separables por medio de la

sustitución z = Ax + By + C .

Ejemplos 2.11.1 Resolver la ED dada, usando una sustitución adecuada:

dy
1) = 1 + ey−x+5
dx
Solución: Sea z = y − x + 5, derivando respecto a x y sustituyendo en la ED:
dz
= ez ED separable
dx

Separando variables e integrando:

−e−z = x + c

40
2.12. APLICACIÓN DE LA ED DE PRIMER ORDEN

Sustituyendo z =y−x+5

x + ex−y−5 = c Solución general


dy
2) = sen (x + y)
dx
Solución: Sea z = x + y, derivando respecto a x y sustituyendo en la ED:
dz
= sen z + 1 ED separable
dx

Separando variables e integrando:

tan z − sec z = x + c

Sustituyendo z =x+y

tan(x + y) + sec(x + y) = x + c Solución general

2.12. Aplicación de la ED de primer orden

2.12.1. Modelo de crecimiento de una población (Analizado en la
unidad I)
Ejemplo 2.12.1.1 Si la población de un pais se duplica en 10 años, ¾En cuántos años
se habrá triplicado, suponiendo que la velocidad de aumento de la población es proporcional
al número de habitantes?

Solución: Sean y = población y t =tiempo. La velocidad de aumento de la población

es proporcional al número de habitantes, así:

 dy = ky

dt
y(0) = y , y(10) = 2y
0 0

41
2.12. APLICACIÓN DE LA ED DE PRIMER ORDEN

Separando variables e integrando:

ln |y| = kt

y = cekt (∗)

Al aplicar las condiciones y(0) = y0 y y(10) = 2y0 para determinar las constantes c y
k se tiene que c = y0 y k = ln 2/10. Sustituyendo estos valores en (∗), se obtiene la ED
que modela el aumento de la población en un tiempo dado.

y = y0 e(ln 2/10)t

Sea t3 el tiempo en triplicarse la población, entonces y (t3 ) = 3y0 . Luego,

3y0 = y0 e(ln 2/10)t3 . Despejando t3 :

10 ln 2
t3 = ≈ 15,85 años.
ln 2


2.12.2. Cambio de temperatura de un cuerpo (Ley de enfriamien-

to de Newton)
Tratemos de descubrir la ley según la cual cambia la temperatura de un cuerpo inmerso

en un medio que mantiene su temperatura constante:

Según la experiencia, sabemos que si el cuerpo está más caliente que el medio, su tem-

peratura disminuye; si el cuerpo está más frío que el medio, su temperatura aumenta; si

el cuerpo tiene igual temperatura que el medio, la mantiene constante.

Hipótesis: Supongamos que el cambio de temperatura está en relación directa a la

diferencia de temperatura entre el cuerpo y el medio, ¾Qué ecuación diferencial podrá

modelar este fenómeno?

42
2.12. APLICACIÓN DE LA ED DE PRIMER ORDEN

Sean:

T = Temperatura del cuerpo en cada instante.

Tm = Temperatura del medio (Constante).

t = Tiempo
dT
= Velocidad de cambio de la temperatura.
dt

dT
⇒ = k (T − Tm ) , k < 0 (∗)
dt
Notas:
dT
1. Si T < Tm , entonces es positiva.
dt
dT
2. Si T > Tm , entonces es negativa.
dt
dT
3. Si T = Tm (solución de equilibrio), entonces es nula.
dt
La solución general de la ED (∗) está dado por:

T (t) = Tm + cekt

Ejemplo 2.12.2.1 Suponga que un cuerpo se introduce en un medio, con una tem-
peratura inicial de 100◦ C y 10 minutos después su temperatura se reduce a 80◦ C . ¾Qué
temperatura tendrá el cuerpo después de 30 minutos? Asumir que la temperatura del medio
es Tm = 27◦ C .
Solución:

T (t) = Tm + cekt

T (t) = 27 + cekt

Condiciones iniciales:

T (0) = 100◦ C

T (10) = 80◦ C

T (30) =?

43
2.12. APLICACIÓN DE LA ED DE PRIMER ORDEN

Para determinar la constante c, aplicamos la condición T (0) = 100◦ C , así c = 73.

ln |53/73|
Luego, aplicando la condición T (10) = 80◦ C , se tiene que k = = −0,032.
10
Sustituyendo en la ED:

T (t) = 27 + 73e−0,032t

Para determinar la temperatura después de 30 minutos, evaluamos la ED anterior en

◦
t = 30, así T (30) = 54,94 C . 

2.12.3. Trayectorias ortogonales

Una trayectoria ortogonal de una familia de curvas es una curva que corta cada uno

de los miembros de una familia de manera ortogonal, es decir, formando ángulos rectos.

Figura 2.1: Trayectorias ortogonales

Por ejemplo, cada miembro de la familia y = mx de rectas que pasan por el origen es
una trayectoria ortogonal de la familia x2 + y 2 = r2 de círculos concéntricos que tienen
centro en el origen. Decimos que ambas familia son trayectorias ortogonales una de la otra.

44
2.12. APLICACIÓN DE LA ED DE PRIMER ORDEN

Ejemplo 2.12.3.1 Encuentre la trayectorias ortogonales de las familias de curvas x =

ky , donde k es una constante arbitraria. Las curvas x = ky 2 forman una familia de
2

parábolas cuyo eje de simetría es el eje X .

Solución: El primer paso consiste en hallar una ED que satisfaga todos los miembros
de la familia.

x = ky 2 (a)

Al derivar respecto de x se obtiene:

dy
1 = 2ky (b)
dx
x
De (a) : k = . Sustituyendo en (b):
y2

dy y
=
dx 2x
Esto signica que la pendiente de la recta tangente en cada punto (x, y) de una de las
dy y
parábolas es = .
dx 2x

En una trayectoria ortogonal la pendiente de la recta tangente debe ser reciproca

negativa de esta pendiente. Por tanto, las trayectorias ortogonales deben satisfacer la ED:

dy 2x
=−
dx y

Al resolver esta ED se obtiene:

2 y2
x + = c Trayectoria ortogonal
2


y
Ejemplo 2.12.3.2 Hallar el tipo más general de curvas
tal que la normal a cualquier punto coincida con la recta que
Tangente

une el origen de coordenadas con dicho punto.

N ormal
dy
Solución: Sean m = pendiente; mT = ; mN = y

dx
α
x
x
45
2.12. APLICACIÓN DE LA ED DE PRIMER ORDEN

dx y
− ; tan α = m = , entonces
dy x

mN = m
dx y
− = ⇒ −xdx = ydy
dy x

Integrando y multiplicando por 2:

−x2 = y 2 + c1

x2 + y 2 = c Solución general

Ejemplo 2.12.3.3 Encontrar curva que tiene las siguientes propiedades:

a) Su ordenada es igual a 5 veces el ln de su pendiente.
b) Cruza a la recta x = 10 con un angulo de α = 45◦ .
Solución:
 
dy
a) Sean y= ordenada; y = 5 ln , entonces
dx

dy
ey/5 = . Separando variables e integrando:
dx
x = −5e−y/5 + c Solución general
y

b) α = 45◦ ⇒ tan α = tan 45◦ = 1.

y
tan 45◦ = y
x
y
1 = ⇒ x = y.Si x = 10 ⇒ y = 10 α
x x
x

Al sustituir en la solución general, c = 10 + 5e−2 , así:

x = −5e−y/5 + 10 + 5e−2

46
2.12. APLICACIÓN DE LA ED DE PRIMER ORDEN

2.12.4. Circuitos eléctricos

Una aplicación importante de las ED lineales de primer orden (segundo orden) surge

del análisis de circuitos eléctricos simples, formados por una fuente de voltaje, un resistor

y un inductor o un condensador (capacitor).

Para los circuitos que consideraremos, el comportamiento de la corriente en el circuito

se rige por la ley de voltaje de Kichho :

La suma algebraica de los cambios instantáneos del potencial (caída de voltaje) en torno
de cualquier lazo cerrado debe anularse.

Caída de voltaje
a) Resistor:

ER = RI (1)

donde:

ER : Caída voltaje resistor.

I : Corriente. Unidad de medida amperio (A).

R : Constante de proporcionalidad llamada resistencia. Unidad medida Ohm(Ω).

Ecuación (1) recibe el nombre de ley de Ohm.

b) Inductor:
dI
EL = L (2)
dt
donde:

EL : Caída voltaje inductor.

L : Constante de proporcionalidad llamada inductancia. Unidad medida henrio(H ).

dI
: Razón de cambio instantánea de la corriente.
dt

47
2.12. APLICACIÓN DE LA ED DE PRIMER ORDEN

c) Capacitor o condensador:
1
Ec = q (2)
c
donde:

Ec : Caída voltaje Capacitor o condensador.

c : Constante llamada capacitancia. Unidad medida faradio (F ).

q : Carga eléctrica. Unidad de medida coulumb (C)

Nota: La corriente en el condensador es la razón de cambio de su carga:

dq
I = (∗)
dt
Fuente de voltaje:
a) Batería

b) Generador

Produce una fuerza electromotriz (fem), E, que impulsa a las cargas eléctricas y pro-

duce una corriente. Dependiendo de la naturaleza de la fuente E, puede ser:

Constante.

Función del tiempo.

Nota: Las dos formas más importantes para E(t) son:

1. E(t) = E0 , E0 constante. Corresponde a una fuerza de fem como la batería. La

corriente resultante es llamada corriente continua (DC).

2. E(t) = E0 sen ωt, E0 , ω constantes. La fem oscila entre ±E0 . La corriente resultante
en el circuito es llamada corriente alterna (AC).

2.12.5. Circuitos RL
R
I(t)
E L

48
2.12. APLICACIÓN DE LA ED DE PRIMER ORDEN

Para determinar la corriente I(t) en el circuito, supongamos que una fuente de voltaje

suma energía potencial al circuito. Sea E(t) el voltaje que proporciona al circuito en el

instante de tiempo t.

Aplicando la ley de voltaje de Kichho:

EL + ER = E(t) (∗∗)

Sustituyendo (1) y (2) en (∗∗)

dI
L + RI = E(t). Al dividir por L, se tiene:
dt
dI R E(t)
+ I= Ecuación diferencial lineal en I
dt L L
El factor integrante para esta ED lineal es:

Z 
R
µ(t) = exp = eRt/L
L

y su solución general es:

Z
−Rt/L E(t)
I(t) = e eRt/L dt
L

Nota: Por lo general, se da I(0) = 0 como condición inicial.

Ejemplo 2.12.5.1 Un circuito RL es controlado por un voltaje E(t) = 10 sen 4t V . Si
2
R = 2Ω, L = H . Determine la corriente en el circuito para t ≥ 0, así como el voltaje en
3
la resistencia y en el inductor. La corriente inicial es nula.
Solución: La corriente en el circuito esta dada por:
Z
−Rt/L E(t) R 2
I(t) = e eRt/L dt + c; donde: = = 3; E(t) = 10 sen 4t. Así:
L L 2
3
Z
−3t 10 sen 4t
I(t) = e e3t dt + c
2
3
Z
−3t
= 15e e3t sen 4t dt + c. Resolviendo esta integral se tiene:

3
I(t) = (3 sen 4t − 4 cos 4t) + ce−3t
5

49
2.12. APLICACIÓN DE LA ED DE PRIMER ORDEN

12
Aplicando la condición inicial I(0) = 0, se tiene que c= , entonces:
5
3
3 sen 4t − 4 cos 4t + 4e−3t


I(t) = Corriente circuito
5
Voltaje:
Resistencia:
6
3 sen 4t − 4 cos 4t + 4e−3t


ER (t) = RL =
5

Inductor:
dI
EL (t) = L . Derivando I(t):
dt
dI 12
3 cos 4t + 4 sen 4t − 3e−3t .


= Así, se tiene que:
dt 5
8
3 cos 4t + 4 sen 4t − 3e−3t


EL (t) =
5


2.12.6. Circuitos RC
R
I(t)
E C

Aplicando la ley de voltaje de Kichho:

E(t) − ER − EC = 0 ⇒ ER + EC = E(t) (∗ ∗ ∗)

Sustituyendo (1) y (2) en (∗ ∗ ∗)

1
RI + q = E(t)
c
La corriente en el capacitor es:

dq
I =
dt

50
2.12. APLICACIÓN DE LA ED DE PRIMER ORDEN

Luego,

dq 1
R + q = E(t). Al dividir por R, se tiene:
dt c
dq 1 E(t)
+ q= Ecuación diferencial lineal en q
dt cR R
Nota: La condición para un condensador es su carga q en t = 0.
Ejemplo 2.12.6.1 Determinar la corriente que uye en un circuito RC, si el capacitor
está inicialmente descargado y la fem de conducción está dada por E(t) = E0 e−at , donde
E0 y a son constantes.

Solución: La corriente en el circuito esta dada por:

dq
I(t) = (1)
dt

 dq + 1 q = E(t) ;

E(t) = E0 e−at
dt cR R
q(0) = 0


dq 1 E0 −at
+ q= e (∗)
dt cR R
Factor integrante:

Z 
1
µ(t) = exp = et/(cR)
cR

Multiplicando (∗) por el factor integrante:

d  t/(cR)  E0 (1/(cR)−a)t
e q = e
dt R
Z
d  t/(cR)  E0
e q = e(1/(cR)−a)t dt
dt R
cE0
q(t) = e−at + ke−t/(cR)
1 − acR

51
2.12. APLICACIÓN DE LA ED DE PRIMER ORDEN

cE0
Aplicando la condición inicial q(0) = 0, se tiene que k= , entonces:
1 − acR
!
cE0
q(t) = e−at − e−t/(cR) Carga en el capacitor
1 − acR

Derivando q respecto a t:
!
dq cE0 1 −t/(cR)
= e − ae−at
dt 1 − acR cR

Sustituyendo en (1):
!
cE0 1 −t/(cR)
I(t) = e − ae−at Corriente en el circuito
1 − acR cR

52