Posgrado en Ciencia e Ingenierı́a de Materiales 2 de octubre de 2013

Fundamentos de Matemáticas
Ecuaciones diferenciales

Solución: Tarea 4 (Total: 18 puntos)

II.2. Ecuaciones diferenciales de primer orden
La ecuación de Ricatti es una ecuación no-lineal
dy
= P (x) + Q(x)y + R(x)y 2 (II.1)
dx
llamada ası́ en honor del matemático y filósofo italiano, conde Jacobo Francesco Ricatti (1676-1754). En muchos
casos, dependiendo de lo que sean P (x), Q(x) y R(x) la solución de II.1 no puede ser expresada en términos de
funciones elementales.

1. Si y1 es una solución particular conocida de la ecuación de Ricatti, demuestre que y = y1 + u es una

familia de soluciones de II.1, en donde u es la solución de
du
− (Q + 2y1 R)u = Ru2 .
dx
(2 puntos)
0
Si y = y1 + u, entonces y = y10 0 2
+ u , y por lo tanto dy/dx = P (x) + Q(x)y + R(x)y se transforma en

y10 + u0 = P + Q(y1 + u) + R(y1 + u)2

⇔ y10 + u0 = P + Qy1 + Ry12 + Qu + 2y1 Ru + Ru2 .

Puesto que y1 es solución de la ecuación de Ricatti, obtenemos

u0 − (Q + 2y1 R)u = Ru2 .

2. Encuentre la solución del problema con valor inicial

dy
= (y 2 + 1)t
dt
y(0) = 1.

(2 puntos)
Por separación de variable se obtiene
Z Z
1
dy = tdt + c (II.2)
y2 + 1
tan−1 y = t2 /2 + c (II.3)
y(t) = tan(t2 /2 + c). (II.4)

La condición inicial implica la solución al problema

y(t) = tan(t2 /2 + π/4).

3. Determine la solución de la ecuación diferencial

dy 2
+ y = t − 1.
dt t
(2 puntos)

V. Moock: [email protected] UNAM 2014-1

Posgrado en Ciencia e Ingenierı́a de Materiales 2 de octubre de 2013

Esta ecuación es lineal con a(t) = 2/t y r(t) = t − 1. Primero calculamos el factor de integración
2
R R
a(t)dt 2/tdt
µ(t) = e =e = e2 ln t = eln(t )
= t2 .

Recordamos que la idea detrás de este método es multiplicar ambos lados por µ(t), de manera que el lado
izquierdo de la nueva ecuación sea el resultado de la regla del producto. En este caso, al multiplicar por
µ(t) = t2 , resulta
dy
t2 + 2ty = t2 (t − 1).
dt
Observe que el lado izquierdo es la derivada del producto de t2 y y(t). Entonces, esta ecuación es la misma
que
d 2
(t y) = t2 (t − 1) = t3 − t2 .
dt
Integrando ambos lados con respecto a t, obtenemos
t4 t3
t2 y = − + c,
4 3
donde c es una constante arbitraria. La solución general es
t2 t c
y(t) = − + 2.
4 3 t
Por supuesto, podemos verificar que esas funciones satisfacen la ecuación diferencial sustituyéndolas de
nuevo en la ecuación.
4. Decaimiento radioactivo. El isótopo radioactivo torio 234 se desintegra con una rapidez proporcional
a la cantidad presente del mismo. Sea Q(t) la cantidad de torio 234 presente en cualquier instante t, en
donde Q se mide en miligramos y t en dı́as. La observación fı́sica de que el torio 234 se desintegra con
una rapidez proporcional a la cantidad presente significa que la razón de cambio con el tiempo dQ/dt es
proporcional a Q; por lo tanto satisface la ecuación diferencial
dQ
= −rQ (II.5)
dt
en donde la constante r > 0 se conoce como la razón de decaimiento. Se busca la solución de II.2 que
también satisfaga la condición que 100 mg de este material se reduce a 82.04 mg en una semana. Halle el
intervalo de tiempo que debe transcurrir para que la masa decaiga hasta la mitad de su valor original.
(4 puntos)
Se busca la solución al problema II.2 que también satisfaga la condición inicial

Q(0) = 100

ası́ como la condición de la frontera

Q(7) = 82,04.
La ecuación II.2 es lineal y también separable; su solución general es

Q(t) = ce−rt ,

donde c es una constante arbitraria. La condición inicial requiere que c = 100, por lo tanto

Q(t) = 100e−rt .

A fin de satisfacer la condición de la frontera se da

82,04 = 100e−7r ,

entonces
ln 0,8204
r=− = 0,02828 dı́as−1 .
7
Por lo tanto, se ha determinado la razón de decaimiento r. Con este valor se obtiene que

Q(t) = 100e−0,02828t mg.

V. Moock: [email protected] UNAM 2014-1

Posgrado en Ciencia e Ingenierı́a de Materiales 2 de octubre de 2013

El lapso durante el que la masa se reduce a la mitad de su valor original se domina vida media del
material. Sea τ el tiempo en que Q(t) es igual a 50 mg. Entonces, por

50 = 100e−rt

o bien
rτ = ln 2.
Usando el valor de r se encuentra que, para el torio 234,
ln 2
τ= ≈ 24,5 dı́as.
0,02828

5. Resuelve la ecuación diferencial

(y cos x + 2xey ) + (senx + x2 ey − 1)y 0 = 0.

(2 puntos)
Es fácil ver que
My (x, y) = cos x + 2xey = Nx (x, y)
de modo que la ecuación dada es exacta. Por lo tanto existe una función f (x, y) tal que

fx (x, y) = y cos x + 2xey

fy (x, y) = senx + x2 ey − 1.

Al integrar la primera de estas ecuaciones se obtiene

f (x, y) = ysen + x2 ey + h(y). (II.6)

Si se hace fy = N da
fy (x, y) = senx + x2 ey + h0 (y) = senx + x2 ey − 1.
Por lo tanto, h0 (y) = −1 y h(y) = −y. Puede omitirse la constante de integración ya que cualquier solución
de la ecuación precedente es satisfactoria; no se requiere la más general. Al sustituir h(y) en la ecuación
II.3 da
f (x, y) = ysenx + x2 ey − y.
De donde la solución al problema queda dada implı́citamente por

ysenx + x2 ey − y = c.

6. Resuelve la ecuación diferencial

dy y 2 + 2xy
= .
dx x2
(2 puntos)
Al escribir esta ecuación como
dy y y
= ( )2 + 2
dx x x
se demuestra que es homogénea. Las variables no pueden separarse, la ecuación no es exacta ni tiene factor
integrante obvio. Por tanto, se ve uno conducido a la sustitución y = xv, que transforma la educación dada
en
dv
x + v = v 2 + 2v.
dx
De donde,
dv
x = v2 + v
dx
o bien, al separar las variables,
dx dv
= .
x v(v + 1)

V. Moock: [email protected] UNAM 2014-1

Posgrado en Ciencia e Ingenierı́a de Materiales 2 de octubre de 2013

Si se desarrolla el segundo miembro por fracciones parciales, se obtiene

dx 1 1
=( − )dv.
x v v+1
Al integrar los dos miembros, se llega a

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

en donde c es una constante arbitraria. Ası́, al combinar los logaritmos y tomar la exponencial de ambos
miembros, se obtiene
v
cx = .
v+1
Por último, al sustituir v en términos de y de la solución del problema en la forma
y/x y
cx = = .
(y/x) + 1 y+x
Al despejar y se obtiene
cx2
y= .
1 − cx
7. Resuelve la ecuación diferencial
dy y 5
− = − x2 y 3 . (II.7)
dx x 2
(2 puntos)
−2 0 −3 0
Aquı́ n = 3, ası́ que hacemos z = y , z = −2y y y multiplicamos la ecuación II.4 por −2y para obtener
una ecuación lineal de primer orden z 0 + a(x)z = f (x)
2
z0 + z = 5x2 , (II.8)
x
con a(x) = 2/x, y ası́ a(x)dx = 2 1/xdx = 2 ln x = ln x2 . Multiplicando ambos lados de la ecuación por
R R
2
eln x = x2 e integrando se llega a
Z
x2 y = 5x2 · x2 dx + c = x5 + c.

Entonces, la solución al problema II.5 es

y −2 = z = x3 + cx−2 ,

y por lo tanto la solución al problema original vale

y = (x3 + cx−2 )−1/2 .

8. Considere el problema de valor inicial

x0 (t) = x(t)
x(0) = 1.

Emplee el método de Picard para encontrar las aproximaciones x1 , x2 , x3 y xn . (2 puntos)

Como sabemos, la ecuación arriba tiene la solución única x(t) = et . En este caso, la función f (t, x) de la
forma general x0 (t) = f (t, x(t)), x(t0 ) = x0 está dada por f (t, x(t)) = x(t), ası́ que las iteradas de Picard
dan sucesivamente

x0 (t) = x0 = 1,
Z t
x1 (t) = 1 + 1ds = 1 + t,
0
Z t
t2
x2 (t) = 1 + 1 + sds = 1 + t + ,
0 2
Z t 2
s t2 t3
x3 (t) = 1 + 1 + s + ds = 1 + t + + ,
0 2 2! 3!

V. Moock: [email protected] UNAM 2014-1

Posgrado en Ciencia e Ingenierı́a de Materiales 2 de octubre de 2013

y claramente,
n
t2 tn X tk
xn (t) = 1 + t + + ... + = .
2! n! k!
k=0

Luego
∞ k
X t
lı́m xn (t) = = et .
n→∞ k!
k=0

V. Moock: [email protected] UNAM 2014-1