SOLUCION PES ECUACIONES DIFERENCIALES

Ingeniera Civil Codigo 10008-11039-19003

Segundo Semestre 2005

Problema 1. Un tanque de 100 litros de capacidad contiene inicial-

mente 25 litros de una soluci on salina con concentraci
on de 1 gramo
de sal por litro. En el instante t = 0 , comienza a entrar solucion con
concentracion 43 gramos por litro y a una velocidad de 3 litros por mi-
nuto. Al mismo tiempo empieza a salir soluci on del tanque a 1 litro
por minuto. Cuando el tanque llega a los tres cuartos de su capaci-
dad se abre una segunda llave que permite la entrada de agua pura al
tanque a 1 litro por minuto. Determine la cantidad de sal que hay en
el tanque en todo tiempo inferior a 33 minutos.

Soluci on. Sean V (t) y x(t) , respectivamente, la cantidad de litros de solucion

y la cantidad gramos de sal que hay en el tanque en el instante t . Entonces
V (0) = 25 y como la concentracion inicial de sal es 1 gramo por litro, tenemos
que x(0) = 25.
El valor de V (t) desde t = 0 hasta el instante t0 en que V (t0 ) = 34 100 = 75
esta dado por la formula
V (t) = 25 + (3 1)t = 25 + 2t .
Entonces
75 = V (t0 ) = 25 + 2t0 = t0 = 25
Por lo tanto la ecuacion diferencial para 0 t 25 es
dx 3 x(t)
= 3 ,
dt 4 V (t)
es decir
dx 9 x(t)
=
dt 4 25 + 2t
La solucion general de la ecuacion es
3 1
x(t) = (25 + 2t) + K (25 + 2t) 2 .
4
Imponiendo la condicion inicial, obtenemos
3 1 75 K 125
25 = x(0) = 25 + + K (25) 2 = + = K = .
4 4 5 4
Luego para 0 t 25 tenemos
3 125 1
x(t) = (25 + 2t) + (25 + 2t) 2 .
4 4
La ecuacion diferencial para t 25 es
dx 9 x(t)
= ,
dt 4 V (t)
pero ahora
V (t) = 75 + 3 (t 25) = 3 t .
Luego nuestra ecuacion es
dx 9 x(t)
= ,
dt 4 3t
1
con condicion inicial
3 125 1 35
x(25) = (25 + 50) + (25 + 50) 2 =
4 4 2 3
La solucion general de la ecuacion es
27 1
x(t) = t + K t 3 .
16
Imponiendo la condicion inicial, obtenemos
 
35 27 13 1 35 675
= x(25) = 25 + + K (25) = K = (25) 3 .
2 3 16 2 3 16
Luego para t 25
 
27 1 35 675 1
x(t) = t + 3 (25) 3 t 3.
16 2 3 16
Entonces como V (33) = 3 33 = 99 es menor que la capacidad maxima (100
litros), la solucion anterior vale para 25 t 33.
Problema 2. a) Usando la definicion de convolucion o aplicando transfor-
mada de Laplace demuestre las tres igualdades siguientes
a1) eat eat = t eat ,
a2) 1 eat eat = a12 [1 eat + at eat ] ,
a3) t eat eat = a13 [2 + at 2eat + at eat ].

b) Usando transformada de Laplace y convolucion, resuelva el problema de

valores iniciales
ty 00 (t) + (4t 2) y 0 (t) 4 y(t) = f (t) , y(0) = 1 ,

0
para 0 t < 2,
donde f (t) = 1 4(t2)
4
1 e para t 2.

Solucion. a) Tenemos
   
1  at at
 1 1 1 a
L 1 e + at e (s) = 2 +
a2 a s sa (s a)2
1 (s a)2 s(s a) + as
 
= 2
a s(s a)2
1
= L 1 eat eat (s) .


=
s(s a)(s a)
lo que prueba la primera igualdad.
Tambien
   
1  at at
 1 2 a 2 a
L 2 + at 2e + at e (s) = 3 + 2 +
a3 a s s sa (s a)2
1 2s(s a)2 + a(s a)2 2s2 (s a) + as2
 
= 3
a s2 (s a)2
1
= L t eat eat (s) .


= 2
s (s a)(s a)

b) Poniendo Y (s) = L(y(t))(s) tenemos

 
1 4(t2) 1 1 1
u(t 2) (s) = e2s



L(f (t))(s) = L 1 e
4 4 s s+4
1
= e2s ,
s(s + 4)
L(4y(t))(s) = 4Y (s) ,
L(2y 0 (t))(s) = 2 [sY (s) 1] = 2sY (s) + 2 ,

L(4ty 0 (t))(s) = 4 (sY (s) 1) = 4sY 0 (s) 4Y (s) ,
s

L(ty 00 (t))(s) = s2 Y (s) s y 0 (0) = 2sY (s) s2 Y 0 (s) + 1 .



s
Luego aplicando transformada de Laplace a la ecuacion y dividiendo por s(s+4)
se obtiene la ecuacion de primer orden
4(s + 2) 3 1
Y 0 (s) + ; , Y (s) = e2s 2
s(s + 4) s(s + 4) s (s + 4)2
Resolviendo la ecuacion se obtiene
s2 + 6 C 1 1
Y (s) = 2
+ 2 2
+ e2 s 2 ,
s(s + 4) s (s + 4) 2 s (s + 4)2
1 4 6 C 1 1
= 2
+ 2
+ 2 2
+ e2 s 2
s+4 (s + 4) s(s + 4) s (s + 4) 2 s (s + 4)2
= L(e4t )(s) 4 L(e4t e4t )(s) + 6 L(1 e4t e4t )(s)
1
+ C L(t e4t e4t )(s) + e2 s L(t e4t e4t )(s)
2
= L(e4t )(s) 4 L(e4t e4t )(s) + 6 L(1 e4t e4t )(s)
1
+ C L(t e4t e4t )(s) + e2 s L(t e4t e4t )(s)
2
y nuestra solucion es
y(t) = e4t 4 e4t e4t + 6 (1 e4t e4t ) + C (t e4t e4t )
1 
+ [t e4t e4t ] u(t 2)

2 t=(t2)
6  C 
= e4t 4te4t + 1 e4t 4te4t 2 4t 2e4t 4te4t
 
16 64
1 
2 4(t 2) 2e4(t2) 4(t 2)e4(t2) u(t 2)


128
1 
12 C + 4 C t + (20 + C) e4t 2(88 C)te4t

=
32
1 
5 2t + 3e4(t2) 2te4(t2) u(t 2) .


64
Problema 3. Para x > 0 y constante fijas A y , considere la ecuaci
on
de Euler
x2 y 00 + (2A + 1) x y 0 + (A2 + ) y = 0 .
a) Escriba la soluci
on general para < 0, para = 0 y para > 0.
b) Usando lo anterior encuentre los autovalores y las autofunciones del
problema de contorno
 2 00
x y + (2A + 1) x y 0 + (A2 + ) y = 0
y(1) = y(e3 ) = 0
c) Usando m
etodo de separaci
on de variables resuelva el problema
x2 uxx + 5 x ux + 4 u = utt , 1 x e3 , t 0
u(1, t) = u(e3 , t) = 0 , t 0
u(x, 0) = 0 , ut (x, 0) = x2 , 1 x e3 .

Soluci
on. a) La ecuacion caracterstica
k(k 1) + (2A + 1) k + (A2 + ) = 0 ,
se reduce a
(k + A)2 + = 0 .
Luego la solucion general para < 0 es
h i
A
y(x) = x c1 x + c2 x ,
para = 0 es
y(x) = xA [c1 + c2 ln(x)] ,
y para > 0 es
h    i
y(x) = xA c1 cos ln(x) + c2 sen ln(x) .

b) Imponiendo las condiciones de frontera en la solucion general para 0,

obtenemos como u nica solucion y(x) 0. Luego debemos tomar > 0 y la
solucion general
h    i
y(x) = xA c1 cos ln(x) + c2 sen ln(x) .
Entonces imponiendo las condiciones de frontera se obtienen los autovalores y
autofunciones son
n2 2 A
 n 
n = , yn (x) = x sin ln(x) ,
9 3
para todo n 1.

c) Pongamos u(x, t) = M (x)N (t). Reemplazando en la ecuacion, separando

variables e imponiendo las condiciones de frontera se obtienen las ecuaciones
 2 00
x M (x) + 5 x M 0 (x) + (4 + )M (x) = 0
N 00 (t) + N (t) = 0 .
M (1) = M (e3 ) = 0
Luego los autovalores y autofunciones son
n2 2  n 
n = , Mn (x) = x2 sin ln(x) ,
9 3
para todo n 1.
Reemplazando en la ecuacion para la variable t obtenemos
n2 2
Nn00 (t) + N (t) = 0 ,
9
cuya solucion general es
n  n 
Nn (t) = an cos t + bn sin t , n 1.
3 3
Tenemos entonces la solucion formal
X h n   n i  n 
2
u(x, t) = x an cos t + bn sin t sin ln(x) .
n=1
3 3 3

Consideremos ahora la condicion inicial

X  n 
2
0 = u(x, 0) = x an sin ln(x) .
n=1
3
Luego
an = 0 , n 1,
y

X n   n 
2
u(x, t) = x bn sin t sin ln(x) .
n=1
3 3
Finalmente la condicion

2 2
X n  n 
x = ut (x, 0) = x bn sin ln(x) ,
n=1
3 3
implica Z 3
n 2  n  2
bn = sin s ds = [1 (1)n ] .
3 3 0 3 n
Por lo tanto
6
bn = [1 (1)n ] ,
n22
y nuestra solucion es

6 2 X 1 n
n   n 
u(x, t) = 2 x 2
[1 (1) ] sin t sin ln(x) .
n=1
n 3 3