Profesor: Eric Padilla Mora.

Nociones elementales sobre ecuaciones diferenciales.

Muchas de las leyes de la fı́sica, la electrónica y la mecánica, entre otras áreas, se rigen por la igualdad
entre una función desconocida y sus derivadas, lo cual define a una ecuación diferencial. Por ejemplo, de
acuerdo con la ley de Kirchhoff, un circuito eléctrico simple que contiene una resistencia de R ohmios y un
inductor con inductancia de L henrios, en serie con una fuente de fuerza electromotriz (baterı́a o generador)
que proporciona un voltaje de E(t) voltios en el tiempo t segundos (llamado circuito RL), satisface que
dI R E(t)
L + RI = E(t) ⇒ I 0 + I = , donde I(t) es la corriente medida en amperios. Esta y otras relaciones
dt L L
que pueden expresarse mediante una o varias razones de cambio entre las variables implicadas, generan lo
que llamaremos una ecuación diferencial.

Definición: si en una ecuación, la incógnita es una función de una o varias variables y si en la ecuación
aparece una o más derivadas de la función (derivadas parciales si se trata de una función de varias variables),
entonces se dice que ella es una ecuación diferencial. Son ejemplos de ecuaciones diferenciales:
dI
1. + 3I = 3 ⇒ I 0 + 3I = 3. (Modelo para la corriente en un circuito RL).
dt
dy
2. = 3x + x2 + 1 ⇒ y 0 = 3x + x2 + 1
dx
d2 y dy
3. 2
+ (2 sen x) + y = ex ⇒ y 00 + (2 sen x)y 0 + y = ex
dx dx
4. y (4) + xy 0 + xy = x − 1

5. En 1798 el economista inglés Thomas Malthus modeló el crecimiento poblacional al suponer que la
rapidez con que crece una población en un paı́s en un cierto tiempo es proporcional a la misma población
del paı́s en ese momento. Si P = P (t) representa la población de ese paı́s en el instante de tiempo t, se
dP
tiene que = kP ⇒ P 0 = kP , con k constante.
dt
6. En el modelo básico un núcleo atómico consiste en una combinaciones de protones y neutrones. Muchas
sustancias poseen núcleos inestables; es decir, núcleos que se desintegran convirtiendo una sustancia en
otra. Tales sustancias inestables se llaman radioactivas. En el modelo matemático para esto se supone
que la rapidez con la cual se desintegran átomos radioactivos es proporcional a la cantidad de átomos
radioactivos presentes. Si M = M (t) representa la cantidad de sustancia radioactiva en el instante de
dM
tiempo t, se tiene que = kM ⇒ M 0 = kM , con k constante.
dt
7. La ley de Torricelli establece que la tasa de cambio con respecto al tiempo del volumen V de agua
dV p
en un tanque que se vacı́a a través de un orificio de área a, está dada por, = −a 2gy, con g la
dt
aceleración debida a la gravedad y y la altura del agua contenida en el recipiente en el tiempo t.

1
 8. Newton propuso una ley experimental sobre el enfriamiento o calentamiento de un cuerpo: “la rapidez
con la que cambia la temperatura de un cuerpo es proporcional a la diferencia de temperaturas entre
el cuerpo y el medio ambiente que lo rodea”. Si T = T (t) representa la temperatura del cuerpo en el
dT
instante t y Ta representa la temperatura del medio ambiente, entonces se satisface que = k(T − Ta )
dt
dy x+y
9. (x + y) dx − 4y dy = 0 ⇒ (x + y) dx = 4y dy ⇒ (x + y) = 4y ⇒ (x + y) = 4yy 0 ⇒ = y0
dx 4y
∂ 2u ∂ 2u
10. = ⇒ uxx = utt . (función u de dos variables)
∂x2 ∂t2
∂u ∂u
11. x +y = u. (función u de dos variables)
∂x ∂y

Ecuación diferencial ordinaria (EDO): es aquella ecuación diferencial en la cual sólo existe una variable
independiente y las derivadas que aparecen son derivadas ordinarias.
Ejemplos:
1. y 0 = 4x + 3
d3 y dy 2 2
2. 3
+ = ey + sen t ⇒ y 000 + y 0 = ey + sen t
dt dt
3. y 0 − (y 2 + 5)y 00 = 0

4. 4y 2 z 00 − 3yz 0 − 2z = 0
3
5. 3 (y 000 ) = −2y 0
dI
6. + 3I = 3
dt
7. y IV + y 000 + 4y 00 − sen x = cos x
dy
8. x2 dy + y dx = 0 ⇒ x2 dy = −y dx ⇒ x2 = −y ⇒ x2 y 0 + y = 0
dx

Ecuación diferencial parcial (EDP): es aquella ecuación diferencial en la cual hay dos o más variables
independientes y las derivadas que aparecen son derivadas parciales.
Ejemplos:
∂z ∂z
1. − = 0. (función z de dos variables)
∂x ∂t
∂ 2z ∂z
2. = . (función z de dos variables)
∂x2 ∂t
∂ 3v ∂ 3v
3. + = v − 2. (función v de dos variables)
∂x3 ∂y 3
∂ 2y ∂y
4. 2
= a2 . Con a constante. (Ecuación de onda)
∂t ∂x

2
Orden de una ecuación diferencial: esta dado por el orden de la derivada más alta que aparece en la
ecuación. Ejemplos:
1. y 0 = 4x + 3. Ecuación diferencial ordinaria (EDO) de primer orden.

2. 4y 2 z 00 − 3yz 0 − 2z = 0. EDO de segundo orden.

dI
3. + 3I = 3. EDO de primer orden.
dt
4. y IV + y 000 + 4y 00 − sen x = cos x. EDO de cuarto orden.
∂z ∂z
5. = . EDP de primer orden.
∂x ∂t
∂ 2z ∂z
6. 2
= . EDP de segundo orden.
∂x ∂t
∂ 4u ∂ 2u
7. a2 + 2 = 0. EDP de cuarto orden.
∂x4 ∂t

Grado de una ecuación diferencial: es el grado de la derivada de mayor orden que aparece en la ecuación.
Ejemplos:
1. y 0 = 4x + 3. EDO de primer orden de primer grado.

2. 5wy 0 + y = 3w2 . EDO de primer orden de primer grado.

2
3. y 000 + 2 (y 00 ) + y 0 = cos x. EDO de tercer orden de primer grado.

4. y 000 + (4 sen x)y 00 + 2y 0 = x2 . EDO de tercer orden de primer grado.

3 2
5. (y 00 ) + (y 0 ) + 3y = 3x. EDO de segundo orden de tercer grado.
∂z ∂z
6. = z + x . EDP de primer orden de primer grado.
∂x ∂y
∂ 2w ∂ 2w
7. 2
+ 2 = x2 + y. EDP de segundo orden de primer grado.
∂x ∂t
Reconocer que estamos ante la presencia de una ecuación diferencial y poder clasificarla es importante; sin
embargo, determinar de si una función es o no solución de una ecuación diferencial se torna necesario para
nuestros intereses.

Solución de una ecuación diferencial ordinaria

Una solución de ecuación diferencial ordinaria, sobre un intervalo I, es cualquier función y tal que y y sus
derivadas existen en I y además, satisfacen la ecuación; es decir, al sustituirlas en la ecuación se obtiene una
identidad.
Ejemplo 1. Verifique si la función dada por y = ex + e−3x , es solución de la ecuación diferencial dada por
y 000 + 2y 00 = 3y 0 .
Solución:
Primero calculamos las derivadas y 0 = ex − 3e−3x ; y 00 = ex + 9e−3x y y 000 = ex − 27e−3x .

3
Sustituyendo en la ecuación diferencial, se

tiene:
000 00 −3x −3x
x x
= ex − 27e−3x + 2ex + 18e−3x = 3ex − 9e−3x = 3 ex − 3e−3x = 3y 0 .



y + 2y = e − 27e + 2 e + 9e
De esta manera se verifica que dicha función sı́ es solución de la ecuación diferencial dada.
1
Ejemplo 2. Verifique si la función definida por y = x 2 , con x ∈]0, +∞[, es o no solución de la ecuación
diferencial dada por 2x2 y 00 + 3xy 0 − y = 0.
Solución: −1 −3
0 x2 00 −x 2
Cálculo de las derivadas y = ;y = .
2 4
Sustituyendo en la ecuación diferencial,
! se tiene:!
−3 −1 1 1 1 1
2 00 0 2 −x 2 x2 1 −x 2 3x 2 1 −3x 2 3x 2
2x y + 3xy − y = 2x · + 3x · −x =
2 + −x =
2 + =0
4 2 2 2 2 2
Por tanto, sı́ es solución de la ecuación diferencial.

Ejemplo 3. Verifique si la función definida por u = x tan x, con tan x bien definida, es o no solución de la
ecuación diferencial dada por xu0 = u + x2 + u2 .
Solución:
Cálculo de las derivadas y 0 = tan x + x sec2 x.
Sustituyendo en la
ecuación diferencial, se tiene:
x · tan x + x sec2 x = x tan x + x2 sec2 x = x tan x + x2 (1 + tan2 x) = x tan x + x2 + x2 tan2 x = u + x2 + u2 .
Por tanto, sı́ es solución de la ecuación diferencial.

Ejercicios
1. Determine si la función dada es o no solución de la ecuación diferencial correspondiente.

(a) y = e3x + 10e2x de la ecuación diferencial dada por y 0 − 2y = e3x

6 6 dy
(b) y = − e−20x de la ecuación diferencial dada por + 20y = 24
5 5 dx
(c) y = x + 1 de la ecuación diferencial dada por (y 0 )3 + xy 0 = y
d2 y 2
(d) y = x2 − x−1 , x 6= 0, de la ecuación diferencial dada por 2
= 2 · y.
dx x
(e) y = c1 e−x + c2 e2x 00 0
de la ecuación diferencial dada por y − y − 2y = 0.
(f) y = x2 + x3 de la ecuación diferencial dada por x2 y 00 − 4xy 0 = 6y.
(g) y = xex de la ecuación diferencial dada por y 00 + y = 2y 0 .
(h) y = c1 sen x + c2 cos x, con c1 y c2 constantes, de la ecuación diferencial dada por y 00 + y = 0
(i) y = c1 e3x + c2 e−4x , con c1 y c2 constantes, de la ecuación diferencial dada por y 00 + y 0 = 12y
(j) y = e3x cos(2x) de la ecuación diferencial dada por y 00 − 6y 0 + 13y = 0
d2 y dy
(k) y = e2x + xe2x de la ecuación diferencial dada por 2
− 4 + 4y = 0
dx dx
(l) y = x cos(ln x), con x > 0, de la ecuación diferencial dada por x2 y 00 + 2y = xy 0
(m) y = x2 ex de la ecuación diferencial dada por y 000 − 3y 00 + 3y 0 = y
(n) y = cos(x) de la ecuación diferencial dada por y 0 + y tan(x) = 0
(o) y = e− sen(x) + sen(x) + 1 de la ecuación diferencial y 0 + y cos(x) = cos(x) sen(x)

4
 2. Determine para cuál valor o valores de r la función f (x) = erx es solución de la ecuación diferencial
10y − 11y 0 + y 00 = 0

3. Determine para cuál valor o valores de k la función f (x) = xk es solución de la ecuación diferencial
x2 y 00 + 4xy 0 + 2y = 0

4. Determine para qué valor o valores de k la ecuación diferencial y 00 + y 0 − 6y = 0 tiene una solución de
la forma y = ekx

5. Determine si la función dada por P = cekt , con c y k constantes, es o no solución de la ecuación

dP
diferencial dada por = kP . (Modelo de crecimiento de la población)
dt

Soluciones implı́citas de una ecuación diferencial ordinaria

En algunas ocasiones la función solución queda definida implı́citamente por una ecuación que involucra a las
variables dependientes e independientes. Por ello se le denomina solución implı́cita.

Ejemplo 4. Suponga que la ecuación y 3 + 3x + 7 = 6y define a y como una función de x. Verifique si esta
función ası́ definida es o no solución implı́cita de la ecuación diferencial ordinaria dada por y 00 − 2y (y 0 )3 = 0.
Solución:
Calculando las derivadas se tiene:
−3 −1
3y 2 y 0 + 3 = 6y 0 ⇒ 3y 2 y 0 − 6y 0 = −3 ⇒ 3y 0 (y 2 − 2) = −3 ⇒ y 0 = 2
⇒ y0 = 2 .
3(y − 2) y −2
2y 2y −1 −2y
y 00 = 0
2 ·y = 2 · 2 = .
2
(y − 2) (y − 2) y − 2
2 (y − 2)3
2
Sustituyendo en la ecuación diferencial,
 3se tiene:
3 −2y −1 −2y −2y −2y 2y
y 00 − 2y (y 0 ) = 3 − 2y = 3 − 3 = 3 + =0
2
(y − 2) 2
y −2 2
(y − 2) 2
(y − 2) 2
(y − 2) (y − 2)3
2
Por tanto se tiene que la función dada sı́ es solución implı́cita de la ecuación diferencial.
−x
Ejemplo 5.Verificar si x2 + y 2 = 5 es o no solución implı́cita de la ecuación diferencial dada por y 0 = ;
y
sujeta a la condición y(−1) = 2.
Solución:
−2x −x
Calculando las derivadas se tiene: 2x + 2yy 0 = 0 ⇒ y 0 = ⇒ y0 = . Con lo cual sı́ es solución de la
2y y
ecuación diferencial. Analizando la condición se tiene que si x = −1, y = 2; ası́ al sustituir en la ecuación
original (ecuación solución) se tiene que (−1)2 + (2)2 = 1 + 4 = 5. Por tanto también se satisface la condición.
xy
Ejemplo 6. Verificar si x2 = 2y 2 ln y es o no solución implı́cita de la ecuación y 0 =
x2 + y2
Solución:
2x
Calculando las derivadas se tiene: 2x = 4yy 0 ln y + 2yy 0 ⇒ 2x = y 0 (4y ln y + 2y) ⇒ y 0 =
2(2y ln y + y)
0 xy 0 xy
⇒y = ⇒y = 2 . Con lo cual sı́ es solución de la ecuación diferencial.
(2y 2 ln y + y 2 ) x + y2

5
Ejercicios

1. Determine si la función dada es o no solución implı́cita de la ecuación diferencial correspondiente.

(a) x2 y + y 2 = c1 , con c1 constante, de la ecuación diferencial dada por 2xy dx + (x2 + 2y) dy = 0
dy 3x2
(b) y 2 − x3 + 8 = 0 de la ecuación diferencial dada por =
dx 2y
dy
(c) 4x2 = y 2 de la ecuación diferencial dada por y − 4x = 0
dx
3x2 − 4x − 2
(d) y 2 − 2y − x3 + 2x2 + 2x = c de la ecuación diferencial dada por y 0 =
2(y − 1)
y y2
(e) y = ce x de la ecuación diferencial dada por y 0 =
xy − x2
(f) xy 2 − y 3 = 7 de la ecuación diferencial dada por (2x − 3y)y 0 + y = 0
y
(g) (x + y)2 = xe x de la ecuación diferencial dada por (x2 + y 2 ) dx + (x2 − xy) dy = 0

Consideraciones sobre la solución de una ecuación diferencial

Solución general: normalmente una ecuación diferencial de orden n tiene una solución que involucra n
constantes; tal solución se llama solución general de la ecuación diferencial.
Solución particular: es la que se obtiene de una solución general al sustituir las constantes por valores
especı́ficos, los cuales puede provenir de ciertas condiciones dadas.
Solución singular: estas son soluciones que no se pueden obtener a partir de la solución general.
Problema con valores iniciales: es un problema que busca determinar la solución de una ecuación dife-
rencial sujeta a condiciones sobre la función desconocida y sus derivadas en un cierto valor de la variable
independiente. Tales condiciones se llaman condiciones iniciales.

Ejemplos
3
a. La función definida por u(t) = c1 et + c2 e−t + 2t es solución de la ecuación dada por u00 − u = e2t ,
e
además, se le llama solución general dado que tiene dos constantes arbitraria que son c1 y c2 .
 2 2
x 1
b. La función definida por y = + c es solución de la ecuación dada por y 0 − xy 2 = 0. Además, es
4
x4
solución general. Si c = 0; la solución y = se llama solución particular, la función y = 0 se le llama
16
solución singular.

Ejercicio.

a. Verificar que la función dada por y = Aex + Bxex es solución de la ecuación diferencial y 00 − 2y 0 + y = 0.
Luego determine la solución particular que satisface las condiciones iniciales y(0) = 1 y y 0 (0) = −2.

Como propósito del curso, interesa, determinar la solución de una ecuación diferencial. Por ello se discutirán
algunos de los métodos que permiten resolverlas.

6
 Algunos métodos para resolver ecuaciones diferenciales no lineales de primer orden.

A pesar que no existe un método general que permita determinar la solución de las ecuaciones diferenciales,
algunas por sus caracterı́sticas ofrecen la posibilidad de aplicar ciertos procedimientos y ası́ resolverlas.
Además, en muchos de los casos, la solución encontrada corresponde a una solución implı́cita (principalmente
en las ecuaciones de primer orden no lineales). Corresponde en este curso, analizar y emplear, en la resolución
de ecuaciones diferenciales, algunos de los métodos más comunes.
dy
• Método de variables separables. Si una ecuación diferencial de la forma M (x, y) + N (x, y) =0
dx
(no lineales de primer orden) puede ser escrita de la forma M (x)M (y) dx + N (x)N (y) dy = 0, se
dice que dicha ecuación diferencial es de variables separables. Esto por que puede expresarse como
M (x) N (y)
dx = − dy, lo cual puede resolverse al integrar a ambos lados. Ası́;
N (x) M (y)
Z Z
M (x) N (y)
dx = − dy
N (x) M (y)

dy
La ecuación diferencial = y −2 xe3x+4y es de variables separables, mientras que la ecuación diferencial
dx
dy
= y + sen(x) no es de variables separables.
dx

Aplicando el método de variables separables, se resolverá algunos ejemplos.

Ejemplo 1. Determine la solución de la ecuación 2(y + 3) dx − xy dy = 0.
Solución: separando variables: Z Z Z Z  
2 y 2 y 2 3
2(y + 3) dx = xy dy ⇒ dx = dy ⇒ dx = dy ⇒ dx = 1− dy ⇒
x y+3 x y+3 x y+3
2 ln x + c = y − 3 ln(y + 3). Que corresponde a la solución general.

x + xy 2
Ejemplo 2. Determine la solución de la ecuación y 0 = , sujeta a y(1) = 0.
4y
Solución: primero debe resolverse la ecuación y al final seZanalizará la condición. Separando variables:
x + xy 2 x(1 + y 2 )
Z
dy dy 4y 4y
= ⇒ = ⇒ dy = x dx ⇒ dy = x dx
dx 4y dx 4y 1 + y2 1 + y2
x2
⇒ 2 ln |1 + y 2 | = + k. Que corresponde a la solución general.
2
x2 (1)2
Para la condición si x = 1 entonces y = 0. Ası́; 2 ln(1 + y 2 ) = + k ⇒ 2 ln(1 + (0)2 ) = +k ⇒
2 2
1 −1
2 ln(1) = + k ⇒ = k.
2 2
x2 1
Por tanto la solución particular corresponde a 2 ln(1 + y 2 ) = − ⇒ 4 ln(1 + y 2 ) = x2 − 1
2 2

Ejemplo 3. Determine la solución de la ecuación xy 0 = y 2 + 1.

Solución: separando las variables: Z Z
dy 2 1 1 1 1
x =y +1⇒ 2 dy = dx ⇒ 2
dy = dx ⇒ arctan(y) = ln |x| + c. (Solución).
dx y +1 x y +1 x

7
 π  π
Ejemplo 4. Determine la solución de la ecuación sen2 (y) dx+ 1 − sen2 (x) dy = 0, sujeta a y


= .
4 4
Solución: separando las variables:
sen2 (y) dx + 1 − sen2 (x) dy = 0 ⇒ sen2 (y) dx + cosZ 2
(x) dy = 0 ⇒Zsen2 (y) dx = − cos2 (x) dy ⇒


1 −1
2
dx = dy ⇒ sec2 (x) dx = − csc2 (y) dy ⇒ sec2 (x) dx = − csc2 (y) dy
cos (x) sen2 (y)
1
⇒ tan(x) + k = cot(y) ⇒ tan(x) + k = ⇒ tan(y) [tan(x) + k] = 1 (solución general).
tan(y)
π π π  h π  i
Para la condición x = y y = , con lo cual tan · tan + k = 1 ⇒ 1 · (1 + k) = 1 ⇒ k = 0.
4 4 4 4
La solución particular corresponde a tan(y) tan(x) = 1

−y
Ejemplo 5. Determine la solución de la ecuación y 0 = , sujeta a y(1) = 3.
x
Solución: separando las variables:
−y −1 −1
Z Z
dy 1 1
= ⇒ dy = dx ⇒ dy = dx ⇒ ln |y| = − ln |x| + k (solución general).
dx x y x y x
Para la condición x = 1 y y = 3. Ası́; ln |3| = ln |1| + k ⇒ ln(3) = 0 + c ⇒ ln(3) = c. Solución particular
ln |y| = − ln |x| + ln(3) ⇒ ln |y| + ln |x| = ln(3) ⇒ yx = 3.
π 
0
Ejemplo 6. Determine la solución de la ecuación y sen(x) = y ln(y), sujeta a y = e.
2
Solución: separando las variables: Z Z
dy 1 1 1
sen(x) = y ln(y) ⇒ dy = dx ⇒ dy = csc(x) dx.
dx Z y ln(y) sen(x) y ln(y) Z
1 1
La integral dy se resuelve por el método de sustitución ası́, dy = ln(ln(y)), por tanto:
Z y Zln(y) y ln(y)
1 1 cos(x)
dy = csc(x) dx ⇒ ln(ln(y)) = ln | csc(x)−cot(x)|+c ⇒ ln(ln(y)) = ln − +c.
y ln(y) sen(x) sen(x)
cos π2


π 1
Para condición se tiene que x = y y = e. Por tanto, ln(ln(e)) = ln
−
+ c ⇒ ln(1) =
2 sen π2 sen π2
ln(1 − 0) + c ⇒ 0 − 0 = c ⇒ 0 = c.
1 cos(x)
Ası́ la solución particular corresponde a: ln(ln |y|) = ln − ⇒ ln(y) = csc(x) − cot(x).
sen(x) sen(x)
x2
Ejemplo 7. Determine la solución de la ecuación y 0 = , sujeta a y(0) = 2.
y (1 +Zx3 )
x2 x2 x2 y2 ln |1 + x3 |
Z
dy
Solución: = ⇒ y dy = dx ⇒ y dy = dx ⇒ = + c.
dx y (1 + x3 ) (1 + x3 ) (1 + x3 ) 2 3
(2)2 ln |1 + (0)3 | ln(1)
Para la condición x = 0 y y = 2. = +c ⇒ 2 = + c ⇒ 2 = c, con lo cual la
2 3 3
2 3
y ln |1 + x |
solución particular corresponde a = + 2.
2 3

Ejemplo 8. Determine la solución de la ecuación y 0 = 1 + x + y 2 + xy 2 , sujeta a y(0) = −1.

dy dy dy
Solución: = (1 + x) + (y 2 + xy 2 ) ⇒ = (1 + x) + y 2 (1 + x) ⇒ = (1 + x)(1 + y 2 )
dx dx dx
x2
Z Z
1 1
⇒ dy = (1 + x) dx ⇒ dy = (1 + x) dx ⇒ arctan(y) = x + + c. (solución general).
1 + y2 1 + y2 2

8
 (0)2 −π
Para la condición x = 0 y y = −1. Ası́, arctan(−1) = 0 + +c ⇒ = c, con lo cual la solución
2 4
x2 π
particular corresponde a arctan(y) = x + − .
2 4

Ejercicios
1. Resolver cada una de las siguientes ecuaciones.
(a) (4 + x)y 0 = y 3
(b) e−y (1 + y 0 ) = 1
dQ
(c) = k(Q − 70), con k constante.
dt
(d) 2y cos(x) dx + cos2 (x) dy = 0
(e) cos x cos y dx + sen x sen y dy = 0
(f) 3y dx = 2x dy
(g) ex (y − 1) dx + 2(ex + 4) dy = 0
(h) y 0 = xy 2
(i) x2 yy 0 = ey
(j) (e2x + 4)y 0 = y
(k) (1 + ln x) dx + (1 + ln y) dy = 0
(l) dx + e3x dy = 0
dx x2 y 2
(m) =
dy 1+x
0 5x−3y
(n) y = e
(o) (1 + x2 + y 2 + x2 y 2 ) dy = y 2 dx
 2
dy y+1
(p) y ln(x) =
dx x
(q) (e + 1) e dx + (ex + 1)3 e−x dy = 0
y 2 −y

(r) y ln y dx − x dy = 0
2
(s) ex+y sen(x) dx + (2y + 1) e−y dy = 0
dN
(t) + N = N tet+2
dt
dy
(u) e2y − y cos(x) = ex sen(2x). Sugerencia sen(2x) = 2 sen(x) cos(x).


dx
(v) (1 + x2 + y 2 + x2 y 2 ) dy = (1 + y 2 )dx
2. Resolver cada una de las siguientes ecuaciones sujetas a las condiciones indicadas
(a) (e−y + 1) sen(x) dx = (1 + cos(x)) dy, sujeta a y(0) = 0
dy
(b) + xy = y, sujeta a y(1) = 3
dx
dy
(c) = y 2 + 1, sujeta a y(1) = 0
dx
dy
(d) = ex−y , sujeta a y(0) = 1
dx
dy xy + 3x
(e) = 2 , sujeta a y(2) = 2
dx x +1

9
 −x
(f) y 0 = , sujeta a y(4) = −3
y
(g) y 0 + y 2 − y =, sujeta a y(2) = 4
(h) ey cos(x) + cos(x) + (ey sen(x) + ey ) y 0 = 0, sujeta a y(π) = 0
dy y − 4x
3. Verifique que la ecuación dada por = no es de variables separables, sin embargo al
dx x−y
y
reemplazar la variable y por v con v = , entonces la ecuación es de variables separables en
x
términos de x y v.
4. Newton propuso una ley experimental sobre el enfriamiento o calentamiento de un cuerpo: “la
rapidez con la que cambia la temperatura de un cuerpo es proporcional a la diferencia de tempe-
raturas entre el cuerpo y el medio ambiente que lo rodea”. Si T = T (t) representa la temperatura
del cuerpo en el instante t y Ta representa la temperatura del medio ambiente, entonces se satisface
dT
que = k(T − Ta ). Haciendo uso de lo anterior resolver los siguientes ejercicios.
dt
(a) Una barra metálica a una temperatura de 100◦ C se pone en un cuarto a una temperatura
constante de 0◦ C. Si después de 20 minutos la temperatura de la barra es de 50◦ C ¿Cuánto
tiempo tardará para que la barra tenga una temperatura de 25◦ C?
(b) Un cuerpo a una temperatura desconocida se coloca en un refrigerador a una temperatura
constante de 1◦ C. Si después de 20 minutos la temperatura del cuerpo es de 40◦ C y después
de 40 minutos la temperatura del cuerpo es de 20◦ C. Determine la temperatura inicial del
cuerpo.

Estudiemos ahora otro método para resolver ecuaciones diferenciales.

• Método de ecuaciones homogéneas. Algunas ecuaciones diferenciales que no son de variables se-
parables, pueden transformarse en una de variables separables, mediante una sustitución adecuada. Tal
es el caso que se presenta en las ecuaciones diferenciales homogéneas.
Definición: una ecuación diferencial de primer orden se llama homogénea, si puede escribirse de la
dy y
forma =F , para alguna función apropiada F .
dx x
x2 + 3y 2
0
Ejemplo 1. La ecuación diferencial dada por y = es homogénea, dado que si se toma a
2xy
2
1 + 3 xy


y 0 x2 + 3y 2 0 0 1 + 3z 2
z = , se tiene que y = ⇒y = ⇒ y = .
2 xy


x 2xy 2z
dy x2 + xy + y 2
Ejemplo 2. La ecuación = es homogénea, dado que se puede escribir de la forma
  dx x2 + y 2
2 y y2
dy x 1 + x
+ x2 y dy 1 + z + z2
=   , si se hace z = , se tiene que = .
dx x2 1 + y
2
x dx 1 + z2
x2

dy x+y
Ejemplo 3. La ecuación = ln(x) − ln(y) + es homogénea, dado que se puede escribir de la
dx x − 2y
dy y x 1 + y

y dy 1+z
x
forma = − ln + 2y , si se hace z =

, se tiene que = − ln(z) + .
dx x x 1− x x dx 1 − 2z

10
Si una ecuación diferencial es homogénea, al hacer las sustituciones correspondientes, se puede con-
vertir en una ecuación de variables separables; por tanto, se resuelve aplicando el método de variables
separables. Algunos ejemplos.

x+y
Ejemplo 1. Resolver la ecuación dada por y 0 = .
x
x+y x y y
Solución: veamos si es homogénea y 0 = ⇒ y 0 = + ⇒ y 0 = 1 + . Al hacer el cambio de
x x x x
y 0
variable z = , se tiene y = 1 + z. Sin embargo, como y depende de x y z depende de x y de y, faltarı́a
x
0 y dz
escribir a y en términos de z y de x. Como z = ⇒ z · x = y ⇒ x + z = y 0 . De esta manera:
x dx
dz dz dz
y0 = 1 + z ⇒ x+z = 1+z ⇒ x = 1+z−z ⇒ x = 1 que es una ecuación de variables
dx dx dx
separables. Z Z
dz 1 1 y y
x = 1 ⇒ dz = dx ⇒ dz = dx ⇒ z = ln |x| + c. Sustituyendo z = ⇒ = ln |x| + c. Que
dx x x x x
corresponde a la solución general.
 y  y
Ejemplo 2. Resolver la ecuación dada por y ln + 1 dx − x ln dy = 0.
x x
Solución:
  y  veamos
 si es homogénea
y  y  y
y ln + 1 dx − x ln dy = 0 ⇒ y ln + 1 dx = x ln dy
 x y  xy  dy y   y
x   y  x
⇒ y ln + 1 = x ln ⇒ ln + 1 = ln y 0 . Al hacer el cambio de variable
x x dx x x x
y z (ln(z) + 1) y dz
z = , se tiene z (ln(z) + 1) = ln(z) y 0 ⇒ = y 0 . Como z = ⇒ z · x = y ⇒ x + z = y0.
x ln(z) x dx
De esta manera:
z (ln(z) + 1) dz z ln(z) z dz z dz z dz
= x+z ⇒ + = x+z ⇒ z+ = x+z ⇒ = x⇒
ln(z) dx ln(z) ln(z) dx ln(z) dx ln(z) dx
ln2 xy


ln2 z
Z Z
1 ln(z) 1 ln(z)
dx = dz ⇒ dx = dz ⇒ ln(x) + c = ⇒ ln(x) + c =
x z x z 2 2

Ejemplo 3. Resolver la ecuación dada por 2xyy 0 = 4x2 + 3y 2 .

4x2 3y 2 x 3y y
Solución: veamos si es homogénea, y 0 = + ⇒ y0 = 2 + . Sea z = .
2xy 2xy y 2x x
dz 2 3z dz 2 3z dz 4 + 3z 2 − 2z 2 dz z2 + 4
x+z = + ⇒ x= + −z ⇒ x= ⇒ x=
dx z 2 dx z 2 Z dx 2z dx 2z
Z 2
2z 1 2z 1 y
⇒ 2 dz = dx ⇒ dz = dx ⇒ ln |z 2 + 4| = ln |x| + c ⇒ ln 2 + 4 = ln |x| + c
z +4 x z2 + 4 x x
h  y i
Ejemplo 4. Resolver la ecuación dada por x y 0 − csc = y.
 y i x
h y  y  y
Solución: x y 0 − csc = y ⇒ y 0 = + csc . Sea z = .
x x x xZ Z
dz dz 1 1 1
x + z = z + csc(z) ⇒ x = csc(z) ⇒ dz = dx ⇒ sen(z)dz = dx
dx y dx csc(z) x x
⇒ − cos = ln |x| + c.
x
 
dx x 2 x
Ejemplo 5. Resolver la ecuación dada por = + sen .
dy y y

11
    
dx x 2 x 0 x 2 x x
Solución: = + sen ⇒ x = + sen . Resulta conveniente tomar z = ; donde se
dy y y y y y
dz
tiene que zy = x, ası́ y + z = x0 .
dy Z Z
dz 2 dz 2 1 1 2 1
Al sustituir y + z = z + sen (z) ⇒ y = sen (z) ⇒ 2
dz = dy ⇒ csc (z) dz = dy ⇒
  dy dy sen (z) y y
x
− cot = ln |y| + c
y
y y2
Ejemplo 6. Resolver la ecuación dada por y 0 = + 2 , sujeta a y(1) = 1.
x x
0 y y2 0 y  y z y
Solución: y = + 2 ⇒ y = + . Sea z = .
x x x x xZ Z Z Z
dz 2 dz 2 1 1 1 1 −2 1
x+z =z+z ⇒ x = z ⇒ 2 dz = dx ⇒ 2
dz = dx ⇒ z dz = dx
dx dx z x z x x
−1 x
⇒ = ln |x| + c ⇒ − = ln |x| + c.
z y
1
Para la condición se tiene que si x = 1 entonces y = 1. Al sustituir − = ln |1| + c ⇒ −1 = c, ası́ la
1
x
solución particular corresponde a ⇒ − = ln |x| − 1 ⇒ −x = y (ln |x| − 1)
y

Tal como se indicó al inicio, para ciertas ecuaciones diferenciales es necesario realizar los cambios por
expresiones que contribuyan a lograr separar las variables. Por ejemplo:

dy
Ejemplo 7. Resolver la ecuación dada por = x + y, sujeta a y(2) = 2.
dx
dz dz
Solución: al hacer el cambio de variable z = x + y, se tiene además que = 1 + y0 ⇒ − 1 = y0,
Z Zdx dx
dz dz 1 1
con lo cual: −1 = z ⇒ = z+1 ⇒ dz = dx ⇒ dz = dx ⇒ ln(z + 1) = x + c ⇒
dx dx z+1 z+1
ln (x + y + 1) = x + c. Que corresponde a la solución general.

Para la condición se tiene que si x = 2, entonces y = 2. Ası́ ln(2 + 2 + 1) = 2 + c ⇒ ln(5) − 2 = c. Por

tanto la solución particular corresponde a ln (x + y + 1) = x + ln(5) − 2 ⇒ y = ex+ln(5)−2 − x − 1.

Ejercicios

1. Resolver cada una de las siguientes ecuaciones diferenciales.

x+y
(a) y 0 =
x
0 y
(b) y = − 1
x
dy −(x + 2y)
(c) =
dx 2x + y
(d) (x − 2y) dx + (2x + y) dy = 0
(e) 2(2x2 + y 2 ) dx − xy dy = 0
(f) xy dx − (x2 + 3y 2 ) dy = 0

12
 dy x2 + y 2
(g) =
dx xy
y
0
(h) xy = y + xe x
(i) (x3 + y 3 ) dx − xy 2 dy = 0
2xy dx + y 2 − 3x2 dy = 0


(j)
(k) y (ln(x) − ln(y)) dx = (x ln(x) − x ln(y) − y) dy
−y 2
dy e x + xy 2
(l) = y
dx x
dy x+y y−1
(m) = . Sugerencia: use z =
dx x+1   x+1
dx x x x
(n) = − sen2 . Sugerencia: use z = .
dy y y y
dy p
(o) = 3 3 y − 3x. Sugerencia: use z 3 = y − 3x
dx
2. Resuelva cada una de las siguientes ecuaciones diferenciales, sujetas a las condiciones indicadas
(a) 2xy + 3y 2 dx = 2xy + x2 dy, sujeta a y(1) = 1

y + x cos2 xy


dy π
(b) = , sujeta a y(1) =
dx x 4
3 3 2
(c) (x + y ) dx − (xy ) dy = 0, sujeta a y(1) = 0
(d) y (ln(x) − ln(y)) dx = (x ln(x) − x ln(y) − y) dy, sujeta a y(0) = 1
−y
dy x2 e x + y 2
(e) = , sujeta a y(1) = 0
dx xy
y + x cos2 xy


dy π
(f) = , sujeta a y(1) =
dx x 4
Otros métodos, empleados para resolver ecuaciones diferenciales, requieren de conceptos no desarrollados
hasta entonces ¿Cuáles son?

Derivadas parciales y diferenciales totales

Derivadas parciales: si una función presenta más de una variable, al calcular sus derivadas si estas existen
(llamadas derivadas parciales), es importante indicar respecto a cual de ellas se está derivando. Veamos
algunos ejemplos (en este curso solo se emplearán funciones en dos variables):
∂ ∂
Ejemplo 1. Considere la función Z(x, y) = x2 + 2y 5 . Calcular Z(x, y) y Z(x, y)
∂x ∂y
∂
Solución: al calcular Z(x, y) se debe considerar a x como variable y a y como si fuese constante. Ası́,
∂x
∂
Z(x, y) = 2x + 0 = 2x
∂x
∂
Al calcular Z(x, y) se debe considerar a y como variable y a x como si fuese constante. Ası́,
∂y
∂
Z(x, y) = 0 + 10y 4 = 10y 4
∂y

13
 ∂ ∂
Ejemplo 2. Considere la función Z(x, y) = 2x2 cos(y) + x4 sen(y). Calcular Z(x, y) y Z(x, y)
∂x ∂y
∂ ∂
Solución: Z(x, y) = 4x cos(y) + 4x3 sen(y) y Z(x, y) = −2x2 sen(y) + x4 cos(y)
∂x ∂y
∂Z ∂Z
Ejemplo 3. Considere la función Z(x, y) = ex+y sen(x) + ln(x) tan(y). Calcular y
∂x ∂y
∂Z x+y x+y tan(y) ∂Z x+y
Solución: =e sen(x) + e cos(x) + y =e sen(x) + ln(x) sec2 (y)
∂x x ∂y
∂Z ∂Z
Ejemplo 4. Considere la función Z(x, y) = x3 y + y 4 x2 + sen(x). Calcular y
∂x ∂y
∂Z 2 4 ∂Z 3 3 2
Solución: = 3x y + 2y x + cos(x) y = x + 4y x
∂x ∂y
∂Z ∂Z
Ejemplo 5. Considere la función Z(x, y) = [y + ln(x)]5 . Calcular y
∂x ∂y
4
∂Z 1 5 [y + ln(x)] ∂Z
Solución: = 5 [y + ln(x)]4 · = y = 5 [y + ln(x)]4
∂x x x ∂y
x2 ∂Z ∂Z
Ejemplo 6. Considere la función Z(x, y) = 4 . Calcular y
y ∂x ∂y
x2 ∂Z 2x ∂Z −4x2
Solución: observe que Z(x, y) = 4 = x2 y −4 . Ası́, = 2xy −4 = 4 y = −4x2 y −5 =
y ∂x y ∂y y5
Ejemplo 7. Considere la función Z(x, y) = cos(x) − x sen(x) + y 2 . Calcular las derivadas parciales.
Solución:
∂Z
= − sen(x) − [sen(x) + x cos(x)] + 0 = − sen(x) − sen(x) − x cos(x) = −2 sen(x) − x cos(x)
∂x
∂Z
= 2y
∂y
Ejercicios.

1. Para cada una de las siguientes funciones, determine las derivadas parciales

(a) Z(x, y) = cos(xy) + x3 y 2

(b) Z(x, y) = −5x4 + 9x2 y 2 − 5y 4
sen(x)
(c) Z(x, y) =
cos(y)
(d) Z(t, m) = ln(tm)
(e) I(x, y) = x cos(y) + 2xey
(f) I(t, p) = t sen(p) + t2 ep + 2p
(g) I(x, y) = exy + 2x3 y 2
(h) F (x, t) = arctan(xt) − x + t
(i) F (x, t) = ex sen(t) − 2t sen(x)
(j) F (x, y) = ln(x) + ln(y) − yex
x
(k) F (x, t) = cos(x) + ln(t) + − tan(t)
t

14
 (l) Z(x, y) = x ln(y) + xy + y ln(x) + x2 y 4

Diferencial total para una función de dos variables: sea Z una función de dos variables, si sus derivadas
∂Z ∂Z
parciales existen, entonces se define su diferencial total como dZ = dx + dy
∂x ∂y
Ejemplo 1. Determine el diferencial total de Z(x, y) = cos(xy) + 3x2 y 2 .
∂Z ∂Z
Solución: dado que = − sen(xy) y + 6xy 2 y además, = − sen(xy) x + 6x2 y, se tiene que
∂x ∂y
∂Z ∂Z
dy = −y sen(xy) + 6xy 2 dx + −x sen(xy) + 6x2 y dy



dZ = dx +
∂x ∂y
Ejemplo 2. Determine el diferencial total de F (x, y) = exy + 2 cos(x).
∂F ∂F
Solución: dado que = exy · y − 2 sen(x) = yexy − 2 sen(x) y además, = exy · x = xexy , se tiene
∂x ∂y
∂F ∂F xy xy
dF = dx + dy = (y e − 2 sen(x)) dx + (x e ) dy
∂x ∂y
Ejemplo 3. Determine el diferencial total de Z(x, t) = sen(x + t) + ln(xt).
∂Z t 1 ∂F x 1
Solución: dado que = cos(x + t) + = cos(x + t) + y además, = cos(x + t) + = cos(x + t) + ,
∂x  xt  x  ∂t xt t
1 1
se tiene que: dZ = cos(x + t) + dx + cos(x + t) + dt
x t
Ejercicios.
1. Para cada una de las siguientes funciones, determine el diferencial total.

(a) Z(x, y) = tan(xy) + x5 y 3

(b) Z(x, y) = cos(x) y 4 − 6y 2
cos(xy)
(c) Z(x, y) =
sec(y)
(d) Z(t, m) = ln(tm) + m cos(t)
(e) I(x, y) = y cos(y) − 3x ex
t
(f) I(t, p) = + t3 e2p
p
(g) I(x, y) = e2xy − 3x4 y 3
(h) F (x, t) = arctan(xt) − x2 + t3
(i) F (x, t) = e3x cos(t) − 3t sen(x)
(j) F (x, y) = ln(y) − ln(x2 ) + yex
(k) F (x, t) = cos(x) + ln(xt) − tan(t)
(l) Z(x, y) = y ln(y) + x2 y + y ln(x) + x3 y 2
p
(m) Z(x, y) = (4x + 5y) + sen(x) + y

Otro concepto importante y que emplearemos en la solución de ecuaciones diferenciales es el de integración

parcial.

15
Integración parcial: si una función presenta varias variables y se necesita integrar, debe indicarse respecto
a cual de la variables se hará y estaremos realizando una integración parcial.
Z
x3 + y ∂x. En este caso integraremos respecto a la variable x, por tanto conside-


Ejemplo 1. Resolver
x4
Z
x3 + y ∂x =


raremos a y como una constante. Ası́, + xy + C(y), es importante destacar que como se
4
integra respecto a x, al agregar la constante de integración esta puede ser cualquier función que no dependa
de x y ello se indica al escribir C(y).
Z
x3 + y ∂y. En este caso integraremos respecto a la variable y, por tanto conside-


Ejemplo 2. Resolver
y2
Z
x3 + y ∂y = x3 y + + C(x), se destaca una vez más que como se


raremos a x como una constante. Ası́,
2
integra respecto a y, al agregar la constante de integración esta puede ser cualquier función que no dependa
de y y ello se indica al escribir C(x).
Z
Ejemplo 3. Resolver (sen(x) + cos(y)) ∂x.
Z
Solución: (sen(x) + cos(y)) ∂x = − cos(x) + x cos(y) + C(y)
Z
Ejemplo 4. Resolver (sen(x) + cos(y)) ∂y.
Z
Solución: (sen(x) + cos(y)) ∂y = y sen(x) + sen(y) + C(x)
Z
x2 ey ∂x.


Ejemplo 5. Resolver
x3
Z
x2 ey ∂x = ey + C(y)


Solución:
3
Z
x2 ey ∂y.


Ejemplo 6. Resolver
Z
x2 ey ∂y = x2 ey + C(x)


Solución:
Z
5x4 − 9x2 y 2 + 5y 4 ∂x.


Ejemplo 7. Resolver
Z
5x4 − 9x2 y 2 + 5y 4 ∂x = x5 − 3x3 y 2 + 5xy 4 + C(y)


Solución:
Z
Ejemplo 8. Resolver (y cos(x) + 2xey ) ∂y.
y2
Z
Solución: (y cos(x) + 2xey ) ∂y = cos(x) + 2xey + C(x)
2
Z 
y 
Ejemplo 9. Resolver + ln(x) + tan(y) ∂x.
Z  x
y 
Solución: + ln(x) + tan(y) ∂x = y ln(x) + x ln(x) − x + x tan(y) + C(y)
x

16
Ejercicios.
1. Resolver cada una de las siguientes integrales parciales
Z
exy + x2 y 2 ∂x


(a)

1 y2
Z  
(b) + ∂y
y x2
Z
5 − 3x2 y 3 ∂x


(c)
Z
x3 y 2 + sen(x) cos(y) ∂y


(d)
Z  
2y 2x
(e) + ∂x
x y
Z  
− cos(y)
(f) + ln(x) ∂y
sen(x)
Z  
1
(g) + sen(2x) ∂x
x2
Z
x2 + y 3 tan(x) ∂y


(h)
Z
x ln(x) + xy 2 ∂x


(i)
Z  
x 2xy
(j) +e ∂y
y
Z  
1
(k) + x3 y ∂x
1 + x2
Z
(l) (y cos(xy)) ∂y

Los conocimientos anteriores, pueden ser empleados en la resolución de algunas ecuaciones diferenciales, por
tanto veamos de cuales se tratan.

• Ecuaciones diferenciales exactas: al calcular las derivadas parciales de una función de dos varia-
bles, se nota que ambas derivadas son, generalmente, a su vez funciones de dos variables, por tanto
∂ ∂
se pueden escribir como Z(x, y) = M (x, y) y Z(x, y) = N (x, y). Con ello el diferencial total de
∂x ∂y
dicha función se puede escribir como: dF = M (x, y) dx + N (x, y) dy.

Además, toda ecuación diferencial de primer orden puede escribirse como M (x, y) dx + N (x, y) dy = 0,
para ciertas funciones M y N . La relación entre el diferencial total y la expresión anterior dice que
dada una ecuación diferencial de la forma M (x, y) dx + N (x, y) dy = 0, si pudiéramos encontrar una
∂Z ∂Z
función Z tal que = M (x, y) y = N (x, y) obtendrı́amos la solución general de dicha ecuación
∂x ∂y
y corresponde a Z(x, y) = c.

17
Definición: una ecuación diferencial de la forma M (x, y) dx + N (x, y) dy = 0 se dice que es exacta si
existe una función Z(x, y) tal que

∂Z ∂Z
= M (x, y) y = N (x, y)
∂x ∂y

x2 y2
Ejemplo 1. La ecuación (x + y) dx + (x − y) dy = 0 es exacta por que existe Z(x, y) = + xy − ,
2 2
∂Z ∂Z
tal que = x+y y = x − y. Además, se tiene que la solución de dicha ecuación corresponde a
∂x ∂y
x2 y2
c= + xy −
2 2

Ejemplo 2. La ecuación (y cos(x) + 2xey ) dx + sen(x) + x2 ey + 2 dy = 0 es exacta por que existe

∂Z ∂Z
Z(x, y) = y sen(x) + x2 ey + 2y, tal que = y cos(x) + 2xey y = sen(x) + x2 ey + 2. Además, se
∂x ∂y
tiene que la solución de dicha ecuación corresponde a c = y sen(x) + x2 ey + 2y.

∂Z ∂Z
En estos ejemplos se logró determinar la función Z tal que = M (x, y) y = N (x, y); sin
∂x ∂y
embargo, en muchos de los casos, eso no es fácil de lograr. Además, generalmente se presenta una
ecuación diferencial y debe analizarse si esta es o no exacta, sin tener la función Z correspondiente;
por tanto, lo que debe hacerse es preguntar primero sobre la existencia de dicha función y en caso que
exista, el cómo poder determinarla; el siguiente resultado nos ayuda a resolver dicha situación.
∂M ∂N
Sean las funciones M , N , y , continuas en una región R. Entonces la ecuación diferencial dada
∂y ∂x
∂M ∂N
por M (x, y) dx + N (x, y) dy = 0 es exacta si y solo si = , en todo punto de R.
∂y ∂x
Mediante la resolución de algunos ejemplos, se empleará lo comentado y ası́ determinar la solución de
algunas ecuaciones diferenciales.

Ejemplo 1. Resolver la ecuación diferencial dada por (x − y cos(x)) dx − (y + sen(x)) dy = 0.

Solución: dado que está escrita de la forma M (x, y) dx + N (x, y) dy = 0, con M (x, y) = (x − y cos(x))
y N (x, y) = −(y + sen(x)), veamos si la ecuación diferencial es exacta.

∂M ∂(x − y cos(x))
= = − cos(x)
∂y ∂y
∂N ∂(−(y + sen(x)))
= = − cos(x)
∂x ∂x

∂M ∂N
Dado que = la ecuación es exacta.
∂y ∂x
Si una ecuación de la forma M (x, y) dx + N (x, y) dy = 0 es exacta, significa que existe una función
∂Z ∂Z
Z en dos variables, tal que su diferencial total dZ = dx + dy = M (x, y) dx + N (x, y) dy = 0.
∂x ∂y
Se observa además que para que Z cumpla dicha condición Z = c, con c constante. Ası́, Z = c es la
solución general de la ecuación dada.

18
 ∂Z
Para determinar la función Z, como la ecuación es exacta, se cumple entonces que = M (x, y) y
∂x
∂Z
= N (x, y). Lo cual nos permite encontrar fácilmente la solución. En el ejemplo:
∂y

∂Z ∂Z
= (x − y cos(x)) (∗) y = −(y + sen(x)) (∗∗)
∂x ∂y

x2
Z Z
∂Z
Usando (∗), = x − y cos(x) ⇒ ∂Z = (x − y cos(x)) ∂x ⇒ Z = − y sen(x) + C(y), para
∂x 2
alguna función C(y). (Se ha integrado de forma parcial)

∂Z
Para determinar C(y) se emplea (∗∗), ası́ derivando Z respecto y, se tiene = − sen(x) + C 0 (y) ⇒
∂y
−(y + sen(x)) = − sen(x) + C 0 (y) ⇒ −y − sen(x) + sen(x) = C 0 (y) ⇒ −y = C 0 (y).
−y 2
Z Z
0 dC
Como −y = C (y) ⇒ −y = ⇒ −y dy = dC ⇒ = C(y).
dy 2
x2 y2
Finalmente la solución general corresponde a c = − y sen(x) −
2 2

Lo anterior estable un procedimiento general para resolver una ecuación diferencial homogénea. Este
se resume a continuación.

Si una ecuación diferencial de la forma M (x, y) dx+N (x, y) dy = 0 es exacta, entonces existe una función
∂Z ∂Z
de dos variables, la cual podrı́a denotarse Z(x, y), tal que = M (x, y) y además, = N (x, y).
∂x ∂y
Z Z
∂Z
Usando que = M (x, y) ⇒ ∂Z = M (x, y) ∂x, ası́ al mantener a y como constante e integrar
∂x Z
respecto a x, se logra Z(x, y) = M (x, y) dx + C(y). Cabe indicar que esta función C(y) es arbitraria
∂Z
y juega el papel de constante de integración. (Pudo haberse elegido = N (x, y) e integrar respecto
∂y
a y).
Para determinar
Z la función C(y), puede derivarse
Z Z(x, y) respecto a y y lograr:
∂Z ∂ ∂Z ∂ ∂Z
= M (x, y) dx + C 0 (y) ⇒ = M (x, y) dx + C 0 (y). Pero como = N (x, y), al
∂y ∂y Z ∂y ∂y Z ∂y
∂M (x, y) 0 0 ∂
sustituir N (x, y) = dx + C (y) ⇒ C (y) = N (x, y) − M (x, y) dx. Ası́, al integrar
∂y Z  Z  ∂y
∂
respecto a y se tiene que C(y) = N (x, y) − M (x, y) dx dy.
∂y
De esta manera la solución general de la ecuación diferencial corresponde a:
Z Z  Z 
∂
k = M (x, y) dx + N (x, y) − M (x, y) dx dy
∂y

Resolvamos algunos ejemplos:

19
 1 y2
 
2y
Ejemplo 2. Resolver + dx − dy = 0.
x x2 x
Solución:
1 y2
 
2y
Como está en la forma M (x, y) dx + N (x, y) dy = 0, con M (x, y) = + y N (x, y) = − =
x x2 x
−2yx−1 .
∂ 2y ∂ 2y
Verifiquemos si es o no exacta: M (x, y) = 2 y N (x, y) = 2 . Por tanto es exacta.
∂y x ∂x x
∂Z 1 y2 ∂Z −2y
= + 2 y
Como es exacta debe existir una función Z en dos variables tal que = .
∂x x x ∂y x
−y 2
  Z  
−2y −2y −2y
Z
∂Z
Usando que = ⇒ ∂Z = ∂y ⇒ ∂Z = ∂y ⇒ Z(x, y) = + C(x) =
∂y x x x x
−y 2 x−1 + C(x).

∂Z y2 1 y2 y2 1
Derivando Z respecto a x, se tiene = 2 + C 0 (x) ⇒ + 2 = 2 + C 0 (x) ⇒ = C 0 (x)
Z Z ∂x x x x x x
1 dC 1
⇒ = ⇒ dx = dC ⇒ ln |x| + c = C(x).
x dx x
−y 2
Por tanto la solución general corresponde a k = + ln |x|.
x
Z   Z  
−2y
Z
1 2 −2 2y
Si se usa la fórmula general se tiene: k = +y x dx + − dx dy
x x x2
y2 y2
 
−2y 2y
Z
⇒ k = ln |x| − − + dy ⇒ k = ln |x| − .
x x x x
dy 2 + yexy
Ejemplo 3. Resolver = .
dx 2y − xexy
Solución: al escribirla de la forma M (x, y) dx + N (x, y) dy = 0 ⇒ (2y − xexy ) dy = (2 + yexy ) dx ⇒
(−2 − yexy ) dx + (2y − xexy ) dy = 0. Ası́ M (x, y) = (−2 − yexy ) y N (x, y) = (2y − xexy ).
Verifiquemos si es o no exacta.
∂ ∂
M (x, y) = −exy − yxexy y N (x, y) = −exy − xyexy = −exy − yxexy . Por tanto es exacta.
∂y ∂x
∂Z ∂Z
Como es exacta existe una función Z(x, y) tal que = (−2 − yexy ) y = (2y − xexy ).
∂x ∂y
Z Z
∂Z
Usando que = (−2 − ye ) ⇒ ∂Z = (−2 − ye ) ∂x ⇒ ∂Z = (−2 − yexy ) ∂x
xy xy
∂x
yexy
⇒ Z(x, y) = −2x − + C(y) = −2x − exy + C(y).
y
∂Z
Derivando Z respecto a y, se tiene = −xexy + C 0 (y) ⇒ (2y − xexy ) = −xexy + C 0 (y)
∂y
⇒ 2y = C 0 (y) ⇒ y 2 = C(y).

Por tanto la solución general corresponde a k = −2x − exy + y 2 .

20
Ejemplo 4. Resolver cos(x) − x sen(x) + y 2 dx + 2xy dy = 0, sujeta a la condición y(π) = 1

Solución: para este ejemplo se tiene que M (x, y) = cos(x)−x sen(x)+y 2 y N (x, y) = 2xy. Verifiquemos
si la ecuación diferencial es exacta
∂ ∂
M (x, y) = 2y y N (x, y) = 2y. Por tanto es exacta.
∂y ∂x
∂Z ∂Z
= cos(x) − x sen(x) + y 2 y


Como es exacta existe una función Z(x, y) tal que = 2xy).
∂x ∂y
Z Z
∂Z
Usando que = 2xy ⇒ ∂Z = (2xy) ∂y ⇒ ∂Z = (2xy) ∂y ⇒ Z(x, y) = xy 2 + C(x).
∂y
∂Z
= y 2 + C 0 (x) ⇒ cos(x) − x sen(x) + y 2 = y 2 + C 0 (x) ⇒


Derivando Z respecto a x, se tiene que
∂x
2y = C 0 (x) ⇒ (cos(x) − x sen(x)) = C 0 (x) ⇒ sen(x) + x cos(x) − sen(x) = C(x) ⇒ x cos(x) = C(x).

∴ la solución general corresponde a k = xy 2 + x cos(x)

Para la condición se tiene que y(π) = 1, ası́ k = (π) (1)2 + π cos(π) ⇒ k = 0. De esta manera, la
solución particular corresponde a 0 = xy 2 + x cos(x)
 
y
Ejemplo 5. Resolver dx + (ln(x − 1) + 2y) dy = 0, sujeta a la condición y(2) = 4
x−1  
y
Solución: para este ejemplo se tiene que M (x, y) = y N (x, y) = (ln(x − 1) + 2y). Verifique-
x−1
mos si la ecuación diferencial es exacta:
∂ 1 ∂ 1
M (x, y) = y N (x, y) = . Por tanto es exacta.
∂y x−1 ∂x x−1
 
∂Z y ∂Z
Como es exacta existe una función Z(x, y) tal que = y = ln(x − 1) + 2y.
∂x x−1 ∂y
Z Z
∂Z
Usando que = (ln(x − 1) + 2y) ⇒ ∂Z = (ln(x − 1) + 2y) ∂y ⇒ ∂Z = (ln(x − 1) + 2y) ∂y ⇒
∂y
Z(x, y) = y ln(x − 1) + y 2 + C(x).

∂Z y y y
Derivando Z respecto a x, se tiene que = + C 0 (x) ⇒ = + C 0 (x) ⇒ 0 = C 0 (x) ⇒
∂x x−1 x−1 x−1
c = C(x).

∴ la solución general corresponde a k = y ln(x − 1) + y 2

Para la condición se tiene que y(2) = 4, ası́ k = 4 ln(2 − 1) + (4)2 ⇒ k = 16. De esta manera, la solución
particular corresponde a 16 = y ln(x − 1) + y 2

dy 6x2 + y 3
Ejemplo 6. Resolver = , sujeta a la condición y(1) = 1
dx 5 − 3xy 2
Solución: expresada como M (x, y)dx + N (x, y)dy = 0, se tiene: −6x2 − y 3 dx + 5 − 3xy 2 dy = 0.

21
 ∂ ∂
−6x2 − y 3 = −3y 2 5 − 3xy 2 = −3y 2 .



Verifiquemos si la ecuación diferencial es exacta: y
∂y ∂x
Por tanto es exacta.
Z Z  Z 
∂M
Como es exacta la solución corresponde a: k = M (x, y) dx + N (x, y) − (x, y) dx dy
Z Z  Z  ∂y
−6x2 − y 3 dx + 5 − 3xy 2 − −3y 2 dx dy



⇒k=
Z Z
3 3
5 − 3xy + 3xy dy ⇒ k = −2x − xy + 5 dy ⇒ k = −2x3 − xy 3 + 5y.
2 2 3 3


⇒ k = −2x − xy +

Para la condición y(1) = 1, k = −2 (1)3 − (1) (1)3 + 5 (1) ⇒ k = −2 − 1 + 5 = 2. Por tanto la solución
particular es 2 = −2x3 − xy 3 + 5y

Ejemplo 7. Resolver x3 + xy 2
dx + x2 y + y 3 dy = 0

Solución: M (x, y) = x3 + xy 2 y N (x, y) = x2 y + y 3 . Verifiquemos si la ecuación diferencial es

∂ ∂
x3 + xy 2 = 2xy y x2 y + y 3 = 2xy. Por tanto es exacta.



exacta:
∂y ∂x
∂Z ∂Z
= x3 + xy 2 y = x2 y + y 3 .



Como es exacta existe una función Z(x, y) tal que
∂x ∂y
Z Z
∂Z
= x3 + xy 2 ⇒ ∂Z = x3 + xy 2 ∂x ⇒ ∂Z = x3 + xy 2 ∂x




Usando que
∂x
x4 x2 y 2
⇒ Z(x, y) = + + C(y).
4 2
Z Z
∂Z 2 0 2 3 2 0 3
Derivando Z respecto a y, se tiene = x y + C (y) ⇒ x y + y = x y + C (y) ⇒ y dy = dC
∂y
y4
⇒ = C(y).
4
x4 x2 y 2 y 4
2
∴ la solución general corresponde a k = + + ⇒ 4k = x2 + y 2
4 2 4

Ejemplo 8. Resolver 2x + y + (x + 2y)y 0 = 0

dy
Solución: al escribirla de la forma M (x, y) dx + N (x, y) dy = 0 se tiene 2x + y + (x + 2y) =0⇒
dx
(x + 2y) dy = (−2x − y) dx ⇒ (2x + y) dx + (x + 2y) dy = 0, ası́ M (x, y) = (2x + y) y N (x, y) = (x + 2y).
∂ ∂
Verifiquemos si la ecuación diferencial es exacta: (2x + y) = 1 y (x + 2y) = 1. Por tanto es
∂y ∂x
exacta.

∂Z ∂Z
Como es exacta existe una función Z(x, y) tal que
= (2x + y) y = x + 2y.
∂x ∂y
Z Z
∂Z
Usando que = (x + 2y) ⇒ ∂Z = (x + 2y) ∂y ⇒ ∂Z = (x + 2y) ∂y ⇒ Z(x, y) = xy + y 2 + C(x).
∂y

22
 Z Z
∂Z
Derivando Z respecto a x, se tiene = y + C 0 (x) ⇒ 2x + y = y + C 0 (y) ⇒ (2x) dx = dC
∂x
⇒ x2 = C(x).
∴ la solución general corresponde a k = x2 + xy + y 2

Ejercicios.

1. Determinar si la ecuación diferencial es exacta. Si lo es, hallar su solución general.

(a) yex dx + ex dy = 0
(b) (4x3 − 6xy 2 ) dx + (4y 3 − 6xy) dy = 0
 
2
(c) y dx + x + dy = 0
y
 y 0
(d) 1 + ln(x) + y = 1 − ln(x)
x
 dy  2
(e) 3y 2 + x cos(xy)
 
+ 3x + y cos(xy) = 0
dx

(f) e y − 2x dx + e2xy dy = 0
2x 2

(g) (2x + y + 1) dx + (x − 3y 2 + 2) dy = 0
(h) (ex−y + x) dx = ex−y + y dy

(i) (−2x − xy 2 ) dx + (4y + yx2 ) dy = 0

(j) 2xyy 0 = x2 − y 2
dy x − y cos(x)
(k) =
dx sen(x) + y
dy y(y − ex )
(l) = x
dx e − 2xy
(m) (1 + cos(x + y)) dx + cos(x + y) dy = 0
(n) 2xy dx + (x2 − sen(y)) dy = 0
(o) e2y − y cos(xy) dx + 2xe2y − x cos(xy) + 2y dy = 0

(p) x(6xy + 5) dx + (2x3 + 3y) dy = 0

(q) x ln(y) dx − (x + y ln(x)) dy = 0
(r) (2u + t + 2ut2 ) du + (u + 2u2 t) dt = 0
(s) u2 t du + (1 + ut2 ) dt = 0
2. Determine la solución particular de cada una de las siguientes ecuaciones diferenciales.
(a) 2xy dx + (x2 + 1) dy = 0, sujeta a y(1) = −3
yex + 2ex + y 2 dx + (ex + 2xy) dy = 0, sujeta a y(0) = 6


(b)
y 2 cos(x) − 3x2 y − 2x dx + 2y sen(x) − x3 + ln(x) dy = 0, sujeta a y(0) = e



(c)
(d) (3x2 y + 7xy 2 + 1) dx + (x3 + 7x2 y + 2y 2 ) dy = 0, sujeta a y(2) = 1
dy y 2 − 2x
(e) = , sujeta a y(1) = 2
dx 2y − 2xy
(f) (x2 + 2ye2x )y 0 + 2xy + 2y 2 e2x = 0, sujeta a y(1) = 1
(g) 3x2 + 2xy + 3y 2 + (x2 + 6xy)y 0 = 0, sujeta a y(1) = 2
(h) 1 + y cos(xy) + x cos(xy)y 0 = 0, sujeta a y(1) = 0

23
 3. Determine el valor de k, de manera tal que ecuación diferencial sea exacta
(a) y 3 + kxy 4 − 2x dx + 3xy 2 + 20x2 y 3 dy = 0

(b) 2x − y sen(xy) + ky 4 − 20xy 3 + x sen(xy) dy = 0

(c) (xy 2 + kx2 y) + (x + y)x2 y 0 = 0

(d) 2xyy 0 + y k = 0

• Ecuaciones diferenciales hechas exacta: algunas ecuaciones diferenciales no son exactas, sin em-
bargo, pueden multiplicarse por alguna función adecuada y hacerse exactas (dicha función recibe el
nombre de factor integrante). Por ejemplo, si se requiere resolver la ecuación dada por y 2 +(1+xy)y 0 = 0
∂ 2 ∂
se tiene que: y 2 dx + (1 + xy) dy = 0. Al verificar si es o no exacta (y ) = 2y y (1 + xy) = y.
∂y ∂x
Por tanto no es exacta.

1 1 2  1
Al multiplicar la ecuación por el factor integrante , se logra: y dx + (1 + xy) dy = · 0
  y y   y
1 1
⇒ y dx + + x dy = 0. En esta ecuación M (x, y) = y y N (x, y) = + x , al verificar si es o no
y   y
∂ ∂ 1
exacta se tiene: (y) = 1 y + x = 1. Que es exacta.
∂y ∂x y

La solución
Z esta dadaZ por:
 Z  Z Z   Z 
∂M 1
k = M (x, y) dx + N (x, y) − (x, y) dx dy ⇒ k = (y) dx + + x − 1 dx dy
Z    ∂y Z y
1 1
⇒ k = xy + + x − x dy ⇒ k = xy + dy ⇒ k = xy + ln |y|.
y y

Al igual que en el ejercicio anterior existen algunas ecuaciones diferenciales que no son exactas y que,
en algunos casos, pueden ser multiplicadas por un factor integrante, el cual se denotará µ, y lograr
hacerlas exactas. Pero ¿Cómo calcular dicho factor integrante?
∂M
Si una ecuación diferencial de la forma M (x, y) dx + N (x, y) dy = 0 no es exacta, considere =Ay
∂y
∂N
a = B y veamos algunos casos.
∂x
A−B
Caso 1. Si después de simplificar el cociente se obtiene una función constante o que dependa
Z  N Z ∂M − ∂N
A−B ∂y ∂x
dx dx
exclusivamente de x, entonces µ = e N =e N .

B−A
Caso 2. Si después de simplificar el cociente se obtiene una función constante o que dependa
M Z
∂N
Z
B−A ∂x
− ∂M
∂y
dy dy
exclusivamente de y, entonces µ = e M =e M .

24
 N (x, y) M (x, y)
Caso 3. Si existen m y n, números reales, tales que la expresión A − B = m −n
x y
entonces µ(x, y) = xm y n .

Caso 4. Si existen funciones

R R P (x) y Q(y) tales que la expresión A − B = N (x, y) P (x) − M (x, y) Q(y)
P (x) dx
entonces µ(x, y) = e e Q(y) dy .

Como en todos los casos el interés es solo determinar el factor integrante, al resolver las integrales para
dicho factor se asigna como valor de la constante de integración un cero por simplicidad, es decir c = 0.
Ejemplos. Para cada una de las siguientes ecuaciones diferenciales, compruebe que no son exactas y
determine un posible factor integrante.

x + y2
Ejemplo 1. y 0 = .
2xy
dy x + y2
Solución: = ⇒ 2xy dy = (x + y 2 ) dx ⇒ (−x − y 2 ) dx + 2xy dy = 0. Al verificar si es o no
dx 2xy
∂ ∂N
exacta se tiene: A = (−x − y 2 ) = −2y y B = (2xy) = 2y, que evidentemente no es exacta.
∂y ∂x
A−B −2y − 2y −4y −2
Para el cálculo de µ, haciendo = = = .
N 2xy 2xy x
Como esta expresión depende solo de x, un posible factor integrante corresponde a
R −2 −2 1
µ = e x dx ⇒ µ = e−2 ln |x| ⇒ µ = eln |x| ⇒ µ = x−2 = 2 .
x
2
Ejemplo 2. sen(y) dx + (sen (y) − x cos(y)) dy = 0.
∂ ∂N
Solución: ¿Será exacta? A = (sen(y)) = cos(y) y B = (sen2 (y) − x cos(y)) = − cos(y), que no es
∂y ∂x
exacta.
B−A − cos(y) − cos(y) −2 cos(y)
Para el cálculo de µ, haciendo = = = −2 cot(y). Como esta
M sen(y) sen(y)
expresión depende solo de y, un posible factor integrante corresponde a
R −2 1
µ = e −2 cot(y) dy ⇒ µ = e−2 ln(sen(y)) ⇒ µ = eln(sen(y)) ⇒ µ = (sen(y))−2 = .
sen2 (y)
Ejemplo 3. (2xy + y 4 ) dx + (3x2 + 6xy 3 ) dy = 0.
∂ ∂N
Solución: ¿Será exacta? A = (2xy + y 4 ) = 2x + 4y 3 y B = (3x2 + 6xy 3 ) = 6x + 6y 3 , que no es
∂y ∂x
exacta.

B−A (6x + 6y 3 ) − (2x + 4y 3 ) 4x + 2y 3 2(2x + y 3 ) 2

Para el cálculo de µ, haciendo = = = = . Dicha
M 2xy + y 4 R 2
2xy + y 4 y(2x + y 3) y
2
expresión solo depende de y, ası́ el factor integrante es µ = e y dy ⇒ µ = e2 ln y ⇒ µ = eln(y) ⇒ µ = y 2 .

Ejemplo 4. (6y − 24xy 5 ) dx + (9x − 56x2 y 4 ) dy = 0.

∂ ∂N
Solución: ¿Será exacta? A = (6y − 24xy 5 ) = 6 − 120xy 4 y B = (9x − 56x2 y 4 ) = 9 − 112xy 4 , que
∂y ∂x
no es exacta.

25
Para el cálculo de µ, consideremos:
9x − 56x2 y 4 6y − 24xy 5
   
4 4 4



A − B = 6 − 120xy − 9 − 112xy = −3 − 8xy = 1 · −2· =
x y
N (x, y) M (x, y)
1· −2· , ası́ el factor integrante corresponde a µ(x, y) = xy 2
x y

Además, dicho factor integrante pudo resolverse al considerar

∂M ∂N 1
2
= 6 − 120xy 4 − 9 − 112xy 4 = −3−8xy 4 = · 9x − 56x2 y 4 − · 6y − 24xy 5 =




A−B = −
∂y ∂x R 1 R 2
x y
dx dy ln(x) ln(y) 2
P (x) N (x, y) − Q(y) M (x, y), ası́ µ(x, y) = e x e y = e e = xy .

Ejemplo 5. (2x2 + e−y ) dx + (x3 + xy) dy = 0.

∂ ∂N 3
Solución: ¿Será exacta? A = (2x2 + e−y ) = −e−y y B = (x + xy) = 3x2 + y, que no es exacta.
∂y ∂x

Para el cálculo de µ, consideremos:

∂M ∂N −1
A−B = − = −e−y − (3x2 + y) = (x3 + xy) · − (2x2 + e−y ) · 1 = N (x, y) P (x) − M (x, y) Q(y).
∂y ∂x x
R −1 R ey
Por tanto el factor integrante corresponde a µ(x, y) = e x dx e 1 dy = e− ln(x) ey =
x

Dado que al multiplicar la ecuación diferencial por dicho factor integrante, ésta se convierte en exacta,
se emplea el método de ecuaciones exactas para resolver el ejercicio. Veamos algunos ejemplos:

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

∂M ∂N
Solución: verificar si es o no exacta, A = =1yB= = 2xy − 1. Que no es exacta.
∂y ∂x
A−B 1 − (2xy − 1) −2xy + 2 −2(xy − 1) −2
Cálculo del factor integrante: = 2
= 2 = = . Como dicha
N x y−x x y−x x(xy − 1) x
R −2 −2 1
expresión depende solo de x, entonces µ = e x dx = e−2 ln(x) = eln(x) = x−2 = 2 .
x

Multiplicando la ecuación diferencial por dicho factor integral se tiene

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

∂M 1 ∂N 1
¿Será exacta? = 2 y = 2 . Que es exacta.
∂y x ∂x x
∂Z y ∂Z 1
Como es exacta existe una función Z(x, y) tal que =2+ 2 y =y− .
∂x x ∂y x
y2
Z Z  
∂Z 1 1 y
Usando que = y − ⇒ ∂Z = y− dy ⇒ Z(x, y) = − + C(x). Derivando Z(x, y)
∂y x x 2 x
∂Z y y y
respecto a x, se tiene = 2 + C 0 (x) ⇒ 2 + 2 = 2 + C 0 (x) ⇒ 2 = C 0 (x) ⇒ 2x = C(x).
∂x x x x

26
 y2 y
Por tanto la solución general de la ecuación corresponde a: k = − + 2x
2 x
π
Ejemplo 2. Resolver sen(y) dx + (sen2 (y) − x cos(y)) dy = 0, sujeta a y(0) =
2
∂M ∂N
Solución: verificar si es o no exacta, A = = cos(y) y B = = − cos(y). Que no es exacta.
∂y ∂x
A−B cos(y) + cos(y) 2 cos(y)
Cálculo del factor integrante: = 2
= 2
. Como dicha expresión
N sen (y) − x cos(y) sen (y) − x cos(y)
no depende x, entonces utilizaremos el otro caso.

B−A − cos(y) − cos(y) −2 cos(y)

= = = −2 cot(y). Como dicha expresión depende solo de y, en-
M sen(y) sen(y)
R 1
tonces µ = e −2 cot(y) dy = e−2 ln(sen(y)) = sen−2 (y) = .
sen2 (y)

Multiplicando la ecuación diferencial por dicho factor integral se tiene

   
1  2
 1 x cos(y)
sen(y) dx + (sen (y) − x cos(y)) dy = 0 ⇒ dx + 1 − dy = 0
sen2 (y) sen(y) sen2 (y)

∂M − cos(y) ∂N − cos(y)
¿Será exacta? = 2
y = . Que es exacta.
∂y sen (y) ∂x sen2 (y)
∂Z 1 ∂Z x cos(y)
Como es exacta existe una función Z(x, y) tal que = y =1− .
∂x sen(y) ∂y sen2 (y)
Z Z  
∂Z 1 1 x
Usando que = ⇒ ∂Z = dx ⇒ Z(x, y) = + C(y). Derivando Z(x, y)
∂x sen(y) sen(y) sen(y)
∂Z −x cos(y) x cos(y) −x cos(y)
respecto a y, se tiene = 2
+ C 0 (y) ⇒ 1 − 2
= 2
+ C 0 (y) ⇒ 1 = C 0 (y) ⇒
∂y sen (y) sen (y) sen (y)
y = C(y).

x
Por tanto la solución general de la ecuación corresponde a: k = + y. Al utilizar la condición
sen(y)
π 0 π π π x
x=0yy = ;k = π

+ ⇒ k = . Por tanto la solución particular corresponde a = +y
2 sen 2 2 2 2 sen(y)
2y
Ejemplo 3. Resolver y 0 =
x3
x + sen(2y)
2y dy 2y
Solución: y 0 = 3
⇒ = 3
⇒ −2y dx + (x + x3 sen(2y)) dy = 0, al verificar
x + x sen(2y) dx x + x sen(2y)
∂M ∂N
si es o no exacta, A = = −2 y B = = 1 + 3x2 sen(2y). Que no es exacta.
∂y ∂x

Cálculo del factor integrante:

A−B −2 − (1 + 3x2 sen(2y)) −3 − 3x2 sen(2y) −3(1 + x2 sen(2y)) −3
= 3
= 3
= 2
= . Como dicha ex-
N x + x sen(2y) x + x sen(2y) x(1 + x sen(2y)) x
R −3 −3 1
presión depende solo de x, entonces µ = e x dx = e−3 ln(x) = eln(x) = 3 .
x

27
Multiplicando la ecuación diferencial por dicho factor integral se tiene
   
1  3
 −2y 1
−2y dx + (x + x sen(2y)) dy =0⇒ dx + + sen(2y) dy = 0
x3 x3 x2

∂M −2 ∂N −2
¿Será exacta? = 3 y = 3 . Que es exacta.
∂y x ∂x x
∂Z −2y ∂Z 1
Como es exacta existe una función Z(x, y) tal que = 3 y = 2 + sen(2y).
∂x x ∂y x
Z  
−2y −2y
Z
∂Z y
Usando que = 3 ⇒ ∂Z = 3
dx ⇒ Z(x, y) = 2 + C(y). Derivando Z(x, y) respecto
∂x x x x
∂Z 1 1 1 − cos(2y)
a y, se tiene = 2 + C 0 (y) ⇒ 2 + sen(2y) = 2 + C 0 (y) ⇒ sen(2y) = C 0 (y) ⇒ = C(y).
∂y x x x 2
y cos(2y)
Por tanto la solución general de la ecuación corresponde a: k = 2
− .
x 2
 y  dx
Ejemplo 4. Resolver = ln(x) − y 3
x dy
y y
Solución: dx = (ln(x) − y 3 ) dy ⇒ dx + (− ln(x) + y 3 ) dy = 0, al verificar si es o no exacta,
x x
∂M 1 ∂N −1
A= = yB= = . Que no es exacta.
∂y x ∂x x
−1 1 −2
B−A x
− x x −2x −2
Cálculo del factor integrante: = y = y = = . Como dicha expresión depende
M x x
xy y
R −2
dy −2 1
solo de y, entonces µ = e y = e−2 ln(y) = eln(y) = .
y2

Multiplicando la ecuación diferencial por dicho factor integral se tiene

   
1 hy 3
i 1 − ln(x)
dx + (− ln(x) + y ) dy = 0 ⇒ dx + + y dy = 0
y2 x xy y2

∂M −1 ∂N −1
¿Será exacta? = 2 y = 2 . Que es exacta.
∂y xy ∂x xy
∂Z 1 ∂Z − ln(x)
Como es exacta existe una función Z(x, y) tal que = y = + y.
∂x xy ∂y y2
Z Z  
∂Z 1 1 ln(x)
Usando que = ⇒ ∂Z = dx ⇒ Z(x, y) = + C(y). Derivando Z(x, y) respecto
∂x xy xy y
∂Z − ln(x) 0 − ln(x) − ln(x) 0 0 y2
a y, se tiene = + C (y) ⇒ + y = + C (y) ⇒ y = C (y) ⇒ = C(y).
∂y y2 y2 y2 2
ln(x) y 2
Por tanto la solución general de la ecuación corresponde a: k = + .
y 2

28
Ejercicios.

1. Determine un factor integrante, si depende de solo una variable, para cada una de las siguientes
ecuaciones diferenciales.
(a) (x + 2y 2 ) dx + xy dy = 0
(b) 2y dx + (x + y) dy = 0
(c) 2y dy + (x + y) dx = 0
(d) (5y − x2 )dx + x dy = 0
(e) (x + y − xy) dx + x dy = 0
(f) (u + 5t) du + 3u dt = 0
(g) (u + 5t) du + 3t dt = 0
(h) xy 2 dx + (1 + x2 y + x2 y 2 ) dy = 0
(i) y dx − (x + 6y 2 ) dy = 0
(j) (2x3 + y) dx − x dy = 0
(k) (x2 + 2x + y) dx + 2 dy = 0
(l) y 2 dx + (xy − 1) dy = 0
2. Usar el factor integrante dado para resolver cada una de las siguientes ecuaciones diferenciales.
(a) (4x2 y + 2y 2 ) dx + (3x3 + 4xy) dy = 0, con µ(x, y) = xy 2
(b) (3y 2 + 5x2 y) dx + (3xy + 2x3 ) dy = 0, con µ(x, y) = x2 y
1
(c) y dx − x dy = 0, µ(x) = 2
x
1
(d) y dx − x dy = 0, µ(x, y) =
xy
3. Determinar la solución general de cada una de las siguientes ecuaciones diferenciales.
(a) (y 2 e2x + 1) dx + (2ye2x − ex ) dy = 0
(b) (4x3 − 6xy 2 ) dx + (4y 3 − 6xy) dy = 0
(c) (x2 + y 2 + x) dx + xy dy = 0
y
(d) y 0 =
x − y2 − 1
(e) (x4 ln(x) − 2xy 3 ) dx + (3x2 y 2 ) dy = 0
(f) (x3 + x + y) dx − x dy = 0
(g) (y ln(y) + yex ) dx + (x + y cos(y)) dy = 0
(h) ex dx + (ex cot(y) + 2y csc(y)) dy = 0
(i) (x + 3y 2 ) dx + 2xy dy = 0
(j) (4xy + 3y 2 − x) dx + x(x + 2y) dy = 0
(k) y(x + y + 1) dx + x(x + 3y + 2) dy = 0
4. Determinar la solución general de cada una de las siguientes ecuaciones diferenciales, sujetas a las
condiciones iniciales dadas.
(a) x dx + x2 y + 4y dy = 0, con y(4) = 0

(b) x2 + y 2 − 5 dx = (y + xy) dy, con y(0) = 0

y − 2x3
(c) y 0 = , con y(2) = −4
x
29
 −y 3 − 2yex
(d) y 0 = , con y(0) = 1
ex + 3y 2
(e) y dx + (3x − y 3 ) dy = 0, con y(1) = 1

• Ecuaciones diferenciales lineales de primer orden: en términos generales, una ecuación diferen-
cial lineal de primer orden se puede expresar como a0 (x)y 0 + a1 (x)y = f (x) con a0 (x) 6= 0. Si se divide
a1 (x) f (x)
por el coeficiente de la primera derivada, es decir por a0 (x), se tiene y 0 + y= , lo que nos
a0 (x) a0 (x)
lleva a la siguiente definición.

Definición: una ecuación diferencial lineal de primer orden es aquella que puede escribirse de la forma
y 0 + P (x)y = Q(x); donde P (x) y Q(x) son funciones continuas en la variable x en la región donde esté
definida.

Solución de una ecuación diferencial lineal de primer orden: se puede verificar que la ecuación
diferencial lineal de primer orden dada Rpor y 0 + P (x)y = Q(x), generalmente no es exacta y que su
factor integrante corresponde a µ(x) = e P (x) dx .

Ası́ al multiplicar la ecuación diferencial por dicho factor integrante se tiene:

R R R R R
e P (x) dx
[y 0 + P (x) y] = Q(x) e P (x) dx
⇒ y0e P (x) dx
+e P (x) dx
P (x)y = Q(x) e P (x) dx
.

Sin embargo, se puede observar que el término de la izquierda corresponde a la derivada respecto a x,
del producto de la variable dependiente y por el factor integrante. Es decir;

d h R P (x) dx i R R R R
ye = y 0 e P (x) dx + e P (x) dx P (x) y = y 0 e P (x) dx + y e P (x) dx P (x)
dx
Z
d h R P (x) dx i R
P (x) dx
R
P (x) dx
R
Con lo cual se tiene ye = Q(x) e ⇒ ye = Q(x)e P (x) dx dx
R
Zdx R

− P (x) dx P (x) dx
⇒y=e · Q(x)e dx Que corresponde a la solución general.

Nótese que si Q(x) = 0, la fórmula anterior

R es válida y en este caso se tiene que la solución general de la
ecuación y 0 + P (x)y = 0 es y = c e− P (x) dx . Además, a dicha ecuación se le llama ecuación diferencial
lineal homogénea.
Ejemplo 1. Resolver la ecuación x dy = (x sen(x) − y) dx.
Solución: al escribirla de la forma y 0 + P (x)y = Q(x), se tiene
dy dy x sen(x) − y y 1 1
x = (x sen(x) − y) ⇒ = ⇒ y 0 = sen(x) − ⇒ y 0 + y = sen(x), donde P (x) =
dx dx x x x x
y Q(x) = sen(x). Por tanto la solución corresponde a:
Z  Z  Z 
−
R 1
dx
R 1
dx − ln(x) ln(x) 1
y=e sen(x) e
x dx ⇒ y = e sen(x) e
x x sen(x) dx ⇒ dx ⇒ y =
 Z  x
1 1 sen(x) c
y= −x cos(x) + cos(x) dx ⇒ y = [−x cos(x) + sen(x) + c] ⇒ y = − cos(x) + + .
x x x x

30
 dy
Ejemplo 2. Resolver la ecuación − 3y = 0.
dx 0
Solución: al escribirla de la forma y + P (x)y = Q(x), se tiene y 0 − 3y = 0. Por tanto la solución
corresponde a:
R
y = ce− −3 dx = ce3x . Nótese además
Z que pudoZ haberse resuelto por el método de variables separables
dy 1
ası́ − 3y = 0 ⇒ dy = 3y dx ⇒ dy = 3 dx ⇒ ln(y) = 3x + k ⇒ y = e3x+k = ce3x , con c = ek .
dx y
dy
Ejemplo 3. Resolver la ecuación x − 4y = x6 sen(x).
dx
4 x6 sen(x) 4
Solución: y 0 +P (x)y = Q(x), se tiene xy 0 −4y = x6 sen(x) ⇒ y 0 − y = ⇒ y 0 − y = x5 sen(x),
x x x
−4 5
donde P (x) = y Q(x) = x sen(x). Por tanto la solución corresponde a:
x
Z R  Z 
− −4 −4
R
dx dx 5 4 ln(x) −4 ln(x) 5
y=e x e x x sen(x) dx ⇒ y = e e x sen(x) dx
Z  Z  Z 
ln(x)4 ln(x)−4 5 4 −4 5 4
⇒y=e e x sen(x) dx ⇒ y = x x x sen(x) dx ⇒ y = x x sen(x) dx
⇒ y = x4 [−x cos(x) + sen(x) + c] ⇒ y = −x5 cos(x) + x4 sen(x) + c x4 .

dy
Ejemplo 4. Resolver la ecuación + y = x, con y(0) = 4.
dx
Solución: y + P (x)y = Q(x), se tiene y 0 + y = x, con lo cual P (x) = 1 y Q(x) = x. Por tanto la
0

solución corresponde a:
Z  Z 
−
R R
−x c
y=e 1 dx
e 1 dx
x dx ⇒ y = e e x dx ⇒ y = e−x [xex − ex + c] ⇒ y = x − 1 + x .
x
e
c c
Para la condición se tiene que si x = 0, entonces y = 4. Ası́ y = x − 1 + x ⇒ 4 = 0 − 1 + 0
e e
5
⇒ 4 + 1 = c ⇒ 5 = c. Por tanto la solución particular corresponde a y = x − 1 + x .
e

dy 1, 0 ≤ x ≤ 1
Ejemplo 5. Resolver la ecuación + 3y = f (x), con y(0) = 0, donde f (x) = .
dx 0, x>1
Solución: como f (x) es una función discontinua, se analizará cada una de las partes en la cual es con-
tinua y ası́ resolver la ecuación diferencial.
Si 0 ≤ x ≤ 1, entonces f (x) = 1, ası́ debe resolverse y 0 + 3y = 1, con lo cual P (x) = 3 y Q(x) = 1. Por
tanto la solución corresponde a:

e3x
Z  Z   
−
R
3 dx
R
3 dx −3x 3x −3x 1 c
y = e e · 1 dx ⇒ y = e e dx ⇒ y = e + c ⇒ y = + 3x . Como
3 3 e
1 c −1 1 1
y(0) = 0, entonces 0 = + 3·0 ⇒ c = . Ası́ y = − 3x
3 e 3 3 3e

31
Si x > 1, entonces f (x) = 0, ası́ debe resolverse y 0 + 3y = 0, y la solución corresponde a:
R
Z R 
− 3 dx
y=e e 3 dx
· 0 dx ⇒ y = c e−3x .

1
− 3e13x , 0 ≤ x ≤ 1

De esta manera la solución de dicha ecuación diferencial corresponde a: y = 3 .
ce−3x , x>1

Ejemplo 6. Resolver la ecuación y dx = (yey − 2x) dy.

Solución: para este caso resulta más práctico escribir la ecuación de la forma x0 + P (y)x = Q(y), ası́
dx yey − 2x dy 2x 2 2
= ⇒ = ey − ⇒ x0 + x = ey , con lo cual P (y) = y Q(y) = ey . Por tanto la
dy y dx y y y
solución corresponde a:
Z  Z  Z 
2 2
R R
− dy dy y −2 ln(y) 2 ln(y) y −2 2 y
x=e y e y e dy ⇒ x = e e e dy ⇒ x = y y e dy
2ey 2ey c
⇒ x = y −2 y 2 ey − 2yey + 2ey + c ⇒ x = ey −
 
+ 2 + 2.
y y y

Ejemplo 7. Una barra con una temperatura de 100◦ C se pone en un cuarto cuya temperatura constante
es de 0◦ C. Si después de 20 minutos la temperatura de la barra es de 50◦ C. Determine:

1. El tiempo que empleará la barra para llegar a una temperatura de 25◦ C.

2. La temperatura de la barra después de 10 minutos.

Solución:
dT
De acuerdo con la ley de Newton = −k(T − Tm ) ⇒ T 0 + kT = kTm . Que corresponde a una
dt
ecuación diferencial lineal de primer orden.
0
Como Tm = 0, se Z debe resolver
 T + kT = 0 cuya solución esta dada por
R R
T (t) = e− k dt
e k dt
· 0 dt = ce−kt .
Además, si T = 100 cuando t = 0, al sustituir se logra T = 100e−kt . Luego, dado que para t = 20
entonces T = 50, con ello se obtiene que 50 = 100e−20t ⇒ k = 0, 035.

Ası́ en cualquier instante t, la temperatura de la barra estará determinada por T (t) = 100e−0,035t

Sustituyendo T = 25, se tiene que la barra obtiene una temperatura de 25◦ C a los 39,6 minutos; además,
cuando t = 10, la temperatura es de 70,5◦ C.

Existen ecuaciones diferenciales que no son lineales pero se pueden transformar en lineales. Una de
estas es la denominada Ecuación de Bernoulli, la cual tiene la forma y 0 + P (x)y = Q(x)y n , con n 6= 0
y n 6= 1. Dicha ecuación puede transformarse en lineal al hacer el cambio v = y 1−n .

32
Ejercicios.

1. Determinar la solución general de cada una de las siguientes ecuaciones diferenciales.

dy
(a) + y = e3x
dx
dy
(b) x − 4y = x6 ex
dx
dy
(c) + 5y = 50
dx
dy 10y
(d) + = 10
dx 2x + 5
(e) y 0 + y cos(x) = sen(x) cos(x)
(f) y 0 + ex y = e2x
(g) xy 0 + (2 + 3x)y = xe−3x
dy
(h) = x2 e−4x − 4y
dx
(i) xy 0 + 2y = x−3
(j) xy 0 + 2y = sen(x)
dx
(k) y + 2x = 5y 3
dy
dy
(l) x + 3y + 2x2 = x3 + 4x
dx
dy 1
(m) = 2y
dx e + 2x
dy 1
(n) = 2
⇒ x + y 2 = x0 ⇒ x0 − x = y 2
dx x+y
2 0
(o) x y + 2xy = senh(3x)
(p) y 0 + 2xy = x3
(q) xy 0 + 2y = ex + ln(x)
dI
(r) L + RI = E, con L, R y E constantes.
dt
(s) 2xy 0 − 3y = 4x2
1
(t) y 0 = y
e +x
2. Determinar la solución general de cada una de las siguientes ecuaciones diferenciales, sujetas a las
condiciones indicadas.
dy
(a) x + y = 2x, con y(1) = 0
dx
(b) y 0 = 2y + x e3x − e2x , con y(0) = 2

(c) xy 0 + y = ex , con y(1) = 2

dy
(d) (x + 1) + y = ln(x), con y(1) = 0
dx
dx
(e) y − x = 2y 2 , con y(1) = 5
dy
dT
(f) = k(T − tm ); con T (0) = T0 ; k, Tm y T0 constantes.
dt
33
 dT
(g) = k(T − 50), con k constante y T (0) = 200
dt
(h) xy 0 + y = ex , con y(−1) = 4
(i) xy 0 + 2y = 4x2 , con y(1) = 2
(j) y 0 − y = 2xe2x , con y(0) = 1

dy 1, 0 ≤ x ≤ 3
(k) + 2y = f (x), con y(0) = 0, donde f (x) = .
dx 0, x>3

dy 1, 0 ≤ x ≤ 1
(l) + y = f (x), con y(0) = 1, donde f (x) = .
dx −1, x>1

dy x, 0 ≤ x ≤ 1
(m) + 2xy = f (x), con y(0) = 2, donde f (x) = .
dx 0, x>1

3. Resuelva los siguientes problemas.

(a) Un marca pasos consiste en un interruptor, una baterı́a de voltaje constante E0 , un capacitor
con capacitancia constante C y el corazón como resistor con resistencia constante R. Cuando
se cierra el interruptor, se carga el capacitor; cuando se abre el interruptor, se descarga el
capacitor, enviando un estı́mulo eléctrico al corazón. Durante el tiempo que el corazón está
dE −1
siendo estimulado, el voltaje en el corazón satisface la ecuación diferencial lineal = E.
dt RC
Resolver dicha ecuación diferencial sujeta a E(4) = E0 . (Ejercicio tomado del libro Ecuaciones
diferenciales con problemas de valores en la frontera de Dennis Zill, página 66).
(b) Según la ley de enfriamiento (o calentamiento) de Newton: “la rapidez con la que cambia
la temperatura de un cuerpo es proporcional a la diferencia de temperaturas entre el cuerpo
y el medio ambiente que lo rodea”. Si T = T (t) representa la temperatura del cuerpo en
el instante t y Tm representa la temperatura del medio ambiente, entonces se satisface que
dT
= k(T − Tm ), con k constante de proporcionalidad. Usando lo anterior resuelva los
dt
siguientes problemas.
(a) Se toma un termómetro de una habitación donde la temperatura es de 70◦ C y se lleva al
exterior, donde la temperatura del aire es de 10◦ C. Después de medio minuto el termómetro
marca 50◦ C ¿Cuál es la lectura del termómetro en t = 1 minuto? ¿Cuánto tarda el
termómetro en alcanzar 15◦ C?
(b) Suponga que la temperatura T de un pequeño objeto caliente situado en un medio de
temperatura Ta decrece a razón proporcional T − Ta . Si la temperatura del objeto pasa de
100◦ C a 80◦ C en 20 minutos, cuando la temperatura del medio es de 2◦ C ¿Cuánto tardará
en enfriarse de 100◦ C a 60◦ C?
(c) Un cuerpo que tiene una temperatura de 70◦ C es depositado en un lugar donde la tempe-
ratura se mantiene a 40◦ C. Después de 3 minutos la temperatura del cuerpo ha disminuido
a 60◦ C ¿Cuál es la temperatura del cuerpo después de 5 min? ¿Cuánto tiempo pasará
para que el cuerpo tenga 50◦ C?
(d) Un objeto que tiene una temperatura 50◦ F se coloca a las 10 : 00 horas en un horno
que se mantiene a 375 ◦ F. A las 11 : 15 horas su temperatura era 125◦ F ¿A qué hora la
temperatura del objeto será de 150◦ F?
(e) Una taza de café cuya temperatura es 40◦ C se coloca en un cuarto cuya temperatura es
18◦ C. Dos minutos más tarde la temperatura del café es 20◦ C ¿Después de cuánto tiempo
la temperatura del café será de 10◦ C?

34
 (f) La temperatura de un motor en el momento en que se apaga es de 200◦ C y la temperatura
del aire que lo rodea es de 30◦ C. Después de 10 min la temperatura del motor ha bajado
a 180◦ C ¿Cuánto tiempo transcurrirá para que la temperatura del motor disminuya hasta
40◦ C?
(g) Un cuerpo caliente se enfrı́a a un ritmo proporcional a la diferencia de temperatura respecto
al ambiente. Si el cuerpo se calienta a 110◦ C y se coloca en aire con temperatura de 10◦ C.
Tras una hora la temperatura del cuerpo es de 60◦ C ¿Cuánto tiempo tarda en llegar a
30◦ C?
(h) Si la temperatura del aire es de 10◦ C y un cuerpo se enfrı́a desde 80◦ C hasta 30◦ C en 30
minutos ¿Cuánto tiempo tardará en descender su temperatura hasta 20◦ C?
(c) Para un circuito en serie que contiene sólo un resistor y un inductor, establece la segunda ley
di
de Kirchhoff que la suma de las caı́das de voltaje en el inductor L y la caı́da de voltaje en el
dt
resistor iR es la misma que el voltaje impreso E(t) en el circuito. Ası́, se obtiene la ecuación
di
diferencial lineal para la corriente corresponde a L + Ri = E(t), donde L y R son constantes
dt
que se conocen como la inductancia y la resistencia respectivamente.

q(t)
La caı́da de voltaje en un capacitor con capacitancia C está dado por , donde q es la carga
C
1
en el capacitor. Por tanto, para el circuito en serie se satisface que Ri + q = E(t).
C
dq
Además, la corriente i y la carga q se relacionan mediante i = , con lo cual se obtiene la
dt
dq 1
ecuación diferencial dada por R + q = E(t).
dt C
Usando los elementos anteriores resuelva los siguientes ejercicios
1
(a) Una baterı́a de 12 volts se conecta a un circuito en serie en el que la inductancia es de
2
henry y la resistencia es de 10 ohms. Determine la corriente i si la corriente inicial es cero.
(b) Una baterı́a de 30 volts se conecta a un circuito en serie en el que la inductancia es de
0, 1 henry y la resistencia es de 50 ohms. Determine la corriente i si la corriente inicial es
cero. Además, calcule la corriente cuando el tiempo tiende a infinito.
500t

3e
Respuesta: I(t) = e−500t +c
5
(c) Una baterı́a de 5 volts se conecta a un circuito en serie en el que la inductancia es de 1
henry y la resistencia es de 50 ohms. Determine la corriente i si la corriente inicial es cero.
Además, calcule la corriente
 50t cuando
 el tiempo tiende a infinito.
e
Respuesta: I(t) = e−50t +c
10
(d) Una baterı́a de 60 volts se conecta a un circuito en serie en el que la inductancia es de 4
henry y la resistencia es de 12 ohms. Determine la corriente i si la corriente inicial es cero.
Además, calcule la corriente cuando
 el tiempo tiende a infinito.
Respuesta: I(t) = e−3t 20e3t + c


35
 (e) Una baterı́a de 4 sen(t) volts se conecta a un circuito en serie en el que la inductancia es
de 4 henry y la resistencia es de 100 ohms. Determine la corriente i si la corriente inicial
es cero.  25t 
−25t e (25 sen(t) − cos(t))
Respuesta: I(t) = e +c
626
(f) Una baterı́a de 3 sen(2t) volts se conecta a un circuito en serie en el que la inductancia es
de 0, 5 henry y la resistencia es de 10 ohms. Determine la corriente i si la corriente inicial
es cero.  20t 
−20t 3e (20 sen(2t) − 2 cos(2t))
Respuesta: I(t) = e +c
404
(g) Si se aplica una fuerza electromotriz de 30 voltios a un circuito LR en serie en el que
la inductancia es 1 henry y la resistencia es de 50 ohms. Determine la corriente i(t) si
i(0) = 0.
(h) Se aplica una fuerza electromotriz a un circuito en serie en el que la resistencia es de
200 ohms y la capacitancia es de 10−4 farads. Determine la carga q(t) en el capacitor se
q(0) = 0. Además, encuentre la corriente i(t).

Ecuaciones diferenciales lineales de orden superior.

Ecuación diferencial lineal de enésimo orden.

Definición: una ecuación diferencial lineal de orden n es toda aquella de la forma

a0 (x)y n + a1 (x)y n−1 + a2 (x)y n−2 + ... + an−2 (x)y 00 + an−1 (x)y 0 + an (x)y = F (x)

con a0 (x) 6= 0. Además, todos los ai (x) y F (x) dependen solo de x y son continuas en un intervalo I.

Ejemplos

1. Para n = 1, a0 (x)y 0 + a1 (x)y = F (x), que corresponde a una ecuación diferencial lineal de primer orden,
equivalente a y 0 + P (x)y = Q(x).

2. Para n = 2, a0 (x)y 00 + a1 (x)y 0 + a2 (x)y = F (x), que corresponde a una ecuación diferencial lineal de
segundo orden.

3. y 000 + 4y 00 + 3y 0 + 5 = sen(x).

4. y IV + 3x2 y 00 − 2xy = x2 + cos(3x).

Además, a una ecuación diferencial lineal de orden n se le denomina:

(a) De coeficientes constantes, si todos los ai (x) con i = 0, 1, 2, 3, ..., n son constantes reales.
Ejemplos:

(a) 5y 00 − 3y 0 + 2y = 5 cos(2x).
(b) y 000 + 2y 00 − 4y = 7.
(c) 3y 00 − 6y 0 + y = 0.

(b) De coeficientes variables, si al menos existe un ai (x) con i = 0, 1, 2, 3, ..., n que depende de x.

36
 Ejemplos:

(a) xy 00 − 3y 0 + 2y = 5 cos(2x) + ln(x).

(b) y 000 + 3xy 00 − 4y = 7x3 .
(c) sen(x)y 00 − 6xy 0 + y = 3x2 .
(d) 2y IV + 2y 00 − 4xy = 3

Solución de una ecuación lineal de orden n con coeficientes constantes

Para resolver una ecuación de la forma a0 y n + a1 y n−1 + a2 y n−2 + ... +
an−1 y 0 + an y = F (x), primero se expresa
en la forma φ(D)y = a0 Dn + a1 Dn−1 + a2 Dn−2 + ... + an−1 D + an y = F (x) con Dn la enésima derivada
de y.

La solución general de la ecuación φ(D)y = F (x) está dada por y = yc + yp , donde yc es la solución
complementaria de la ecuación homogénea φ(D)y = 0, mientras que yp es una solución particular de la
ecuación no homogénea φ(D)y = F (x).

Solución complementaria de una ecuación homogénea

Para determinar la solución de una ecuación de la forma φ(D)y = 0, analicemos primero la solución de la
ecuación diferencial y 0 − 2y = 0, se puede verificar fácilmente que su solución general es y = ce2x . Lo anterior
permite suponer que toda ecuación homogénea de la forma φ(D)y = 0, con coeficientes constantes, tiene por
solución particular una función exponencial de la forma y = emx , donde m es un parámetro por determinar.

Pero, puede observarse que y 0 = memx y si y = emx es una solución de y 0 − 2y = 0 ⇒ m emx − 2emx = 0 ⇒
emx (m − 2) = 0 ⇒ m = 2, con lo cual y = e2x es una solución particular de la ecuación y 0 − 2y = 0 y su
respectiva solución es y = ce2x .

Luego, al usar la notación (D − 2)y = 0, se puede observar que la solución principal de la ecuación propuesta,
y por tanto su solución general, se puede obtener del cero de la ecuación φ(m) = 0 (llamada ecuación auxiliar)
donde φ(m) = 0 tiene la misma forma que φ(D) = 0.

Ecuación diferencial lineal homogénea de segundo orden

Considere la ecuación diferencial lineal homogénea de segundo orden dada por a0 y 00 + a1 y 0 + a2 = 0, con
a0 , a1 , a2 constantes, a0 6= 0.

Dicha ecuación en notación de operadores corresponde a (a0 D2 +a1 D +a2 )y = 0, ası́ φ(D) = a0 D2 +a1 D +a2 .
Si se sabe que φ(m) tiene la misma forma que φ(D), entonces la respectiva ecuación auxiliar o ecuación ca-
racterı́stica está dada por φ(m) = a0 m2 + a1 m + a2 = 0. Si y = emx , con m número real, es una solución
particular de la ecuación, entonces su solución complementaria yc se obtiene de las raı́ces de la ecuación
auxiliar φ(m) = 0.

Dado que φ(m) = 0 es una ecuación de segundo grado y los ceros de ella dependen del valor del discriminante
4 = a21 − 4a0 a2 , se debe considerar tres casos:

37
 (1) Si 4 > 0 entonces φ(m) = 0 tiene dos raı́ces reales distintas, sean m1 y m2 . Con lo cual y1 = em1 x
y y2 = em2 x son soluciones linealmente independientes de la solución original, por tanto la solución
complementaria o general es y = yc = c1 y1 + c2 y2 = c1 em1 x + c2 em2 x

(2) Si 4 = 0 entonces φ(m) = 0 tiene dos raı́ces reales iguales, sean m1 . Ası́ solo se puede obtener una
−a1 −a1
solución exponencial y1 = em1 x . Por la fórmula cuadrática se sabe que m1 = o 2m1 = . Ası́ las
2a0 a0
Z − a R a1 −a 1 x
−a1
Z
e 0 e a0
soluciones de la ecuación son y1 y y2 donde y2 = y1 2
dx = y1 2
pero como 2m1 =
Z 2m1 x Z y1 y1 a0
e
se tiene que em1 x dx = em1 x dx = xem1 x . Por lo anterior y1 = em1 x y y2 = xem1 x son
e2m1 x
soluciones linealmente independientes de la ecuación, de esta manera la solución complementaria es
y = yc = c1 em1 x + c2 xem1 x .

(3) Si 4 < 0, las soluciones corresponden a números complejos no reales de la forma a ± bi, (con i la unidad
imaginaria y a, b ∈ R).
Ası́ la solución complementaria está dada por yc = c1 eax sen(bx) + c2 eax cos(bx).

Ejemplos

1. Resolver la ecuación y 00 + 2y 0 − 8y = 0.
Solución: La ecuación en notación de operadores corresponde a D2 + 2D − 8 y = 0, la ecuación

auxiliar es m2 + 2m − 8 = 0 ⇒ (m − 2)(m + 4) = 0, cuyas soluciones o raı́ces son m = 2 y m = −4.

La solución general de la ecuación corresponde a yc = c1 e2x + c2 e−4x .

2. Resolver la ecuación y 00 − 3y 0 + 2y = 0.
Solución: La ecuación en notación de operadores corresponde a D2 − 3D + 2 y = 0, la ecuación

auxiliar es m2 − 3m + 2 = 0 ⇒ (m − 1)(m − 2) = 0, cuyas soluciones o raı́ces son m = 1 y m = 2.

La solución general de la ecuación corresponde a yc = c1 ex + c2 e2x .

3. Resolver la ecuación 4y 00 − y = 0.
Solución: La ecuación en notación de operadores corresponde a 4D2 − 1 y = 0, la ecuación auxiliar


1 −1
es 4m2 − 1 = 0 ⇒ (2m − 1)(2m + 1) = 0, cuyas soluciones o raı́ces son m = y m = .
−x x
2 2
La solución general de la ecuación corresponde a yc = c1 e 2 + c2 e 2 .

4. Resolver la ecuación 9y 00 − 4y = 0.
Solución: La ecuación en notación de operadores corresponde a 9D2 − 4 y = 0, la ecuación auxiliar


2 −2
es 9m2 − 4 = 0 ⇒ (3m − 2)(3m + 2) = 0, cuyas soluciones o raı́ces son m = y m = .
−2 2
3 3
La solución general de la ecuación corresponde a yc = c1 e 3 x + c2 e 3 x .

5. Resolver la ecuación y 00 + 4y 0 + 4y = 0.
Solución: La ecuación en notación de operadores corresponde a D2 + 4D + 4 y = 0, la ecuación

auxiliar es m2 + 4m + 4 = 0 ⇒ (m + 2)(m + 2) = 0, cuya solución es m = 2 (raı́z doble).

La solución general de la ecuación corresponde a yc = c1 e2x + c2 xe2x .

38
6. Resolver la ecuación y 00 − 6y 0 + 9y = 0.
Solución: La ecuación en notación de operadores corresponde a D2 − 6D + 9 y = 0, la ecuación

auxiliar es m2 − 6m + 9 = 0 ⇒ (m − 3)(m − 3) = 0, cuyas soluciones o raı́ces son m = 3 (Multiplicidad

2). La solución general de la ecuación corresponde a yc = c1 e3x + c2 xe3x .
7. Resolver la ecuación y 00 + 4y 0 + 5y = 0.
Solución: La ecuación en notación de operadores corresponde a D2 + 4D + 5 y = 0, la ecuación


√
2 −4 ± −4 −4 ± 2i
auxiliar es m + 4m + 5 = 0, cuyas soluciones o raı́ces son m = = = −2 ± i.
2 2
−2x −2x
La solución general de la ecuación corresponde a yc = c1 e cos(x) + c2 e sen(x).
8. Resolver la ecuación y 00 − 2y 0 + 2y = 0.
Solución: La ecuación en notación de operadores corresponde a D2 − 2D + 2 y = 0, la ecuación


√
2 ± −4 2 ± 2i
auxiliar es m2 − 2m + 2 = 0, cuyas soluciones o raı́ces son m = = = 1 ± i.
x
2 x 2
La solución general de la ecuación corresponde a yc = c1 e cos(x) + c2 e sen(x).

Ecuación diferencial lineal homogénea de orden superior a 2

El método empleado para resolver una ecuación diferencial lineal homogénea de orden n, con n ≥ 3
de la forma φ(D)y = 0, es una generalización de la técnica desarrollada para las ecuaciones lineales de
orden dos.

De la misma manera, los ceros o raı́ces de su respectiva ecuación auxiliar φ(m) = 0, determinan la
solución general de φ(D)y = 0. Sabiendo que una solución principal o particular de ésta es una función
exponencial de la forma y = emx , donde m es un número real o complejo por determinar. Por tanto se
tiene las siguientes condiciones:
(a) La ecuación auxiliar φ(m) = 0, tiene n raı́ces reales distintas, sean m1 ,m2 , m3 ,..., mn . Ası́ las n
soluciones son linealmente independientes, por lo que la solución corresponde a:
y = c1 em1 x + c2 em2 x + c3 em3 x + ... + cn emn x .
Ejemplos:
(1) Resolver y 000 + y 00 − 6y 0 = 0.
La ecuación en notación de operadores corresponde a (D3 + D2 − 6D)y = 0.
La auxiliar φ(m) = m3 + m2 − 6m = 0 ⇒ m(m − 2)(m + 3) = 0, ası́ las soluciones son m = 0,
m = 2 y m = −3.
Por tanto la solución general corresponde a: y = c1 e0x + c2 e2x + c3 e−3x = c1 + c2 e2x + c3 e−3x

(2) Resolver (D3 + 5D2 + 2D − 12)y = 0.

3
La auxiliar φ(m) = m√ + 5m2 + 2m − 12 = 0 ⇒ (m + 3)(m2 + 2m − 4) = 0, ası́ las soluciones son
m = −3 y m = −1 ± 5. √ √
Por tanto la solución general es: y = c1 e−3x + c2 e(−1− 5)x + c3 e(−1+ 5)x

(3) Resolver (D3 + 2D2 − 5D − 6)y = 0.

La auxiliar φ(m) = m3 + 2m2 − 5m − 6 = 0 ⇒ (m − 2)(m + 3)(m + 1) = 0, ası́ las soluciones son
m = 2, m = −3 y m = −1.
Por tanto la solución general es: y = c1 e2x + c2 e−3x + c3 e−x

39
 (b) La ecuación auxiliar tiene a m1 como raı́z repetida k veces (multiplicidad k), por lo que (m − m1 )k
es un factor de dicha ecuación.

Por tanto la solución complementaria es y = c1 em1 x + c2 xem1 x + c3 x2 em1 x + ... + ck xk−1 em1 x .
Si la multiplicidad coincide con el orden de la ecuación diferencial, es decir k = n, entonces la
anterior es la solución buscada. En caso que k < n, entonces debe combinarse con las otras solu-
ciones de la ecuación auxiliar para obtener la solución general de φ(D)y = 0

Ejemplos:
(1) Resolver y 000 + 3y 00 − 4y = 0. En notación de operadores se tiene (D3 + 3D2 − 4)y = 0, la
auxiliar m3 + 3m2 − 4 = 0 ⇒ (m − 1)(m + 2)2 = 0, ası́ las soluciones son m = 1 y m = −2 (con
multiplicidad 2). Por tanto la solución general corresponde a y = c1 ex + c2 e−2x + c3 xe−2x

(2) Resolver 2y 000 + y 00 − 18y 0 − 9y = 0. En notación de operadores (2D3 + D2 − 18D − 9)y = 0, la

auxiliar 2m3 + m2 − 18m − 9 = 0 ⇒ (m − 3)(m + 3)(2m + 1) = 0, ası́ las soluciones son m = 3,
−1 −x
m = −3, m = . Por tanto, la solución general corresponde a y = c1 e3x + c2 e−3x + c3 e 2
2
(3) Resolver y 000 − 3y 00 + 3y 0 − y = 0. En notación de operadores se tiene (D3 − 3D2 + 3D − 1)y = 0,
la auxiliar m3 − 3m2 + 3m − 1 = 0 ⇒ (m − 1)3 = 0, ası́ las soluciones son m = 1 (con multiplicidad
3). Por tanto la solución general corresponde a y = c1 ex + c2 xex + c3 x2 ex

Práctica
1. Resuelva cada una de las siguientes ecuaciones diferenciales.

(a) y 00 + y 0 − 6y = 0
(b) y 00 + 6y 0 + 12y = 0
(c) 4y 00 + 12y 0 + 9y = 0
(d) y 00 − y 0 = 0
d2 y dy
(e) 2 2
+5 +y =0
dx dx
00 0
(f) y − 3y + 2y = 0
(g) 3y 00 − 8y 0 − 3y = 0
(h) 2y 00 − 5y 0 − 3y = 0
(i) y 00 − 10y 0 + 25y = 0
(j) y 00 + y 0 + y = 0
d2 y dy
(k) 2
+ 2 + 5y = 0
dx dx
2
dy dy
(l) 2
− 4 + 4y = 0
dx dx
2
dy dy
(m) 2
+2 −y =0
dx dx
00 0
(n) y − 3y + 2y = 0

40
 (ñ) y 00 + 16y = 0
(o) y 00 + 2y 0 − 8y = 0
(p) 2y 00 − 5y 0 = 0

2. Resuelva cada una de las siguientes ecuaciones diferenciales

(a) y 000 − 3y 00 + 7y 0 − 5y = 0
(b) y 000 + 3y 00 + 3y 0 + y = 0
(c) y 000 − 6y 00 + 11y 0 − 6y = 0
(d) y 000 − y 00 − y 0 + y = 0
(e) y 000 − 3y 00 + 3y 0 − y = 0
(f) y 000 − 3y 00 − 10y 0 = 0
(g) y 000 + 3y 00 − 4y 0 = 0
(h) y (4) + 2y 00 + y = 0
(i) 4y (4) − 16y 000 + 7y 00 + 4y 0 − 2y = 0
(j) 4y (4) − 4y 000 − 13y 00 + 7y 0 + 6y = 0
(k) 4y (4) − 7y 000 + 18y 00 − 20y 0 + 8y = 0
(l) 4y (4) − 4y 0 + 1 = 0
(m) 4y (4) − 4y 000 − 23y 00 + 12y 0 + 36y = 0
(n) y (4) − y = 0
(ñ) y (5) + y (3) = 0

3. Resuelva cada una de las siguientes ecuaciones diferenciales, sujetas a las condiciones iniciales dadas.

(a) y 00 − 4y 0 + 13y = 0, con y(0) = −1 y y 0 (0) = 2

(b) y 00 − 10y 0 + 25y = 0, con y(0) = 1 y y 0 (1) = 0
(c) y 00 − 2y 0 + y = 0, con y(2) = 0 y y 0 (2) = 1
π  π 
(d) y 00 + 4y = 0, con y =1yy =0
6 6
(d) 3y 000 + 5y 00 + y 0 + 7y = 0, con y(1) = 0, y 0 (1) = 0 y y 00 (1) = 0

Solución particular de la ecuación φ(D)y = F (x), F (x) 6= 0.

Se estudiará el procedimiento para resolver una ecuación diferencial lineal de coeficientes constantes, también
conocida como ecuación no homogénea, que como ya se habı́a indicado la solución es de la forma y = yc + yp .

Como primer paso se debe determinar la solución complementaria yc de la ecuación homogénea asociada
a la ecuación, es decir, la solución de φ(D)y = 0, posteriormente se determina la solución yp , lo cual se
puede hacer por diversos métodos, en particular estudiaremos el denominado Método de los coeficientes
indeterminados.

41
 Método de los coeficientes indeterminados.

Este método se aplica para obtener la solución particular yp de la ecuación lineal homogénea de orden n,
φ(D)y = F (x), F (x) 6= 0; siempre y cuando se cumplan simultáneamente las siguientes dos condiciones:

(a) La ecuación diferencial sea solo de coeficientes constantes.

(b) La función F (x) contenga: expresiones polinomiales (constantes, lineales, de grado dos, ...); o expre-
siones de la forma: eax , sen(bx), cos(bx) con a, b ∈ R; o combinaciones lineales de sumas y productos
finitos de estas expresiones.

Esto obedece a que las expresiones enunciadas tiene la propiedad que las derivadas de sus sumas y productos
son sumas y productos de constantes, polinomios, exponenciales, senos y cosenos.

Si como ejemplo se tiene la ecuación de segundo orden a0 y 00 + a1 y 0 + a2 y = F (x), si yp es una solución

particular de esta, entonces debe satisfacerla, es decir, a0 yp00 + a1 yp0 + a2 yp = F (x) es una identidad.

Con base en lo anterior, se puede suponer que una propuesta correcta para yp es aquella que tiene la misma
forma que F (x). Para comprender la forma de seleccionar yp consideremos dos casos.

Caso 1 . Ninguna de las funciones que conforman la propuesta de yp se repite con respecto a las funciones que
componen la solución complementaria yc .

Ejemplo 1. Resolver la ecuación y 00 + 3y 0 + 2y = 3x.

Solución:
Paso 1. Se resuelve la ecuación homogénea asociada y 00 + 3y 0 + 2y = 0. Lo cual equivale a resolver
m2 + 3m + 2 = 0 ⇒ (m + 2)(m + 1) = 0 ⇒ m = −2 y m = −1, ası́ yc = c1 e−2x + c2 e−x .
Paso 2. Se determina la forma de yp con base en la forma de F (x). Ası́, para el ejemplo como F (x) = 3x
una adecuada selección para yp es yp = Ax + B.
Paso 3. Partiendo del supuesto que yp es solución, se sustituye yp y sus derivadas en φ(D)y = F (x). Por
tanto yp = Ax+B; yp0 = A, yp00 = 0 sustituyendo en y 00 +3y 0 +2y = 3x se tiene que 0+3 (A)+2 (Ax+B) =
3x ⇒ 2Ax + (3A + 2B) = 3x.
Paso 4. Se resuelve el sistema de ecuaciones (o la ecuación) que se genera con la igual, para ası́
determinar los valores de los parámetros o coeficientes indeterminados. En el ejemplo se tiene que
3 3
2Ax + (3A + 2B) = 3x, lo cual se satisface si 2A = 3 ⇒ A = y 3A + 2B = 0 ⇒ 3 · + 2B = 0
2 2
−9 −9 3x 9
⇒ 2B = ⇒B= . Ası́ yp = − .
2 4 2 4
3x 9
Paso 5. Se escribe la solución general de la ecuación es decir y = yc + yp = c1 e−2x + c2 e−x + −
2 4

Ejemplo 2 Resolver la ecuación y 00 + 3y 0 + 2y = 10e3x .

Solución:
Paso 1. Se resuelve la ecuación homogénea y 00 + 3y 0 + 2y = 0. Cuya solución es: yc = c1 e−2x + c2 e−x .
Paso 2. Se determina la forma de yp con base en la forma de F (x), para ello puede usar coeficientes inde-
terminados A, B, C, D, entre otros. Ası́, para el ejemplo, una adecuada selección para yp es yp = Ae3x .
Paso 3. Partiendo del supuesto que yp = es solución, se sustituye yp y sus derivadas en φ(D)y = F (x).
Ası́ para el ejemplo, se tiene que yp = Ae3x ; yp0 = 3Ae3x , yp00 = 9Ae3x sustituyendo en y 00 + 3y 0 + 2y = 3x
se tiene que 9Ae3x + 3 · 3Ae3x + 2 Ae3x = 10e3x ⇒ 20Ae3x = 10e3x .

42
 Paso 4. Se resuelve el sistema de ecuaciones (o la ecuación) que se genera con la igual para ası́ deter-
minar los valores de los parámetros. En el ejemplo se tiene que 20Ae3x = 10e3x , lo cual se satisface si
1 e3x
20A = 10 ⇒ A = . Ası́ yp = .
2 2
e3x
Paso 5. Se escribe la solución general de la ecuación es decir y = yc + yp = c1 e−2x + c2 e−x +
2

Ejemplo 3 Resolver la ecuación y 00 + 3y 0 + 2y = sen(x).

Solución:
Paso 1. Se resuelve la ecuación homogénea y 00 + 3y 0 + 2y = 0. Cuya solución es: yc = c1 e−2x + c2 e−x .
Paso 2. Se determina la forma de yp con base en la forma de F (x). Ası́, para el ejemplo, una adecuada
selección para yp es yp = A sen(x) + B cos(x).
Paso 3. Partiendo del supuesto que yp es solución, se sustituye yp y sus derivadas en φ(D)y = F (x).
Ası́ para el ejemplo, se tiene que:
yp = A sen(x) + B cos(x); yp0 = A cos(x) − B sen(x), yp00 = −A sen(x) − B cos(x) sustituyendo en
y 00 + 3y 0 + 2y = sen(x) se logra:
[−A sen(x) − B cos(x)] + 3 [A cos(x) − B sen(x)] + 2 [A sen(x) + B cos(x)] = sen(x)
⇒ (A − 3B) sen(x) + (3A + B) cos(x) = sen(x).
Paso 4. Se resuelve el sistema de ecuaciones (o la ecuación) que se genera con la igual para ası́ deter-
minar los valores de los parámetros. En el ejemplo se tiene: 
A − 3B = 1
(A − 3B) sen(x) + (3A + B) cos(x) = sen(x), lo cual se satisface si
3A + B = 0
1 −3 1 3
Lo cual tiene por solución A = yB= . Ası́ yp = sen(x) − cos(x).
10 10 10 10
sen(x) 3 cos(x)
Paso 5. Se escribe la solución general de la ecuación y = yc + yp = c1 e−2x + c2 e−x + −
10 10

Como es notorio la dificultad del método radica en la elección correcta de yp . Algunas propuestas:

Función Propuesta
F (x) = 4x + 2 yp = Ax + B
F (x) = 7x2 + 2x − 3 yp = Ax2 + Bx + C
F (x) = 4e−2x yp = Ae−2x
F (x) = (2x − 5)e3x yp = (Ax + B)e3x
F (x) = (2x2 + 9x + 2)e7x yp = (Ax2 + Bx + C)e7x
F (x) = sen(20x) yp = A sen(20x) + B cos(20x)
F (x) = cos(−3x) yp = A sen(−3x) + B cos(−3x)
F (x) = e−x cos(5x) yp = e−x [A sen(5x) +
B cos(5x)]
yp = Ax2 + Bx + C sen(−3x) + Dx2 + Ex + F cos(−3x)


F (x) = 4x2 cos(−3x)
F (x) = xex cos(x) yp = ex [(Ax + B) sen(x) + (Cx + D) cos(x)]

Caso 2 . Una o varias de las funciones que conforman la propuesta de yp también forma parte de algún o
algunos de los componentes de yc .

La forma de resolver este tipo de ecuaciones es muy similar a la propuesta en el caso 1. La diferencia
radica en que la propuesta de yp debe considerar cuántas veces se presenta en la solución complemen-
taria, el o los componentes que también están en la función F (x) y que, por tanto, forman parte de yp .

43
 Esto conlleva a considerar que la propuesta de yp debe aplicar lo de las raı́ces repetidas de la ecuación
auxiliar asociada a la ecuación homogénea. Veamos algunos ejemplo.

Ejemplo 1. Resolver y 00 − 4y 0 + 3y = 27 − 32e3x .

Paso 1. Al resolver la ecuación homogénea y 00 − 4y 0 + 3y = 0, la solución es yc = c1 ex + c2 e3x .
Paso 2. Al plantear la solución particular, según la forma de F (x) = 27 − 32e3x parece evidente que
yp = A + Be3x , sin embargo, Be3x está contenida en la solución yc , por tanto aplicando la técnica de
las raı́ces repetidas la propuesta para yp es yp = A + Bxe3x .
Paso 3. Al derivar y sustituir se tiene que yp0 = Be3x + 3Bxe3x , yp00 = 3Be3x + 3Be3x + 9Bxe3x ⇒ yp00 =
6Be3x + 9Bxe3x . Ası́:
y 00 − 4y 0 + 3y = 27 − 32e3x ⇒ 6Be3x + 9Bxe3x − 4 Be3x + 3Bxe3x + 3 A + Bxe3x = 27 − 32e3x ⇒

3A + 2Be3x = 27 − 32e3x
Paso 4. Al resolver el sistema se tiene que 3A = 27 ⇒ A = 9 y 2B = −32 ⇒ B = −16. Por tanto
yp = 9 − 16xe3x
Paso 5. Se escribe la solución y = yc + yp = c1 ex + c2 e3x + 9 − 16xe3x .

Ejemplo 2. Resolver y 000 + 2y 00 = −8x + 3e−2x .

Paso 1. Al resolver la homogénea se tiene que y 000 + 2y 00 = 0. yc = c1 + c2 x + c3 e−2x .
Paso 2. Al plantear la solución particular, según la forma de F (x) = −8x + 3e−2x , parece evidente
que yp = A + Bx + Ce−2x , sin embargo A + Bx está contenida en la solución yc con c1 + c2 x, además,
Ce−2x está contenida también en yc con c3 e−2x , por tanto aplicando la técnica de las raı́ces repetidas
la propuesta para yp es yp = Ax2 + Bx3 + Cxe−2x .
Paso 3. Al derivar y sustituir se tiene que yp0 = 2Ax + 3Bx2 + Ce−2x − 2Cxe−2x , yp00 = 2A + 6Bx −
4Ce−2x + 4Cxe−2x y yp000 = 6B + 12Ce−2x − 8Cxe−2x . Ası́:

y 000 + 2y 00 = −8x + 3e−2x

⇒ 6B + 12Ce−2x − 8Cxe−2x + 2 2A + 6Bx − 4Ce−2x + 4Cxe−2x = −8x + 3e−2x

⇒ (4A + 6B) + 12Bx + 4Ce−2x = −8x + 3e−2x

−2 3
Paso 4. Al plantear y resolver el sistema se tiene que 12B = −8 ⇒ B = , 4C = 3 ⇒ C = y
3 4
4A + 6B = 0 ⇒ A = 1.
2 3
Paso 5. Se escribe la solución y = yc + yp = c1 + c2 x + c3 e−2x + x2 − x3 + xe−2x .
3 4
Ejercicios.

1. Resuelva cada una de las siguientes ecuaciones diferenciales

12x x 62ex
(a) y 00 + 4y 0 = 5x2 ex . Solución:y = c1 + c2 e−4x + x2 ex − e +
5 25
√ ! √ !!
00 0 2 x2 5x 11 3x 7x 7x
(b) y − 3y + 4y = x + x. Solución: + + +e2 c1 sen + c2 cos
4 8 2 2 2

(c) y 00 + 4y 0 + 4y = e3x
(d) y 00 + y 0 + y = x ex
(e) y 00 − 2y 0 + y = 2x

44
 (f) y 00 − y = ex sen(x)
(g) y 00 − 9y 0 = 81x2 + 5
(h) y 00 − y 0 − 2y = −2x3 + x2 + 4x + 1
(i) y 00 − y 0 − 2y = 16xe3x
(j) y 00 + 25y = 6 sen(x)
(k) y 00 − y 0 = (3x2 + 6x + 1) ex
(l) y 00 − y = x2 ex
(m) y 00 + 6y 0 + 8y = 3e5x + x2 − 1
(n) y 00 + 3y 0 + 2y = 3x2 − x + 1
(ñ) y 00 − 2y 0 − 8y = 16e2x − 21e−3x

2. Proponga una solución particular de y 00 + 3y 0 + 2y = x2 cos(πx).

Solución: yp = (Ax2 + Bx + C) sen(πx) + (Ax2 + Bx + C) cos(πx)

Otro método que permite resolver cierto tipo de ecuaciones diferenciales requiere del empleo de la transfor-
mada de Laplace.

Nociones básicas de la transformada de Laplace.

La transformada de Laplace, nos puede conducir a un método para resolver ecuaciones diferenciales lineales.
La ventaja del método es que se puede transformar una ecuación diferencial, en otra de tipo algebraico. De
esta manera, al resolver la ecuación algebraica, por métodos de matemática básica, paralelamente resolvemos
la ecuación diferencial.

Definiremos en primer instancia el concepto de transformada de Laplace.

Definición Sea F (t); t ≥ 0, la transformada de Laplace de F (t) denotada por L{F (t)}; está definida por la
Z ∞ Z b Z ∞
−st −st
integral impropia L{F (t)} = e F (t) dt = lim e F (t) dt, siempre y cuando la e−st F (t) dt sea
0 b→∞ 0 0
convergente. En cuyo caso conduce a una función que depende de la variable s, es decir L{F (t)} = F (s).

Ejemplos. Calcular
Z ∞ b
−e−st −e−sb e−s·0
 
1 1
1. L{1} = e−st · 1 dt = lim = lim + =0+ = , con s > 0.
0 b→∞ s 0
b→∞ s s s s
∞ b b
e(−s+a)t e(−s+a) b e(−s+a)·0
Z Z  
at −st at (−s+a)t
2. L{e } = e · e dt = lim e dt = lim = lim −
0 b→∞ 0 b→∞ −s + a b→∞ −s + a −s + a
0
1 1
=0− = , con s > a.
−s + a s−a
Z ∞ b b
−e−st ae−st cos(at) a2
Z
3. L{sen(at)} = e−st · sen(at) dt = lim sen(at) − − 2 e−st sen(at) dt
0 b→∞ s s2 0 s 0
De esta manera:

45
 a2
 
− sen(ba) a cos(ba) a a
L{sen(at)} = lim sb
− bs
+ 2 − 2
L{sen(at)} ⇒ L{sen(at)} = 2 . Con s > 0
b→∞ se se s s s + a2
sen(ba) sen(ba) ba sen(ba) ba a
Observación: lim sb
= lim sb
· = lim · lim sb = 0 · lim sb = 0. De igual
b→∞ e b→∞ e ba b→∞ ba b→∞ e b→∞ se
cos(ba) cos(ba) ba
manera se tiene que lim = lim · sb = 0
b→∞ esb b→∞ ba e
s
4. De forma similar se puede comprobar que L{cos(at)} = . Con s > 0
a2 + s 2
Además se tiene que

(a) L{cF (t)} = c L{F (t)}

(b) L{F (t) ± G(t)} = L{F (t)} ± L{G(t)}

Transformada de Laplace de algunas funciones básicas

1
(a) L{1} = , con s > 0
s
1
(b) L{t} = , con s > 0
s2
n!
(c) L{tn } = , con n = 1, 2, 3, 4, 5, ... y s > 0
sn+1
1
(d) L{eat } = , con s > a
s−a
1
(e) L{t eat } = , con s > a
(s − a)2
n!
(f) L{tn eat } = , con s > a
(s − a)n+1
a
(g) L{sen(at)} = , con s > 0
s2
+ a2
a
(h) L{ekt sen(at)} = , con s > a
(s − k)2 + a2
2as
(i) L{t sen(at)} = , con s > 0
(s2 + a2 )2
s
(j) L{cos(at)} = , con s > 0
s2 + a2
s−k
(k) L{ekt cos(at)} = , con s > k
(s − k)2 + a2
s 2 − a2
(l) L{t cos(at)} = , con s > 0
(s2 + a2 )2

46
 a
(m) L{senh(at)} = , con s > |a|
s2
− a2
s
(n) L{cosh(at)} = 2 , con s > |a|
s − a2

(ñ) L{eat F (t)} = f (s − a) con f (s) = L{F (t)}

La transformada de Laplace de derivadas

(a) L{y 0 (t)} = sY (s) − y(0); Y (s) = L{y(t)}

(b) L{y 00 (t)} = s2 Y (s) − sy(0) − y 0 (0)

(c) L{y (n) (t)} = sn Y (s) − sn−1 y(0) − ... − y (n−1) (0)

Ejemplos. Resuelva cada uno de los siguientes ejercicios.

1 e7
1. L{e−3t + et+7 } = +
s+3 s−1
2. L{6 sen(3t) cos(5t) + e2t cos(3t)} = L {3 [sen(8t) + sen(−2t)] + e2t cos(3t)}
24 6 s−2
= 2 − 2 +
s + 64 s + 4 (s − 2)2 + 9
4! 24
3. L{e−3t t4 } = 5
=
(s + 3) (s + 3)5
3! 10 1 6 10 1
4. L{t3 − 10t + 1} = 4
− 2 + = 4− 2 +
s s s s s s
11 5 12
5. L{11 + 5e4t − 6 sen(2t)} = + − 2
s s−4 s +4
7 1 3! 7 1 6
6. L{7 + e−3t + t2 } = + + 3 = + + 3
s s+3 s s s+3 s
2 120
7. L{t2 + t5 } = +
s3 s6
2s 5! 2s 120
8. L{t sen(t) + t5 } = + 6 = 2 + 6
(s2 + 1) 2 s (s + 1) 2 s
3 5 4 7s
9. L{3 − 5e2t + 4 sen(t) − 7 cos(3t)} = − + 2 − 2
s s−2 s +1 s +9
s−3
10. L{e3t · cos(5t)} =
(s − 3)2 + 25
2 s
11. L{senh(2t) + cos(3t)} = + 2
s2
−4 s +9
  
3 4 −t 2 3 4 −t 1 cos(2t) 12 24 2 2s
12. L{2t − t e + 4 cos (t)} = L 2t − t e + 4 + = 4 − 5
+ + 2
2 2 s (s + 1) s s +4

47
 Para algunos ejercicios es importante recordar que:
(a) 2 sen(A) cos(B) = sen(A + B) + sen(A − B) o que 2 cos(A) sen(B) = sen(A + B) − sen(A − B)
(b) 2 cos(A) cos(B) = cos(A + B) + cos(A − B) o que 2 sen(A) sen(B) = cos(A − B) − cos(A + B)
 
3 1 3 2 1 3 3 1 1 3
13. L sen(2t) − sen(3t) = · 2 − · 2 = · 2 − · 2
10 5 10 s + 4 5 s + 9 5 s +4 5 s +9
−3 2 1
14. L −3et + 2e2t + e3t =

+ +
s−1 s−2 s−3
4s 1
15. L t sen(2t) + e3t = 2

2
+
(s + 4) s−3
√ o √
n 2 1 7
16. L t2 e3t + te4t + sen 7t = 2
+ 2
+ 2
(s − 3) (s − 4) s +7
Práctica. Determine la transformada de Laplace para cada una de las funciones dadas.
1. L 2t4 − 4t + 10


2. L 1 + e4t − 3e−5t


3. L 4t2 − 5 sen(3t))


4. L t2 − e−9t + 5t


5. L 5e3t − 2e−t + 7


6. L t3 + 3t2 + 3t + 1


7. L {sen(2t) + cos(2t)}
8. L 1 + 2e2t − e−4t


9. L e3t + 2e−2t


10. L {2 sen(t) − 3 cos(3t)}

11. L 3e−t − sen(6t)


12. L t3 − 3t + cos(2t)


13. L {−3 cos(2t) + 5 sen(4t)}

14. L e−3t (t − 2)


La transformada inversa de Laplace

Si F (s) representa la transformada de Laplace de una función f (t), es decir L {f (t)} = F (s), se dice que f (t)
es la transformada inversa de Laplace de F (s) y se escribe L−1 {F (s)} = f (t).

Práctica. Usando la tabla de la transformada de Laplace, determine las siguientes transformadas inversas
básicas.

48
  
−1 1
1. L =1
s
 
−1 1
2. L =t
s2
 
−1 5 3
3. L + = 5 + 3t
s s2
t4
     
−1 1 −1 4! 1 −1 4!
4. L =L = L =
s5 4!s5 4! s5 24
  ( √ )
√ 
−1 1 1 −1 7 1
5. L = √ L = √ sen 7t
s2 + 7 7 s2 + 7 7
 
−1 n!
6. L = tn , n = 1, 2, 3, 4, 5, ...
sn+1
 
−1 2 4 5
7. L + − = 2e−t + 4e2t − 5t
s+1 s−2 s−1
 
−1 1
8. L = eat
s−a
 
−1 4 3!
9. L + = 4e3t + t3
s − 3 s4
 
−1 n!
10. L = tn eat
(s − a)n+1
 
−1 4!
11. L = t4 e2t
(s − 2)5
 
−1 k
12. L = senh(t)
s2 + k 2
t5 t5
   
−1 1 2 −1 5! 2 3t
13. L 6
+ = L 6
+ = + 2e = + 2e3t
s s−3 5! s s−3 5! 120
   
−12 −12 −1 3! −12 3 −5t
14. L−1 4
= L 4
= t e = −2t3 e−5t
(s + 5) 3! (s + 5) 3!
    t
√
1 1 e sen( 2 t)
15. L−1 2
= L−1 2
= √
s − 2s + 3 (s − 1) + 2 2
       
−1 s −1 s+2−2 −1 s+2 −1 1
16. L =L =L − 2L
s2 + 4s + 5 (s + 2)2 + 1 (s + 2)2 + 1 (s + 2)2 + 1
−2t −2t
= e cos(t) − 2e sen(t)

49
 s2 + 6s + 9
 
−1
17. L (usando fracciones parciales)
(s − 1)(s
 − 2)(s + 4)    
−16 −1 1 25 −1 1 1 −1 1
= L + L + L
5 s−1 6 s−2 30 s+4
t 2t −4t
−16e 25e e
= + +
5 6 30
 
−1 8 −2s + 6
18. L + 2 = 8e−3t − 2 cos(2t) + 3 sen(2t)
s+3 s +4
 
16s
19. L−1 = 8t sen(t)
(s2 + 1)2
Práctica. Determine la inversa de la Transformada de Laplace para una de las siguiente funciones.
 
−1 7
1. L
s3
 
−1 1 48
2. L − 5
s2 s
 
−1 1 1 1
3. L − +
s2 s s − 2
 
−1 1 5
4. L +
4s + 1 s + 7
       
−1 2 11 2 −1 1 −1 1 2 −1 4!
5. L + + = 2L + 11 L + L
s − 3 (s − 3)2 (s − 3)5 s−3 (s − 3)2 4! (s − 3)5
 
−1 7s − 1
6. L
(s − 3)(s + 2)(s − 1)
 
−1 1
7. L
s(s2 + 4)
 
−1 2s + 3
8. L
s2 − 4s + 20
 
−1 s 1
9. L −
s2 + 4 s2 + 4
 
−1 3s + 5
10. L
s2 + 7
 
−1 1 6 6
11. L + +
s7 s3 s2 + 9
Utilizando las propiedades de la transformada de Laplace se puede reducir una ecuación diferencial lineal con
coeficientes constantes, a una ecuación algebraica donde la incógnita de esta ecuación es la transformada de
Laplace de la solución de la ecuación diferencial, desde este punto de vista, para determinar la solución es
necesario calcular la transformada inversa de Laplace del resultado de la ecuación algebraica.

50
Ejemplos. Resolver.

(a) y 0 + 3y = 13 sen(2t); con y(0) = 6

Solución:
y 0 + 3y = 13 sen(2t) ⇒ L {y 0 + 3y} = L {13 sen(2t)}
26
⇒ sy(s) − y(0) + 3y(s) = 2
s +4
26
⇒ sy(s) − 6 + 3y(s) = 2
s +4
26
⇒ (s + 3) y(s) − 6 = 2
s +4
26
⇒ (s + 3) y(s) = 2 +6
s +4
6s2 + 50
⇒ (s + 3) y(s) = 2
s +4
6s2 + 50 8 −2s + 6 8 s 2
⇒ y(s) = 2
= + 2 = −2 2 +3 2
(s + 3) (s + 4) s+3 s+ 4 s + 3  s + 4  s +  4  
−1 −1 8 −2s + 6 −1 8 −1 s −1 2
⇒ L {y(s)} = L + 2 =L − 2L + 3L
s+3 s +4 s+3 s2 + 4 s2 + 4
−3t
⇒ y(t) = 8e − 2 cos(2t) + 3 sen(2t).

(b) y 00 + y = 16 cos(t); con y(0) = (0) y y 0 (0) = (0)

Solución:
y 00 + y = 16 cos(t) ⇒ L {y 00 + y} = L {16 cos(t)}
16
⇒ s2 y(s) − sy(0) − y 0 (0) + y(s) = 2
s +1
2

16s
⇒ s + 1 y(s) = 2
s +1
8 · 2s
⇒ y(s) = 2
(s + 1)2  
−1 −1 8 · 2s
⇒ L {y(s)} = L
(s2 + 1)2
⇒ y(t) = 8t sen(t).

(c) y 00 + 2y 0 + 5y = 0; con y(0) = 3; y 0 (0) = −7

Solución:
L {y 00 + 2y 0 + 5y} = L {0}
⇒ s2 y(s) − sy(0) − y 0 (0) + 2[sy(s) − y(0)] + 5y(s) = 0
⇒ s2 y(s) − 3s + 7 + 2sy(s) − 6 + 5y(s) = 0
⇒ (s2 + 2s + 5) y(s) = 3s − 1
3s − 1 3s − 1 3(s + 1) − 4 3(s + 1) 2·2
⇒ y(s) = 2 = = = −
s + 2s + 5 (s + 1)2 + 4 (s + 1)2 + 4 (s + 1)2 + 4 (s + 1)2 + 4
 
−1 −1 3(s + 1) 2·2
Ası́ L {y(s)} = L − ⇒ y(t) = 3e−t cos(2t) − 2e−t sen(2t)
(s + 1)2 + 4 (s + 1)2 + 4
(d) y 00 − 4y 0 + 3y = 0; con y(0) = 3 y y 0 (0) = 5
Solución:
L {y 00 − 4y 0 + 3y} = L {0}
⇒ s2 y(s) − sy(0) − y 0 (0) − 4[sy(s) − y(0)] + 3y(s) = 0

51
 ⇒ s2 y(s) − 3s − 5 − 4[sy(s) − 3] + 3y(s) = 0
⇒ s2 y(s) − 3s − 5 − 4sy(s) + 12 + 3y(s) = 0
⇒ y(s) (s2 − 4s + 3) = 3s − 7
3s − 7
⇒ y(s) = 2
s − 4s + 3
3s − 7 1 2
⇒ y(s) = = +
(s − 3) (s −1) s − 3 s− 1
1 2
⇒ L−1 {y(s)} = L−1 +
s−3 s−1
⇒ y(t) = e3t + 2et

(e) y 00 + 2y 0 = 4; con y(0) = 1 y y 0 (0) = 4

Solución:
4
L{y 00 + 2y 0 } = L{4} ⇒ s2 y(s) − sy(0) − y 0 (0) + 2[sy(s) − y(0)] =
s
2 4
⇒ s y(s) − s + 4 + 2[sy(s) − y(0)] =
s
2 4
⇒ s y(s) − s + 4 + 2sy(s) − 2 =
s
2 4
⇒ y(s)(s + 2s) + 2 − s =
s
2 4
⇒ y(s)(s + 2s) = + s − 2
s
2 s2 − 2s + 4
⇒ y(s)(s + 2s) =
s
s2 − 2s + 4 s2 − 2s + 4 3 2 2
⇒ y(s) = = = − +
s(s2 + 2s) s2 (s + 2)  s + 2 s s2
3 2 2
⇒ L−1 {y(s)} = L−1 − + 2
s+2 s s
⇒ y(t) = 3e−2t + 2t − 2

(f) y 00 − 3y 0 + 2y = 12e4t ; con y(0) = 1 y y 0 (0) = 0

Solución:
12
L{y 00 − 3y 0 + 2y} = L{12e4t } ⇒ s2 y(s) − sy(0) − y 0 (0) − 3[sy(s) − y(0)] + 2y(s) =
s−4
12
⇒ s2 y(s) − s − 3[sy(s) − 1] + 2y(s) =
s−4
2 12
⇒ s y(s) − s − 3sy(s) + 3 + 2y(s) =
s−4
12
⇒ s2 y(s) − 3sy(s) + 2y(s) = +s−3
s−4
12 + s2 − 4s − 3s + 12
⇒ y(s)(s2 − 3s + 2) =
s−4
s2 − 7s + 24 2 7 6
⇒ y(s) = = − +
(s − 4)(s −2)(s − 1) s−4 s− 2 s−1
2 7 6
⇒ L−1 {y(s)} = L−1 − +
s−4 s−2 s−1
⇒ y(t) = 2e4t − 7e2t + 6et

52
 (g) y 00 − 3y 0 + 2y = e−4t ; con y(0) = 1 y y 0 (0) = 5
Solución:
1
L{y 00 − 3y 0 + 2y} = L{e−4t } ⇒ s2 y(s) − sy(0) − y 0 (0) − 3[sy(s) − y(0)] + 2y(s) =
s+4
2 1
⇒ s y(s) − s − 5 − 3[sy(s) − 1] + 2y(s) =
s+4
2 1
⇒ s y(s) − s − 5 − 3sy(s) + 3 + 2y(s) =
s+4
2 1
⇒ (s − 3s + 2)y(s) = s + 2 +
s+4
2 s2 + 6s + 9
⇒ (s − 3s + 2)y(s) =
s+4
s2 + 6s + 9 s2 − 3s + 2 −16 1 25 1 1 1
⇒ y(s) = 2 = = + +
(s − 3s +2)(s + 4) (s − 1)(s − 2)(s + 4)  5 s − 1 6 s − 2 30 s + 4
−16 1 25 1 1 1
⇒ L−1 {y(s)} = L−1 + +
5 s−1 6 s − 2 30 s + 4
−16 t 25 2t 1
⇒ y(t) = e + e + e−4t
5 6 30
Práctica Usando la transformada de Laplace resuelva las siguientes ecuaciones diferenciales.

1. y 0 − y = 1; con y(0) = 0

2. 2y 0 + y = 0; con y(0) = −3

3. y 0 + 6y = e4t ; con y(0) = 2

4. y 00 + 5y 0 + 4y = 0; con y(0) = 1 y y 0 (0) = 0

5. y 00 + 9y = et ; con y(0) = 0 y y 0 (0) = 0

1
6. 4y 00 − y = 1; con y(0) = 0 y y 0 (0) =
2
7. y 00 − 3y 0 + 2y = 0; con y(0) = 0 y y 0 (0) = 4

Material de referencia.

(a) Acuña, L y Calderón C. (2011). Sucesiones y Series. Editorial Tecnológica de Costa Rica.

(b) Ávila, J. (2003). Ejercicios de Cálculo. Editorial Tecnológica de Costa Rica.

(c) Barrantes, H. (2008). Introducción a las Ecuaciones Diferenciales. EUNED.

(d) Becerrel, J y Elizarraraz, D.(2004). Ecuaciones diferenciales. Técnicas de solución y aplicaciones. Uni-
versidad Autónoma Metropolitana. México.

(e) González, J y De León B. (1994). Conceptos básicos de ecuaciones diferenciales. Material complemen-
tario curso Ecuaciones Diferenciales. Universidad Nacional. Costa Rica.

(f) Larson, R; Hostetler, R y Edwards, B. (1999). Cálculo. McGrawHill.

(g) Meneses, S. (2005). Ecuaciones diferenciales. Editorial Instituto Tecnológico de Costa Rica.

53
 (h) Meneses, S. (2007). Aplicaciones de las ecuaciones diferenciales de primer orden a la.... Editorial Ins-
tituto Tecnológico de Costa Rica.

(i) Murillo, M. (2003). Sucesiones y Series. Editorial Tecnológica de Costa Rica.

(j) Nagle, Saff y Snider. (2005). Ecuaciones diferenciales y problemas con valores en la frontera. PEAR SON.

(k) Stewart, J. (1994). Cálculo. Editorial Iberoamericana.

(l) Zill, D.(2006).Ecuaciones diferenciales con problemas de valores en la frontera.Editorial THOMSON.

Direcciones consultadas

(a) Sucesiones y series.

Disponible en: http://www.fic.umich.mx/ lcastro/sucesiones.pdf

(b) Sucesiones de números reales.

Disponible en: http://www.uoc.edu/in3/emath/docs/Sucesiones.pdf

(c) Introducción a las sucesiones y series numéricas.

Disponible en: http://www.matematica.ciens.ucv.ve/labfg/mat3/ssmat3.pdf

(d) http://analisisycalculo.blogspot.com/2012/01/ecuaciones-diferenciales-exactas.html

(e) http://pcmap.unizar.es/∼manuel calvo/edo1.pdf

(f) http://analisisycalculo.blogspot.com/2012/01/ecuaciones-diferenciales-lineales-de.html

54