Aplicaciones de las

Ecuaciones Diferenciales

Gloria Aguilar
Natalia Boal
Carmelo Clavero
Francisco Gaspar

Departamento de Matemática Aplicada

Universidad de Zaragoza

1
Selección de ejemplos obtenidos de los libros de Nagle-Saff, Braun, Finzio-Ladas, Campbell,
entre otros.

Trayectorias ortogonales

1. En ingenierı́a se presenta a menudo el problema geométrico de encontrar una familia de curvas

(trayectorias ortogonales) que intersequen ortogonalmente en cada punto a una familia dada
de curvas. Por ejemplo, es posible que se den las lı́neas de fuerza y se pida obtener la ecuación de
las lı́neas equipotenciales. Consideremos la familia de curvas descrita por la ecuación F (x, y) = k,
donde k es un parámetro real.
i) Usando diferenciación implı́cita, demuestra que, para cada curva de la familia, la pendiente
está dada por
∂F
dy
= − ∂x .

dx ∂F
∂y
ii) Usando que la pendiente de una curva ortogonal (perpendicular) a una curva es la inversa
de la pendiente de la curva dada, demuestra que las curvas ortogonales a la familia F (x, y) = k
satisfacen la ecuación diferencial
∂F ∂F
(x, y)dx − (x, y)dy = 0.
∂y ∂x
iii) Utilizando la ecuación diferencial precedente, demuestra que las trayectorias ortogonales
de la familia de circunferencias x2 + y2 = k son rectas que pasan por el origen.
2. Usa el método del problema anterior para encontrar las trayectorias ortogonales de cada
familia de curvas dada. Dibuja conjuntamente la familia de curvas y sus trayectorias ortogonales.

i) xy = k, ii) 2x2 + y2 = k, iii) x2 − y2 = k.

Problemas de mezclas

3. En los problemas de mezclas se quiere calcular la cantidad de una sustancia, x(t), que hay
en un tanque en cada instante de tiempo t. Usando que la derivada de x respecto a t expresa la
razón de cambio de la sustancia presente en el tanque, se cumple la relación
dx
= velocidad de entrada - velocidad de salida .
dt
Dada la velocidad a la que un fluido que contiene la sustancia entra en el tanque y la concen-
tración de la sustancia, se cumple la relación

velocidad de flujo entrante × concentración = velocidad de entrada .

Suponiendo que la concentración de la sustancia es uniforme, para calcular la concentración se

divide x(t) por el volumen de la mezcla que hay en el instante t. As´ı,

velocidad de flujo saliente × concentración = velocidad de salida .

Por ejemplo, en un tanque de 1000 l (litros) una solución salada de salmuera empieza a fluir a
velocidad constante de 6 l/min. La solución se mantiene bien agitada y fluye hacia el exterior
con velocidad de 6 l/min. Sabiendo que la concentración de salmuera que entra es de 1 Kg/l,
2
calcula cuándo la concentración de sal será de 1/2 Kg/l. Si la velocidad de salida es 5 l/min,
determina la concentración de sal en función del tiempo.
4. La corriente sanguı́nea lleva un medicamento hacia el interior de un órgano a razón de
3cm3 /seg y sale de él a la misma velocidad. Se sabe que el volumen del órgano es 125 cm3
y la concentración del medicamento que entra es de 0.2 g/cm3 . ¿Cuál es la concentración del
medicamento en el órgano si inicialmente no habı́a vestigio alguno del medicamento? ¿Cuándo
la concentración será de 0.1 g/cm3 ?
5. El agua del rı́o Aguadulce fluye hacia el lago Magdalena a razón de 300 gal/min (1 galón
es una medida inglesa de capacidad que equivale a 4.5459 litros; en Estados Unidos equivale a
3.7853 litros). El lago Magdalena contiene aproximadamente 100 millones de galones de agua.
La fumigación de los naranjales cercanos ha ocasionado que la concentración de plaguicidas en
el lago llegue a ser de 35 partes por millón. Si se suspende la aplicación de plaguicidas, ¿cuánto
tiempo transcurrirá antes de que la concentración de los mismos en el lago esté por debajo de
10 partes por millón? (Se supone que el rı́o no contiene plaguicidas y que el volumen del lago
es constante).

Problemas de crecimiento de población

6. El modelo malthusiano de crecimiento de una población p(t) supone que la tasa de creci-
miento es proporcional a la población presente. Sabiendo que la población de Estados Unidos
en 1790 era de 3.93 millones y en 1800 de 5.31 millones, usa el modelo anterior para conocer
la población en función del tiempo. Este modelo supone que la tasa de mortalidad es nula,
que desde luego es errónea. Parece natural pensar que la tasa de mortalidad natural también
es proporcional al tamaño de la población. No obstante, debido a otros factores de mortalidad
(desnutrición, enfermedades, crı́menes violentos, etc), se puede suponer que la tasa de morta-
lidad es proporcional al número de interacciones bipartitas. Para una población de tamaño p,
existen p(p − 1)/2 interacciones de este tipo. Prueba que con esta hipótesis, el PVI que rige el
modelo tiene la forma
dp
= ap − bp2, p(0) = p 0 .
dt
Esta ecuación se conoce con el nombre de ecuación logı́stica. Calcula la población en función
del tiempo. En el caso particular de que la población en 1790 sea de 3.93 millones, en 1840 de
17.07 millones y en 1890 de 62.95 millones, determina la solución usando el modelo logı́stico.
7. Un modelo de crecimiento de población que se utiliza en predicciones actuariales se basa en
la ecuación de Gompertz
dp
= p(a − b ln p),
dt
con a y b constantes reales.
i) Resuelve la ecuación diferencial para calcular p(t).
ii) Si p(0) = p0 > 0, proporciona una fórmula para p(t) en términos de a, b, p0 y t.
iii) Describe el comportamiento de p(t) cuando t → ∞, considerando los casos b > 0 y b < 0.
8. Una bola de nieve se derrite de forma que la razón de cambio de su volumen es proporcional
al área de su superficie. Si el diámetro inicial de la bola es de 4 pulgadas y al cabo de 30
min. es de 3 pulgadas, ¿cuándo será de 2 pulgadas?, ¿cuando desaparecerá la bola de nieve?
Si suponemos que la bola se derrite tal que la razón de cambio de su diámetro es proporcional
al área de su superficie, con los mismos datos, ¿cuándo será su diámetro de 2 pulgadas?; en
términos matemáticos ¿cuándo desaparecerá la bola de nieve?.
9. Uno de los métodos más precisos para determinar la edad de restos arqueológicos es el
método del carbono 14 (C 14 ), basado en que para cualquier organismo vivo una proporción

3
constante de átomos de carbono está formada por el isótopo radiactivo C 14 . La proporción
permanece prácticamente constante durante toda la vida y cuando el organismo muere el C 14
sigue su proceso de desintegración, con lo cual la proporción disminuye. Un modelo simple para
describir el fenómeno supone que la cantidad de átomos que se desintegran es proporcional a la
cantidad de átomos presentes, es decir, se cumple la ecuación diferencial
dN
= −λN (t),
dt
siendo N (t) la cantidad de C14 en una muestra en el tiempo t. Suponiendo que N (0) = N0
calcula la solución del PVI. A R(t) = λN (t) se le llama tasa de desintegración y R(0) es
la tasa original de desintegración que coincide con la tasa de desintegración de la materia viva.
Sabiendo que la edad media del C14 (la edad media es el tiempo que debe transcurrir para que
se desintegre la mitad de la sustancia) es aproximadamente de 5600 años, resuelve los siguientes
casos:
i) El nivel de carbón vegetal extraı́do en las grutas de Lascaux (Francia) en 1950 dio una
medida de 0.91 desintegraciones por minuto y gramo, mientras que la materia viva da 6.68
desintegraciones. Calcula la época en que las grutas estuvieron habitadas.
ii) En la excavación en Ninpur (Babilonia), en 1950 el carbón vegetal de una viga dio
4.09 desintegraciones por minuto y gramo. Suponiendo que este carbón se formó en la época
Hamurabi, calcula una fecha probable de la sucesión de Hamurabi.
iii) En una cueva de Sudáfrica se encontró un cráneo humano junto a los restos de una fogata.
Los arqueólogos creen que la edad del cráneo es igual a la de la fogata. Sabiendo que solo un
2% de la cantidad original de C 14 queda en la madera, calcula la edad aproximada del cráneo.

Problemas de calefacción y de enfriamiento de edificios

10. Sea T (t) la temperatura interior de un edificio para el que la razón de cambio de temperatura
es la diferencia entre la razón a la que aumenta y la razón a la que disminuye. Suponemos
que afectan tres factores en la temperatura. El primero (calor producido por las personas,
luces y máquinas) incrementa la temperatura a razón de H(t). El segundo es el calentamiento
(enfriamiento) producido por la calefacción (aire acondicionado) a razón de U (t). El tercer factor
es el efecto de la temperatura exterior M (t), para el que la razón de cambio es proporcional a
la diferencia entre la temperatura exterior y la interior (ley de Newton del enfriamiento).
Recogiendo todo, la ecuación diferencial que modela el fenómeno es
dT
= k(M (t) − T (t)) + H(t) + U (t),
dt
con k constante real que depende de las propiedades del edificio (número de puertas y ventanas,
aislamiento, material, etc); al valor 1/k se le llama constante de tiempo del ediftcio. Resuelve
la ecuación diferencial para calcular la temperatura en función del tiempo.
i) Si la temperatura exterior es constante M = 0 y se cumple H = U = 0, escribe la solución
sabiendo que T (t0) = T0 (esto refleja como var´ıa la temperatura del edificio).
ii) Supongamos que H es constante H0 , U = 0 (no hay calefacción o enfriamiento) y la
temperatura exterior está dada por M (t) = M0 − B cos ωt, donde B es una constante real
positiva y ω = π/12 (onda senoidal de perı́odo 24 horas, con mı́nimo en t = 0 (medianoche) y
máximo en t = 12 (mediodı́a)). Calcula la solución sabiendo que a medianoche la temperatura
es T0 (esto refleja como varı́a la temperatura en primavera u otoño cuando no hay calefacción
ni aire acondicionado).
iii) En el apartado anterior ii) supongamos que hay un termostato para controlar la tem-
peratura del edificio en relación con la temperatura deseada Td . Suponiendo que el calen-
tamiento o enfriamiento suministrado es proporcional a la diferencia de temperatura, es decir,

4
U (t) = ku(Td T− (t)), calcula T (t) siendo de nuevo T (0) = T0 (esto refleja como var´ıa la tem-
peratura en verano (aire acondicionado) y en invierno (calefacción)).
11. En una mañana, mientras las personas trabajan, la calefacción mantiene la temperatura
interior a 210 C. A mediodı́a se apaga y se sabe que la temperatura exterior es de 120 C durante
la tarde. Sabiendo que la constante de tiempo del edificio es de 3 horas, ¿cuándo la temperatura
será de 160 C? Si se dejan abiertas algunas ventanas, la constante de tiempo se reduce a 2 horas;
¿cuándo se alcanzarán ahora los 160 C?.
12. La ley de Stefan de la radiación establece que la razón de cambio de la temperatura de
un cuerpo que se encuentra a T grados Kelvin en un medio a M grados Kelvin, es proporcional
a M 4 − T 4 . Plantea y resuelve la ecuación diferencial. Explica por qué las leyes de Newton
y Stefan son aproximadamente iguales cuando T se acerca a M y M es constante. (Ayuda:
factoriza M 4 − T 4).
13. Un vino tinto se saca de la bodega, que es un lugar fr´ıo a 100C, y se deja reposar en un
cuarto con temperatura de 230 C. ¿En qué momento la temperatura del vino será de 180 C si
transcurren 10 minutos para alcanzar los 150C?

Mecánica Newtoniana
14. A un objeto de masa m se le aplica una velocidad inicial dirigida hacia abajo v0 y cae bajo
la acción de la fuerza de la gravedad, que se supone constante, y la fuerza debida a la resistencia
del aire que es proporcional a la velocidad. Calcula la ecuación del movimiento del objeto.
15. Un paracaidista cuya peso (masa) es de 100 Kg se deja caer desde un helicóptero situado a
3000 metros. Suponiendo que la resistencia del aire es proporcional a la velocidad con constante
k1 = 20 kg/seg si el paracaı́das está cerrado y k1 = 100 kg/seg si está abierto y que el paracaı́das
se abre transcurridos 30 segundos, ¿cuándo llegará a la superficie?. Si el paracaı́das se abre al
cabo de 1 minuto, ¿cuándo llegará el paracaidista a la superficie?.
16. Vuelo de un cohete. Un cohete con masa inicial de m0 Kg se lanza verticalmente desde
la superficie de la Tierra. El cohete expele gas a razón de α Kg/seg y a una velocidad constante
de β m/seg relativa al cohete. Suponiendo que el campo gravitacional es constante de g kg/seg2,
la segunda ley de Newton da lugar a la ecuación
dv
(m0 − αt) − αβ = −g(m0 − αt),
dt
donde v = dx/dt es la velocidad del cohete, x es su altura respecto de la superficie de la Tierra
y m0 αt− es la masa del cohete a los t segundos del lanzamiento. Sabiendo que la velocidad
inicial es cero, resuelve la ecuación anterior para calcular la velocidad y altura del cohete para
0 ≤ t ≤ m0/α.

Problemas de reacciones quı́micas

17. Se considera una reacción quı́mica que ocurre en una solución bien mezclada. Se supone
que la reacción es irreversible, mientras que la temperatura y el volumen son constantes. Al
inicio, la reacción incluye dos reactivos A, B y un producto E; se supone que una molécula
de A y una de B producen una de E. Sean a, b las concentraciones en el tiempo t de A y B
respectivamente. Sabiendo que se satisface la ley de acción de masas que establece que la tasa
instantánea de producción de un producto por unidad de volumen es proporcional al producto
de las concentraciones de los reactivos, prueba que el PVI es
dx
= k(a(0) − x)(b(0) − x), x(0) = 0,
dt

5
donde k es la constante de reacción y a(0), b(0) son conocidos. Resuelve en el caso que
a(0) = 3, b(0) = 1 mol/l y k = 1 l/mol.s. Interpreta el resultado.
18. Supongamos que dos tanques A, B de volúmenes V1 , V2 están separados por una frontera,
a través de la cual una sustancia fluye proporcionalmente (con constante k) de la diferencia de
concentraciones C1 − C2 . Calcula la concentración en función del tiempo en cada tanques.

Aplicaciones de las EDO de segundo orden

19. Consideremos un sistema mecánico que consiste de un resorte (muelle) en espiral suspendido
de un soporte rı́gido con una masa sujeta al extremo (ver figura 1). Para analizar el fenómeno
usamos dos leyes: ley de Hooke y la segunda ley de Newton. La ley de Hooke establece
que el resorte ejerce una fuerza de restitución opuesta a la dirección de alargamiento del resorte
(peso=kl donde k es la constante de restitución y l es el alargamiento hasta la posición de
equilibrio). Para aplicar la ley de Newton, hay que tener en cuenta que sobre la masa m actúa
la fuerza de la gravedad (mg, g vale 9.8m/seg 2 ), la fuerza de restitución (− kx − mg donde x
es el alargamiento), la fuerza de amortiguación que es proporcional a la velocidad de la masa
(− b(dx/dt) con b la constante de amortiguación) y las fuerzas externas (f (t)). Ası́, la ecuación
diferencial que resulta es
d 2x dx
m 2 + b + kx = f (t).
dt dt
Si b = 0 el sistema es no amortiguado; en otro caso se dice amortiguado. Si f (t) = 0 se dice
libre; en otro caso se dice forzado. Sabiendo esto, resuelve los siguientes casos.
i) Una masa de 5 Kg se sujeta de un resorte y ocasiona un estiramiento de 0.5 m. A
continuación se tira de la masa 0.1 m abajo del punto de equilibrio y se le aplica una velocidad
dirigida hacia arriba de 0.1 m/seg. Calcula la ecuación del movimiento armónico simple de la
masa. ¿Cuando alcanzará la masa por primera vez su mı́nima altura después de haberse puesto
en movimiento?
ii) Una masa de 2 kg estira un resorte 49 cm al llegar al reposo en equilibrio. La constante
√
de amortiguación es 8 5N ·seg/m. Si la masa se tira 10 cm hacia abajo del punto de equilibrio
y se le aplica una velocidad de 2 m/seg dirigida hacia abajo, ¿cuál es el desplazamiento máximo
que alcanzará a partir de la posición de equilibrio?
iii) Una masa de 2 kg estira un resorte 20 cm al llegar al reposo en equilibrio. La constante
de amortiguación es 5N · seg/m. En t = 0 la masa se desplaza 5 cm hacia abajo y se suelta. En
el mismo instante se aplica una fuerza externa f (t) = 0.3 cos t N al sistema. Calcula la ecuación
de movimiento de la masa.
Sabiendo que la frecuencia de resonancia del sistema es γr /2π donde
.
k b2
γr = − ,
m 2m2

determina la frecuencia de resonancia del sistema. Se dice que un sistema está en resonancia
cuando es estimulado por una fuerza externa cuya frecuencia coincide con la frecuencia de
resonancia del sistema. En el caso de que la constante de amortiguamiento sea pequeña, si la
fuerza externa tiene una frecuencia cercana a la frecuencia del sistema, el sistema está sujeto a
grandes oscilaciones. Estas grandes vibraciones son las que preocupan a los ingenieros, puesto
que pueden ocasionar que las alas de los aviones se rompan, los puentes se desplomen o (menos
importante) que las copas de vidrio se hagan añicos.
20. Se considera un servomecanismo que modela a un piloto automático. Si y(t) es la dirección
real de la nave en el instante t y g(t) la dirección deseada, entonces el error o desviación es

6
e(t) = y(t) − g(t). La segunda ley de Newton, expresada en términos de momentos de las fuerzas
proporciona la ecuación
Iyjj(t) = −ke(t),
donde I es el momento de inercia del eje del timón y k es una constante positiva. Calcula el
error cuando g(t) = at con a constante, sabiendo que el eje del timón se encuentra inicialmente
en reposo. Supongamos ahora que para controlar las oscilaciones se suministra al eje del timón
una fuerza giratoria proporcional a ej (t), pero de signo opuesto. Prueba que ahora la ecuación
diferencial es
Iyjj(t) = −ke(t) − µej(t),
donde µ es una constante positiva. Determina ahora el error sabiendo que g(t) = a con a
constante y de nuevo el eje del timón se encuentra inicialmente en reposo.
21. Consideremos un circuito eléctrico elemental que consta de una fuerza electromotriz (baterı́a
o generador), una resistencia, una bobina y un condensador en serie (ver figura 2). Estos circuitos
se rigen por dos principios fı́sicos: la conservación de la carga eléctrica y la conservación de la
energ´ıa. Las leyes que formulan estos principios son las leyes de Kirchoff que establecen:
i) La corriente I que pasa a través de cada elemento de un circuito en serie debe ser la misma.
ii) La suma algebraica de los cambios instantáneos de potencial (caı́da de voltaje) en un
circuito cerrado debe ser cero.
En cada elemento la ca´ıda de voltaje es la siguiente:
a) Según la ley de Ohm, la caı́da de voltaje en una resistencia es proporcional a la corriente
I que pasa por ella, es decir, ER = RI donde R es la resistencia.
b) Según las leyes de Faraday y Lenz, la caı́da de voltaje en una bobina es proporcional
dI
a la razón de cambio instantánea de la corriente I que pasa por ella, es decir, EL = L donde
dt
L es la inductancia.
c) La caı́da de voltaje en un condensador es proporcional a la carga eléctrica q del conden-
1 1
sador, es decir, E = q donde C es la capacitancia y es la elastancia.
C
C C
dq
Se sabe que la carga y la corriente están relacionadas por I = . Ası́, la ecuación diferencial
dt
que rige el modelo es
d2q dq 1
L 2 +R + q = E(t),
dt dt C
donde E(t) es el voltaje aplicado por la fuerza electromotriz., o en términos de la corriente
d 2I dI + 1 I = dE(t)
L 2 + R .
dt dt C dt
A partir de esto calcula la corriente (amperios) y la carga (culombios) en los siguientes casos;
en todos ellos analiza la solución.
i) La resistencia es de 0.02 ohmios, la inductancia es de 0.001 henrios, la capacitancia es de
2 faradios, la fuerza electromotriz es E(t) = sen(100t) voltios y la corriente y carga inicial son
cero.
ii) La resistencia es de 120 ohmios, la inductancia es de 4 henrios, la capacitancia es de
2200−1 faradios, la fuerza electromotriz es E(t) = 10cos(2t) voltios y la corriente y carga inicial
son cero.
iii) La resistencia es de 10 ohmios, la inductancia es de 4 henrios, la capacitancia es de 0.01
faradios, la fuerza electromotriz es E(t) = E0cosγt voltios con E0, γ constantes y la corriente y
carga inicial son cero.
22. En el diseño de una planta de tratamiento de aguas residuales, se propuso el siguiente PVI
60 − H = 77.8Hjj + 19.42(Hj)2, H(0) = Hj(0) = 0,

7
donde H(t) es el nivel del fluido en una cámara de eyección y t es el tiempo en segundos. Usa
el método de Euler para aproximar H(t) en [0, 5].

Aplicaciones de los sistemas de EDO

23. Un material radiactivo A se descompone de acuerdo a la reacción A → B C con

constantes k1 y k2 respectivamente, siendo B, C el producto intermedio y final respectivamente.
Prueba que el sistema diferencial es →
dCA
= −k C ,
1 A
dt
dCB = k C − k C ,
1 A 2 B
dt
dCC = k C ,
2 B
dt
con CA, CB, CC las concentraciones de las tres sustancias. Sabiendo que k1 = 3, k2 = 1, CA(0) =
1, CB (0) = CC (0) = 0, calcula las concentraciones en función del tiempo. Representa la concen-
tración de las 3 sustancias en [0, 10].
24. En una reacción quı́mica intervienen 3 sustancias A, B, C. Se sabe que la constante de paso
de A a B es k1 , de B a A es k2 , de B a C es k3 y de C a B es k4 . Calcula la concentración de
cada sustancia, sabiendo que k1 = 1, k2 = 0, k3 = 2, k4 = 3, CA(0) = 1, CB(0) = CC (0) = 0, y
representa la concentración de las 3 sustancias en [0, 10].
25. En el estudio de la fermentación cinética que permite generar penicilina, la ley logı́stica
dy1
= k1 y1(1 − y 1/k 2),
dt
ha sido usada para describir la dinámica de crecimiento. El término (1 − y1 /k2 ) incluye los
efectos del cese de crecimiento por falta de nutriente. Se sabe que el grado de producción de
penicilina satisface la ecuación
dy2
= k3 y1 − k 4y 2,
dt
es decir, es proporcional al grado de concentración de células y se degrada por hidrólisis. Resuelve
el modelo en el intervalo [0, 12] usando el método de Euler e interpreta los resultados sabiendo
que k1 = 0.0312, k2 = 47.7, k3 = 3.374, k4 = 0.01268, y1(0) = 5, y2(0) = 0.
26. La conversión de glucosa a ácido glucónico es una simple oxidación producida por un
microorganismo en un proceso de fermentación. El mecanismo del proceso de fermentación es

glucosa + células → células (crecimiento celular) ,

glucosa + O2 → gluco... + H2 O2 (oxidación) ,
gluco... + H2 O → ácido glucónico (hidrólisis) ,
1
H2 O2 → H2 O + 2O2 (descomposición) .

Un modelo está dado por

dy1
= b 1y (1 1 ),2
1 − y /b
dt2
dy byy
= 3 1 4 − 0.9082b y ,
5 2
dt
dy b4 + y4
3
=b y,
5 2
dt . Σ
dy4 b3y 1y 4
= −1.011 ,
dt b4 + y4

8
donde y1 es la concentración de células, y2 es la concentración de gluco..., y3 es la concentracion de
ácido glucónico, y4 es la concentración de glucosa y bi , i = 1, · · ·, 5 son parámetros. Sabiendo que
b1 = 0.949, b2 = 3.439, b3 = 18.72, b4 = 37.51, b5 = 1.169, y1(0) = y4(0) = 5, y2(0) = y3(0) = 0,
resuelve el modelo en el intervalo [0, 9] usando el método de Euler e interpreta los resultados.
27. Sobre una superficie lisa horizontal se sujeta una masa de 2 Kg a una superficie vertical por
un medio de un resorte cuya constante es de 4 N/m. Otra masa de 1 kg se conecta al primer
objeto mediante un resorte con constante 2 N/m (ver figura 3). Sabiendo que los objetos se
desplazan 3 m hacia la derecha de sus posiciones de equilibrio y luego se sueltan, comprueba
que el sistema diferencial que rige el modelo es
d 2x
2 + 6 x − 2y = 0,
2
2dt
dy
+2 2=0
y− x ,
dt2
dx dy
x(0) = y(0) = 3, (0) = = 0.
dt dt
Calcula un sistema diferencial de primer orden equivalente al dado. Resuelve el problema,
representa la solución e interpreta los resultados.
28. Se considera el circuito eléctrico dado en la figura 4. Demuestra que el sistema diferencial
que rige el modelo es
1 t
I + 64I2jj = −2 sen ,
9 1 24
1 jj jj
64 I + 9I − 64I =
=3 I2 + I3,
I1 0,
3 2
donde I1.I2, I3 son las intensidades que circulan por las distintas ramas de la red.
i) Elimina I1, I3 del sistema para probar que se cumple
(92)(642)I(4) + (82)(64)Ijj + I2 = 0.
2 2

ii) Resuelve el problema de ii). A partir de la solución de ii) calcula I1 , I3 .

29. Calcula el sistema diferencial y las condiciones iniciales asociadas a los circuitos de las
figuras 5 y 6, suponiendo que las corrientes iniciales son todas cero. Resuelve el problema para
obtener las corrientes que circulan por cada rama.
30. Dos grandes tanques de 100 galones cada uno se encuentran interconectados por medio de
tubos. El lı́quido fluye del tanque A hacia el tanque B a razón de 3 gal/min. y de B a A a razón
de 1 gal/min. El lı́quido contenido en cada tanque permanece bien agitado. Una solución de
salmuera con una concentración de 2 lb/gal de sal fluye hacia el tanque A a razón de 6 gal/min.
La solución (diluida) fluye hacia afuera del sistema del tanque A a razón de 4 gal/min y del
tanque B a razón de 2 gal/min. Si inicialmente el tanque A contiene sólo agua y el tanque B
contiene 200 lb de sal (1 libra equivale a 456.3 gramos), demuestra que el sistema diferencial que
rige el modelo es
dx 1 7
= 12 + y− x,
dt 100 100
dy 3 3
= x− y,
dt 100 100
x(0) = 0, y(0) = 200.
Resuelve el problema, representa las soluciones e interpreta los resultados.
Supongamos ahora que la solución no fluye de B hacia afuera, que solamente 1 gal/min fluye
de A a B y que sólo 4 gal/min de salmuera fluyen hacia el interior del sistema por el tanque

9
A. Manteniendo los restantes datos iguales, resuelve el problema, representa las soluciones e
interpreta los resultados.
31. Un famoso experimento de fı́sica llevado a cabo por Thomson para determinar la proporción
de carga respecto a la masa m para un electrón, usa las ecuaciones diferenciales

d2x + He dy
m = Ee,
dt2 2 dt
dx
dy
m 2 − He = 0,
dt dt

donde H es la intensidad del campo magnético y E es la intensidad del campo eléctrico que se
han aplicado.
i) Transforma el sistema anterior en un sistema equivalente de primer orden.
ii) Resuelve para el caso particular H = 4, e = 1, m = 2, E = 6.
iii) Resuelve usando transformada de Laplace, suponiendo que H = 1/4, e = 4, m = 1, E = 2
con condiciones iniciales x(0) = y(0) = 0, xj(0) = yj(0) = 1

Proyectos asociados a ecuaciones diferenciales

32. (Linealización de problemas no lineales). Un planteamiento útil para analizar una

ecuación no lineal consiste en estudiar su ecuación linealizada, la cual se obtiene reemplazando
los términos no lineales mediante aproximaciones lineales.
i) Por ejemplo, la ecuación
d2θ
+ sen θ = 0,
dt2
que rige el movimiento de un péndulo simple, se aproxima por
d2θ
+ θ = 0.
dt2
Resuelve el problema lineal con θ(0) = π/12, θ j (0) = 0. Esta solución es una aproximación a la
solución del problema no lineal
ii) Da una linealización del PVI
. Σ
xjj(t) + 0.1 1 − x2(t) xj(t) + x(t) = 0, x(0) = 0.5, x‘(0) = 0.

Resuelve el problema lineal planteado.

33. (Limpieza de los grandes lagos). Supongamos que cada uno de los 5 grandes lagos de
Estados Unidos se considera como un depósito que contiene un lı́quido en el que se encuentra
disuelto un contaminante (DDT, fósforo, mercurio). Suponemos que el volumen de cada lago
es constante, las velocidades de flujo de agua son constantes durante todo el año, cuando un
lı́quido entra en un lago se produce una mezcla perfecta y que los contaminantes están disueltos
en el agua y entran o salen de los lagos por medio de afluencias o corrientes.
i) Usa las velocidades de las corrientes salientes dadas en la figura 7 para calcular el tiempo
que llevar´ıa “vaciar” cada lago. Esto proporciona una cota inferior para el tiempo que llevar´ıa
eliminar todos los contaminantes.
b) Para mejorar el resultado, suponemos que cada lago es un depósito separado al que fluye
solamente agua limpia. Con este planteamiento calcula el tiempo necesario para que el nivel de
contaminación en cada lago se reduzca a la mitad y al 5 por ciento del nivel original.
c) Par tener en cuenta que la contaminación fluye de un lago al siguiente, usa el modelo
compartimental de la figura 7 para calcular cuando el nivel de contaminación en cada lago se

10
reduce a la mitad suponiendo que la contaminación ha cesado, es decir, las afluencias que no
provienen de ninguno de los lagos son de agua limpia. Suponiendo que todos los lagos tienen la
misma concentración original de contaminación p, ¿cuánto tardará en reducirse la contaminación
al 5 por ciento del nivel original?.

Problemas asociados a EDP

34. La conducción de calor a lo largo de una tuberı́a completamente aislada de longitud 1

metro, por la que circula agua a una velocidad de 1 m/s y cuyos extremos izquierdo y derecho se
mantiene a temperatura constante de 00 C y 10 C respectivamente, está modelada por el problema
∂2u ∂u
∂u + =0 0 1 0
−k 2 , < x< , t > ,
∂t ∂x ∂x
u(0, t) = 0, t ≥ 0, (1)
u(1, t) = 1, t ≥ 0,
u(x, 0) = f (x), 0 ≤ x ≤ 1,

donde u(x, t) es la temperatura en el instante t, x es la posición, k es la constante de transmi-

sividad calorı́fica del agua y f (x) es la distribución de temperaturas en t = 0.
i) Calcula φ(x) tal que v(x, t) = u(x, t) − φ(x) sea solución de
∂2v ∂v
∂v + =0 0 1 0
−k 2 , < x< , t > ,
∂t ∂x ∂x
v(0, t) = 0, t ≥ 0, (2)
v(1, t) = 0, t ≥ 0,
v(x, 0) = g(x), 0 ≤ x ≤ 1,

donde g(x) = f (x) − φ(x) y u(x, t) es la solución de (1)

ii) Resuelve el problema (2) cuando k = 1, g(x) = (1 − x)ex/2.
35. Se desea estudiar el movimiento vibratorio de una cuerda de longitud 1 fijada por sus
extremos, amortiguado por la acción del rozamiento del aire. Llamando u(x, t) al desplazamiento
de la cuerda de cada punto x en el instante t respecto de su posición de equilibrio, el fenómeno
está modelado por
∂u ∂ 2u
∂2u = 2 =0 0 1 0
+2
c k 2 , < x< , t > ,
∂t2 ∂t ∂x
u(0, t) = u(1, t) = 0, t ≥ 0,
(3)
u(x, 0) = f (x), 0 ≤ x ≤ 1,
∂u
(x, 0) = 0, 0 ≤ x ≤ 1,
∂t
≤ ≤

donde c es una constante que depende del rozamiento y k es una constante que depende de la
elasticidad y tensión de la cuerda.
i) Resuelve (3) cuando c < k y f (x) = sen πx + sen 3πx cos πx.
ii) Si c = 1/2, k = 1 calcula T tal que
1
"u(x, t)"L2(0,1) ≤ "f (x)"L2(0,1), ∀ t ≥ T.
2
36. La transmisión de impulsos eléctricos en un cable largo con una distribución de capacitancia,
inductancia y resistencia, viene modelada por la ecuación de los telegraftstas:

11
∂ 2u ∂u 2
2∂ u
+2
a + bu − c = 0, a, b, c ∈ R, c > 0, 0 < x < π, t > 0.
∂t2 ∂t ∂x2

12
Calcula la solución formal del problema anterior cuando las condiciones iniciales y de contorno
son
u(0, t) = u(π, t) = 0, t ≥ 0,
u(x, 0) = 0, 0 ≤ x ≤ π,
∂u
(x, 0) = g(x), 0 ≤ x ≤ π.
∂t

13