Universidad Nacional Autónoma de México

Ecuaciones Diferenciales
Tarea 1
Cruz Reyes Mirya Lizbeth
Merced González Yadira Josseline
Zárate Casas Axel Vladimir
Ciudad de México a 2 de marzo de 2024

Ejercicios
1.Encuentre la solución general a las siguientes ecuaciones diferenciales.

a) dy
dt + ycos(t) = 0
Resolveremos por el metodo de variables separables: Separando las variables tenemos:

dy
dy = (−ycos(t))dt −→ = cos(t)dt
−y
Integrando ambos lados de la ecuacion tenemos:
Z Z Z Z
dy 1
= cos(t)dt −→ − dy = cos(t)dt −→ −ln(y) = sin(t) −→ ln(y) = −sin(t)
−y y

−→ y = exp(−sin(t) + c)
dy 2x
b) dx+ 1+x 2 = 0
Resolveremos por el método de variables separables: Separando las variables tenemos:
−2x
dy = dx
1 + x2
Integrando ambas partes de la ecuación tenemos:
−2x
Z Z Z
x
dy = dx −→ y = −2 dx...(⋄)
1 + x2 1 + x2

Resolviendo la integral del lado derecho: declaramos t = x2 + 1 po lo tanto nuestra integral queda:
Z Z Z
1 1 1 1
−2 dt −→ −2 dt = −1 dt = −ln(t)
2t 2 t t
Sustituyendo el valor de t tenemos:
Z
x
−2 dx = −ln(x2 + 1)
1 + x2

sustituyendo este valor en (⋄) obtenemos:

y = −ln(x2 + 1) + c

1
c) dy 2
dt + t y = t
2

Resolveremos por el método de variables separables: Separando las variables tenemos:

dy dy dy
= t2 − t2 y −→ = t2 (1 − y) −→ = t2 dt
dt dt 1−y
Integrando ambos lados de la ecuación tenemos:
Z Z
dy
= t2 dt...(⊛)
1−y
Resolviendo la integral del lado izquierdo: declaramos t = 1 − y por lo tanto nuestra integral queda:
Z Z
1 1
− dt −→ − dt = −ln(t)
t t
Sustituyendo el valor de t: Z
dy
= −ln(1 − y)
1−y
sustituyendo en ⊛ tenemos:

t3 t3 t3
Z
−ln(1 − y) = t2 dt −→ −ln(1 − y) = −→ 1 − y = exp(− ) −→ −y = −exp(− )
3 3 3

t3
−→ y = exp(− + c)
3
d) dy
dt + y = 2t
dy
Notemos que se trata de una ecuación diferencial lineal de la forma: dx + p(x)y = q(x), resolveremos por
factor integrante: iniciemos declarando lo siguiente.

p(x) = 1; q(x) = 2t

Veamos ahora que: Z Z

µ = exp( p(t)dt) −→ µ = exp( 1dt) = exp(t)

Ahora aplicando: Z
yµ = qµdx
Z Z Z
−→ yexp(t) = 2texp(t)dt = 2 texp(t)dt = 2(texp(t) − exp(t)dt) = 2(texp(t) − exp(t))

2texp(t) − 2exp(t) exp(t)(2t − 2)

−→ yexp(t) = 2texp(t) − 2exp(t) −→ y = = = 2t − 2 + c
exp(t) exp(t)
−→ y = 2t − 2 + c
f)y dy
dt+t=0
Resolveremos por el método de variables separables: Separando las variables tenemos:
dy
y = −t −→ ydy = −tdt
dt
Integrando ambos lados de la ecuación tenemos:

y2 t2
Z Z
ydy = −tdt −→ = − −→ y 2 = −t2
2 2

2
 p
−→ y = ± −t2 + c
dy dy
g) dx + y = 2exp(x) Notemos que se trata de una ecuación diferencial lineal de la forma: dx + p(x)y = q(x),
resolveremos por factor integrante: iniciemos declarando lo siguiente.

p(x) = 1; q(x) = 2exp(x)

Veamos ahora que: Z Z

µ = exp( p(t)dt) −→ µ = exp( 1dt) = exp(x)

Ahora aplicando: Z
yµ = qµdx
Z Z
yexp(x) = 2exp(x)exp(x)dx = 2exp(x)dx

declaramos t = 2x obteniendo: Z
exp(t)dt = exp(t)

remplazando a t: Z
2exp(x)dx = exp(2x)

exp(2x)
−→ yexp(x) = exp(2x) −→ y = = exp(x) + c
exp(x)
−→ y = exp(x) + c
2.Encuentre la solución a las siguientes ecuaciones diferenciales dado el valor inicial.

a) dy
dt − 2yt = 1; y(0) = 1
dy
Notemos que se trata de una ecuación diferencial lineal de la forma: dx + p(x)y = q(x), resolveremos por
factor integrante: iniciemos declarando lo siguiente.

p(x) = −2t; q(x) = 1

Veamos ahora que: Z Z

µ = exp( p(t)dt) −→ µ = exp( −2tdt) = exp(−t2 )

Ahora aplicando: Z
yµ = qµdx

1√
Z Z
2 2
yexp(−t ) = (1)(exp(−t ))dt = exp(−t2 )dt = πerf t(t) + C
2
Por lo tanto nuestra familia de soluciones es:
1√
y = (exp(t2 ))( πerf (t) + C)
2
Con erf (t) como la función error
Z x
2
erf (x) := √ exp(−t2 ) dt
π 0

3
Ahora resolvamos para la solución particular: sabemos que

y(0) = 1 −→ t = 0 : y = 1

Sustituyendo en la familia de soluciones:

1√ 1√
1 = (exp(02 ))( πerf (0) + C) −→ 1 = (1)( π(0) + C)
2 2
−→ C = 1
Sustituyendo C tenemos la solución particular:
1√
y = (exp(t2 ))( πerf (t) + 1)
2
b) dy
dt − 2ty = t; y(0) = 1
Resolveremos por el método de variables separables: Separando las variables tenemos:
dy
= t + 2ty = t(1 + 2y)
dt
dy 1 1
−→ = t −→ dy = tdt
dt 1 + 2y 1 + 2y
Integrando ambos lados tenemos: Z Z
1
dy = tdt...(⊡)
1 + 2y
resolviendo la integral del lado izquierdo: declaramos t = 1 + 2y obteniendo ası́:
Z Z
1 1 1 1
dt = dt = ln(t)
2t 2 t 2
Sustituyendo el valor de t obtenemos:
Z
1 1
dy = ln(1 + 2y)
1 + 2y 2
Sustituyendo en ⊡:
t2
Z
1
ln(1 + 2y) = tdt =
2 2
1 t2
−→ ln(1 + 2y) = −→ ln(1 + 2y) = t2 −→ 1 + 2y = exp(t2 + c)
2 2
Por lo tanto nuestra familia de soluciones es:
exp(t2 + c) 1
1 + 2y = exp(t2 + c) −→ y = − ...(†)
2 2
Ahora resolvamos para la solución particular: sabemos que

y(0) = 1 −→ t = 0; y = 1

Sustituyendo en la familia de soluciones:

exp(02 + c) 1 exp(c) 1
1= − = − −→ 2 = exp(c) − 1
2 2 2 2
−→ 2 − exp(c) = 1 −→ −exp(c) = −1 − 2 = −3 −→ exp(c) = 3

4
 −→ C = ln(3)
Sustituyendo C en † tenmos:

exp(t2 + ln(3)) 1 exp(t2 )exp(ln(3)) 1

y= − = −
2 2 2 2
Por lo tanto nuestra solución particular es:

3exp(t3 ) 1
y= −
2 2
dy 1
c) dx = exp(y)+x ; y(1) = 0
Primero veamos que:
1 dy dy
− + = 0 −→ (exp(y) + x) − 1 = 0...(⋆)
exp(y) + x dx dx

Ahora: sea M = −1 y N = exp(y) + x

∂M ∂N
−→ = 0; =1
∂y ∂x
Como las parciales no son iguales entonces se trata de una ecuación inexacta. Por lo tanto debemos de
hallar un factor de integración u que permita convertir a la E.D en Exacta.
∂ ∂ dU (y)
−→ (U (y)M ) = (U (y)N ) −→ − = U (y)
∂y ∂x dy
∂U (y) Z ∂U (y) Z
∂y ∂y
−→ = −1 −→ = −1 −→ ln(u(y)) = −y −→ U (y) = exp(−y)
U (y) U (y)
−→ U (y) = exp(−y)
Ahora multiplicando U (y) en (⋆) tenemos:
dy dy
(exp(−y)x + 1) − exp(−y) = 0 −→ −exp(−y) + (exp(−y)x + 1) =0
dx dx
Donde: M = −exp(−y) y N = exp(−y)x + 1 Veamos lo siguiente
∂M ∂N
= exp(−y); = exp(−y)
∂y ∂x
Como las parciales son iguales, entonces se trata de una E.D Exacta, procedemos a resolver, integraremos
para obtener f (x, y)
Z
x
f (x, y) = −exp(−y)dx = −exp(−y)x = − + g(y)
exp(y)

Ahora queremos saber cuanto vale g(y)

∂f ∂ x dg
−→ = (− + g(y)) = exp(−y)x +
∂y ∂y exp(y) dy
Igualando el anterior resultado tenemos:
dg dg
exp(−y)x + = exp(−y)x + 1 −→ =1
dy dy

5
Como queremos saber el valor de g(y) entonces integramos lo anterior
Z
1dy = y

Ahora sustituyendo g(y) en f (x, y)

f (x, y) = −exp(−y)x + y
De esta forma expresaremos la solución general:

−→ −exp(−y)x + y = C

−→ y = W (exp(−C)X) + C
Ahora aplicaremos los valores iniciales:

y(1) = 0 −→ x = 1; y = 0

Sustituyendo en la familia de soluciones:

−exp(−0)1 + 1 = C −→ −1 + 1 = C −→ C = 0

Sustituyendo C tenemos la solución particular:

y = W (exp(−0)X) + 0 = W (1) = W
dy
d) dx − y = 2xexp(2x); y(0) = 1
dy
Notemos que se trata de una ecuación diferencial lineal de la forma: dx + p(x)y = q(x), resolveremos por
factor integrante: iniciemos declarando lo siguiente.

p(x) = −1; q(x) = 2xexp(2x)

Veamos ahora que: Z Z

µ = exp( p(t)dx) −→ µ = exp( −1dx) = exp(−x)

Ahora aplicando: Z
yµ = qµdx
Z Z
−→ yexp(−x) = exexp(2x)exp(−x)dx = 2 xexp(2x)exp(−x)dx
Z Z
= 2 xexp(x)dx = 2(xexp(x) − exp(x)dx)

= 2(xexp(x) − exp(x)) = 2xexp(x) − 2exp(x) + c

−→ yexp(−x) = 2xexp(x) − 2exp(x) + c −→ y = 2xexp(2x) − 2exp(2x) + Cexp(x)
Por lo que nuestra familia de soluciones es:

y = 2exp(2x)(x − 1) + Cexp(x)...(⊙)

Ahora resolvamos para la solución particular: sabemos que

y(0) = 1 −→ x = 0; y = 1

Sustituyendo en la familia de soluciones:

1 = 2exp(2(0))(0 − 1) + Cexp(0) = 2exp(0)(−1) + C = 2(1)(−1) + C = −2 + C

6
 −→ C = 3
Sustituyendo C en (⊙)

y = 2exp(2x)(x − 1) + 3exp(x) = 2xexp(2x) − 2exp(2x) + 3exp(x)

Por lo tanto nuestra solución particular es:

y = 2xexp(2x) − 2exp(2x) + 3exp(x)

dy
f) dx − 2xy = 1; y(0) = 21
dy
Notemos que se trata de una ecuación diferencial lineal de la forma: dx + p(x)y = q(x), resolveremos por
factor integrante: iniciemos declarando lo siguiente.

p(x) = −2x; q(x) = 1

Veamos ahora que: Z Z

µ = exp( p(x)dx) −→ µ = exp( −2xdx) = exp(−x2 )

Ahora aplicando: Z
yµ = qµdx

1√
Z
yexp(−x2 ) = (1)exp(−x2 )dx = πerf (x) + c
2
Por lo tanto nuestra familia de soluciones es:
1√
y = exp(x2 )( πerf (x) + C)
2
Ahora resolvamos para la solución particular: sabemos que
1 1
y(0) = −→ x = 0; y =
2 2
Sustituyendo en la familia de soluciones:
1 1√ 1 1√
= exp(02 )( πerf (0) + C) −→ = πerf (0) + C
2 2 2 2
1 1√ 1
−→ = π(0) + C −→ = C
2 2 2
√
Sustituyendo C en la familia de soluciones tenemos: y = exp(x2 )( 12 πerf (x) + 12 )
dy
g) (1 − x2 ) dx − xy = x(1 − x2 ); y(0) = 2
Veamos que:
dy dy xy x − x3
(1 − x2 ) − xy = x − x3 −→ − =
dx dx 1 − x2 1 − x2
dy
Notemos que se trata de una ecuación diferencial lineal de la forma: dx + p(x)y = q(x), resolveremos por
factor integrante: iniciemos declarando lo siguiente.

x x − x3
p(x) = − ; q(x) =
1 − x2 1 − x2
Veamos ahora que: Z Z
x
µ = exp( p(t)dx) −→ µ = exp( − dx)
1 − x2

7
Resolviendo la integral: Z
x
− dx
1 − x2
aplicando t = 1 − x2 tenemos:
Z Z Z
1 1 1 1 1 1
− − dt = −1(− ) dt = dt = ln(t)
2t 2 d 2 t 2
remplazando el valor de t Z
x 1
− 2
dx = ln(1 − x2 )
1−x 2
De esta forma veamos que:
1 1 1
p
µ = exp( ln((1 − x2 )) = exp(ln(1 − x2 ) 2 ) = (1 − x2 ) 2 = 1 − x2
2
Ahora aplicando: Z
yµ = qµdx

x − x3 p x(1 − x2 ) p
p Z Z Z p
2
−→ y( 1 − x ) = 2
( 1 − x )dx = 1 − x dx = x 1 − x2 dx
2
1 − x2 1 − x2
Tomado t = 1 − x2 √
1√ 1 √
Z Z Z
1 1 1 2t t
− tdt = − tdt = t 2 dt = −
2 2 2 2 3
√ √
1 2(1 − x2 ) 1 − x2 (1 − x2 ) 1 − x2
=− =− +C
2 3 3
√
p (1 − x2 ) 1 − x2 (−1 + x2 ) C
−→ y( 1 − x2 ) = − + C −→ y = +√
3 3 1 − x2
Por lo tanto nuestra familia de soluciones es:
(−1 + x2 ) C
y= +√
3 1 − x2
Ahora resolvamos para la solución particular: sabemos que

y(0) = 2 −→ x = 0; y = 2

Sustituyendo en la familia de soluciones:

(−1 + 02 ) C −1 1 1 7
2= +√ = + C = − + C −→ −C = − − 2 = −
3 1 − 02 3 3 3 3
7
−→ C =
3
Sustituyendo C en la familia de soluciones tenemos:
7
(−1 + x2 )
y= +√ 3
3 1 − x2
Por lo tanto nuestra solución particular es:
7
(−1 + x2 )
y= +√ 3
3 1 − x2

8
3. Una partı́cula se mueve en lı́nea recta con x, v y a representando la posición, velocidad, y la aceleración,
respectivamente. Considerando que la velocidad v(t) es una función compuesta por la posición, es decir,
v(t) = v(x(t)). Use la regla de la cadena para encontrar que se cumple la relación v 2 = vo2 + 2a (x − x0 )
cuando dv dt = a y se considera que la aceleración es constante y v0 y x0 son la velocidad y la posición
iniciales.

Si tenemos una función compuesta y = f (g(x)), donde g(x) es la función interna y f (u) es la función
externa, la derivada de y con respecto a x se expresa como:

dy dy du
= ·
dx du dx
En nuestro caso, consideramos v(t) = v(x(t)), donde v(t) es la velocidad y x(t) es la posición. Queremos
encontrar dv dv
dt cuando dt = a.
Aplicando la regla de la cadena, tenemos:

dv dv dx
= ·
dt dx dt
dv
Dado que dt = a, podemos expresar la derivada de la velocidad con respecto a la posición como:

dv
a= ·v
dx
dv
Resolviendo para dx :

dv a
=
dx v
Ahora, resolviendo por el método de variables separables: iniciamos separando las variables obteniendo:

dv
−→ v= a −→ vdv = adx
dx
Integrando ambos lados de la igualdad tenemos:

v2
Z Z
vdv = adx −→ + c1 = ax + c2
2
Por lo tanto nuestra familia de soluciones es:
v2
= ax + C
2
Ahora apliquemos nuestros valores iniciales V(0) = v0 (con v0 nuestra velocidad inicial) sustituyendo en
nuestra familia de soluciones tenemos:
(v0 )2 (v0 )2 (v0 )2
= a(0) + C −→ = C −→ C =
2 2 2
Sustituyendo en nuestra familia de soluciones para ası́ obtener nuestra solución particular obtenemos:

v2 (v0 )2
= ax + −→ v 2 = 2ax + (v0 )2 −→ v 2 = (v0 )2 + 2ax
2 2
Tomando x = x − x0 ya que se trata de nuestra posición obtenemos:

v 2 = (v0 )2 + 2a(x − x0 )

9
4. Use el factor integrante para encontrar la solución general de la ecuación lineal

y ′ + y cos(x) = cos(x)

dy
+ P (x)y = Q(x)
dx
Factor Integrante R R
P0 dx cos(x)dx
µ(x) = e =e = esin(x)

Entonces tenemos:
 
sin(x) dy
e + y cos(x) = cos(x)
dx
dy
esin(x) + esin(x) y cos(x) = esin(x) cos(x)
dx
Z Z
sin(x) ′
(e y) = esin(x) cos(x)dx

Se realiza un cambio de variables:

u : sin x =⇒ du : cos x
Z
eu du = eu + c = esin(x) + c

y = esin(x) + c
−sin(x)
y = [esin(x) + c]e
−sin(x)
y = 1 + ce
5. La resistencia del aire provoca que algunos objetos se desaceleren cuando caen y en algunos casos al-
cancen una velocidad terminal. Ası́, considerando que la resistencia del aire es proporcional a la velocidad,
la ecuación diferencial que rige el movimiento es
dv r
= −g − v,
dt m
donde v es la velocidad, g es la aceleración de la gravedad y es constante, r es una constante de resistencia
y m es la masa del objeto. Use la solución homógenea de la ecuación diferencial para encontrar la solución
general para v(t) y determine el valor de la velocidad terminal en el lı́mite t → ∞
Tenemos

dv r
= −g − v
dt m
Donde:
- v es la velocidad del objeto.
- g es la aceleración debida a la gravedad (constante).
- r es la constante de resistencia del aire.

10
- m es la masa del objeto.
rv(t)
dividimos ambos lados por −g − m

dv(t)
dt
=1
−g − rv(t)
m
Integramos ambos lados respecto a t:
Z dv(t) Z
dt
dt = 1dt
−g − rv(t)
m
Simplificamos y reescribimos:
Z
v
− mg+rvt dt
m
simplificamos:
Z
vm
− dt
mg + rvt
Z
1
−vm dt
mg + rvt
usamos sustitución u = mg + rvt
Z
1
−vm du
rvu
Z
1 1
−vm du
rv u
simplificamos:
Z
m 1 m
− du = − log(|u|)
r u r
Recordamos que u = mg + rvt
m
− log(|mg + rv(t)|)
r
simplificamos

m log(|gm + rv(t)|)
− +c
r
Resolvemos para v(t):
gm −(rt)
v(t) = ce m −
r
Cuando (t) tiende a infinito, el exponente en el primer término se vuelve muy negativo. Esto significa que
−(rt)
(e m ) se acerca a cero. Por lo tanto, el primer término se vuelve (c · 0 = 0).

El segundo término (− gmr ) no depende de (t), por lo que permanece constante.

El lı́mite de (v(t)) cuando (t → ∞) es simplemente el segundo término:

11
 gm
lı́m v(t) = −
t→∞ r
6. Demostrar que la ecuación

dy x2
=
dx 1 − y2
es separable, y encontrar la ecuación de sus curvas solución. Dibuje el diagrama de las curvas solución

♢ Se debe de encontrar (x0 , y 0 ) = (3, 2)

x2 32
f (x, y) = 2
f (3, 2) = = −3
1−y 1 − (2)2

♢ Se calculara f (x, y)f (x0 , y 0 ) = f (3, y)f (x, 2)

32 x2 −3x2
( )( ) =
1 − (2)2 1 − 22 1 − y2

♢ Comprobación si f (x0 , y 0 )(x, y) = f (x0 , y)f (x, y 0 )

x2 −3x2
f (3, 2)( ) =
1 − y2 1 − y2
Entonces es separable

dy x2
= =⇒ (1 − y 2 )dy = x2 dx
dx 1 − y2
Z Z Z Z Z
1 − y = x dx =⇒ 1dy − y dx = x2 dx
2 2 2

y3 x3
=⇒ y − = +C
3 3

12
Figura 1: la curva esta dada por la constante C que tiene un valor de -6.1

Figura 2: la curva esta dada por la constante C que tiene un valor de 6.1

Figura 3: la curva esta dada por la constante C que tiene un valor de 0

13