INTRODUCCIÓN AL

CÁLCULO
DIFERENCIAL E
INTEGRAL

CURSADA 2018
 PROGRAMA INTRODUCCIÓN AL CÁLCULO DIFERENCIAL E INTEGRAL

Unidad I:
Naturaleza de las ecuaciones diferenciales, definición, grado, diferencia entre
ecuaciones diferenciales ordinarias y ecuaciones en derivadas parciales.
Interpretación gráfica: Curvas Integrales, Teorema de Picard, monotonía, concavidad,
simetría, singularidad e isoclinas.

Unidad II:
Ecuaciones diferenciales lineales de primer orden, factor integrante, Ecuaciones a
variables separables Ecuaciones homogéneas. Ecuaciones exactas. Aplicaciones:
Trayectoria Ortogonal, crecimiento, desintegración, mezclas etc.

Unidad III
Métodos de aproximación numérica de Euler, Euler mejorado y Runge Kuta.

Unidad IV:
Espacio solución: dimensión del espacio solución, Wronskiano. Ecuaciones
diferenciales lineales de segundo orden. Teorema de Abel. Reducción de Orden.
Resolución de la ecuación homogénea a coeficientes constantes. Método de
coeficientes indeterminados. Método de variación de parámetros. Aplicaciones.

Unidad V:
Sistemas de ecuaciones de primer orden: generalidades sobre sistemas, sistemas
lineales, sistemas lineales homogéneos a coeficientes constantes. Métodos por
operadores y por autovectores y autovalores. Sistemas lineales inhomogeneos.
Aplicaciones

Unidad VI:
Sistemas autónomos lineales. Tipos de puntos críticos. Estabilidad. Trayectoria.
Sistemas autónomos no lineales. Aplicación: Sistema presa predador de Volterra.

Bibliografía:

 Braun, M ; Differential Equations and their Applications. Springer, New York 1991.
 Lomen, D & Lovelock, D; Exploring Differential Equations via Graphics and Data.
John Wiley, New York, 1996.
 Boyce, W, Di Prima, R; Introducción a las Ecuaciones Diferenciales, Jhon Wiley,
New York, 1970.
 Simmons, G; Robertson, J Ecuaciones Diferenciales con aplicaciones, McGraw-
Hill, New York, 1993.
 Práctico No 1- Introducción a las ecuaciones diferenciales ordinarias

1) Clasificar las siguientes ecuaciones diferenciales según sean ordinarias o parciales. Determinar orden
y grado de las mismas, como así también la función desconocida y la variable independiente.

u  2 u u
b) y ´´ 4 y ´ 5 y  e 4 2 
3x
a) y´ = x 2 + 5 y c)
t x y

2 3 dr
 d 3u   d 2 s  = rq
d)  3    2   s  3t e) d q
 dt   dt 
2) Verificar que 𝑥 3 + 3𝑥𝑦 2 = 1 es una solución implícita de la ecuación diferencial,
𝑑𝑦
2𝑥𝑦 + 𝑥2 + 𝑦2 = 0 en el intervalo 0 < 𝑥 < 1.
𝑑𝑥

3) En cada uno de los siguientes ítems, se describe una función y = g(x) por medio de alguna propiedad
de su gráfico. Escribir una ecuación de la forma y’ = f(x, y) que tenga a g(x) como solución (o como una de
sus soluciones).
a. La pendiente del gráfico de g en el punto (x, y) es la suma de x e y.
b. La tangente al gráfico de g en el punto (x, y) interseca al eje x en el punto (x/2, 0).
c. Toda línea recta perpendicular al gráfico de g pasa por el punto (0, 1).
4) Escribir una ecuación diferencial que sea un modelo matemático de la situación descripta.
a. La aceleración de un auto deportivo es proporcional a la diferencia entre 250 km/h y la velocidad
del auto.
b. En una población fija de P habitantes, la razón de cambio del número N de personas que han
escuchado cierto rumor es proporcional al número de personas que aún no lo han escuchado.
c. En una población fija de P habitantes, la razón de cambio del número N de personas infectadas
con una enfermedad es proporcional al número de personas infectadas y al número de personas
no infectadas.
5) Los cambios en el número de individuos de una población dependen de un balance entre el número de
nacimientos y muertes. Si b es una constante que representa el número de nacimientos por individuo y por
unidad de tiempo y m es lo mismo, pero en términos de mortalidad, entonces 𝑑𝑃/𝑑𝑡 = 𝑏𝑃 − 𝑚𝑃 . Luego
dP/dt = (b-m) P , donde r = b-m es la tasa intrínseca de crecimiento. Este es el modelo de crecimiento
exponencial de poblaciones. a) Analizar qué ocurre con la población cuando 𝑏 < 𝑚, 𝑏 = 𝑚 y 𝑏 > 𝑚. b)
Bosquejar el comportamiento de P(t) a los largo del tiempo, asumiendo 𝑃(0) = 𝑃0 , en los tres casos. c)
Encontrar una expresión para P(t). d) Graficar luego la solución para tres valores diferentes de la constante C.
6) Durante décadas, las ballenas han sido cazadas, llegando a valores mínimos extremos de 5000
individuos en el año 1978. En condiciones óptimas, la tasa intrínseca de crecimiento fue de 0.047 año-1.
Asumiendo un crecimiento exponencial. a) Plantear la ecuación diferencial. b) Graficar el campo de
direcciones e isoclinas. c) ¿Cuántos años tardaría la población en alcanzar 50.000 individuos? ¿Y 150.000?
7) Los procesos de decaimiento de un material radiactivo o de eliminación de drogas en sangre pueden
modelarse siguiendo el mismo proceso que el descripto en el crecimiento exponencial de poblaciones, pero
considerando que no existe el término de natalidad, sino solo el término de mortandad negativo denominado,
en estos casos, tasa específica de desintegración o de decaimiento, que indica una destrucción o de
eliminación de materia en cada caso.
 a. Plantear la ecuación diferencial para estos casos indicando las unidades de la variable
dependiente y del parámetro r en la ecuación.
b. Graficar el campo de direcciones y algunas soluciones para condiciones iniciales dadas.
8) La vida media VM, de un isótopo radiactivo se define como el tiempo necesario para que la mitad de
su masa se desintegre.
a. Hallar VM y verificar que es independiente del tiempo inicial.
b. Calcular la VM del C14, sabiendo que r = 0.000120.
9) Todos los organismos vivos poseen C14 en pequeñas cantidades, dado que toman la parte del mismo
que se encuentra en la atmósfera (6x10 10 átomos de C14 por gramo de C12) a través del CO2 utilizado en el
proceso de fotosíntesis de los productores primarios. A partir de allí, continúa el proceso a través de la cadena
trófica llegando a consumidores. Cuando el organismo muere, deja de incorporar C 14 y comienza el proceso de
decaimiento del mismo. Si consideramos que A0 era la cantidad presente en el organismo al momento de su
muerte y A1 la cantidad medida en el momento en el cual se quiere realizar la determinación, la solución de la
ecuación que modela el proceso será: A1  A0e kt . Demuestre que la expresión de la edad, despejando el
A0
tiempo de la solución dada y considerando la vida media y la tasa de decaimiento del C 14, es t  8.310 ln .
A1

10) La variación de la temperatura de un objeto, o tasa de cambio de la temperatura, es directamente

proporcional a la diferencia entre la temperatura del objeto en cada tiempo y la del ambiente exterior al cual se
lo somete. Ésta es la ley de enfriamiento de Newton. a) Exprese la ecuación diferencial que modela este
fenómeno. b) Encuentre la solución general. c) Aplique el modelo para resolver el siguiente problema: Al
sacar una torta del horno, su temperatura es 150°C y después de 3 minutos de 90°C. Si la temperatura
ambiente es de 20°C, a) Plantee el modelo y grafique el campo de pendientes e isoclinas, b) ¿Cuál será la
temperatura de la torta a los 10 minutos? c) ¿Cuándo la temperatura de la torta alcanza los 20°C?
11) Resolver las ecuaciones de variables separables:
a. 4𝑥𝑦 𝑑𝑥 + (𝑥 2 + 1)𝑑𝑦 = 0
b. (𝑥𝑦 + 2𝑥 + 𝑦 + 2)𝑑𝑥 + (𝑥 2 + 2𝑥)𝑑𝑦 = 0
c. 2𝑟(𝑠2 + 1)𝑑𝑟 + (𝑟 4 + 1)𝑑𝑠 = 0
d. 𝑐𝑜𝑠𝑒𝑐 𝑦 𝑑𝑥 + sec 𝑥 𝑑𝑦 = 0
e. (𝑒 𝑥 + 1) cos(𝑢) 𝑑𝑢 + 𝑒 𝑥 (𝑠𝑒𝑛(𝑢) + 1) 𝑑𝑦 = 0
f. 𝑡𝑔 𝜃 𝑑𝑟 + 2𝑟 𝑑𝜃 = 0
12) Resolver los problemas de valor inicial
a. (y + 2)dx + y(x + 4)dy = 0; y(−3) = −1
b. 8 cos 2 y dx + cosec 2 x dy = 0; y(π⁄12) = π⁄4
13) Usar los conceptos de monotonía, concavidad, simetría, singularidad, isoclinas y unicidad para dibujar
las soluciones de las siguientes ecuaciones diferenciales:
dy dy dy
a) = x 3 , P ( 1 ,1 ) b ) = x 4 , P ( 1 ,1 ) c ) = cos x, P ( 0 , 0 )
dx dx dx
dx dy dy 1
d ) = sent, P ( p , 2 ) e ) = e - t , P ( 0 ,1 ) f ) = 2 , P ( 1 ,1 )
dt dt dx x
dy 1 dx 1
g) = , P ( 2 ,1 ) h ) = , P ( 1 ,p / 4 )
dx x( x - 1 ) dt 1 + t 2
dy dy
i ) = ln t, P ( 1 ,1 ) j ) = x 2e - x , P ( 0 ,1 )
dt dx
 Práctico No 2: Ecuaciones diferenciales lineales de orden 1

1) Resolver las ecuaciones diferenciales lineales:

a. 𝑦 ′ − 3𝑥 2 𝑦 = 𝑥 2 b. 𝑦 ′ + 𝑦 𝑡𝑔(𝑥) = 𝑠𝑒𝑐(𝑥) + 2𝑥 𝑐𝑜𝑠(𝑥) c. 𝑦 ′ − 3𝑦 = 2
d. 𝑥𝑦 ′ + (1 + 𝑥)𝑦 = 5 e. 𝑥 2 𝑑𝑦 + (2𝑥𝑦 − 𝑒 𝑥 )𝑑𝑥 = 0 f. 𝑦 ′ + 𝑦 𝑡𝑔(𝑥) = 𝑠𝑒𝑛(𝑥)
g. 𝑦 ′ + 2𝑦 = 𝑒 2𝑥 h. 𝑦 ′ + 𝑦 𝑐𝑜𝑡𝑔(𝑥) = 4𝑥 2 𝑐𝑜𝑠𝑒𝑐(𝑥) i. 𝑥𝑦 ′ + (2 + 3𝑥)𝑦 = 𝑥𝑒 −3𝑥
j. 𝑥𝑦 ′ + 𝑦 + 𝑥 = 𝑒 𝑥 k. 𝑥 2 𝑑𝑦 + (𝑥 − 3𝑥𝑦 + 1)𝑑𝑥 = 0 l. 𝑥 −1 𝑦 ′ + 2𝑦 = 3
2) Encontrar la solución particular de la ecuación que satisface la condición dada:
a. 𝑥𝑦 ′ − 𝑦 = 𝑥 2 + 𝑥; 𝑦(1) = 2
′ −3𝑥
b. 𝑦 + 2𝑦 = 𝑒 ; 𝑦(0) = 2
′ −𝑥
c. 𝑥𝑦 + 𝑦 + 𝑥𝑦 = 𝑒 ; 𝑦(1) = 0
3) Una ciudad tenía una población de 25000 habitantes en 1970 y una población de 30000 habitantes en
1980. Suponiendo que la población continuará creciendo exponencialmente a una razón constante, ¿Qué
población se puede esperar para el año 2020?
4) Las poblaciones descriptas en los ejercicios anteriores crecen indefinidamente a tasa constante,
dP/dt = rP, r > 0. En la realidad, el ambiente le impone un límite al crecimiento de la población
disminuyendo su tasa de crecimiento. Esta situación se puede reflejar proponiendo una función lineal, f(P),
dP
decreciente tal que = f(P)P como se observa en la figura siguiente, donde K representa la capacidad
dt
máxima de carga del ambiente.
a) Plantear la ecuación diferencial para esta función. (Ecuación
logística, denso dependencia)
b) Encontrar la expresión de P(t).
c) Determinar el comportamiento de la solución para valores
iniciales de población 𝑃0 < 𝐾, 𝑃0 = 𝐾 y 𝑃0 > 𝐾 y bosquejar las
soluciones que ilustren las condiciones iniciales anteriores.

5) En 1934, Gause realizó experimentos de competición involucrando dos especies de protozoos

(Paramecium caudatum y Paramecium aurelia). Cuando cada población se cultivó de manera separada, las
105 64
expresiones logísticas que ajustó fueron: N1 (t )  1.1244t
y N1 (t )  . Determine, para
1  34 e 1  15 e  0.794 t
ambas poblaciones: a) población inicial, b) capacidad de carga, c) ecuación diferencial correspondiente d)
tiempo donde el crecimiento es mayor
6) Dadas dos curvas C1 y C2 que se cortan en un punto ( x, y) diremos que se cortan ortogonalmente si sus
rectas tangentes en dicho punto son ortogonales. Dada una familia de curvas que depende de un parámetro
h( x, y, c)  0 diremos que una curva C es una trayectoria ortogonal a dicha familia si en cada punto en que C
corta a una curva de la familia lo hace ortogonalmente, es decir, las dos rectas tangentes en el punto de corte
son perpendiculares. Asimismo, dadas dos familias de curvas Y 1 y Y2 , diremos que la familia de curvas Y2
son las trayectorias ortogonales a la familia Y 1 si cada curva de la familia Y2 corta ortogonalmente a la familia
Y1 y toda curva ortogonal a esta última familia pertenece a Y 2.
Dada una familia de curvas uniparamétrica h( x, y, c)  0 nos planteamos cómo obtener la familia de sus
trayectorias ortogonales. Supongamos que conocemos la ecuación diferencial y y  f ( x, y, c) de la familia
h( x, y, c)  0 . Fijamos un punto ( x, y) de una de las curvas de la familia h( x, y, c)  0 . La pendiente de su
recta tangente es f ( x, y) . Recordemos que dos rectas de pendientes m1 y m2 , respectivamente, son
ortogonales si, y sólo si m1.m2  1 . Por tanto, la pendiente de cualquier curva ortogonal a ella en el punto
1 1
( x, y) es  . Y esta nueva curva debe ser solución de la ecuación diferencial y   . En
f ( x, y ) f ( x, y )
resumen, si disponemos de la ecuación diferencial de una familia de curvas, la ecuación diferencial de la
1
familia de curvas ortogonales se obtiene sustituyendo en la ecuación anterior y por  .
y
a. Dada la familia de curvas integrales 𝑥 2 + 𝑦 2 = 𝑐. Determinar la familia de trayectorias ortogonales.
b. Calcular el valor de K para que las parábolas 𝑦 = 𝑐1 𝑥 2 + 𝐾; sean trayectorias ortogonales de la familia
de elipses 𝑥 2 + 2𝑦 2 − 𝑦 = 𝑐2 .
2
c. Determinar las trayectorias ortogonales de la familia de parábolas y = cx .

7) Dado un compartimento como el de la figura con un volumen V fijo [capacidad] conteniendo una
sustancia que se renueva a través de una entrada y salida que permite un flujo de caudal expresado por F(t)
[capacidad/tiempo]. Una nueva sustancia, en general denominada traza, y(t) [masa], es agregada y mezclada
dentro del tanque a una tasa I(t) [masa/tiempo]. Plantee la ecuación diferencial que modela la variación de la
cantidad de trazador dentro del compartimento si ésta es la diferencia de sustancia entrante y saliente del
volumen considerado en cada instante de tiempo

I(t)

F (t)
8) En un tanque de 500 litros hay un flujo de agua de 10 litros/segundo.
a. Si se arrojan al tanque 20 kg de sal, 1) Plantee el modelo y grafique el campo direcciones. 2)
¿Cuánto tiempo tomará que la cantidad remanente de sal sea de 5 kg?
b. Si inicialmente el tanque no tenía sal y en un determinado momento se comienzan a incorporar al
tanque 3 kg de sal por segundo, 1) Plantee el modelo y grafique el campo direcciones. 2) ¿Cuánto
tiempo debo esperar para que la cantidad de sal sea 50 kg?
c. ¿Y si la cantidad de sal que ingresa al tanque por segundo está dada por 𝑓(𝑡) = 2 𝑡?
d. Analizar cómo varió el modelo en función de los enunciados.
9) La tasa de crecimiento de un pez es directamente proporcional a L  L , donde L es la longitud
máxima que puede alcanzar la especie, y L es la longitud del pez en el instante que es estudiado. a) Encontrar
la expresión que exprese la longitud del pez en cualquier instante de tiempo. b) Sabiendo que en los peces
existe una relación entre peso y longitud dada por W   Lv . Escribir la ecuación que describa la variación de
peso en función del tiempo y la expresión del peso del pez en cualquier instante de tiempo dado.
10) Se sabe que una pelota se infla a razón de 20 pies3/s. ¿A qué velocidad se incrementa el radio de la
pela cuando su diámetro es de 1pie? (Recuerde que 𝑉 = (4/3) 𝜋 𝑟 3 ).
11) Por un agujero circular de área A, en el fondo de un tanque, sale agua.
El agua del tanque se reduce por segundo a cA(2gh)1/2, donde 0<c<1. Deduzca
una ecuación diferencial que exprese la altura h del agua en cualquier
momento t, que hay en el tanque de la figura. El radio del agujero es de 2
pulgadas y g=32 pies/s2

dv
12) La aceleración representa la variación de la velocidad y se expresa como a = . La segunda ley de
dt
Newton expresa que podemos relacionar la masa de un cuerpo con la acelaración que alcanza al aplicarle
dv
fuerzas exteriores mediante la ecuación F = m a = m . Un modelo simple de caída libre considerando la
dt
dv
resistencia del aire está dado por la expresión m = mg − kv, considerando que la resistencia del aire es
dt
proporcional a la velocidad de caída y en sentido contrario a la misma y que el peso del cuerpo se puede
expresar como el producto entre su masa y la aceleración de la gravedad.
a. Encontrar la solución para v(t) de la ecuación planteada considerando v(0) = 0.
b. Determinar la expresión para la velocidad límite, considerando que la misma se define como el límite
de la solución para t → ∞
c. Aplicar el modelo para resolver el siguiente problema: Se arroja un cuerpo de 6 kg desde la terraza de
un edificio. Asumiendo que la magnitud de la resistencia del aire en cada instante es igual al doble de la
magnitud de la velocidad, hallar la velocidad y la distancia recorrida al cabo de t segundos.

Más problemas!!!
13) Un cultivo de bacterias duplica su tamaño cada 6hs. ¿Cuánto tiempo tarda en alcanzar 5 veces su
tamaño original?
14) 5 gr de una sustancia radiactiva decaen a 4.1 gr en 4 días. Encontrar la vida media y la tasa específica
de decaimiento.
15) Se calienta agua para mate, pero el cebador se distrae y el agua alcanza los 110ºC. La pava se deja
entonces a temperatura ambiente, que en ese día otoñal ronda los 10ºC. Luego de media hora, el agua está a
60º ¿En cuánto tiempo alcanzará los 30ºC?
16) El flujo de agua a través de un tanque de 500 litros es de 10 litros/seg. Si se arrojan al tanque 20 Kg de
sal, ¿cuánto tiempo tomará que la cantidad remanente de sal sea de 5 Kg?
17) La media vida del cobalto radiactivo es 5,27 años. Supongamos que un accidente nuclear ha dejado el
nivel de radiación de cobalto en una cierta región 100 veces más alto que el nivel aceptable para la vida
humana. ¿Cuántos años deberán pasar hasta que la región se vuelva habitable?
18) Después de que te aplican un antibiótico inyectable, la concentración de la droga en el cuerpo
disminuye al 50% en 10hs. ¿Cuánto tiempo pasa hasta que la droga alcanza el 10% de la concentración
original?
19) Si la vida media de una sustancia es 1600 años ¿qué proporción quedará al cabo de 2400 años? ¿Y de
8000 años?
20) Se arroja un cuerpo de 6 kg desde la terraza de un edificio. Asumiendo que la magnitud de la
resistencia del aire en cada instante es igual al doble de la magnitud de la velocidad, hallar la velocidad y la
distancia recorrida al cabo de t segundos.
21) Por razones obvias, la sala de disección de un forense se mantiene a una temperatura constante de 5ºC.
Ocurre un crimen. Mientras realiza la autopsia de la víctima, el forense es asesinado, y el cuerpo del delito
robado. A las 7 hs, el ayudante del forense descubre su cadáver a 23ºC; al mediodía, su temperatura es de
16ºC. Si la temperatura normal de un cuerpo puede asumirse en 37ºC ¿a qué hora murió el forense,
aproximadamente?
22) Un modelo de crecimiento de población estacional propone que la tasa de crecimiento puede ser
1 dN
expresada en la forma;  k sen 2t , donde t se mide en años y t  0 corresponde al comienzo de la
N dt
primavera. De esta forma, la población crecerá durante la primavera y el verano, y decrecerá en otoño e
invierno.

a. Plantear la ecuación diferencial correspondiente dado que N (1/ 2)  N 0 .

b. De acuerdo al modelo planteado, ¿cuándo se tiene el máximo nivel de población?

23) Hallar la trayectoria ortogonal a 𝑦 = 𝑐/𝑥.

24) Un depósito de 100 litros de capacidad contiene 20 litros de agua en la que hay disueltos 10 gramos de
sal. A partir de cierto instante se vierten en el depósito 4 litros por minuto de agua que contienen 2 gramos de
sal por litro, mientras que se deja salir la disolución bien mezclada a un ritmo de 2 litro por minuto. (a)
Encontrar la ecuación que modele la cantidad de sal en el depósito en cada instante. (b) ¿Cuál es la
concentración de sal en el momento en que el depósito se llena?
25) Un reactor convierte el uranio 238 relativamente estable en el isótopo plutonio 239. Después de 15
años se determina que se desintegró 0.043% de la cantidad inicial A 0 de plutonio. Calcule la vida media de
este isótopo si la rapidez de desintegración es proporcional a la cantidad presente.
26) Un vaso de agua a 25 oC se introduce en un congelador a -20 oC. En 15’ el agua está ya a 20 oC.
¿Cuánto tiempo tarda en llegar el agua a cero grado?
27) Un paracaidista de masa m se lanza de un avión y abre el paracaídas a una altura ho del suelo, cuando
ha adquirido en caída libre una velocidad vo. Si al lanzarse, el paracaidista tenía velocidad nula, ¿a qué altura
estaba el avión en el momento del lanzamiento? (considerar la fuerza de rozamiento del aire proporcional a la
velocidad)
28) El 2% de los alumnos de los alumnos que cursan matemática difunden a las 12 hs. del lunes el rumor
de que el martes a las a las 11 hs. tendrán un examen sorpresa. Sabiendo que el lunes a las 14hs la noticia es
conocida por el 5% de los alumnos ¿cuál es el porcentaje de alumnos que a la hora de comienzo de la clase del
martes conocen la noticia? (el rumor se transmite proporcional al número de encuentros entre alumnos que
conocen la noticia y aquellos que no la conocen).
29) Se encuentra un hueso fosilizado que contiene una milésima de la concentración de C-14 que se
encontraba en la materia viva. Si se sabe que la vida media del C-14 es de 5700 años, estime la edad del fósil.
30) Una colonia de bacterias crece de acuerdo con la ley logística, con una capacidad de carga de 5.108
individuos y una tasa de crecimiento natural r = 0.01 día-1. a) ¿Cuál será la población después de 3 días si
inicialmente era de 108 individuos? b) Si la población inicial fue de 8.108, la población crecerá o decrecerá?
Justifique.
31) Suponga que una gota de agua se evapora a una velocidad proporcional a su área superficial. Si
originalmente el radio es de 3 mm y una hora después se ha reducido a 2 mm, encontrar una expresión para el
radio de la gota como función del tiempo ¿Cuánto tardará en desaparecer?
32) Un Hombre provisto de un paracaídas se lanza desde una gran altura. El peso conjunto entre el
hombre y el paracaídas es de 98 Kg. Sea v(t)(la velocidad 8 en metros por segundos)t segundos después del
lanzamiento. Durante los 10 primeros segundos, antes de abrirse el paracaídas se supone que la resistencia del
aire es de 0,1v(t) kg. Después una vez abierto el paracaídas, la resistencia del aire es 2v(t). Suponiendo que la
aceleración de la gravedad es 9,8m/s2. Hallar una expresión adecuada para la velocidad v(t) y el espacio
recorrido.
33) El radio 226 tiene una vida media de 1620 años. Encuentre el período de tiempo en que un cuerpo de
este material se reduce a tres cuartas partes de su tamaño original.
34) Un contenedor es arrastrado por el hielo en un trineo, incluido el trineo, el peso total es de 80 libras.
Suponiendo que es despreciable la resistencia del hielo y que el aire opone una resistencia en libras igual a 5
veces la velocidad (en pies/seg) del trineo. (Sabiendo que la gravedad es de 32pies /seg 2), encontrar: a) la
fuerza constante ejercida por el trineo para obtener una velocidad terminal de 10 millas por hora. b) La
velocidad y distancia recorrida al cabo de 48 segundos
 Práctico No 3: Métodos numéricos
para la resolución de ecuaciones diferenciales ordinarias.

1) Utilice los tres métodos para estimar y( 1 ) si y' = - y; y( 0 ) = 1 considerando tres pasos diferentes.

2) Dado el PVI xdy + 2 ydx = 0; y( 2 ) = 1 ,

a. Aproximar numéricamente la solución por el método de Runge Kutta en X=2.4

b. Comparar el error con dos pasos h=0.2 y h=0.1.
3) Estimar y( 0.2 ) de la ecuación y' - 3x - 3 y = 0, y( 0 ) = 1 utilizando el método de Runge Kutta
tomando h=0.1. Luego comparar con la solución exacta.
4) Estimar y( 0.02 ) de la ecuación y' = 2 x + 2 y, y( 0 ) = 1 utilizando el método de Euler mejorado y
Runge Kutta comparando con la solución exacta.
5) Para la población de ballenas azules del ejercicio 7 del práctico 1, estimar el número de ballenas que
hubo en 1998 mediante una aproximación numérica.
6) Sobre la base del problema 10.a del práctico 2,
a. estimar la cantidad sal en el tanque luego de 50 seg, utilizando Euler mejorado y dos pasos
diferentes.
b. calcular el error cometido al trabajar con cada paso.
7) Estimar 𝑦(−2.5) de la ecuación (𝑦 + 2)𝑑𝑥 + 𝑦(𝑥 + 4)𝑑𝑦 = 0; 𝑦(−3) = −1 utilizando el método
de Euler mejorado y Runge Kutta comparando con la solución exacta que fue obtenida en 3.5.a).
 Práctico No 4: Ecuaciones diferenciales homogeneas y exactas

1) Resolver las ecuaciones homogéneas y problemas de valores iniciales

a. (𝑥 + 𝑦)𝑑𝑥 − 𝑥 𝑑𝑦 = 0
b. (2𝑥𝑦 + 3𝑦 2 )𝑑𝑥 − (2𝑥𝑦 + 𝑥 2 )𝑑𝑦 = 0
c. 𝑣 3 𝑑𝑢 + (𝑢3 + 𝑢𝑣 2 )𝑑𝑣 = 0; 𝑣(3) = 9
d. (𝑥 2 + 3𝑦 2 )𝑑𝑥 − 2𝑥𝑦 𝑑𝑦 = 0; 𝑦(2) = 6
2) Resolver las siguientes ecuaciones reducibles a homogéneas:
a. (𝑥 + 𝑦 − 2)𝑑𝑥 + (𝑥 − 𝑦 + 4)𝑑𝑦 = 0
b. (𝑥 + 𝑦 + 1)𝑑𝑥 + (2𝑥 + 2𝑦 − 1)𝑑𝑦 = 0
c. (2𝑥 − 4𝑦)𝑑𝑥 + (𝑥 + 𝑦 − 3)𝑑𝑦 = 0
3) Verificar si las siguientes ecuaciones diferenciales son exactas y en ese caso, resolverlas.
a) 𝑦 2 𝑑𝑥 + 2𝑥𝑦 𝑑𝑦 = 0
b) (3 𝑥 2 + 4𝑥𝑦)𝑑𝑥 + (2 𝑥 2 + 2𝑦)𝑑𝑦 = 0
c) (2𝑥 𝑐𝑜𝑠 𝑦 + 3𝑥 2 𝑦)𝑑𝑥 + (𝑥 3 − 𝑥 2 𝑠𝑒𝑛 𝑦 − 𝑦)𝑑𝑦 = 0
d) (𝑒 2𝑦 − 𝑦 cos(𝑥𝑦))𝑑𝑥 + (2𝑥𝑒 2𝑦 − 𝑥 cos(𝑥𝑦) + 2𝑦) 𝑑𝑦 = 0
e) (cos 𝑥 𝑠𝑒𝑛𝑥 − 𝑥𝑦 2 )𝑑𝑥 = 𝑦(𝑥 2 − 1)𝑑𝑦; 𝑦(0) = 2
f) (𝑥 + 2𝑦 2 )𝑑𝑥 + 𝑥𝑦 𝑑𝑦 = 0
4) En cada uno de los siguientes problemas determinar la constante a para que la e.d. sea exacta, y
entonces resolverla:
𝑑𝑦
a) 𝑡 + 𝑦 𝑒 2𝑡𝑦 + 𝑎𝑡 𝑒 2𝑡𝑦 =0
𝑑𝑡
1 1 (𝑎𝑡+1) 𝑑𝑦
b) + 𝑦2 + =0
𝑡2 𝑦3 𝑑𝑡
dy
c) eat+y + 3t 2 y 2 +(2yt 3 + eat+y ) =0
dt
5) Para cada una de las ecuaciones siguientes encontrar si es posible un factor de integración que
dependa de una sola variable. En caso de obtenerlo resolver la ecuación.
a) (𝑥 + 𝑦 2 )𝑑𝑥 + 𝑥𝑦 𝑑𝑦 = 0 d) (5𝑦 − 𝑥 2 )𝑑𝑥 + 𝑥𝑑𝑦 = 0
b) 2𝑦𝑑𝑥 + (𝑥 + 𝑦) 𝑑𝑦 = 0 e) (𝑥 + 𝑦 − 𝑥𝑦)𝑑𝑥 + 𝑥𝑑𝑦 = 0
2
c) (𝑥𝑦 − 1)𝑑𝑥 + (𝑥 − 𝑥𝑦)𝑑𝑦 = 0 f) (𝑢 + 5𝑡)𝑑𝑢 + 3𝑢 𝑑𝑡 = 0
6) Mediante un estudio se observó que la variación de la cantidad de pastura por hectárea en función de
la carga ganadera está dada por la siguiente ecuación: dy⁄dx = (2y − x 2 )/(−2x + y). Encontrar una
expresión para la función que permite conocer la cantidad de pastura por hectárea en función de la carga
ganadera del lugar.
7) Cuando dos especies coexisten en un mismo ambiente, la variación en la abundancia de la especie 2
(y) es en función de la abundancia de la especie 1 (x) y viceversa. Consideremos que la abundancia de la
especie 2 varía según la ecuación:
a. dy/dx = (−x + 2y)/(2x − y).
b. 𝑑𝑦/𝑑𝑥 = (−2𝑥 + 4𝑦 − 6)/(𝑥 + 𝑦 − 3).
c. 𝑑𝑦/𝑑𝑥 = (𝑥 − 𝑦 − 1)/(𝑥 − 𝑦 − 2).
En cada caso encontrar una función que describa la abundancia de la especie 2 en función de la abundancia de
la especie 1.
8) Determinar las trayectorias ortogonales de la familia de:
a. circunferencias dadas por: x 2 + y 2 = 2cx,
b. todas las circunferencias que pasan por el origen y tienen centro en el eje x.
 Práctico N° 5: Ecuaciones Lineales de 2° orden

1- Wronskiano – conjunto fundamental de soluciones

1.
a Verificar que y 1 = t y y2 = 1 / t son soluciones de la ecuación diferencial 2t 2 y'' + 3ty' – y = 0
sobre el intervalo 0 < t < ∞ .
b Evaluar W [ y1 , y2 ](t ) . ¿Qué ocurre cuando t → 0 ?
c Mostrar que y1 (t ) e y2 (t ) forman un conjunto fundamental de soluciones de la ecuación dada en el
intervalo 0 < t < ∞ .
d Resolver el problema del valor inicial: 2t 2 y'' + 3ty' – y = 0 , y (1) = 2 , y' (1) = 1

d2y dy
2. Considerar la ecuación diferencial 2
− 4 + 3y = 0 .
dx dx

a Demostrar que cada una de las funciones ex, e3x es una solución de la ecuación diferencial en el intervalo
a<x<b donde a y b son reales positivos arbitrarios tales que a<b.
b ¿Por que podemos concluir que cada una de las funciones 5ex +2 e3x, 6ex-4e3x; -7ex+5 e3x es también una
solución de la ecuación diferencial dada en a<x<b ?
c ¿Podemos afirmar lo mismo de 3ex; -4ex; 5ex; 6ex en a<x<b ?

d2y dy
3. Considerar la ecuación diferencial 2 − 5 + 6y = 0
dx dx
a Demostrar que e2x y e3x son soluciones L.I de esta ecuación en el intervalo - ∞ < t < ∞ .
b Escribir la solución general de la ecuación dada.
c Determinar la solución que satisface las condiciones y (0) = 2 , y' (0) = 3 . Explicar por qué esta solución
es única. ¿En que intervalo está definida?

2- Reducción de Orden
4. Encuentre una segunda solución para cada una de las ecuaciones diferenciales, utilizando el teorema de
Abel:
a) x 2 y'' + 2 xy' − 2 y = 0 , y1 ( x ) = x b) y'' + 2y' + yx = 0, y1 ( x) = e− x c) xy'' + y' = 0 , y1 ( x ) = 1

d) y'' + y = 0 , y1 ( x) = senx e) y'' − y = 0 , y1 ( x) = e x

5. Mostrar que y=x es una solución de la ecuación de Legendre de orden uno:
2
( 1− x )y'' − 2 xy' + 2 y = 0; −1 < x < 1 y encontrar una segunda solución linealmente independiente
(Por el método de reducción de orden). Analizar los problemas de validez de las soluciones.
6. Encuentre una segunda solución para cada una de las ecuaciones diferenciales especificadas:
a) x 2 y'' + 2 xy' − 2 y = 0 , y1 ( x ) = x b) xy'' + y' = 0 , y1 ( x ) = 1
1
c) x 2 y'' + 2 xy' = 0 , y1 ( x) = 1, x > 0 y = 0, y1 ( x ) = x, x > 0
d) xy'' − y' +
x
7. Mostrar que y1( x ) = x−1/ 2 senx es solución de la ecuación de Bessel de orden α = 1 / 2 :
x 2 y'' + xy' + ( x 2 − α 2 )y = 0 , y encontrar una segunda solución linealmente independiente de y1.
3- Ecuaciones Diferenciales de Orden dos no lineales (ausencia de variables).
8. Ausencia de variable independiente o Ausencia de variable dependiente:
2 3
a ) yy’’ + ( y’ ) = 0 b )y’’ + y’ = 0 c ) yy’’ − ( y’ ) = 0 d )xy’’ + y’ = 1, x > 0
2 2 2 3
e )2y 2 y’’ + 2y ( y’ ) = 1 f )y’’ + ( y’ ) = 2e-y g )y’’ + x ( y’ ) = 0 h )y’’ + y ( y’ ) = 0

4- Ecuaciones Diferenciales de Orden dos lineales a coeficientes constantes

9. Resolver los siguientes problemas a valores iniciales de orden dos con coeficientes constantes:

a) y'' − 5 y' + 6 y = 0 , y (1) = e 2 e y' (1) = 3e 2 b) y'' − 6 y' + 5 y = 0 , y (0) = 3 e y' (0) = 11
c) y'' − 6 y' + 9 y = 0 , y (0) = 0 e y' (0) = 5 d) y'' + 4 y' + 5 y = 0 , y (0) = 1 e y' (0) = 0

e) y'' + 4 y' + 2 y = 0 , y (0) = −1 e y' (0) = 2 + 3 2 f) y'' + 8 y' − 9 y = 0 , y (1) = 2 e y' (1) = 0

10. Comportamiento de las soluciones cuando x → ∞ :

a) Si a, b y c son constantes positivas, demuestre que todas las soluciones de ay’’ + by’ + cy = 0 tienden
a cero cuando x → ∞ .
b) Si a>0 y c>0, pero b=0. Demuestre que el resultado anterior deja de ser verdadero, pero que todas las
soluciones son acotadas cuando x → ∞ .
c) Si a>0 y b>0, pero c=0, demuestre que el resultado de a) deja de ser cierto pero que las soluciones tienden a
una constante que depende de las condiciones iniciales cuando x → ∞ . Determine esta constante para las
condiciones iniciales y(0)=y0, y’(0)=y0’.

11. Resolver las ecuaciones lineales no homogéneas mediante el método de Variación de Parámetros:
a) y'' − 4y' + 4y = ( x + 1) e2x b) y'' + y = sen t c) y'' − y = 1 / x

d) y'' + y = cos2 x e) y'' + y = sec x f) y'' − 3 y' + 2y = e3x / (1 + e x )

12. En los siguientes ejercicios resolver los problemas a valores iniciales lineales de orden dos
donde y (0) = 1, y' (0) = 0
a) y'' − y = xe x b) y'' + 2y' − 8y = 2e−2x − e− x

13. Aplicar el método de variación de parámetros, para hallar la solución general de:
x3 y''' + x 2 y'' − 2 xy' + 2 y = 2 x 4 , x > 0 ,sabiendo que x , x 2 y 1 / x son soluciones de la ecuación
homogénea.

14. Resolver las siguientes ecuaciones diferenciales lineales de segundo orden no homogéneas con coeficientes
constantes utilizando el método de coeficientes indeterminados:

a) y'' − 3 y' −10 y = 6e 4 x b) y'' + 10 y' + 25 y = 14e−5 x

c) y'' − 2 y' + 5 y = 25 x 2 + 12 d) y'' + 4 y = 4 cos 2 x + 6 cos x + 8 x 2 − 4 x
e) y'' − 2 y' + 2 y = e x senx f) y'' − 3 y' + 2 y = 14 sen 2 x −18 cos 2 x
 Práctico No 6: Sistemas de Ecuaciones Diferenciales lineales de primer orden

1. Considere el problema con valores iniciales u’’ + p (t ) u’ + q (t ) u = g (t ) , u (0) = uo , u’ (0) = u’ 0 .

Transforme este problema a un problema con valores iniciales para dos ecuaciones lineales de primer orden.

2. Ley de Fick (mecanismo de transporte): La tasa de transporte es directamente proporcional a la diferencia de

concentración en los dos compartimentos.

V1 V2

La figura es una representación compartamental de la membrana celular, el fluido extracelular y el fluido dentro de
la célula. Los componentes que necesita la célula para su crecimiento son transportados y pasan a través de la
membrana mediante un mecanismo llamado difusión pasiva el cual es gobernado por la Ley de Fick.
Sean V1 y V2 los volúmenes correspondientes, A1(t) y A2(t) la cantidad de sustancia en cada compartimento y k el
factor de permeabilidad. La concentración de sustancia en cada uno de los compartimentos queda definida entonces
como x(t) = A1(t) / V1 e y(t) = A2(t) / V2. Siendo el volumen constante, A’1(t) = x’(t) * V1 e A’2(t) = y’(t) *V2 y
aplicando la ley de Fick a la tasa de transporte de cantidad de sustancia obtenemos:

 A1 ' = k (y - x )  x' = k/V1 (y - x )

 ó, equivalentemente, 
 A 2 ' = k (x - y)  y' = k/V2 (x - y)

3. Lleve el sistema de ecuaciones dado a una ecuación lineal de orden dos. Encuentre las soluciones y analice el
comportamiento de las mismas a largo plazo. Bosqueje una gráfica de ambas en un mismo gráfico.

4. Una población de peces migra desde el lago 1 al lago 2. Sea x(t) el número de peces en el lago 1 e y(t) el
número de peces en el lago 2. La tasa de migración entre los dos lagos se asume que está gobernada por la ley
de Fick. Los volúmenes de ambos lagos son V1 = 2 Km3 y V2 = 1 Km3. Al comienzo de la estación, el stock de
peces en el lago 1 es de 30 000. a) Plantee el sistema de ecuaciones correspondiente. b) encuentre el número
final de peces en cada lago asumiendo que no han habido pérdidas por pesca o muerte. (Nota: el valor de k es
uno).

5. Determinar la naturaleza y propiedades de estabilidad del punto crítico (0,0) para cada uno de los siguientes
sistemas autónomos lineales. Hallar la ecuación de las trayectorias, resolver esa ecuación y dibujar algunas de
las trayectorias indicando la dirección de t creciente.

 dx  dx 
 dx 
 dx
 = 2 x  = −x − 2 y 
 dt = 3x − 4 y 
 dt = −4 x − y
 dt dt
a ) b ) c ) d )
 dy  dy 
 dy 
 dy
 = 3 y  = 4 x − 5 y 

= 2x + 3y 

= x−2y
 dt  dt  dt
  dt

 dx  dx  dx


 dt = 4 y  dt = −x 
 dt
=x
e ) f ) g )
 dy  dy 
 dy
 = −x  = −2 y  = −x + 2 y
 dt  dt  dt

6. El modelo de crecimiento poblacional no acotado para una especie se describe con la ecuación: y'= k y , cuya
solución es y= C e kt .Para el caso de interacción entre dos especies, siendo x (t ) la población de la especie 1
x' = a x + by
en el tiempo t e y (t ) la población de la especie 2 en el tiempo t, será: 
 y' = c x + d y

Si las especies compiten por un mismo recurso: b, c < 0 (el término b y c, representa el tipo de interacción). En
una interacción depredador-presa, asumimos que el número de presas consumidas por unidad de tiempo varía
directamente con el número de depredadores. Si la especie 1 consume exclusivamente especie 2, luego b>0 y c>0.

En una relación simbiótica, la presencia de la especie 2 puede incrementar la tasa de supervivencia de la especie 1.
A su vez, un miembro de 2 puede proveer refugio a un miembro de 1, b > 0.

En los siguientes ejercicios determine: a) las soluciones x(t) e y(t); b) el tipo de relación entre las poblaciones y c)
el comportamiento de las poblaciones a largo plazo.

 x' = x - 2.5y x' = x - 2.5 y x' = - x + 2 y x' = 2x − y

 y' = - 1.6x + y  y' = 1.6x - y  y' = x - 2 y  y' = - x + 2y
   
a)  b)  c)  d) 
 x(0) = 1500  x(0) = 1500  x(0) = 450  x(0) = 500
 y(0) = 1000.  y(0) = 1000  y(0) = 150  y(0) = 300

6-a) 6-b)
3000 1500

2000 1000

500
1000
2 4 6 8 10
-500
0.2 0.4 0.6 0.8 1
-1000
-1000
-1500

6-c) 6-d)
450

400 2000

350
1000
300

250
0.2 0.4 0.6 0.8 1
200
-1000

1 2 3 4 5
 I (t)
1 F21 2
7. Dado el sistema
F01
F12
F20
F10

donde: Fij = tasa de flujo del tanque i al tanque j [l /min];I (t) trazador agregado por unidad de tiempo [gr/min]; Vi
volumen del tanque i [l]; x (t) cantidad de trazador en el tanque 1[gr]; y (t) cantidad de trazador en el tanque 2 [gr];
ci concentración en el tanque i (c1 = x(t) / V 1 ; c2 = y(t) / V 2) [gr/l].
Las ecuaciones serán:

 x ' = −( F12 + F10 ) c1(t ) + F21c 2 (t ) + I (t )


 y' = F12 c1(t ) − ( F21 + F20 )c2 (t )

8. Sea un sistema de dos tanques como se muestra en la figura. Si se vierten 50 gr de pintura en el tanque 1,
plantee las ecuaciones para x (t) e y (t), cantidades de pintura en cada uno de los tanques al tiempo t. Si la
concentración de pintura es de 1 gr/gal en el flujo del tanque 1 y consideramos x (0) = y (0) = 0, encuentre y
grafique las soluciones.

5 gal/min
5 gal/min

5 gal/min

10 gal/min

9. En cada uno de los problemas propuestos resuelva el sistema de ecuaciones correspondiente por eliminación

 x1' = x1 + x2  x 1' = 3 x 1 − 2 x 2 − e −t sent  x 1' = 2 x 1 − 5 x 2 − sen 2t x 1 (0) = 0

 '  −t 
a) x2 = 4 x1 − 2 x2 b) x 2' = 4 x 1 − x 2 + 2e cos t c ) x 2' = x1 − 2 x 2 + t x 2 (0) = 1
  
  

10. En cada uno de los problemas siguientes resuelva el sistema de ecuaciones dado, comprobando en cada caso
que aparece el número apropiado de constantes en la solución general.

( D 2 − 3D + 2 )x 1 + ( D − 1)x 2 = 0 ( D 2 − 4D + 4 )x 1 + 3Dx 2 = 1

a ) b )
 ( D − 2 )x 1 + ( D + 1)x 2 = 0  ( D − 2 )x 1 + ( D + 2 )x 2 = 0

( 2D − 1)x 1 + Dx 2 = t ( D 2 − 4D + 4 )x 1 + ( D 2 + 2D )x 2 = 0
c ) d )
 ( D − 1)x 1 + Dx 2 = 2  ( D 2 − 2D )x 1 + ( D 2 + 4D + 4 )x 2 = 0

12. La quema de combustibles fósiles por la industria adiciona dióxido de carbono desde la tierra hacia la atmósfera
a una tasa de aproximadamente 6x109 to /año. Este exceso de carbono debe ser absorbido por el océano para evitar
acumulaciones a largo plazo. En la figura se muestra un modelo simple de dos compartimentos para el ciclo del
dióxido de carbono. Bolin* planteo el siguiente problema de valores iniciales, estimando las condiciones iniciales x
(0) = 7 x 109 to e y (0) = 35000 x 109 to.
 x' = - 0.2 x + 0.0025y + I

 y' = 0.2x - 0.0025y

a) si I = 0, encontrar las soluciones x(t) e y(t). Determine el comportamiento a largo plazo de las soluciones y
analizar la estabilidad del sistema.
b) idem al punto a si I = 6x109 to /año.

I
0.22
Atmósfera Océano
0.002

13. Plantear las ecuaciones que describen el siguiente sistema.

I (t) = I
a12

a10

14. En el flujo sanguíneo, los iones potasio se transportan desde el plasma hacia los glóbulos rojos y viceversa. Este
comportamiento sugiere un sistema simplificado como el que se muestra en la figura

a12
Plasma Glóbulos
rojos
a21

Si se inyecta una cantidad fija de potasio radioactivo ( K42) en el flujo sanguíneo, y no existen pérdidas en el
sistema: Plantee las ecuaciones correspondientes y determine el funcionamiento a largo plazo de las mismas. Desde
x(t )
la muestra, se ajustó la ecuación = 0.06 + 0.94 e −0.3 t a los datos obtenidos, encuentre a12 y a 21.
x0

15. La Creatina (fosfato de creatina) es una componente de la orina que es utilizada para suplir energía a los
músculos. Cuando una dosis de Creatina se inyecta en el flujo sanguíneo, subsecuentes muestras de sangre
muestran que la concentración puede representarse por una función de la forma: c1 e − r1t + c 2 e − r2t .
Esto sugiere el uso de un modelo de dos compartimentos como se muestra en la figura:

a12
Sangre Tejido
a21
a10
Eliminado a través
de la orina

Plantee el correspondiente sistema de ecuaciones. Muestre que x'' + ( a12 + a21 + a10 )x' + a10 a21 x = 0 y que
ambas raíces de la ecuación característica son negativas. Encuentre las soluciones x (t) e y (t).
16. Para cada sistema no lineal hallar sus puntos críticos y resolver la ecuación de la trayectoria asociada
 dx 2  dx


 dx


 dt = y( x −1 )  = y( x 2 −1 )  = −x
 dt  dt
a ) b ) c )
 dy 2 
 dy 2 
 dy
 = 2 xy  = −x( x −1 )  = 2 x2 y2
 dt 

 dt 

 dt
17. Hallar el punto crítico del sistema, llevar mediante cambio de variables el punto crítico al origen y analizar su
estabilidad.

 dx

 = 2 x − 2 y + 10
 dt


 dy
 = 11x − 8 y + 49

 dt


d2x dx
18. La ecuación 2
+ 2b + a 2 x = 0 permite estudiar las vibraciones libres de una masa sujeta a un muelle.
dt dt

a) Convertir la ecuación de orden 2 en un sistema de dos ecuaciones de orden 1.

b) En los cuatro casos siguientes, describir la naturaleza y las propiedades de estabilidad del punto crítico y dar
una breve interpretación del movimiento correspondiente de la masa.
i) b=0 ii) 0<b<a iii) b=a iv) b>a

19. Encontrar el punto crítico de un modelo presa predador de Volterra, analizar su estabilidad y encontrar su
trayectoria.

 dx

 = x( a − by )
 dt


 dy
 = − y( c − dx )

 dt


20. El lago Erie tiene un volumen de 480 km3. Los flujos de entrada, proveniente del lago Hurón y de salida hacia
el lago Ontario son de 350 km3/año. Suponga que en el instante t = 0 (años) la concentración de
contaminación del lago Erie es cinco veces la del lago Hurón. a) plantee y resuelva la ecuación diferencial que
modela la contaminación en el lago Erie. b) ¿Cuánto tiempo le tomara reducir la concentración de
contaminación en el lago Erie al doble de la del lago Hurón?