Facultad de Ciencias Fı́sicas y Matemáticas Universidad de Chile

MA2601 - Ecuaciones Diferenciales Ordinarias. Semestre 2009-03

Profesor: Julio López
Auxiliar: Sebastián Reyes Riffo.

Pauta Control 1, MA2601 EDO (07/01/10)

P1. (i) (3 pts.) Muestre que la ecuación diferencial

y y
y0 = + xm y n f ( )
x x
se transforma en una ecuación de variables separables usando el cambio de variable y = vx.
Use la parte anterior para encontrar la solución de

dy y sec2 (y/x)
= +
dx x y2

(ii) Resuelva las siguientes ecuaciones diferenciales:

p
(a) (1.5 pts.) y 2 dx + (x y 2 − x2 − xy)dy = 0
x − 4y − 9
(b) (1.5 pts.) y 0 = −
4x + y − 2
SOLUCIÓN:
(i) Cambio de variable y = vx. Derivando c/r a x: y 0 = v 0 x + v. Luego, sustituyendo en la EDO

v 0 x + v = v + xm (vx)n f (v)

de donde
v 0 = xm+n−1 v n f (v)

Asumiendo v n f (v) 6= 0, se tiene

dv
= xm+n+1 dx
v n f (v)
.................................................................................................................................................... (1.5 pts.)
Aplicando ahora a la EDO con m = 0, n = −2, y teniendo en cuenta que v −2 sec2 (v) 6= 0 obtenemos:

dv
= x−3 dx,
v −2 sec2 (v)

el cual es equivalente a:
1
v 2 cos2 (v)dv = dx
x3
.................................................................................................................................................... (0.3 pts.)
Integrando Z
(1 + cos(2v)) 1
v2 dv = − 2 + c
2 2x
.................................................................................................................................................... (0.3 pts.)
Luego de integrar por partes llegamos a:

v3 1 1 1 1
+ ( v 2 sen2 (2v) − sen(2v) + v cos(2v)) = − 2 + c
6 4 2 2 2x

1
Facultad de Ciencias Fı́sicas y Matemáticas Universidad de Chile

Reemplazando el valor de v = y/x obtenemos

(y/x)3 1 1 1 1
+ ( (y/x)2 sen2 (2y/x) − sen(2y/x) + y/x cos(2y/x)) = − 2 + c
6 4 2 2 2x
.................................................................................................................................................... (0.7 pts.)

(ii) (a) La EDO se puede escribir como:

dy y 2 /x2
= −p
dx (y/x)2 − 1 − y/x

.................................................................................................................................................... (0.2 pts.)

Cambio de variable: v = y/x, de donde y 0 = v 0 x + v. Luego

−v 2
v0 x + v = √
v2 − 1 − v
el cual es equivalente a: √
0 −v v 2 − 1
vx= √
v2 − 1 − v
.................................................................................................................................................... (0.5 pts.)

√
Soluciones constantes: v v 2 − 1 = 0. Esto implica que v = 0 ó v 2 − 1 = 0, de esto último v = ±1. Por tanto,
las soluciones constantes son:
y = 0, y = −x, y = x.

.................................................................................................................................................... (0.3 pts.)

√
Asumiendo v v 2 − 1 6= 0, tenemos  
1 1 1
− −√ dv = dx
v 2
v −1 x
Integrando
p
− ln(v) + ln(v + v 2 − 1) = ln(cx)

el cual es equivalente a: √
v+ v2 − 1
= cx
v
Sustituyendo el valor de v = y/x se tiene
p
y+ y 2 − x2
= cx
y

.................................................................................................................................................... (0.5 pts.)

(b) Sean L1 : x − 4y − 9 = 0 y L2 : 4x + y − 2 = 0. Tales rectas no son paralelas.

Punto de intersección: x0 = 1, y0 = −2.
.................................................................................................................................................... (0.2 pts.)

2
Facultad de Ciencias Fı́sicas y Matemáticas Universidad de Chile

Cambio de variable: ξ = x − 1, η = y + 2, de donde dξ = dx, dη = dy

.................................................................................................................................................... (0.2 pts.)

Luego
dη ξ − 4η 1 − 4η/ξ
=− =−
dξ 4ξ + η 4 + η/ξ
.................................................................................................................................................... (0.2 pts.)

Sea z = η/ξ, entonces η 0 = z 0 ξ + z. Esto implica

z2 + 1
z0ξ = −
z+4
.................................................................................................................................................... (0.4 pts.)

Como z 2 + 1 6= 0, tenemos:  
z 4 dξ
+ dz = −
z2 + 1 z2 + 1 ξ
Integrando:
1
ln(z 2 + 1) + 4 arctan(z) = − ln(ξ) + c
2
Sustituyendo z = η/ξ = (y + 2)/(x − 1) se obtiene:
 2 !  
1 y+2 y+2
ln + 1 + 4 arctan = − ln(ξ) + c
2 x−1 x−1

.................................................................................................................................................... (0.5 pts.)

3
Facultad de Ciencias Fı́sicas y Matemáticas Universidad de Chile

P2. (a) Considere la ecuación diferencial

(7x4 y − 3y 8 )dx + (2x5 − 9xy 7 )dy = 0

(i) (1 pto.) Determine m y n de modo que µ(x, y) = xm y n sea factor integrante.

(ii) (2 pts.) Encuentre su solución general, usando el factor integrante de la parte (i).
(b) Considere la siguiente ecuación de Ricatti:
y2 4mn
y0 = − y − (m − n)2
(m + n)2 (m + n)2
(i) (1 pto.) Encuentre una solución particular constante.
(ii) (2 pts.) Encuentre la solución general.
SOLUCIÓN:

(a) Se tiene la ecuación diferencial

(7x4 y − 3y 8 )dx + (2x5 − 9xy 7 )dy = 0

(i) multiplicando por el factor integrante, resulta:

(7x4+m y n+1 − 3xm y 8+n ) dx + (2x5+m y n − 9x1+m y 7+n ) = 0

| {z } | {z }
M (x,y) N (x,y)

Imponemos que la nueva ecuación diferencial sea exacta:

∂M ! ∂N
=
∂y ∂x
...................................................................................................................................... (0.3 pts.)
∂M ∂N
Luego calculamos , :
∂y ∂x
∂M
= 7(1 + n)x4+m y n − 3(8 + n)xm y 7+n
∂y
∂N
= 2(5 + m)x4+m y n − 9(1 + m)xm y 7+n
∂x
...................................................................................................................................... (0.3 pts.)
A partir de lo anterior:
∂M ∂N
=
∂y ∂x
4+m n m 7+n
7(1 + n)x y − 3(8 + n)x y = 2(5 + m)x4+m y n − 9(1 + m)xm y 7+n

[7(1 + n) − 2(5 + m)]x4+m y n + [9(1 + m) − 3(8 + n)]xm y 7+n = 0

Se concluye el sistema de ecuaciones:

7n − 2m = 3
n − 3m = −5

...................................................................................................................................... (0.3 pts.)

El cual tiene por solución m = 2, n = 1.

Ası́, el factor integrante es µ(x, y) = x2 y

...................................................................................................................................... (0.1 pts.)

4
Facultad de Ciencias Fı́sicas y Matemáticas Universidad de Chile

(ii) La nueva ecuación diferencial es

(7x6 y 2 − 3x2 y 9 ) dx + (2x7 y − 9x3 y 8 ) = 0

| {z } | {z }
M (x,y) N (x,y)

Ya que es exacta, ∃E : Ω ⊂ R2 → R tal que:

∂E
(x, y) = M (x, y) (1)
∂x
∂E
(x, y) = N (x, y) (2)
∂y
E(x, y) = k (3)

para y fijo, integramos (1) sobre x:

Z
E(x, y) = (7x6 y 2 − 3x2 y 9 )dx + g(y) = x7 y 2 − x3 y 9 + g(y)

...................................................................................................................................... (0.5 pts.)

y derivando respecto a y, para luego usar (2) (1.0):

∂E ∂ 7 2
(x, y) = (x y − x3 y 9 ) + g 0 (y)
∂y ∂y
N (x, y) = 2x7 y − 9x3 y 8 + g 0 (y)
2x7 y − 9x3 y 8 = 2x7 y − 9x3 y 8 + g 0 (y)

...................................................................................................................................... (1.2 pts.)

De lo que se determina g(y) = k1 .
Usando (3), se concluye:

x7 y 2 − x3 y 9 + k1 = k
7 2 3 9
x y −x y = k

...................................................................................................................................... (0.5 pts.)

5
Facultad de Ciencias Fı́sicas y Matemáticas Universidad de Chile

(b) Se tiene la ecuación de Ricatti

y2 4mn
y0 = − y − (m − n)2
(m + n)2 (m + n)2

(i) Consideremos yp = k, constante. Reemplazando en la ecuación la determinamos:

k2 4mn
0 = − k − (m − n)2
(m + n)2 (m + n)2
0 = k 2 − 4mnk − (m2 − n2 )2
p
4mn ± 16m2 n2 + 4(m2 − n2 )2
k =
p 2
k = 2mn ± 4m n + m4 − 2m2 n2 + n4
2 2
p
k = 2mn ± m4 + 2m2 n2 + n4
p
k = 2mn ± (m2 + n2 )2
k = 2mn ± (m2 + n2 )

Ası́, se obtienen dos soluciones: k1 = (m + n)2 , k2 = −(m − n)2 .

...................................................................................................................................... (1.0 pts.)
(ii) Ya que es una EDO de Ricatti, usando como solución particular yp = (m + n)2 y el cambio de
1
variable y(x) = (m + n)2 + , se obtiene:
z(x)
!
0 4mn 2 2 1
z − 2
− 2
(m + n) z = −
(m + n) (m + n) (m + n)2
2(m2 + n2 ) 1
z0 + z = −
(m + n)2 (m + n)2
| {z }
A

...................................................................................................................................... (1.2 pts.)

Esta última ecuación se puede resolver con factor integrante, el cual es eAx :

eAx
Z
d
(zeAx ) = − / dx
dx (m + n)2
eAx
zeAx = − +k
A(m + n)2
1
z(x) = − + ke−Ax
2(m + n2 )
2

...................................................................................................................................... (0.6 pts.)

Finalmente, se encuentra y(x)
1
y(x) = (m + n)2 +
ke−Ax − 2(m21+n2 )

...................................................................................................................................... (0.2 pts.)

6
Facultad de Ciencias Fı́sicas y Matemáticas Universidad de Chile

P3. (a) (2 pts.) Muestre que la solución de problema y 00 + 2y 0 = 0, y(0) = a, y 0 (0) = b tiende a una constante
cuando x → ∞. Encuentre dicha constante.
(b) (4 pts.) Un tanque de forma de un cilindro dispuesto en forma horizontal está inicialmente lleno de agua.
La altura del cilindro es H0 y el radio r. El tanque tiene un orificio en el fondo cuyo dı́ametro es ρ. Se
abre el orificio y el lı́quido cae libremente.
(i) (2 pts.) Encuentre la ecuación diferencial que describe el modelo.
(ii) (1 pto.) Solucione la ecuación diferencial.
(iii) (1 pto.) Sabiendo que h(0) = 2r, encontrar el valor de la constante de integración.

SOLUCIÓN:

(a) Usando el cambio de variable u(x) = y 0 (x), resulta

u0 + 2u = 0

ecuación cuyo factor integrante es e2x . Al integrarla, se obtiene

u(x) = k1 e−2x

................................................................................................................................................. (0.6 pts.)

La condición inicial y 0 (0) = b en términos de u(x) es u(0) = b. Al imponerla, da como resultado k1 = b.
Reemplazando esto en la ecuación y volviendo a la variable original:

y 0 (x) = be−2x

Integrando:
b
y(x) = − e−2x + k2
2
................................................................................................................................................. (0.6 pts.)
Imponiendo la condición inicial y(0) = a, se obtiene k2 = a + 2b . Finalmente, la solución es:

b b
y(x) = − e−2x + a +
2 2
................................................................................................................................................. (0.4 pts.)

Ya que lı́m e−2x = 0 se tiene

x→∞
b
lı́m y(x) = a +
x→∞ 2
................................................................................................................................................. (0.4 pts.)

7
Facultad de Ciencias Fı́sicas y Matemáticas Universidad de Chile

(b) (i) Sea A0 y AT las areas transversales del orificio y del tanque respectivamente. Supongamos que el
tanque esta lleno hasta una altura h. Por la ley de Torricelli tenemos:
p
v(t) = 2gh(t),

con g la aceleración de la gravedad. La disminución del nivel del agua h(t) con el flujo de salida viene
dado por:
dh A0 A0 p
=− v(t) = − 2gh(t).
dt AT AT
Con los datos expuestos, tenemos las siguientes areas transversales:
ρ
A0 = π( )2 ,
2
AT = 2xH0 ,

donde H0 el largo del cilindro horizontal y x el ancho el cual es variable en cada tiempo. Como r es
el radio del cilindro, entonces

(x − 0)2 + (h − r)2 = r2 , ⇒ x2 + h2 − 2rh + r2 = r2 ,

de donde
p
x= 2rh − h2 .

Con esto
p
AT = 2 2rh − h2 H0 ,

y ası́
dh πρ2 p πρ2 p
=− √ 2gh = − √ 2g.
dt 2
4 2rh − h H0 4 2r − hH0
...................................................................................................................................... (2 pts.)
(ii) Luego,
√ πρ2 p
2r − h dh = − 2gdt.
4H0
Esta ecuación tiene por solución:

2 πρ2 p
(2r − h)3/2 = 2g t + c.
3 4H0

...................................................................................................................................... (1 pto.)
(iii) Teniendo en cuenta que h(0) = 2r, obtenemos el valor de c = 0.
...................................................................................................................................... (1 pto.)