Iedo 1

CAPITULO I
Georg Friedrich Riemann (1826-1866)
Georg Friedrich Riemann

Ejemplo vivo de la timidez y de la fragilidad física, Riemann impactó, sin embargo, el mundo de las
matemáticas como pocos lo han hecho en la historia.
Hijo del pastor de un pequeño pueblo en Alemania, recibió no obstante una buena educación que
lo llevó a presentar su tesis doctoral delante de Gauss en Gottingen.
Este último, reconocido como difícil de sorprender, quedó entusiasmado por el desarrollo que hizo
Riemann sobre la teoría de la función de una variable compleja. Este episodio se recuerda como la
única vez en la que Gauss haya expresado admiración por un trabajo ajeno.
Ahí aparecen las famosas superficies de Riemann, las cuales generarían el enfoque topológico del
análisis. Un poco más tarde clarificó la noción de integral mediante una nueva definición conocida
como la "Integral de Riemann".
Sus trabajos sobre los fundamentos de la geometría le permitieron generalizar la noción de
espacio y son precursores de las teorías del siglo XX sobre los espacios abstractos.
Pero su complexión débil lo hizo presa de la tuberculosis, un mal entonces incurable, y Riemann
murió en 1866 a los 40 años. Sus obras, que caben en pocas páginas, son de una densidad tal que
dejan trabajo e ideas incluso para los matemáticos de hoy en día.
Las ecuaciones diferenciales aparecen con frecuencia en muchas ramas de la matemática, y sirven
para plantear y resolver problemas provenientes de la física, la ingeniería, la economía, la biología
y de las ciencias sociales, entre otras disciplinas.
Julio Flores Dionicio Página 1

En esta notas de clase solamente se estudiara ecuaciones ordinarias, para empezar nuestro
estudio es útil introducir ciertas definiciones básicas y la noción de solución de una ecuación
diferencial.
DEFINICIONES Y TERMINOLOGÍA
DEFINICIÓN 1: ECUACION DIFERENCIAL
Una ecuación que contiene derivadas de una o más variables respecto a una o más
variables independientes, se dice que es una ecuación diferencial (ED).
En todo este trabajo las derivadas ordinarias se escribirán usando la notación de Leibniz
dy d 2 y dny
, 2 ,..., n , o la notación prima y, y, y,..., y  n 
dx dx dx
ECUACION DIFERENCIAL ORDINARIA

Si una ecuación contiene sólo derivadas de una o más variables dependientes respecto a una sola
variable independiente se dice que es una ecuación diferencial ordinaria (EDO)
dx dy
Ejemplos de ecuaciones diferenciales ordinarias   y  3 x , y  xy y
dt dx
xy  csen( x) y  tan( x)
ECUACION DIFERENCIAL PARCIAL

Una ecuación que involucra derivadas parciales de una o más variables dependientes de dos o más
variables independientes se llama ecuación diferencial parcial (EDP).
u u  3u  2u u
Las x  y u , y  2 4 son ecuaciones diferenciales parciales
x y x 3
t t
Las ecuaciones diferenciales ordinarias son identidades que vinculan una función con sus
derivadas. Por ejemplo, si y  t  denota el número de bacterias en una colonia en función del
tiempo, la ecuación diferencial y t   ky t  , donde k es una constante positiva, expresa que el
aumento de la población bacteriana, representada por la derivada y , es proporcional a la propia
población y, esto es, mientras más bacterias hay, más rápido ellas se multiplican.
DEFINICIÓN 2: ORDEN DE UNA ECUACION DIFERENCIAL

El orden de una ecuación diferencial (ya sea EDO o EDP) es el orden de la mayor
derivada en la ecuación.
Por ejemplo y  xy es de orden 1 ; y   y   xy  0 es de orden 2.

3

Simbólicamente podemos expresar una ecuación diferencial ordinaria de n-ésimo orden con una
variable dependiente por la forma general F ( x, y( x), y( x), y( x)...., y ( n ) ( x))  0 , donde F es una
función con valores reales de n + 2 variables: x, y ( x), y( x), y( x)...., y ( n ) ( x) .
NOTA
Por razones tanto prácticas como teóricas, de ahora en adelante supondremos que
es posible resolver una ecuación diferencial ordinaria en la forma de la ecuación
F ( x, y( x), y( x), y( x)...., y ( n ) ( x))  0 únicamente para la mayor derivada y ( n )
en términos de las n + 1 variables restantes.
dny
La ecuación diferencial n
 f ( x, y, y, y...., y ( n 1) ) , donde f es una función continua con valores
dx
reales, se conoce como la forma normal de la ecuación F ( x, y( x), y( x), y( x)...., y ( n ) ( x))  0 . Así
dy
que cuando sea adecuado para nuestros propósitos, usaremos las formas normales  f ( x, y ) y
dx
d2y
 f ( x, y, y) para representar en general las ecuaciones diferenciales ordinarias de primer y
dx 2
segundo orden respectivamente.
DEFINICIÓN 3: DE UNA ECUACION DIFERENCIAL ORDINARIA LINEAL

Una EDO lineal de orden n es de forma
an ( x) y ( n )  an ( x) y ( n 1)  ....a1 ( x) y  a0 ( x) y  f ( x)
Donde las funciones ai ( x)  , i 1,...n son llamadas coeficientes de la EDO lineal
DEFINICION 4. f(x) (llamado comúnmente lado derecho de la EDO) es idénticamente nulo, la EDO
lineal se dice homogénea. Si f ( x)  0 la EDO lineal se dice no homogénea.
DEFINICION 5. Si los coeficientes ai ( x) no dependen de x, se dice que la EDO lineal es a

coeficientes constantes. De lo contrario se dice que ella es a coeficientes variables.
En el caso que an ( x)  0 se puede dividir la EDO (1) por an ( x) . La EDO que así se obtiene queda
an 1 ( x ) ( n 1) a ( x) a ( x) f ( x)
y(n)  y  .... 1 y  0 y  g ( x ) con , g ( x) 
an ( x ) an ( x ) an ( x ) an ( x)
Se dice que está normalizada.
EJEMPLO 1. xy  sen  x   tg  x  . EDO lineal de orden 1 a coeficientes variables, no homogénea

no normalizada.

EJEMPLO 2. y  2 y  0 . EDO lineal de orden 3 de coeficientes constantes, homogénea y
normalizada.
Una EDO no lineal es simplemente una EDO que no es lineal. Gran parte de esta nota de clases
concierne en el estudio de EDO lineales.
EJEMPLO.  
y 1   y  k 2 .EDO No lineal, de orden 1. Esta es la EDO de curva
2
braquistócrona
SOLUCION DE UNA ECUACION DIFERENCIAL ORDINARIA
Cualquier función   x  , defi nida en un intervalo I y que tiene al menos n derivadas

continuas en I , las cuales cuando se sustituyen en una ecuación diferencial ordinaria de n-
ésimo orden reducen la ecuación a una identidad, se dice que es una solución de la ecuación
en el intervalo.
En otras palabras, una solución de una ecuación diferencial ordinaria de n-ésimo orden
a ( x) y
n
(n)
 a ( x) y
n
( n 1)
 ....a ( x) y  a ( x) y  f ( x)
1 0 es una función  que
posee al menos n derivadas para las que F ( x,  ( x),  ( x),  ( x)....,  ( n ) ( x))  0 para toda x en I , y
decimos que  satisface la ecuación diferencial en I . Para nuestros propósitos supondremos que
una solución  es una función con valores reales. En la mayoría de los casos será conveniente
denotar una solución con el símbolo alternativo y(x).
El intervalo I en la definición también se conoce con otros nombres como son intervalo de defi
nición, intervalo de existencia, intervalo de validez, o dominio de la solución y puede ser un
intervalo abierto (a, b), un intervalo cerrado [a, b], etcétera.
FAMILIAS DE SOLUCIONES
El estudio de ecuaciones diferenciales es similar al del cálculo integral. En algunos libros una
solución  es algunas veces llamada integral de la ecuación y su gráfi ca se llama curva integral.
Cuando obtenemos una antiderivada o una integral indefi nida en cálculo, usamos una sola
constante C de integración.
De modo similar, cuando resolvemos una ecuación diferencial de primer orden F  x, y, y  0 ,
normalmente obtenemos una solución que contiene una sola constante arbitraria o parámetro C .
Una solución que contiene una constante arbitraria representa un conjunto G  x, y, C   0 de
soluciones llamado familia de soluciones uniparamétrica.
 
Cuando resolvemos una ecuación diferencial de orden n, F x, y, y, y,..., y    0 , buscamos una
n
familia de soluciones n-paramétrica G  x, y, C1, C2 ,..., Cn   0 . Esto signifi ca que una sola
ecuación diferencial puede tener un número infi nito de soluciones correspondiendo a un número
ilimitado de elecciones de los parámetros. Una solución de una ecuación diferencial que está libre
de la elección de parámetros se llama solución particular. Por ejemplo, la familia uniparamétrica

y  cx  x cos x es una solución explícita de la ecuación lineal de primer orden xy  y  x 2 senx en
el intervalo  .    . La siguiente fi gura muestra las gráfi cas de algunas de las soluciones en
esta familia. La solución y   x cos x , la curva azul en la fi gura, es una solución particular
correspondiente a c = 0.
En lo sucesivo no haremos distinción entre las curvas integrales y las funciones que las definen.
Por otra parte, la solución general de una ecuación diferencial es una expresión que agrupa, de
manera compacta y explicita, la familia de todas las soluciones. Si toda solución de una ecuación
 
de orden n, F x,   x  ,    x  ,    x  ,...,     x   0 , en un intervalo I, se puede obtener partiendo
n
de una familia n-paramétrica G  x, y, C1, C2 ,..., Cn   0 , con valores adecuados de los parámetros
Ci (i=1,2, . . ., n), se dice que la familia es la solución general de la ecuación diferencial. Al resolver
las ecuaciones diferenciales lineales vamos a imponer restricciones relativamente sencillas a los
coeficientes de esas ecuaciones.
Con estas restricciones siempre nos aseguraremos no sólo de que exista una solución en un
intervalo, sino también de que una familia de soluciones contenga todas las soluciones posibles.
Las ecuaciones no lineales, a excepción de algunas de primer orden, son difíciles de resolver e
incluso resultan irresolubles, en términos de las funciones elementales comunes (combinaciones
finitas de potencias o raíces enteras de x, de funciones exponenciales y logarítmicas, o funciones
trigonométricas o trigonométricas inversas). Además, si en cierto momento nos encontramos con
una familia de soluciones de una ecuación no lineal, no es obvio cuando la familia es una solución
general. Por lo anterior, y en un nivel práctico, el nombre “solución general” sólo se aplica a las
ecuaciones diferenciales lineales.
EJEMPLO 1. Consideremos la ecuación diferencial y’ = 0 en un intervalo I. recordemos que la

derivada de una función en un intervalo es cero si, y solo si, la función es constante. Por lo tanto,
la solución general de la ecuación es y( x)  C , con C 
En muchas ocasiones las ecuaciones diferenciales revelan una gran información cualitativa acerca
de las curvas integrales, como vemos en el siguiente ejemplo:

EJEMPLO 2. Consideremos la ecuación diferencial 𝑦 ′ = 1 − |𝑦| en ℝ , más adelante estudiaremos
las técnicas que nos permitirán resolver este tipo de ecuación. Sin embargo, al examinar la
ecuación diferencial podemos tener una idea de cómo son las curvas integrales. Vemos que las
funciones constantes 𝑦(𝑥) ≡ 1 y 𝑦(𝑥) ≡ −1 son de clase 𝐶 1 y satisfacen 𝑦 ′ = 1 − |𝑦| en ℝ, por lo
que define curvas integrales de la EDO. Supongamos ahora que 𝑦: ℝ → ℝ define una curva integral
de la EDO. En los puntos donde y toma valores entre -1 y 1, su derivada es positiva; luego y es
creciente. En otro caso, y es decreciente. Además mientas más cerca esté y(x) de los valores 1 y -1,
más cercana a 0 será su derivada. Si al contrario, |𝑦(𝑥)| es muy grande, la pendiente será muy
pronunciada.
Otro dato que obtenemos directamente de la ecuación es que si 𝑦: ℝ → ℝ define una curva
integral, entonces la función 𝑦𝑐 : ℝ → ℝ definida por 𝑦𝑐 (𝑥) = 𝑦(𝑥 + 𝑐) también define una curva
integral (basta sustituir en la ecuación). Esto nos dice que al trasladar una curva integral hacia la
derecha o hacia la izquierda obtenemos mas curvas integrales. Con toda esta información
deducimos que, de existir, las curvas integrales deben tener un aspectos bastantes similar a las
que se muestran a en la gráfica siguiente.

PROBLEMAS DE VALOR INICIAL (PVI)
A menudo nos interesa resolver una ecuación diferencial sujeta a condiciones prescritas, que son
las condiciones que se imponen a y(x) o a sus derivadas. En algún intervalo I que contenga a x0 , el
problema
d n y

 n  F x, y, y, y,... y
 dx
 n 1

, en donde y0 , y1 , y2 ,..., yn 1
 y( x )  y , y( x )  y , y( x )  y ,... y ( n1) ( x )  y
 0 0 0 1 0 2 0 n 1
son constantes reales especificadas arbitrariamente, se llama problema de valor inicial (PVI). Los
valores dados de la función desconocida, y(x), y de sus primeras n - 1 derivadas en un solo punto
x0 : y ( x0 )  y0 , y( x0 )  y1 , y( x0 )  y2 ,... y ( n 1) ( x0 )  yn 1 se llaman condiciones iniciales. EI PVI
emmciado con las ecuaciones iniciales, también se denomina problema de valor inicial de
enésimo orden
 dy d2 y
  F  x, y   2  F  x, y, y 
LOS PVI DE PRIMER Y SEGUNDO ORDEN  dx  dx ,
 y ( x0 )  y0 y
 y ( x0 )  y0 , y( x0 )  y1
respectivamente son fáciles de interpretar en términos geométricos.
 dy
  F  x, y 
Para las ecuaciones  dx estamos buscando una solución de la ecuación diferencial
 y ( x0 )  y0
en un intervalo I que contenga a x0 , tal que la curva de solución pase por el punto prescrito
( x0 , y0 )

d2 y
 2  F  x, y, y 
Para las ecuaciones  dx , deseamos determinar una solución de la ecuación
 y ( x )  y , y( x )  y
 0 0 0 1
diferencial cuya gráfica no sólo pase por ( x0 , y0 ) , sino que también pase por ese punto de tal
manera que la pendiente de la curva en ese lugar sea y1 (ver figura). El término condición inicial
y ( x0 )  y0  y( x0 )  y1 procede de los sistemas físicos en que la variable independiente es el
tiempo t y donde y (t0 )  y0  y(t0 )  y1 representan, respectivamente, la posición y la velocidad
de un objeto en cierto momento o tiempo inicial t0 .
A menudo, la solución de un problema de valor inicial de orden n entraña la aplicación de una
familia n-paramétrica de soluciones de la ecuación diferencial dada para determinar n constantes
especializadas, de tal modo que la solución particular que resulte para la ecuación “se ajuste” (0
satisfaga) a las n condiciones iniciales.
EXISTENCIA Y UNICIDAD DE SOLUCION DE UNA ECUACION DIFERENCIAL

Al resolver un problema de valor inicial surgen dos asuntos fundamentales:
¿Existe una solución al problema? Si la hay, ¿es única?
Para un problema de valor inicial, como

 dy
  F  x, y 
 dx
 y ( x0 )  y0 , lo que se pregunta es:
dy
¿La ecuación diferencial  F  x, y  tiene solución?(EXISTENCIA).¿Alguna de las curvas
dx
solución pasa por el punto ( x0 , y0 ) ?. ¿Cuándo podemos estar seguros de que hay precisamente
una curva solución que pasa por el ( x0 , y0 ) ?(UNICIDAD)

Dentro de los confines seguros de un curso formal de ecuaciones diferenciales, se puede asumir,
que la mayor parte de las ecuaciones diferenciales tienen soluciones y que las soluciones de los
problemas de valor inicial probablemente sean únicas. Sin embargo, en la vida real las cosas no
son tan simples. Por consiguiente, antes de resolver un problema de valor inicial es preferible
conocer, si existe una solución y, cuando exista, si es la única. Puesto que vamos a manejar
ecuaciones diferenciales de primer orden en estos dos primeros capítulos, enunciaremos aquí, sin
demostrarlo, un teorema que define las condiciones suficientes para garantizar la existencia y
unicidad de una solución a un problema de valor inicial de primer orden, para ecuaciones que
 dy
  F  x, y 
tengan la forma de las ecuaciones  dx
 y ( x0 )  y0
TEOREMA: Existencia de una solución única

Sea R una región rectangular del plano xy, definida por a  x  b, c  y  d , que contiene
F
al punto ( x0 , y0 ) en su interior. Si F  x, y   son continuas en R , entonces existe algún
y
intervalo I  :  x0  h, x0  h  , h  0 contenido en  a, b (centrado en x0 ), y una función única,
 dy
  F  x, y 
y(x) definida en I  , que satisface el problema de valor inicial  dx .
 y ( x0 )  y0
El resultado anterior es uno de los teoremas más comunes de existencia y unicidad para
ecuaciones de primer orden, ya que es bastante fácil comprobar los criterios de continuidad de
F
F  x, y   .
y
En la figura podemos ver la interpretacidn geometrica del teorema.

EJEMPLO 1
La función y  Ce x es una familia uniparamétrica de soluciones de la ecuación de primer orden

y  y . Todas las soluciones en esta familia están defi nidas en el intervalo  ,  . Si
 y  y
imponemos una condición inicial, digamos, y  0  3 , tendríamos el PVI:  . Ahora al
 y  0   3
sustituir x  0  y  3 en la familia se determina la constante 3  Ce0  C por lo tanto y  3e
x
es una solución del PVI planteado

Ahora si hacemos que la curva solución pase por el punto 1, 2 en lugar de  0,3 , entonces
y 1  2 se obtendrá 2  y 1  Ce1  C  2e1 . En este caso y  2e es una solución del
x 1
 y  y

PVI 
 y 1  2

EJEMPLO 2
Se puede demostrar fasilmenta que que una familia uniparamétrica de soluciones de la ecuación
1
diferencial de primer orden y  2 xy 2  0 , es y  2 . Si establecemos la condición inicial
x C
y  0  1 , entonces al sustituir x  0  y  1 en la familia de soluciones, se obtiene
1 1
1   C  1 . Así y  2 . Ahora enfatizamos las siguientes tres diferencias:
C x 1
1
• Considerada como una función, el dominio de y  x   2 es el conjunto de todos los números
x 1
reales x para los cuales y  x  está defi nida, excepto x  1 y en x  1 . Véase la fi gura

• Considerada como una solución de la ecuación diferencial y  2 xy 2  0 , el intervalo I de defi
1
nición de y  x   2 podría tomarse como cualquier intervalo en el cual y(x) está defi nida y es
x 1
derivable. Como se puede ver en la fi gura anterior, los intervalos más largos en los que
1
y  x  2 es una solución son  , 1 ,  1,1 y 1,  .
x 1
 y  2 xy  0
 2
• Considerada como una solución del problema con valores iniciales  , el intervalo I

 y  0    1
1
de defi nición de y  x   2 podría ser cualquier intervalo en el cual y(x) está defi nida, es
x 1
derivable y contiene al punto inicial x  0 ; el intervalo más largo para el cual esto es válido es
 1,1 . Véase la curva roja en la fi gura
INTERVALO DE EXISTENCIA Y UNICIDAD

Suponga que y  x  representa una solución del problema con valores iniciales
 dy
  F  x, y 
 dx .
 y ( x0 )  y0
Los siguientes tres conjuntos de números reales en el eje x pueden no ser iguales: el dominio de la
función y  x  , el intervalo I en el cual la solución y(x) está definida o existe, y el intervalo I  de
existencia y unicidad.
El ejemplo 2 muestra la diferencia entre el dominio de una función y el intervalo I de definición.
Ahora suponga que  x0 , y0  es un punto en el interior de la región rectangular R en el teorema
de existencia y unicidad. Esto da como resultado que la continuidad de la función f  x, y  en R

por sí misma es suficiente para garantizar la existencia de al menos una solución de
 dy
  F  x, y 
 dx , definida en algún intervalo
 y ( x0 )  y0
I . El intervalo I de definición para este problema con valores iniciales normalmente se toma
como el intervalo más grande que contiene x0 en el cual la solución y  x  está definida y es
derivable.
El intervalo I depende tanto de F  x, y  como de la condición inicial y  x0   y0 . La condición

extra de continuidad de la primera derivada parcial F / y en R nos permite decir que no sólo
existe una solución en algún intervalo I  que contiene x0 , sino que esta es la única solución que
satisface y  x0   y0 . Sin embargo, el teorema de existencia de una solución única no da ninguna
indicación de los tamaños de los intervalos I e I  ; el intervalo de definición I no necesita ser tan
amplio como la región R y el intervalo de existencia y unicidad I  puede no ser tan amplio como
I . El número h  0 que define el intervalo I  :  x0  h, x0  h  podría ser muy pequeño, por lo que
es mejor considerar que la solución y(x) es única en un sentido local, esto es, una solución definida
cerca del punto  x0 , y0  ). Véase el siguiente ejemplo.
EJEMPLO
0, x  0
1 4 
Las funciones y  x   x ,   x   y y  x    1 4 , x  0 , tienen el mismo dominio pero
16 16 x
son obviamente diferentes.
Véanse las figuras a) y b), respectivamente.
 dy 3
  xy 2
Ambas funciones son soluciones del problema con valores iniciales  dx en el intervalo
 y  2  1

 ,  .

COMENTARIOS
(i) Las condiciones del teorema de existencia de una solución única son suficientes pero no
necesarias. Esto significa que cuando F  x, y   F / y son continuas en una región rectangular
 dy
  F  x, y 
R , debe siempre seguir que existe una solución de la ecuación  dx y es única siempre
 y  x0   y0

que  x0 , y0  sea un punto interior a R . Sin embargo si las condiciones establecidas en la hipótesis
del teorema de existencia de una solución única no son válidas, entonces puede ocurrir
 dy
  F  x, y 
cualquier cosa: el problema de la ecuación  dx puede tener una solución y esta
 y  x0   y0

 dy
  F  x, y 
solución puede ser única o la ecuación  dx puede tener varias soluciones o puede no
 y  x0   y0

tener ninguna solución.
dy
(ii) Suponga que la ecuación diferencial de primer orden  F  x, y  tiene una familia
dx
uniparamétrica de soluciones y que F  x, y  satisface la hipótesis del teorema de existencia de
una solución única en alguna región rectangular R del plano xy. Dos curvas solución diferentes
no se pueden interceptar o ser tangentes entre sí en un punto  x0 , y0  en R .
Ya podemos resolver algunas ecuaciones diferenciales. Comenzaremos con las de primer orden
elemental y veremos cómo hacerlo; el método dependerá del tipo de ecuación.
ECUACIONES DIFERENCIALES ORDINARIAS

DE PRIMER ORDEN
Una ecuación diferencial ordinaria de primer orden se escribe: y  F  x, y  , o
F  x, y, y  0 y se dice que y   ( x) es solución en el intervalo I de esta ecuación si se
verifica ( x)  F ( x, ( x)) , x  I , es decir, si cuando se sustituye en la ecuación y por su
expresión e y por la expresión de la derivada, lo que se obtiene es una identidad, algo que es
cierto para todo x  I .

EJEMPLO La función y  e
x
es solución de la ecuación y   y en todo , ya que
y( x)  e x   y( x) x  I
x
Pero también es solución cualquier función de la forma y  Ce siendo C una constante
arbitraria, puesto
Que y( x)  Ce
x
  y ( x) x  I
Así pues, la ecuación del Ejemplo tiene infinitas soluciones, lo que no es una particularidad de
esta ecuación concreta. La ecuación diferencial ordinaria y  F  x, y  posee, en general, una
«familia» de infinitas soluciones dependientes de una constante arbitraria, a la que se llama
solución general de y  F  x, y  . Para cada valor de dicha constante arbitraria se obtiene una
solución particular.
Se llama resolver una ecuación diferencial a encontrar su solución general. En realidad, esto sólo
es posible para unas cuantas (pocas) ecuaciones sencillas. Para la inmensa mayoría de las
ecuaciones diferenciales es necesario recurrir a métodos numéricos y calcular soluciones
aproximadas con ayuda de un ordenador.
Con frecuencia lo que interesa en las aplicaciones es encontrar una solución particular que
verifique alguna condición eadicional. Por ejemplo, que toma un valor dado para un valor,
también dado, de la variable independiente.
ECUACIONES DIFERENCIALES ORDINARIAS

ELEMENTALES DE PRIMER ORDEN
A continuación estudiaremos algunas técnicas que nos permitirán determinar las soluciones de
una gran cantidad de EDO. Comenzaremos por analizar tres tipos de ecuaciones de primer orden
que llamaremos elementales.
 PRIMER TIPO: INTEGRACION DIRECTA (ID): y  f ( x) (1)

 SEGUNDO TIPO. VARIABLE SEPARABLE (vs): y  f ( x) g ( y) (2)
 TERCER TIPO: LINEALES DE PRIMER ORDEN (LI): y  a t  y  b t  (3)
Posteriormente estudiaremos algunas otras ecuaciones de primer y segundo orden que pueden
reducirse a estos elementales.
Para resolver los casos elementales, será útil el cálculo de primitivas y recordar el conocido
Teorema Fundamental del Cálculo.
TEOREMA (TFC). Sea f integrable en  a, b , entonces, dado x0  a, b e y0  , la función y,
 
para x  a, b ,
x
definida por y( x)  y0   f (s)ds ,
x0
Es continua en  a, b y se tiene y ( x0 )  y0 . Si además f es continua en  a, b entonces la función

y(x) es también derivable en  a, b con derivada continua igual a f  x  , esto es, se tiene que
x d x
y( x)  y0   y( s)ds , x a, b f ( s )ds  f ( x) x a, b
ds x0
y
x0
Nos enfocaremos en la resolución de las ecuaciones y no siempre seremos demasiados rigurosos

en los aspectos teóricos que justifican los cálculos. Estos aspectos (existencia, unicidad,
regularidad de la solución) serán tratados con más profundidad en otros cursos.
PRIMER TIPO: INTEGRACION DIRECTA (ID)
Es una ecuación de la forma

y  f ( x)
Si la función f es integrable (por ejemplo si es continua o continua por pedazos), las soluciones
existen y están dadas por: y   f ( x)dx  C Donde
C  R es una constante arbitraria.
EJEMPLO 1. Las soluciones de la ecuación y  sen( x) son de la forma
y   sen( x)dx  C   cos( x)  C , con C  R
EJEMPLO 2. La ecuación y  x tiene como soluciones a las funciones de la forma

x2
y   xdx  C  C , con C  R
2
Denotaremos por C genéricamente a una constante arbitraria sin importar las eventualidades
transformaciones biyectivas que la mantienen arbitraria (ponderaciones por un escalar no nulo,
cambios de signo, suma de otra constante). Por ejemplo C , 2C , C 2 , C  4 pueden ser
representados por una misma constante genérica. Si la transformación no es biyectiva, por
ejemplo C 2 , entonces se pierde la arbitrariedad y es mejor escribir explícitamente C 2 o precisar
de algún modo que la constante no puede ser negativa.
1
EJEMPLO 3. Estudiaremos la ecuación y  para x  0 .
x
dx
SOLUCIÓN y   x
 C  ln( x )  C  ln( x )  ln( k )  ln(k x ) , k  0 ,
En el cálculo anterior hemos reemplazado la constante arbitraria C por ln(k) (donde k  e  0 )

C
para escribir la solución de manera mas compacta y sin perder arbitrariedad. Observemos también
1
que el modulo en el logaritmo nos da la primitiva correcta de cuando x <0. Como x ni es
x
derivable en x  0 , se considera la resolución de la EDO separadamente en cada intervalo  ,0 
y  ,0  .
PROBLEMA DE CAUCHY (PVI)

Supongamos que ahora que queremos encontrar la solución de y  f ( x) definida sobre un
intervalo I y que además pase por cierto punto  x0 , y0  dado. Esto es, queremos resolver el
siguiente problema (PROBLEMA DE CAUCHY).
 y( x)  f ( x) para todo x  I ,


 y( x0 )  y0
Integrando la ecuación y  f ( x) entre x0 y x  I obtenemos
x x
 y( s)ds   f ( s)ds
x
y del TFC , se tiene que y( x)  y( x0 )   f ( s)ds .
x0 x0 x0

x
Detengámonos aquí para comprobar (7) y (8). La función F ( x)   f ( s)ds en (8) es una
x0
primitiva bien particular de f: aquella que cumple F ( x0 )  0 . Entonces la constante C en (7)

deja de ser arbitraria y se tiene que C  y ( x0 ) .
Esto nos dice que el problema de Cauchy que nos planteamos tiene una única solución para cada
condición inicial.
SEGUNDO TIPO: VARIABLES SEPARABLES (vs)

Una EDO de primer orden de la forma
dy
 f ( x) g ( y )
dx (*)
Se dice que es de variables separables
Observe que al dividir entre la función g  y  , podemos escribir una ecuación separable
1 dy dy
dy
 f ( x ) g  y  como  f ( x) o h  y   f ( x ) , donde, por conveniencia h  y 
dx g  y  dx dx
1
representa a .
g  y
dy
Ahora sí y    x  representa una solución de la ecuación h  y   f ( x ) , se tiene que
dx
h   x     x   f ( x) , y por tanto  h   x     x  dx   f  x dx . Pero dy     x  dx , por lo que
la ecuación  h   x     x  dx   f  x dx es la misma que
 h  y dy   f  x dx o G( y)  F ( x)  C , donde C  R es una constante arbitraria. G( y) y

F ( x) son antiderivadas de h  y   1 y f  x  , respectivamente.
g  y
Si queremos una formula explicita para las soluciones debemos despejar y en función de x en la
relación anterior. Si no se puede despejar y las soluciones quedan expresadas de manera
implícita o paramétrica

EJEMPLO 1.
dy 6x2
Resuelva la ecuación diferencial 
dx 2 y  cos y
SOLUCIÓN
Al escribir la ecuación en forma diferencial e integrar ambos lados, se tiene
 2 y  cos y  dy  6 x2dx    2 y  cos y  dy   6 x2dx  y 2  seny  2x3  C , donde C es una
constante. La ecuación  seny  2 x3  C da la solución general en forma implícita. En este caso, es
imposible resolver la ecuación para expresar y de forma explícita como una función de x.
En la figura se muestran las gráficas de varios miembros de la familia de soluciones de la ecuación
diferencial del ejemplo 1. Como se ve en las curvas de izquierda a derecha, los valores de C son 3,
2, 1, 0, -1, -2 y -3.
EJEMPLO 2. Analicemos la ecuación y  xy .

SOLUCIÓN
Aquí f ( x)  x y( y)  y . Observemos primero que la función nula y( x)  0 define una

y
dy
hacemos 
y 
solución de la EDO. Si y  0  xds  C . Tenemos entonces que
 x2 
2
x
x2
ln( y )   C , de donde ln( y )  exp   C   ke 2 , donde k  eC  0 . Eliminando el
2  2 
modulo y considerando los posibles valores positivos y negativos vemos que todas las soluciones
x2
son de forma y  ke .2
EJEMPLO 3. Estudiaremos ahora la EDO y  cos ( y ) .

2
SOLUCIÓN

En este caso f ( x)  1 y g ( y)  cos 2 ( y) . Las soluciones constantes son de la forma
dy
y ( x)  2  k k  . Para las demás soluciones,  cos ( y)   dx  C
2
 sec ( y)dy  x  C de aquí, tan( y)  x  C donde C  R

2
,
. Para y  ( 2 , 2 ) esto es
 
y  arctan( x  C) .
Para los demás valores de y debemos tener en cuenta la periodicidad de la función tangente.
Tenemos entonces que todas las soluciones no constantes son de la forma
y  k  arctan( x  C) ,con C  R y k  Z
Ejemplo 4
BRAQUISTOGRONA.
Se denomina así a la forma que debe tener un alambre para que una argolla que se desliza por el
sin roce bajo acción de la gravedad de un punto a otro de menor altura y no en la misma vertical,
lo haga en el menor tiempo posible.
El EDO que describe la forma de la curva es  

y 1   y  k 2 , donde k es una constante
2
positiva.
SOLUCIÓN
1
k y
2 2
Utilizando el método de separación de variables, se tiene y   
 y 

𝑦
De aquí ∫ √2 𝑑𝑦 = ∫ 𝑑𝑥 + 𝐶 , haciendo 𝑦 = 𝑘 2 𝑠𝑒𝑛2 𝜃 ⇒ 𝑑𝑦 = 2𝑘 2 𝑠𝑒𝑛𝜃 𝑐𝑜𝑠𝜃𝑑𝜃 obtenemos
√𝑘 −𝑦
𝑘𝑠𝑒𝑛𝜃2𝑘 2 𝑠𝑒𝑛𝜃 𝑐𝑜𝑠𝜃𝑑𝜃 1−cos(2𝜃)
∫ = 𝑥 + 𝐶 entonces 2𝑘 2 ∫ 𝑠𝑒𝑛2 𝜃𝑑𝜃 = 𝑥 + 𝐶 de aquí 2𝑘 2 ∫ 𝑑𝜃 =
𝑘 𝑐𝑜𝑠𝜃 2
𝜃 𝑠𝑒𝑛2𝜃 𝜃 𝑠𝑒𝑛2𝜃
𝑥 + 𝐶, entonces2𝑘 2 (2 − ) = 𝑥 + 𝐶 ,luego 𝑥 = 2𝑘 2 ( 2 − ) − 𝐶 Por lo tanto, se tiene
4 4
que 𝑥 = 𝑥(𝜃) e 𝑦 = 𝑦(𝜃) con
𝑘2 𝑘2
𝑥= (2𝜃 − 𝑠𝑒𝑛2𝜃) − 𝐶 ; 𝑦= (1 − 𝑐𝑜𝑠2𝜃).
2 2
𝑘2
Si ahora hacemos 𝑤 = 2𝜃 , vemos que 𝑥= (𝑤 − 𝑠𝑒𝑛𝑤) − 𝐶
2
𝑘2
𝑦= (1 − 𝑐𝑜𝑠𝑤), Con 𝑤 ∈ ℝ. La solución es una familia de curvas llamadas cicloides. Ver
2
figura .
Y 0
-0.2
-0.4
-0.6
-0.8
-1
0.5 1 1.5 X
Curva Braquistócrona con 𝑥 ∈ [0, 𝜋2] y ∈ [−1,0], parámetro k =1 y constante C = 0
PÉRDIDA DE UNA SOLUCIÓN

Se debe tener cuidado al separar las variables ya que las variables que sean divisores podrían ser
cero en un punto. Concretamente, si r es una raíz de la función g  y  , entonces sustituyendo
dy
y  r en  f ( x ) g  y  , se encuentra que ambos lados son iguales a cero; es decir, y  r es
dx
una solución constante de la ecuación diferencial. Pero después de que las variables se separan, el
dy
lado izquierdo de  f ( x)dx está indefi nido en r. Por tanto, y  r podría no representar a
g  y
la familia de soluciones que se ha obtenido después de la integración
y simplificación. Recuerde que una solución de este tipo se llama solución singular
Ejemplo 5
PÉRDIDA DE UNA SOLUCIÓN

dy
Analicemos la ecuación  y2  4
dx
SOLUCIÓN
dy  1/ 4 1/ 4 
Poniendo la ecuación en la forma  dx     dy  dx (1)
y 4
2
 y  2 y  2 
La segunda ecuación en la ecuación (1) es el resultado de utilizar fracciones parciales en el lado
izquierdo de la primera ecuación. Integrando y utilizando las leyes de los logaritmos se obtiene
1 1 y2 y2
ln y  2  ln y  2  x  C1 o ln  4 x  C2 o  e4 x C2 . Aquí hemos
4 4 y2 y2
sustituido 4C1 por C 2 . Por último, después de sustituir eC2 por C y despejando y de la última
1  Ce4 x
ecuación, obtenemos una familia uniparamétrica de soluciones y  2 (2)
1  Ce4 x
dy
Ahora, si factorizamos el lado derecho de la ecuación diferencial como   y  2  y  2 
dx
dy_dx _
(y _ 2)(y _ 2), sabemos del análisis de puntos críticos que  y  2 y y  y  2 son dos
soluciones constantes (de equilibrio). La solución y2 es un miembro de la familia de
soluciones definida por la ecuación (2) correspondiendo al valor C  0 . Sin embargo, y  2
es una solución singular; ésta no se puede obtener de la ecuación (2) para cualquier elección del
parámetro C . La última solución se perdió al inicio del proceso de solución. El examen de la
ecuación (1) indica claramente que debemos excluir y  2 en estos pasos.
EJEMPLO 6
Un problema con valores iniciales
 2y
 e  y  cos x  e sen2 x
dy y
Resuelva  dx
 y  0  0

SOLUCIÓN
e2 y  y sen2 x
Dividiendo la ecuación entre e y cos x se obtiene y
dy  dx
e cos x
Antes de integrar se realiza la división del lado izquierdo y utilizamos la identidad trigonométrica
sen 2 x  2 senx cos x en el lado derecho. Entonces integrando   e y  ye y dy  2 senxdx , se
obtiene e y  ye y  e y  2cos x  C
La condición inicial y  0  cuando x  0 implica que C  4 . Por tanto una solución del problema
con valores iniciales es e y  ye y  e y  2cos x  4
Antes de dar el ejemplo 7 veamos algo sobre reacciones químicas.

REACCIONES QUÍMICAS
REACCIONES QUÍMICAS La desintegración de una sustancia radiactiva, caracterizada por la
dX
ecuación diferencial  kX , se dice que es una reacción de primer orden. En química hay
dt
algunas reacciones que siguen esta misma ley empírica: si las moléculas de la sustancia A se
descomponen y forman moléculas más pequeñas, es natural suponer que la rapidez con que se
lleva a cabo esa descomposición es proporcional a la cantidad de la primera sustancia que no ha
experimentado la conversión; esto es, si X  t  es la cantidad de la sustancia A que permanece en
dX
cualquier momento, entonces  kX , donde k es una constante negativa ya que X es
dt
decreciente. Un ejemplo de una reacción química de primer orden es la conversión del cloruro de
terbutilo,  CH3 3 CCI en alcohol t-butílico  CH3 3 COH :
CH3 3 CCI  NaOH  CH3 3 COH  NaCI .
Sólo la concentración del cloruro de terbutilo controla la rapidez de la reacción. Pero en la
reacción CH 3CI  N aOH  CH 3OH  N aCI se consume una molécula de hidróxido de sodio,
N a OH , por cada molécula de cloruro de metilo, CH 3CI , por lo que se forma una molécula de
alcohol metílico, CH 3OH y una molécula de cloruro de sodio, N a CI . En este caso, la razón con
que avanza la reacción es proporcional al producto de las concentraciones de CH 3CI y N a OH
que quedan. Para describir en general esta segunda reacción, supongamos una molécula de una
sustancia A que se combina con una molécula de una sustancia B para formar una molécula de
una sustancia C . Si X denota la cantidad de un químico C formado al tiempo t y si  y  son,
respectivamente, las cantidades de los dos químicos A y B en t  0 (cantidades iniciales),
entonces las cantidades instantáneas no convertidas de A y B al químico C son   X y   X ,
respectivamente. Por lo que la razón de formación de C está dada por
dX
 k   X    X  , donde k es una constante de proporcionalidad. Una
dt
reacción cuyo modelo es la ecuación dX  k   X    X  se dice que es una reacción de
dt
segundo orden.
Suponga que a gramos de una sustancia química A se combinan con b gramos de una sustancia
química B . Si hay M partes de A y N partes de B formadas en el compuesto y X  t  es el
número de gramos de la sustancia química C formada, entonces el número de gramos de la
sustancia química A y el número de gramos de la sustancia química B que quedan al tiempo t
son, respectivamente,
M N
a X y b X
M N M N

La ley de acción de masas establece que cuando no hay ningún cambio de temperatura, la razón
con la que reaccionan las dos sustancias es proporcional al producto de las cantidades de A y de B
que aún no se han transformado al tiempo t:
dX  M  N 
 a  X  b  X ( %)
dt  M  N  M N 
M N
Si se saca el factor del primer factor y del segundo y se introduce una constante
M N M N
de proporcionalidad k  0 , la expresión (%) toma la forma
dX
 k   X    X  (+++)
dt
M N M N
donde   a y  b . Una reacción química gobernada por la ecuación diferencial
M N
no lineal (+++) como sabemos es una reacción de segundo orden.
EJEMPLO 7 Reacción química de segundo orden
Cuando se combinan dos sustancias químicas A y B se forma un compuesto C. La reacción

resultante entre las dos sustancias químicas es tal que por cada gramo de A se usan 4 gramos de B.
Se observa que a los 10 minutos se han formado 30 gramos del producto C. Determine la cantidad
de C en el tiempo t si la razón de la reacción es proporcional a las cantidades de A y B que quedan
y si inicialmente hay 50 gramos de A y 32 gramos de B. ¿Qué cantidad de compuesto C hay a los 15
minutos? Interprete la solución cuando t   .
SOLUCIÓN
Sea X  t  la cantidad de gramos del compuesto C presentes en el tiempo t. Es obvio que
X  0  0g y X 10  30 g .
Si, por ejemplo, hay 2 gramos del producto C, hemos debido usar, digamos, a gramos de A y b
2 1
gramos de B, así a + b = 2 y b= 4a. Por tanto, debemos usar a   2   de la sustancia química
5 5
8 4 1
A y b   2   g de B. En general, para obtener X gramos de C debemos usar X gramos de
5 5 5
4
A y X gramos de B.
5
1 4
Entonces las cantidades de A y B que quedan al tiempo t son 50  X y 32  X
5 5
respectivamente. Sabemos que la razón con la que se forma el compuesto C satisface que
dX  1  4 
  50  X  32  X  . Para simplificar las operaciones algebraicas subsecuentes, factorizamos
dt  5  5 
1 4
del primer término y del segundo y después introduciremos la constante de
5 5

dX
proporcionalidad:  k  250  X  40  X  . Separamos variables y por fracciones parciales
dt
1 1
podemos escribir que 210 dX  210 dX  kdt , Integrando se obtiene
250  X 40  X
250  X 250  X
ln  210kt  C1 o  C2 e 210 kt
40  X 40  X
25
Cuando t  0 , X  0 , se tiene que en este punto C2  . Usando X 10  30 g en t = 10
4
1 88
encontramos que 210k  ln  0.1258 . Con esta información se despeja X de la última
10 25
1  e0.1258t
ecuación: X  t   1000
25  4e0.1258t
En la fi gura se presenta el comportamiento de X como una función del tiempo.
Es claro de la tabla adjunta y de la ecuación anterior que X  40 conforme t   . Esto significa

que se forman 40 gramos del compuesto C, quedando
1 4 4
50   40   42 g de A y 32  X A y 32   40   0 g de B
5 5 5
TERCER TIPO
LA ECUACION DIFERENCIAL LINEAL DE PRIMER ORDEN
Una ecuación diferencial de primer orden de la forma
y t   a t  y t   b t 
Se dice que es una ecuación lineal en la variable dependiente y.
Resolveremos la ecuación y t   a t  y t   b t  en cuatro partes:

i) y t   a y t   0 , donde a es una constante.
ii) y t   a y t   b , donde a y b es una constante diferentes de cero.
iii) y t   a y t   b t  , donde b(t) es una función continua en 0, .
iv) y  t   a t  y t   b t  , donde a(t) y b(t)son funciones continuas en 0, .
(i) y t   ay t   0, a 
Supongamos que y es solución de (i). Luego y(t )  ay(t )  0 . Por consiguiente

e ( y(t )  ay (t ))  0  eat y(t )  aeat y (t )  (eat y (t ))  0 existe una constante C tal que
at
e at y (t )  C , de aquí y (t )  Ce at como C es un constante arbitraria es dable aseverar que

 at
las soluciones de (i) son infinitas y todas ellas están comprendidas en y (t )  Ce , esta
y (t )  Ce at es llamada solución general de (i)
(Recordar eat ( y(t )  ay (t )  (e at y (t )) ). Ahora bien, frecuentemente una ecuación
diferencial y’(t) + ay(t) = 0 viene acompañada de una condición adicional y (t0 )   ( número
real arbitrario) lo que da lugar a un problema de valor inicial (PVI).
PROBLEMA DE VALOR INICIAL (PVI)

 y(t )  ay  t   0


 y (t0 )  

Veremos que el PVI tiene una única, solución. Para ello basta evaluar en 𝑡0 a la función
y (t )  Ce at (solución general de la ecuación diferencial).
En efecto   y  t0   Ceat  C   eat  y t    eat t  , solución única de ecuación

0 0 0
diferencial en estudio
EJEMPLO 1
Resolver y  3 y  0
Solución

Multiplicando por e3x obtenemos e3x y  e3x y  0   e3x y   0 , entonces existe una constante C
talque e3 x y  C  y  Ce3 x
EJEMPLO 2
Resolver y t   ky t  , done k es una constante positiva; y  0  12
Solución
Multiplicando por e kt obtenemos ekt y  ekt ky  0   ekt y   0 , entonces existe una
constante C talque ekt y  C  y t   Cekt .Ahora para t  0 , obtenemos
12  y  0  Cek .0  C  12 .Luego y t   12ekt
APLICACIONES
Ensayo sobre el principio de la población – Thomas Malthus (1798)

Publicada en el siglo XVIII, es probablemente el ensayo sobre demografía más conocido del
mundo. Elaborado por el economista inglés Thomas Maltus, en la que desarrollaba un modelo
matemático en el que hacia una previsión del crecimiento de la población (crecimiento
exponencial según su ensayo) y como se regulaba este crecimiento debido a la futura escasez de
alimentos (hambrunas) así como el estallido de guerras y epidemias.
Todo y que la mayoría de sus predicciones nunca llegaron a cumplirse, su obra aparece como una
de las más influyentes, sobretodo en Inglaterra, donde aprovechándose de sus teorías se
degradaron mucho las condiciones de vida de las ya desfavorecidas clases bajas y obreras inglesas.
También influyó a los biólogos Darwin y Wallace o, más recientemente el economista John
Maynard Keynes.
Aún hoy en día sus teorías son objeto de discusión y controversia.
La curiosidad: Malthus era pastor anglicano y provenía de una familia de clase alta. Gracias a su
padre mantuvo relación directa con destacados filósofos de la época como David Hume y Jean-
Jaques Rousseau.

 y   t   ay  t   0

EL PVI
 y  t0    COMO MODELOS MATEMÁTICOS
INTRODUCCIÓN
En esta sección introduciremos la idea de una ecuación diferencial como un modelo matemático y
analizaremos algunos modelos específicos en biología, química y física
Con frecuencia es deseable describir en términos matemáticos el comportamiento de algunos
sistemas o fenómenos de la vida real, sean físicos, sociológicos o hasta económicos. La descripción
matemática de un sistema de fenómenos se llama modelo matemático y se construye con ciertos
objetivos. Por ejemplo, podemos desear entender los mecanismos de cierto ecosistema al estudiar
el crecimiento de la población animal en ese sistema, o podemos desear datar fósiles y analizar el
decaimiento de una sustancia radiactiva ya sea en el fósil o en el estrato en que éste fue
descubierto.
1) DINÁMICA POBLACIONAL
Uno de los primeros intentos para modelar el crecimiento de la población humana por medio de
las matemáticas fue realizado en 1798 por el economista inglés Thomas Malthus. Básicamente la
idea detrás del modelo de Malthus es la suposición de que la razón con la que la población de un
país en un cierto tiempo es proporcional* a la población total del país en ese tiempo. En otras
palabras, entre más personas estén presentes al tiempo t, habrá más en el futuro. En términos
matemáticos, si P(t) denota la población al tiempo t, entonces esta suposición se puede expresar
como
dP
 kP (1)
dt
Donde k es una constante de proporcionalidad. Este modelo simple, falla si se consideran muchos
otros factores que pueden influir en el crecimiento o decrecimiento (por ejemplo, inmigración y
emigración), resultó, sin embargo, bastante exacto en predecir la población de los Estados Unidos,
durante 1790-1860. Las poblaciones que crecen con una razón descrita por la ecuación (1) son
raras; sin embargo, (1) aún se usa para modelar el crecimiento de pequeñas poblaciones en
intervalos de tiempo cortos (por ejemplo, crecimiento de bacterias en una caja de Petri).
Supóngase que el numero P(t) de una población en el tiempo t sigue una ley de crecimiento
exponencial, entonces la razón de cambio de una población es proporcional a la población
 dP
  kP
existente en dicho instante  dt
 P  0   P0


Esta ecuación lineal tiene como solución P t   P0ekt , donde P 0 es la población inicial Por
lo tanto, se concluye que el crecimiento poblacional continuara indefinidamente, es decir:
lím P(t )  
t 
Con crecimiento exponencial una población crecería en forma indefinida al avanzar el tiempo.
Sin embargo, en realidad, cuando la población crece lo suficiente existen factores del medio
ambiente que reducen la tasa de crecimiento. Algunos ejemplos son disponibilidad de
dP
alimentos, depredadores, hacinamiento, etc. Estos factores ocasionan que al final disminuya
dt
dS
El modelo de la ecuación (1) para crecimiento también se puede ver como la ecuación  rS ,
dt
que describe el crecimiento del capital S cuando está a una tasa anual de interés r compuesto
continuamente
2) DECAIMIENTO RADIACTIVO
El núcleo de un átomo está formado por combinaciones de protones y neutrones. Muchas de esas
combinaciones son inestables, esto es, los átomos se desintegran o se convierten en átomos de
otras sustancias. Se dice que estos núcleos son radiactivos. Por ejemplo, con el tiempo, el radio Ra
226, intensamente radiactivo, se transforma en el radiactivo gas radón, Rn-222.
dN
Para modelar el fenómeno del decaimiento radiactivo, se supone que la razón con la que los
dt
núcleos de una sustancia se desintegran es proporcional a la cantidad (más precisamente, el
dN
número de núcleos), N  t  de la sustancia que queda al tiempo t; esto es  kN …(2)
dt
Por supuesto que las ecuaciones (1) y (2) son exactamente iguales; la diferencia radica sólo en la
interpretación de los símbolos y de las constantes de proporcionalidad. En el caso del crecimiento,
como esperamos en la ecuación (l), k > 0, y para la desintegración como en la ecuación (2), k < 0.
El modelo de desintegración de la ecuación (2) también se aplica a sistemas biológicos tales como
la determinación de la “vida media” de un medicamento, es decir, el tiempo que le toma a 50% del
medicamento ser eliminado del cuerpo por excreción o metabolización. En química el modelo del
decaimiento, ecuación (2), se presenta en la descripción matemática de una reacción química de
primer orden. Lo importante aquí es: Una sola ecuación diferencial puede servir como modelo
matemático de muchos fenómenos distintos. Con frecuencia, los modelos matemáticos se
acompañan de condiciones que los definen. Por ejemplo, en las ecuaciones (l) y (2) esperaríamos
conocer una población inicial P0 y por otra parte la cantidad inicial de sustancia radioactiva N 0 . Si
el tiempo inicial se toma en t  0 , sabemos que P  0  P0 y que N  0  N0 . En otras palabras, un
modelo matemático puede consistir en un problema con valores iniciales o, como veremos más
adelante , en un problema con valores en la frontera.

 dN
  kN
SOLUCION AL PVI:  dt
 N  0  N0

Muchas sustancias ( o materiales ) radioactivos se desintegran a una velocidad que es
proporcional a la cantidad de sustancias ( o material ) presente en dicho instante. En el PVI
planteado N  t  denota la cantidad de sustancia o elemento radiactivo presente en el instante “t”
y k es una constante positiva.
La ecuación N   t   kN t  o equivalentemente N  t   kN t   0 . para resolver multiplicamos
a ambos miembros por 𝑒 𝑘𝑡 obtenemos N'(t) 𝑒 𝑘𝑡 + k N(t) 𝑒 𝑘𝑡 = 0 ⟺ (N(t) 𝑒 𝑘𝑡 )' = 0 entonces
existe una constante C tal que N(t) 𝑒 𝑘𝑡 = C ⟺ N(t) = C e-kt , evaluando en t = 0 obtenemos C =
N(0) remplazando este valor obtenemos N  t   N0ekt
En donde N 0 y k son constante positivas.
 N 0 Representa la cantidad del elemento que está presente al tiempo t = 0 y se le denomina
cantidad inicial.
 La constante k depende del elemento particular implicado y se llama constante de
decrecimiento ó decaimiento.
 El concepto de semi-vida (o vida media), que definiremos a continuación, reviste particular
importancia; denotamos con  el tiempo que debe transcurrir para que la cantidad inicial N(0) se
In 2
= N(0) = N(0)e-k luego – ln2 = -k , es decir  
N(0) N(0)
reduzca a esto es .
2 2 k
Observamos que, conocida como la semi-vida (o vida media) de la sustancia(o elemento) en
estudios , es independiente de la cantidad inicial N(0).
Ejemplo 1
Si el 45% de una sustancia radioactiva se desintegra en 200 años .¿Cuál es su vida media? y ¿En
cuánto tiempo se desintegrará 60% de la cantidad original?.
Solución
a) Si el 45% se desintegra, permanece sin desintegrar el 55% luego:

55 11 11
N(0) = N(0)e-200k es decir = e-200k luego ln = -200k
100 20 20
1 ln 2 200 ln 2
 k= (ln20-ln11). Por tanto  = así  = .
200 k ln 20  ln 11
b) Si el 60% se desintegra permanece sin desintegrar 40%.
40 2 2 l 2
Luego: N(0) = N(0)e-kt   e-kt  ln     kt de ahí que t = - ln   ,
100 5 5 k 5
1
en donde k = (ln20 – ln11).
200
Ejemplo 2
Si después de 50 días se tiene el 60% de una sustancia radioactiva, determina la constante de
decrecimiento y la vida media de la sustancia.
Solución
60 ln( 0,6)
N(0) = N(0)e-50k  In(0,6) = -50k luego k = - y
100 50
In(2)
=  67,82 días.
k
Ejemplo3 Un isótopo radioactivo tiene una vida media de 16 días. Ud desea obtener 30g al final de
30 días. ¿Cuál es la cantidad inicial del isótopo necesaria?
Solución
Como el tiempo de vida media está en días trabajaremos en dicha unidad de tiempo. Sea N(t) la
cantidad presente al momento t y N(0) la cantidad que andamos buscando Sabemos que:
𝑵(𝒕) = 𝑵(𝟎)𝒆−𝒌𝒕 , donde k es una constante. Usando el tiempo de vida media  para determinar
k, se tiene que: k  ln 2 
1
16
30
ln  2 
Por consiguiente, N 30  N (0)e 30 K
 30  N (0)  30e 30 K
 30e 16
 110.04 g
Ejemplo 4
Suponer que 10 gramos del isótopo 239Pu se liberaron en el accidente nuclear de Chernobyl.
¿Cuánto tiempo tomará a los 10 gramos disminuir a 1 gramo?
Solución
Considerar que N representa la masa (en gramos) del plutonio. Dado que la tasa de desintegración
es proporcional a y, se sabe que 𝑵(𝒕) = 𝑵(𝟎)𝒆−𝒌𝒕 donde t es el tiempo en años. Para encontrar
los valores de las constantes N(0) y k, aplicar las condiciones iniciales. Con base en que N(0) = 10
cuando t = 0, lo cual implica que C = 10. Luego, con base en el hecho de que la vida media de
ln 2 ln 2
239Pu es de 24 100 años se puede tener 24100  k  0.000028761
k 24100

Así, el modelo es N t   10e0.000028761t
Para encontrar el tiempo en que 10 gramos decrecen a 1 gramo, se puede despejar para t en la
0.000028761t
ecuación 1  10e  t  80059
La solución es aproximadamente 80 059 años.
NOTA DE ENRIQUECIMIENTO (Carbono 14(C14))

METODO DEL CARBONO 14
La técnica llamada del C14, para datar un objeto consiste en medir la cantidad de C-14 que queda
en la actualidad en dicho objeto, y utilizar la forma de las soluciones de la ecuación de
decaimiento radiactivo para calcular el tiempo que ha pasado.
Por ejemplo, la técnica de C-14 se utilizó en el año 1988 para estimar la edad del Sudario de Turin,
tela de lino hallada en 1356 que muestra la imagen de un hombre que presenta marcas y
traumas fisicos (ver la Figura), y de la que se pensaba que podría ser la tela que cubria a Jesus de
Nazaret en el sepulcro, llamada también Sabana Santa. Se observó que las fibras del tejido
contenían entre un 92% y un 93% del nivel inicial de C-14
Se llego a un resultado que indicaba que el Sudario fue fabricado entre 689 y 599 años antes del
momento en que fueron realizadas las pruebas, en el año de 1988. Es decir, mucho despues de la
epoca en que vivio Jesus. Lo que probo que no podria ser la Sabana Santa .
Este método se debe al quimico Willar Libby cuyo descubrimiento le valio el premio Nobel de
Quimica en 1960. La teoría se basa en lo siguiente. La atmosfera terrestre es continuamente
bombardiada por rayos cósmicos, los cuales producen neutrones libres que se conbinan con el
nitrógeno de la atmosfera para producir el isotopo C-14 (Carbono 14 bien radiocarbono). Este C-14
se combina con el bióxido de carbono presente en la atmosfera, el cual es absorbido por las
plantas y estas a su vez son alimento para los animales. Asi es como se incorpora el radiocarbono a
los tejidos de seres vivos.

El cociente de la cantidad de carbono 14 y la cantidad de carbono ordinario presentes en la
atmosfera es constantes, y en consecuencia la proporción de isotopo presente en todos los
organismos vivos es la misma que en la atmosfera. Cuando un organismo muere, la velocidad de
incorporación de radiocarbono a el se hace nula y entonces comienza el proceso de desintegración
radioactiva del C-14, que se encontraba presente en el momento de su muerte. Asi comparando la
proporción de C-14 que hay en un fosil con la proporción constante encontrada en la atmosfera es
posible obtener una estimación razonable de su edad.
Ejemplo 1
1
Se ha encontrado que un hueso antiguo contiene de la cantidad original de C-14 de un hueso al
8
tiempo actual. ¿Cuál es la antigüedad del fosil?.
Solución
Sea N(t) la cantidad presente de C-14 en el hueso t años después que murió el organismo. En este
caso, N  t   N  0 ekt , en donde N(0) es la cantidad de C-14 cuando t = 0
La vida media del C-14 es de 5568 años, es decir   5568 años , entonces
0, 69315
k  0, 00012448 , asi N  t   N  0 e0,00012448t
5569
N  0
Buscando el valor de t para el cual N  t   , tenemos que
8
N  0 1
N  0  e0,00012448t   0,00012448t  ln  t  16705 . Asi, el fosil tiene una antigüedad de
8 8
16705 años.
Ejemplo 2
En 1950 se hicieron excavaciones en Nipur (Babilonia), en las cuales se encontraron muestras de

carbón que reportaron 4.09 desintegraciones por minuto y por gramo. Una muestra actual
reporto 6.68 desintegraciones por minuto y por gramo. Se sabe que la primer muestra se formo
en la época del reinado de Hammurabi. Con estos datos, determine hace cuanto tiempo
Hammurabi reino en Babilonia.

Solución
dN
Sea N(t) la cantidad presente de C-14 en el tiempo t. entonces , esla velocidad de
dt
dN dN
desintegración del C-14 al tiempo t y  4.09 ;  6.68 …(1)
dt dt
dN
Sabemos que   kN  t  , además la vida media del C-14 es de 5568 años, es decir
dt
 0, 00012448 , asi N  t   N  0 e0,00012448t … (2)

0, 69315
  5568 años , entonces k 
5569
dN dN
Sustituyendo (2) en   kN  t  se tiene   kN  0  e 0,00012448 t ….(3)
dt dt
dN
Considerando (1) en (3) se tiene  6.68e 0,00012448t
dt
Ahora bien para determinar hace cuanto tiempo reino Hammurabi en Babilonia, tendremos que
dN
calcular t para el cual se cumple que  4.09 , esto es 4.09  6.68e0,00012448t , entonces
dt
4.09
ln  0.00012448t  t  3940.9786.
6.68
Aproximadamente 3941 años que Hammurabi reino en Babilonia.

DOSIFICACIÓN DE MEDICAMENTOS
Determinar y recetar dosis de fármacos son aspectos sumamente importantes de la profesión
médica. Con frecuencia se debe tener precaución debido al lado posible adverso o a los efectos
tóxicos de las medicinas (o drogas).
Muchos medicamentos son utilizados por el cuerpo humano de manera que la cantidad presente
sigue una ley exponencial de disminución. Es decir, si N(t) es la cantidad de fármaco presente en el
cuerpo al tiempo t, entonces N'(t) = - k N(t) o equivalentemente N'(t) + k N(t) = 0 . para
resolver multiplicamos a ambos miembros por 𝑒 obtenemos N'(t) 𝑒 𝑘𝑡 + k N(t) 𝑒 𝑘𝑡 = 0 ⟺ (N(t)
𝑘𝑡
𝑒 𝑘𝑡 )' = 0 entonces existe una constante C tal que N(t) 𝑒 𝑘𝑡 = C ⟺ N(t) = C e-kt , evaluando en t = 0
obtenemos C = N(0) remplazando este valor obtenemos 𝑵(𝒕) = 𝑵(𝟎)𝒆−𝒌𝒕
En donde N(0) y k son constante positivas .
 N(0) representa la cantidad del elemento que esta presente al tiempo t = 0 y se le denomina
cantidad inicial.
 La constante k depende del elemento particular implicado y se llama constante de
decrecimiento ó decaimiento.
Si H es la semivida de tal medicamento, entonces 𝐻 = (𝑙𝑛2)/𝑘
Supóngase que se desea anal izar el caso en que se administren a un paciente dosis iguales de un
fármaco como ese, cada I unidades de tiempo, hasta que se logre un cierto nivel terapéutico.
La razón de administrar dosis reducidas de mantenimiento se relaciona con frecuencia con los
efectos tóxicos de los fármacos. En particular, supóngase que existen d dosis de P unidades cada
una, que se aplican dosis en los tiempos t = 0, I, 2I,…, y (d - 1)I, y que el nivel terapéutico T, se
alcanza en t = dI, lo cual se presenta un intervalo después de que se administra la ultima dosis. Se
verá ahora como determinar una formula que dé el nivel terapéutico.
En el tiempo t = 0, el paciente recibe las primeras P unidades, de manera que la cantidad de

medicamento en el cuerpo es P. al tiempo t = I, la cantidad presente que proviene de la primera

dosis es [de la ecuación (1)] 𝑃𝑒 −𝑘𝐼 . Además, a t = I, se administran las segundas P unidades. Por
ello, la cantidad total de fármaco presente es 𝑃 + 𝑃𝑒 −𝑘𝐼 .
Al tiempo t = 2 I, la cantidad que permanece, y que proviene de la primera dosis es 𝑃𝑒 −2𝑘𝐼 , de la

segunda dosis, que ha estado en el sistema durante sólo un intervalo, la cantidad presente es
𝑃𝑒 −𝑘𝐼 . También, al mismo tiempo t = 2 I se administrar la tercera dosis de P unidades, de manera
que la cantidad total de fármaco presente es 𝑃 + 𝑃𝑒 −𝑘𝐼 + 𝑃𝑒 −2𝑘𝐼
Continuando de esta manera, la cantidad de fármaco presente en el sistema al tiempo dI, un

intervalo de tiempo después de que se administra la ultima dosis, esta dad por 𝑇 = 𝑃𝑒 −𝑘𝐼 +
𝑃𝑒 −2𝑘𝐼 + ⋯ + 𝑃𝑒 −𝑑𝑘𝐼 . (2)
Se puede expresar el lado derecho de la Ec. (2) en forma distinta. En primer lugar, se multiplican
ambos lados de la (2) por 𝑒 −𝑘𝐼 :
𝑒 −𝑘𝐼 𝑇 = 𝑒 −𝑘𝐼 (𝑃𝑒 −𝑘𝐼 + 𝑃𝑒 −2𝑘𝐼 + … + 𝑃𝑒 −𝑑𝑘𝐼 )  𝑒 −𝑘𝐼 𝑇 = 𝑃𝑒 −2𝑘𝐼 + 𝑃𝑒 −3𝑘𝐼 + … +

𝑃𝑒 −(𝑑+𝐼)𝑘𝐼 ) (3)
Restando los resultados de la Ec. (3) de los correspondientes de la Ec.(2), se tiene
𝑇 − 𝑒 −𝑘𝐼 𝑇 = 𝑃𝑒 −𝑘𝐼 − 𝑃𝑒 −(𝑑+𝐼)𝑘𝐼 Simplificando y despejando T, (1 − 𝑒 −𝑘𝐼 )𝑇 = 𝑃𝑒 −𝑘𝐼 (1 −

𝑃𝑒 −𝑘𝐼 (1−𝑒 −𝑑𝑘𝐼 ) 𝑃(1−𝑒 −𝑑𝑘𝐼 )
𝑒 −𝑑𝑘𝐼 )  𝑇 = 𝑇=
1 − 𝑒 −𝑘𝐼 𝑒 −𝑘𝐼 (1 − 𝑒 −𝑘𝐼 )
𝑃(1−𝑒 −𝑑𝑘𝐼 )
 𝑇= (4)
𝑒 −𝑘𝐼 −1
La ecuación (4) permite determinar el nivel terapéutico T en términos de la dosis P, los intervalos I,
el número de dosis d, y la semivida H, del medicamento [puesto que 𝑘 = (𝑙𝑛2)/𝐻]. Entre oras
posibilidades, puede determinarse la dosis P si se conocen T, H, I y d.
El objetivo ahora es mantener el nivel terapéutico en el sistema del paciente. Para lograr esto, se
administra una dosis reducida R, a los tiempos 𝑡 = 𝑑𝐼, (𝑑 + 2)𝐼, asi sucesivamente. En el tiempo
𝑡 = (𝑑 + 1)𝐼, pero antes de administrar la segunda dosis reducida, la cantidad de fármaco en el
sistema, proviene de la primera dosis reducida es 𝑅𝑒 −𝑘𝐼 . Y la cantidad que permanece, proviene
del nivel terapéutico, es 𝑇𝑒 −𝑘𝐼 Supóngase que se requiere que la suma de esas cantidades esas el
nivel terapéutico, T. es decir, 𝑇 = 𝑅𝑒 −𝑘𝐼 + 𝑇𝑒 −𝑘𝐼
Despejando R, obtenemos 𝑅𝑒 −𝑘𝐼 = 𝑇 − 𝑒 −𝑘𝐼  𝑅 = 𝑇(1 − 𝑒 −𝑘𝐼 )𝑒 𝑘𝐼 .
𝑃𝑒 −𝑘𝐼 (1−𝑒 −𝑑𝑘𝐼 )

De aqui se obtiene 𝑅= (1 − 𝑒 −𝑑𝑘𝐼 )𝑒 𝑘𝐼 ,
1−𝑒 −𝑑𝑘𝐼
O en términos más simples, 𝑅 = 𝑃(1 − 𝑒 −𝑑𝑘𝐼 ). (5)
Continuando con la dosis reducida a intervalos I se asegura que le nivel de fármacos en el sistema
nunca caiga por debajo de T. además, se debe observar que como – 𝑑𝑘𝐼 < 0, entonces 0 <

𝑒 −𝑑𝑘𝐼 < 1. En consecuencia, el factor 1 − 𝑒 −𝑑𝑘𝐼 de la Ec. (5) se encuentra entre 0 y 1. Esto asegura
que R sea menor que P, en donde R es en realidad una dosis reducida.
Los especialistas afirman que “la cantidad terapéutica T debe ser elegida de entre una gama de
valores determinados en forma empírica. Se requieren discreción y experiencias médicas para
seleccionar los intervalos y las duraciones o tiempos apropiados para administrar el fármaco.
Incluso es posible que varíe la semi-vida de un fármaco entre pacientes distintos. Su distribución
en el sistema, las interacciones entre los medicamentos, la edad de los pacientes, su salud en
general y la salud de órganos vitales, como el hígado y los riñones.
Ejemplo.
La teofilina es un fármaco que se utiliza para tratar asma bronquial y tiene una semi-vida de 8
horas en el sistema de un paciente relativamente saludable, que no fuma. Supóngase que un
paciente como este logra el nivel terapéutico deseado de este fármaco en 12 horas. Aquí, d=
3.Debido a su toxicidad, debe reducirse la dosis después. Al miligramo más próximo, determine (a)
el nivel terapéutico y (b) la dosis reducida.
Solución
P(1  e  dkl )
a) T = Aquí P = 100mg ,  =8h , I = 4h
e kl  1
ln 2 ln 2 ln 2
Pero  =  k , entonces kI =  4h
k 8h 8h
ln 2 3
3 
ln 2 100mg(1  e 2
) 100mg(1  2 2 )
Es decir kI = luego T = ln2
 1
2 2 2
1
e 2
1
3
1 2 1
100(1    ) 100(1  )
 
2 8
T=   T = 156 mg
2 1 2 1

3
 ln 2 1
2
b) R = (1 – e-dkl)  R = 100 mg(1 - e )  100mg(1  )
8
2 4 2
 R = 100 mg (1  ) = 100mg    R = 65 mg
4  4 
 
(ii) y t   ay t   b; a  b  ; a, b  0
La experiencia alcanzado en (i) sugiere suponer que y es solución de (ii); luego
y(t )  ay(t )  b .Para todo t  0 , multiplicamos ambos lados por e at , se obtiene

b
eat ( y(t )  ay(t ))  beat  (e y (t ))  ( e ) , teniendo en mente el hacho de que si la
at at
a
derivada de dos funciones son iguales entonces ambas funciones deben diferir por una constante.
b at b
Luego existe una constante C tal que e y (t ) 
at
e  C ,de aquí y (t )   Ce  at como
a a
C es un constante arbitraria , entonces es dable afirmar que las soluciones de (ii) son infinitas y
b
todas ellas están comprendidas en y (t )   Ce  at
a
at
Por el rol que juega la función e en el proceso arriba descrito, se acostumbra decir que es el
factor integrante de la ecuación diferencial. Ahora bien, frecuentemente una ecuación diferencial
viene acompañada de una condición adicional y (t0 )   ( número real arbitrario) lo que da
lugar a un problema de valor inicial (PVI).
EL PROBLEMA DE VALOR INICIAL (PVI)
 y  ay  b

 y(t0 )   (*)

Veremos que el PVI (*) tiene una única , solución; para ello basta evaluar en 𝑡0 a la función
b
y (t )   Ce  at (solución general de la ecuación diferencial), en efecto
a
b b
    Ce
 at0
  y (t0 )   Ce  at 0
a a
b b b
De aquí C  (  )e
at0
de tal manera que y (t )   (  )e  a ( t t0 ) es la única solución del
a a a
PVI
Como ejemplos aplicativos veamos los siguientes modelos
 y   t   ay  t   b

EL PVI
1) LEY DE ENFRIAMIENTO/CALENTAMIENTO DE NEWTON

De acuerdo con la ley empírica de Newton de enfriamiento/calentamiento, la rapidez con la que
cambia la temperatura de un cuerpo es proporcional a la diferencia entre la temperatura del
cuerpo y la del medio que lo rodea, que se llama temperatura ambiente. Si T  t  representa la
dT
temperatura del cuerpo al tiempo t, Tm es la temperatura del medio que lo rodea y es la
dt
rapidez con que cambia la temperatura del cuerpo, entonces la ley de Newton de
enfriamiento/calentamiento traducida en una expresión matemática es
dT
 k T  Tm  (3)
dt
Donde k es una constante de proporcionalidad en ambos casos, enfriamiento o calentamiento.
En el caso de enfriamiento La ley de enfriamiento de Newton establece que la razón a que la
temperatura T (t ) cambia en un cuerpo que se enfría es proporcional a la diferencia entre la
temperatura en el cuerpo y la temperatura constante Tm del medio circundante; es decir,
T  t   k T t   Tm  , en este caso k es una constante positiva.

Sean T (0)  T0 la temperatura inicial, entonces T   t   k T t   Tm  con la condición inicial
(PVI) se escribe de la forma
T   KT  kTm
T (0)  T0 Es un PVI de primer orden
Para resolver el PVI multiplicaremos a toda la ecuación por 𝑒 𝑘𝑡
T ekt  kTekt  kTm ekt de aquí obtenemos (Te kt )  kTm e kt ,entonces
Integrando obtenemos Tekt  Tm ekt  C por lo tanto T (t )  Tm  Ce  kt

De donde, evaluando en t = 0, se obtiene C  T0  Tm . Con esto,
T (t )  Tm  (T0  Tm )e  kt
Ejemplo 1 Asesinato
Se encontró asesinado en su hogar a un rico industrial. La policía llego al lugar del crimen a las
11:00 P.M. en ese momento la temperatura del cuerpos era 31°C, y una hora después era de 30°C.
La temperatura del cuarto en el que se encontró el cuerpo era de 22°C. Determinar la hora en que
ocurrió el asesinato.
Solución
Sea t el número de horas después de las cuales se descubrió el cuerpo, y T(t) la temperatura (en
grados Celsius) del cuerpo en el tiempo t. Se desea calcular el valor de t para el cual T = 37 (la
temperatura normal de cuerpo). Por supuesto, este valor de t será negativo. De acuerdo a la ley de
𝑑𝑇
Newton del enfriamiento, = 𝑘(𝑇 − 22),
𝑑𝑡
En donde k es una constante( negativa ) y la (temperatura ambiente) es 22. En consecuencia

𝑑𝑇
= 𝑘(𝑇 − 22), Separando variables, se tiene
𝑑𝑡
𝑑𝑇 𝑑𝑇
= 𝑘 𝑑𝑡, → ∫ 𝑇−22 = ∫ 𝑘 𝑑𝑡, → 𝑙𝑛|𝑇 − 22| = 𝑘𝑡 + 𝐶.
𝑇−22
Como T -22 > 0, → ln(𝑇 − 22) = 𝑘𝑡 + 𝐶.
Cuando t = 0, T = 31, por lo tanto, ln(31 − 22) = 𝑘 × 0 + 𝐶, De donde

𝑇−22
ln(𝑇 − 22) = 𝑘𝑡 + ln 9,→ ln(𝑇 − 22) − ln 9 = 𝑘𝑡,→ 𝑙𝑛 = 𝑘𝑡.
9
Cuando t = 1, entonces T = 30, por lo que

30−22 8 𝑇−22
𝑙𝑛 = 𝑘 × 1.→ 𝑘 = 𝑙𝑛 9 ≈ −0.11778. Así, → 𝑘 = 𝑙𝑛 ≈ −0.11778𝑡.
9 9
37−22
Ahora, se evalúa t cuando T = 37: 𝑘 = 𝑙𝑛 ≈ −0.11778𝑡.→
9
15
ln( ) 0.51083
9
𝑡 ≈ − 0.11778 ≈ − 0.11778,→ 𝑡 ≈ −4.34.
Puesto que 4.34 horas es (aproximadamente) 4 horas 20minutos, el industrial fue asesinado
aproximadamente a las 6:40 P.M.
Ejemplo 2 Pastel
Cuando un pastel se retira del horno, su temperatura es de. Tres minutos después su temperatura
es de 200 F . Determine la temperatura del pastel en cualquier instante después que se ha sacado
del horno si la temperatura ambiente es de 70 F .
SOLUCION
Aquí T (0)  300 F  T0 y T (3)  200 F
 kt  kt
Sabemos que T (t )  Tm  (T0  Tm )e  T (t )  70 F  (300 F  70 F )e
 kt
Es decir T (t )  70 F  230 Fe
130
200 F  T (3)  70 F  230 Fe 3k   e 3k , de aquí
230
13 1 13
ln( )  3k  k   ln( )  k  0.1902
23 3 23
0.1902 t
Luego T (t )  70 F  230 Fe
MODELAMIENTO DE UN FENOMENO EN BIOLOGIA

ÓSMOSIS: es el pasaje o difusión de un solvente (agua) a través de una membrana
semipermeable mediante un gradiente de concentración. La membrana plasmática permite el
paso del agua de un sitio a otro pero no el de sustancias disueltas en ella (solutos). Toda vez que
la célula tenga en su interior una concentración de solutos mayor que la del medio externo, la
célula está en una solución hipotónica. Por lo tanto, el agua ingresa a la célula y provoca que se
agrande. Por el contrario, si la concentración de solutos es mayor en su ambiente externo la
célula está en un medio hipertónico, hecho que provoca la salida de agua intracelular y la
crenación o arrugamiento de la célula. Cuando la concentración de solutos es igual a ambos
lados de la membrana, la célula está en un medio isotónico (igual tonicidad) y no hay difusión de
agua. En la difusión simple, en la facilitada y en la ósmosis no hay gasto de energía.

DIÁLISIS: cuando una membrana separa una sustancia con diferente
concentración a ambos lados, el soluto (la sal) difunde desde el lugar de mayor
concentración al de menor concentración, mientras que el agua lo hace desde el
sitio donde está en mayor cantidad (solución diluida) hacia la de menor cantidad
(solución concentrada de sal). Este proceso, denominado diálisis, se define como
el pasaje de una sustancia disuelta a través de una membrana semipermeable a
favor de un gradiente de concentración y sin gasto de energía.
LEY DE OSMOSIS
Analicemos el fenómeno de la osmosis, presente en muchos procesos fisiológicos. Consideremos
un experimento en que disponemos dos medios de salmuera A y B separados por una membrana
impermeable en t = 0 con concentraciones iniciales 𝐴(0) y B(0) con A(0) < B(0) .
En un instante t > 0 la membrana que los separa se vuelve semipermeable y permite el paso de
las moléculas de agua, pero no de las moléculas de sal disueltas (ver figura).
A(0) B(0) 𝑨(𝒕) 𝑩(𝒕)
t=0 t>0
FIGURA 1 . Osmosis por una membrana
semipermeable
El problema es modelar la evolución de las concentraciones de sal A(𝑡) y B(𝑡) en función del
tiempo. Se observa experimentalmente que a medida que el tiempo transcurre, el agua se
desplaza a través de la membrana desde la solución de baja concentración A hacia la de alta
concentración B hasta alcanzar asintóticamente un valor de equilibrio como se muestran en el
gráfico de la figura 2
Se constata, también experimentalmente, que este valor de equilibrio corresponde al promedio M
de concentraciones, el cual es conservado a través del tiempo. Esto tiene dos consecuencias: dicho
promedio debe ser igual al promedio de las concentraciones iniciales y es por lo tanto conocido.
Además, como A(𝑡) + 𝐵(𝑡) = 2𝑀 es constante, podemos obtener en cada instante t la
concentración en B conociendo la de A y viceversa. Así es que el problema se reduce a encontrar
solamente la función A(𝑡).

Si registramos en un gráfico el logaritmo natural de la diferencia 𝑀 − 𝐴(𝑡) en función del tiempo,
observemos que la curva experimental se ajusta bien a una recta de pendiente negativa – 𝑘. Una
hipótesis razonable es entonces un ajuste A(𝑡) a la asíntota de ordenada M, en efecto
ln(𝑀 − 𝐴(𝑡)) = −𝑘𝑡 + 𝐶 ⇒ A(𝑡) = 𝑀 − 𝐻𝑒 −𝑘𝑡 ,
Donde H es una constante que se obtiene imponiendo la condición inicial A(0) para t = 0 lo que
nos da 𝐻 = 𝑀 − 𝐴(𝑜) . Reemplazando la constante H en la expresión anterior, estos nos provee
de la formula siguientes para la concentración buscada:
(1) 𝐴(𝑡) = 𝑀 − (𝑀 − 𝐴(0))𝑒 −𝑘𝑡 .
Esta función representa un buen modelo de la realidad, ya que se ajusta razonablemente bien a
las mediciones experimentales, sin embargo, nos resulta todavía misterioso por qué deberíamos
aceptar este modelo de crecimiento exponencial como un modelo razonable y no otro modelo
diferente, por ejemplo, un ajuste polinomial por pedazos.
Estos nos lleva a preguntarnos ¿hay alguna ley o principio que explique el fenómeno de la
osmosis? Una idea, que resulta ser fundamental, consiste en estudiar si existe una relación simple
entre la concentración A(𝑡) y su aumento 𝐴′ (𝑡). Derivado (1) obtenemos:
𝐴′ (𝑡) = 𝑘𝐻𝑒 −𝑘𝑡 = 𝑘𝐻𝑒 −𝑘𝑡 + 𝑘𝑀 − 𝑘𝑀 = 𝑘(𝑀 − (𝑀 − 𝐻𝑒 −𝑘𝑡 ))
esto es, la relación buscada es: 𝐴′ (𝑡) = 𝑘(𝑀 − 𝐴(𝑡)) (2) Entonces la solución (1)
satisface (2). Interpretando (2) encontramos una relación diferencial simple y comprensible que
podemos enunciar como ley siguiente:
LEY DE OSMOSIS
“El aumento de concentración es proporcional en cada
instante a la diferencia de concentración entre el promedio
asintótico de concentraciones y la concentración actual.
La constante de proporcionalidad cuantifica la
permeabilidad de la membrana.”
La ley de la osmosis representada por la ecuación diferencial (2) provee una interpretación más
intuitiva y profunda del proceso de osmosis. Ahora veremos que (2) tiene como solución (1).
La EDO que modela (ley de osmosis) este fenómeno es A t   k  M  A t  , con k>0 y

A  0  B  0
donde  M es la concentración promedio . La ecuación es equivalente a
2

A  t   kA t   kM multiplicando a esta ecuación por e kt se obtiene
ekt A  t   kA  t  ekt  kMekt   A  t  ekt    Mekt  , entonces existe una constante arbitraria
C 𝝐ℝ tal que A t  e  Me  C  A t   M  Ce
kt kt  kt
.Este resultado es una familia de
curvas (indexadas por la constante C ).
Si evaluamos en el tiempo inicial t = 0, encontramos el valor de la constante C , esto es,

A  0  B  0
A(O) =C + M ⟺ C = A(0) – M ; pero M .
2
A  0   B  0  A  0   B  0   kt
Por lo tanto, la solución es At    e
2 2
2) MEZCLAS
La mezcla de dos líquidos a menudo origina una ecuación diferencial de primer orden.
Más concretamente, se considera un recipiente que contiene una cantidad de V litros de cierto
fluido, en el que se encuentra disuelta una cantidad, y  0  , de cierta sustancia. En el recipiente
entra constantemente fluido con una concentración de C e gramos por litro y a una velocidad de
Ve litros por minuto. Se supone que los fluidos en el recipiente se mezclan de forma instantánea y
que la mezcla sale del recipiente a una velocidad de Vs litros por minuto. Lo que se desea es
determinar una función que indique la cantidad de sustancia que hay en el interior del recipiente
en cada instante, t. Llamemos V  t  a la cantidad de fluido (litros) presente en el recipiente en el
instante t, y y t  a la cantidad de sustancia disuelta (gramos) en el instante t, de forma que la
y t 
concentración de sustancia disuelta en el instante t es gramos por litro.
V t 
La variación de la magnitud y t  por unidad de tiempo es y  t  y viene dada por la diferencia
entre la cantidad de sustancia que entra (por unidad de tiempo) y la cantidad de sustancia que
sale (por unidad de tiempo):
Puesto que entran Ve litros por minuto, que contienen una concentración C e gramos de sustancia
por litro, se tiene que entran C e . Ve Gramos por minuto de sustancia.

La concentración de sustancia en el fluido que sale es la del fluido en el interior del recipiente, es
y t  y t 
decir gramos por litro. Puesto que salen Vs litros por minuto, se tiene que salen Vs
V t  V t 
gramos por minuto de la sustancia disuelta. Así pues, la variación de la concentración, y  t  ,
y t 
y  t   CeVe  Vs
verifica:
V t 
La expresión de V  t  , cantidad de fluido en el recipiente en el instante t, deberá ser determinada
en cada caso, ya que depende de la cantidad inicial y de las velocidades de entrada y salida del
fluido en el recipiente. Si, por ejemplo, la velocidad de entrada de fluido es igual a la velocidad de
salida, entonces el volumen en el interior del recipiente permanecerá constante.
Ejemplo 1
Un depósito contiene 100 litros de una disolución salina cuya concentración es de 2.5 gramos de
sal por litro. Una disolución conteniendo 2 gramos de sal por litro entra en el depósito a razón de 5
litros por minuto y la mezcla (que se supone uniforme de forma instantánea) sale del depósito a la
misma velocidad. Encontrar la cantidad de sal que hay en cada instante en el depósito.
Solución
Puesto que la velocidad a la que entra el líquido en el depósito es la misma a la que sale, en el
depósito siempre hay la misma cantidad de líquido: 100 litros.
Sea y t  la cantidad de sal en el depósito en el instante t.
En el depósito entran 5l. por minuto de una disolución con 2gr. por litro, luego entran 10gr. de sal
por minuto.
Puesto que la cantidad de sal en el depósito es y t  y la cantidad de líquido que hay es 100l., la
y t 
concentración de la disolución en el depósito es de gramos por litro. Esta disolución sale a
100
y t 
una velocidad de 5 litros por minuto, por lo tanto la sal sale a una velocidad de  5 gramos
100
por minuto.

5y
Así pues, se tiene: y  t   10 
100
Esta ecuación es de variables separables
5y 1000  5 y y  t  1 dy dt
y  t   10   y  t      
100 100 1000  5 y 1000 1000  5 y 100
1 1 1
  ln 1000  5 y  t  C  ln 1000  5 y   t  K  1000  5 y  Ke 0.05t
5 100 20
0.05t
1000  Ke
 5 y  1000  Ke 0.05t  y  t    200  Ke 0.05t
5
Así pues, la solución general de la ecuación diferencial es y t   200  Ke0.05t

Puesto que, inicialmente, la concentración de sal en el depósito era de 2:5 gramos por litro, la
cantidad de sal inicial era de y  0  2.5 100  250 . Sustituyendo esta condición inicial en la
expresión de la solución general se tiene 250  y  0  200  Ke0.050  200  K  K  50
Luego la función que nos da la cantidad de sal en cualquier instante t es: y t   200  50e0.05t
Ejemplo 2
La corriente sanguínea lleva un medicamento hacia el interior de un órgano a razón de 3cm3 / sg y
sale de él a la misma velocidad. El órgano tiene un volumen de 125cm3 : Si la concentración del
medicamento en la sangre que entra en el órgano es de 0.2 gr / cm3 , se pide:
1) ¿Cuál es la concentración del medicamento en el órgano en cada instante si inicialmente no
había vestigio alguno del medicamento?
2) ¿Cuándo la concentración del medicamento en el órgano será de 0.1gr / cm3 ?
Solución
La cantidad de medicamento que entra en el órgano por segundo es: 0:2 x 3 = 0:6 gramos
Si denotamos por y(t) la cantidad de medicamento presente en el órgano en el instante t se
tendrá, puesto que la sangre abandona el órgano a la misma velocidad a la que entra ( 3cm3 / sg ),
y t 
que la cantidad de medicamento que abandona el órgano por segundo será de 3
125
En consecuencia, puesto que la variación por unidad de tiempo (i.e., por segundo) de la cantidad
de medicamento viene dada por:
y t   cantidad que entra por segundo-cantidad que sale por segundo
3 75  3 y
Se tiene y  t   0.6  y t  
125 125
Esta ecuación es de variables separables:
75  3 y y 1 dy 1 1 1
Z y      dt   ln 75  3 y  t C
125 75  3 y 125 75  3 y 125 3 125
Despejando aquí y se obtiene la solución general de la ecuación: y t   25  Ke3t /125

Puesto que, inicialmente, no había ninguna cantidad de medicamento en el órgano, la condición
inicial para
Y(t) es y(0)=0, lo que conduce, sustituyendo, a: 0  y  0  25  K  K  25
En consecuencia la función que nos da la cantidad de medicamento en el órgano en cada instante
es y t   25  25e3t /125
La concentración es la cantidad de medicamento dividido por el volumen del órgano, es decir
y  t  25
125 125
 1  e3t /125   1  e3t /125 
1
5
Por lo tanto, la contestación a la primera pregunta es que
la consentracion en el instante t es  1  e 3t /125 
1
5
Para contestar a la segunda pregunta hay que calcular para qué valor de t se verifica
1
0.1 
5
1  e 3t /125   0.5  1  e3t /125  
3t
125
 ln 0.5
125
t  ln 0.5  28.88segundos
3
Ejemplo 3
Un tanque contiene 50 galones de una disolución compuesta por 90% agua y 10% alcohol.
Una segunda disolución que contiene 50% agua y 50% alcohol se agrega al tanque a una tasa de 4
galones por minuto. Conforme se añade la segunda, el tanque empieza a drenar
a una tasa de 5 galones por minuto, como se muestra en la figura . Si se supone que la disolución
en el tanque se agita constantemente, ¿cuánto alcohol permanecerá en el tanque después de 10
minutos?
Solución
Sea y el número de galones de alcohol en el tanque en cualquier instante t.
Se sabe qué y = 5 cuando t =0. Dado que el número de galones en el tanque en cualquier
 5 
tiempo es 50- t, y que el tanque pierde 5 galones por minuto, se debe perder  y galones
 50  t 
de alcohol por minuto.

Además, ya que el tanque gana 2 galones de alcohol por minuto, el ritmo o velocidad de cambio
de alcohol en el tanque está dada por
dy  5  dy  5  5
 2  y   y  2 , Para resolver esta ecuación lineal, sea a (t )  y se
dt  50  t  dt  50  t  50  t
5
obtiene  a (t )dt   dt  5ln 50  t Ya que t <50, se puede eliminar el signo del valor
50  t
e  a ( t ) dt
 e5ln50t  
1
 50  t 
absoluto y concluir que 5
y 2 1 50  t
 dt  C  y   C  50  t 
5
Así, la solución general es
 50  t   50  t  2  50  t 
5 5 4
2
50 20
Dado que t = 5 para t = 0, se tiene 5   C  50   C   5 , lo cual significa que la solución
5
2 50
50  t  50  t 
5
particular es y   20  
2  50 
Por último, cuando t = 10, la cantidad de alcohol en el tanque es

50  10  50  10 
5
y  20    13.45 gal , lo cual representa una solución que contiene 33.6% de

2  50 
alcohol.
Ejemplo 4
Inicialmente, 50 lb de sal se disuelven en un gran tanque que contiene 300 gal de agua. Una
solución de salmuera se bombea hacia el tanque a razón de 3 gal/min, y luego la solución bien
mezclada se extrae al mismo ritmo. Vea la FIGURA. Si la concentración de la solución que entra es
2 lb/gal, determine la cantidad de sal en el tanque en el instante t. ¿Cuánta sal hay después de 50
min? ¿Y después de un gran tiempo?

Solución
Sea A (t) la cantidad de sal (en libras) en el tanque en el instante t. Para problemas de esta clase, la
razón neta a la que A (t) cambia está dada por
dA
 (razon de la sustancia que entra)-(razon de la sustancia que sale)
dt
Luego, la razón a la que la sal entra al tanque, en libras por minuto, es
(3gal / min).(2lb / gal )  6lb / min mientras que la razón a la que la sal sale es
A A
(3 gal / min).( lb / gal )  lb / min
300 100
dA A 1
Por tanto, la ecuación es  6  A  A6
dt 100 100
t
Resolviendo se obtiene A(t )  600  550e100 , para t = 50 se tiene A(50) = 266.41 lb
(iii) y t   ay t   b t  donde b(t) es una función

continua en 0, mientras que a .
Supongamos que y t  es solución de (iii). Luego y(t )  a y(t )  b(t )

Guiados por el camino seguido en (i) y(ii), multiplicamos ambos miembros de
y(t )  ay(t )  b(t ) por e at , obtenemos eat ( y(t )  ay t )  eat b(t ) 
(eat y(t ))  eat b(t ) ……..(1)
Observando el segundo miembro de esta ecuación , el teorema fundamental del cálculo afirma
t
que ( t0 e b( s)ds)  e b(t ) ,donde t0 es cualquier punto (fijo arbitrario) de 0, . Por
as at
t
consiguiente reemplazando en (1), se obtiene (e y(t ))  (  easb(s)ds) , entonces existe una
at
t0
t
constante C tal que
eat y (t )   easb(s)ds  C , de donde
t0
t
y (t )  e at  easb(s)ds  Ce at
t0
es solución de (iii). De esta manera hemos obtenido la solución general de la ecuación diferencial
en discusión; de esta se desprenden todas las infinitas soluciones.
En un primer momento puede perecer que si se escoge otro punto t1 como límite interior, la
formula

t
Ce at  e at  easb(s)ds
t1
Conduce a otro conjunto solución. Veremos que ambas fórmulas dan origen al mismo conjunto, es
decir.
 t
Ce at  e at  eas b( s )ds : C 
t0  t
 He  at  e  at  eas b( s )ds : H 
t1 
t t
Para probar la igualdad basta observar que  t1
easb(s)ds   easb(s)  M
t0
t

t
Para una cierta constante M, al ser 
t0
eas b( s)ds y
t1
eas b( s)ds primitivas de una misma
función. Luego, dado H  tenemos
t t
He at  e at  e at b(s)ds  He at  e at (  easb(s)ds  M ) 
t1 t0
t t
( H  M )e at  e at  easb(s)ds  Ce at  e at  easb(s )ds
t0 t0
donde L = H + M. en forma análoga se demuestra que para cualquier C  se cumple
t t
Ce at  e at  easb(s)ds  He at  e at  easb(s )ds ,para una cierta constante H. la no
t0 t1
importancia del límite inferior lleva muchas veces a escribir la formula general como
t
Ce at  e at  easb(s)ds , formula que lo usaremos en la practica

 y  ay  b  t 

 y  t0   
Si y es solución del PVI entonces debe existir una constante C tal que
t
y(t )  Ce at  e at  eas b( s)ds y además y (t0 )      Ce at0  C   eat0 . Luego
t 0
t
y(t )   e a (t t0 )  e at  easb(s)ds . Este único candidato a ser solución del PVI es, como
t0

se puede ver, solución del problema de valor inicial en discusión. Hemos demostrado la existencia
y unicidad de la solución del PVI.
 y  t   ay  t   b  t 

EL PVI
CIRCUITO SIMPLE
Un circuito eléctrico simple consta de una corriente eléctrica I(en amperes), una resistencia R(en
ohms), una inductancia L(en henrys), y una fuerza electromotriz E constante (en volts), como se
muestra en la figura . Con base en la segunda ley de Kirchhoff, si el interruptor S se cierra cuando
t=0, la fuerza electromotriz aplicada (voltaje) es igual a la suma de la caída de voltaje en el resto
del circuito. De hecho, esto significa que la corriente I satisface la ecuación diferencial
dI
L  RI  E
dt
EJEMPLO
Encontrar la corriente I como función del tiempo t(en segundos), dado que I satisface la ecuación
dI
diferencial L  RI  sen 2t , donde R y L son constantes diferentes de cero.
dt
Solución
dI R 1
En forma normal, la ecuación lineal dada es  I  sen 2t , entonces
dt L L
R R R
 t 1  t 1  t
Ie L
  e L  sen2tdt  2 e L
( Rsen 2t  2L cos 2t )  C
L 4L  R2
R
1  t
Así, la solución general es I  2 ( Rsen 2t  2 L cos 2t )  Ce L
4L  R2
R
1  t
Así, la solución general es I  2 ( Rsen 2t  2 L cos 2t )  Ce L
4L  R2

(IV) y t   a t  y t   b t  , donde a(t) y b(t)son
funciones continuas en 0, .
Guiados por el método aplicado en la solución de y  ay  b(t ) nos preguntamos si será
posible hallar  (t ) tal que e (t )  y(t )  a(t ) y(t )   (e (t ) y(t ))
 (t )
En tal caso e y(t )  e  (t ) a(t ) y (t )   (t )e  (t ) y (t )  e  (t ) y(t )
Aceptando que y(t )  0 para todo t (lo que es razonable en vista de nuestra experiencia
anterior) tendremos µ’(t) = a(t)   (t )   a(t ) ,es decir, cualquier primitiva de a(t) aparece
como buen candidato. Hechas estas consideraciones preliminares. Supongamos que y(t ) es
solución de (iv).Luego y(t )  a(t ) y(t )  b(t ) 
t t
t0 a ( s ) ds t0 a ( s ) ds
e ( y (t )  a(t ) y(t ))  b(t )e
Donde t0 es cualquier punto de 0, , fijo arbitrario, entonces
t s
t0 a ( s ) ds t t0 a ( ) d
(e y(t ))  (  b( s)e ds) , entonces existirá una constante C tal que
t0
t s
t0
 a ( s ) ds t t0 a ( ) d  t a ( s ) ds
 b(s)e
t s
t0 a ( s ) ds t0 a ( ) d y(t )  e ds  Ce
t0
t
e y(t )   b(s)e ds  C , es decir t0
t0
Esta función es solución de (iv) por lo tanto es la solución general esta fórmula impresiona por su
complejidad pero, afortunadamente, en los casos practico no se utiliza sino se repite el proceso.
Deseamos también hacer hincapié que no tiene importancia cual punto t 0 se escoja, es decir
tenemos entera libertad para seleccionar t0.
 y(t )  a(t ) y (t )  b(t )



 y  t0   


Partimos de la solución general
t s
t0
 a ( s ) ds t t0 a ( ) d  t a ( s ) ds
y(t )  e  b(s)e ds  Ce
t0
t0 evaluamos y (t) para t = t0

s
t0 a ( ) d
t t
 t0 a ( s ) ds t  t0 a ( s ) ds
Entonces C =  
y (t )  e  b( s ) e
t0
ds   e
es una solución de
la ecuación ¿esta es la única solución del PVI? Si  (t) es cualquier solución de la ecuación
diferencial en cuestión, entonces existe L tal que
s
t0 a ( ) d
t t
 t0 a ( s ) ds  t0 a ( s ) ds
 (t0 )  
t
 (t )  e 
t0
b( s ) e ds  Le
con
s
t0 a ( ) d
t t
 t0 a ( s ) ds t  t0 a ( s ) ds
Luego
   (t0 )  e  t0
b( s ) e ds  Le
Por consiguiente L   , demostrando de esta manera que  t   y(t ) para todo t. hemos
conseguido justificar la existencia y unicidad de solución del PVI en cuestión.
FUNCIONES DEFINIDAS POR INTEGRALES

En el curso de cálculo integral analizamos el hecho de que algunas funciones continuas simples no
tienen antiderivadas que sean funciones elementales y que las integrales de esa clase de funciones
se llaman no elementales.
Por ejemplo, usted puede haber visto en cálculo que  e x dx   senx 2 dx no son integrales
2
elementales. En matemáticas aplicadas algunas funciones importantes están definidas en términos

de las integrales no elementales. Dos de esas funciones especiales son la función error y la
2 2 
función error complementario: erf  x   et dt y erfc  x  
x
  et dt
2 2
 0
 x
   
 et dt   et dt  1
2 2
Del conocido resultado podemos escribir
0 2 2 0
 x 
Entonces de la forma  0

0

x
se ve de la ecuación que la función error complementario,
erfc(x), se relaciona con erf(x) por erfc  x   erf  x   1 . Debido a su importancia en

probabilidad, estadística y en ecuaciones diferenciales parciales aplicadas se cuenta con extensas
tablas de la función error. Observe que erf  0  0 es un valor obvio de la función.
EJEMPLO La función error
dy
Resuelva el problema con valores iniciales.  2 xy  2, y  0   1
dx

SOLUCIÓN
Puesto que la ecuación es lineal multiplicamos a la ecuación por e  x para obtener
2
d   x2  x
y de aquí y  2e x  et dt  Ce x
2 2 2
 2e  x
2
e
dx   0
Aplicando y  0  1 en la última expresión obtenemos C  1 . Por tanto, la solución del problema
0 e dt  e  e 1   erf  x 

x
es y  2e x t
2 2 2 2
x x
En la fi gura se muestra en azul oscuro, la gráfica de esta solución en el intervalo

 ,  Junto con otros miembros de la familia.
EJERCICIOS PROPUESTOS
1. En el diagnostico medico se utilizan isotopos radioactivos como reactores para determinar
anormalidades que puedan existir en un órgano. Por ejemplo, si se ingiere yodo radioactivo,
después de cierto tiempo es absorbido por la glándula tiroides. Con el uso de un detector puede
medirse la tasa a la cual es absorbido y determinar si la absorción es normal o no. Supóngase que
se va a utilizar en un rastrero cerebral tecnesio-99m, que tiene una semivida de 6 horas, y se le va
a detectar dentro de 2 horas. ¿Cuál debe ser su actividad ahora para que la actividad al momento
se su utilización sea 10 unidades? Proporcione la respuesta con una cifra decimal. , considere N
como la actividad dentro de t horas, y N0 la actividad ahora.]
2. se va a implantar en forma temporal una sustancia radioactiva que tiene una semivida de 8 días
a un paciente de hospital hasta que se tengan tres quintos de la cantidad originalmente presente.
¿Qué tanto tiempo debe permanecer la sustancia implantada en el paciente?
3. La razón en que un medicamento se disemina en el torrente sanguíneo está regida por la

dX
ecuación diferencial  A  BX , donde Ay B son constantes positivas. La función X(t) describe la
dt
concentración del medicamento en el torrente sanguíneo en el instante t. Encuentre X(t). ¿Cuál es
el valor limitante de X(t) cuando t   ¿En qué instante la concentración es la mitad de su valor
limitante? Suponga que X(0) =0.

4. Suponga que una célula está suspendida en una solución que contiene un soluto de
concentración constante C S . Además, suponga que el volumen V de la célula es constante y que el
área de su membrana permeable es la constante A. Por la ley de Fick, la razón de cambio de su
masa m es directamente proporcional al área A y a la diferencia CS  C (t ) , donde C (t ) es la
concentración del soluto dentro de la célula en el instante t. Encuentre C (t ) si m  VC(t ) y
C (0)  C0 . Vea la FIGURA
5. Un marcapasos, mostrado en la FIGURA, consta de una batería, un capacitor y el corazón como

resistor. Cuando el interruptor S está en P , el capacitor se carga; cuando S está en Q , el
capacitor se descarga, enviando un estímulo eléctrico al corazón. Durante este tiempo, la tensión
dE 1
E aplicada al corazón está dada por la ecuación diferencial lineal  E, t1  t  t2 ,
dt RC
donde R y C son constantes. Determine E(t) si E (t1 )  E0 . (Por supuesto, la apertura y el cierre
del interruptor son periódicos a fin de simular el latido cardiaco natural.)
Resta

6. En un modelo para el crecimiento de tejido, sea A(t) el área del cultivo de tejido en el instante t.
Vea la FIGURA Puesto que la mayoría de divisiones celulares se llevan a cabo en la porción
periférica de tejido, el número de células en la periferia es proporcional a Si se supone que la
razón de crecimiento del área es conjuntamente proporcional a A(t ) y M  A(t ) , entonces un
dA
modelo matemático para A está dado por  k A(t )( M  A(t )) , donde M es el área final de
dt
tejido cuando se ha completado el crecimiento.
a) Resuelva la ecuación diferencial
b) Encuentre lim A(t ) .

t 
7. En los ejercicios 8 y 9, un investigador médico quiere determinar la concentración C (en moles

por litro) de un medicamento marcador inyectado en un fluido en movimiento. Resolver este
problema al considerar un modelo de dilución de un compartimento simple
(Ver la figura). Suponer que el fluido está siendo mezclado y que el volumen de éste en el
compartimento es constante.
Figura para 8 y 9

8. Si el marcador es inyectado instantáneamente en el tiempo t=0, entonces la concentración del
dC  R 
fluido en el compartimento se empieza a diluir según la ecuación diferencial     C , para
dt  V 
t  0 se tiene C (0)  C0
a) Resolver la ecuación diferencial para encontrar la concentración C como función de t.
b) Encontrar lim C (t )
t 
9. En el ejercicios 8, se supuso que había una inyección simple inicial del medicamento marcador
dentro del compartimento. Ahora considerar el caso en el cual el marcador es continuamente
inyectado (iniciando en t =0) a una tasa de Q moles por minuto. Si considera Q despreciable
dC Q  R 
comparada con R, usar la ecuación diferencial     C , para t  0 se tiene C (0)  C0
dt V  V 
a) Resolver la ecuación diferencial para encontrar la concentración C como función de t.
b) Encontrar lim C (t )
t 
10. En una habitación de una casa donde la temperatura era constante de 70 F se encontró un
cuerpo sin vida. Al medir la temperatura del cuerpo al momento de su descubrimiento, la lectura
fue de 85 F . Una segunda medición, una hora después, mostró que la temperatura del cuerpo era
de 80 F . Use el hecho de que si t = 0 corresponde al instante de fallecimiento, entonces la
temperatura del cuerpo en ese instante era de 98, 6 F . Determine cuántas horas transcurrieron
entre el fallecimiento y el descubrimiento.
11. Suponga que un sujeto se ha hallado muerto en una habitación de un motel a la medianoche y
su temperatura corporal es de 80° F. La temperatura de la habitación ha permanecido constante
en 60° F. Dos horas después la temperatura del occiso es de 75° F. Determine la hora de
defunción.
RSPTA: murió aproximadamente a las 7:26 P.M
CUERPOS QUE CAEN Y RESISTENCIA DEL AIRE

Se ha establecido experimentalmente que cuando un cuerpo se mueve a través de un medio
resistivo como el aire (o agua), la fuerza retardadora debida al medio, denominada fuerza de
arrastre, actúa en dirección opuesta a la del movimiento y es proporcional a la potencia de la
velocidad del cuerpo, es decir, kv . Aquí k es una constante de proporcionalidad y  es
constante en el intervalo 1    2 . En términos generales, para velocidades más lentas se toma
  1 . Ahora, suponga que en un cuerpo de masa m que cae la resistencia del aire es proporcional
a su velocidad instantánea v . Si en estas circunstancias tomamos la dirección positiva orientada
hacia abajo, entonces la fuerza neta que actúa sobre la masa está dada por mg  kv , donde el
peso mg del cuerpo es una fuerza que actúa en la dirección positiva y la resistencia del aire es una
fuerza que actúa en la dirección opuesta o hacia arriba. Vea la figura.

dv
Luego, puesto que v está relacionada con la aceleración a por  a , la segunda ley de Newton
dt
dv
se vuelve F  ma  m Al igualar la fuerza neta con esta forma de la segunda ley de Newton
dt
obtenemos una ecuación diferencial lineal para la velocidad y del cuerpo en el instante t,
dv
m  mg  kv , donde k es una constante de proporcionalidad
dt
Para movimiento a gran velocidad, como un paracaidista que cae libremente antes de abrir su
paracaídas, suele suponerse que la resistencia del aire es proporcional al cuadrado de la velocidad
instantánea; en otras palabras,   2 . Si de nuevo la dirección positiva se toma hacia abajo,
entonces un modelo para la velocidad v de un cuerpo que cae está dado por la ecuación
dv
diferencial no lineal m  mg  kv 2 , donde k es una constante de proporcionalidad. Vea la
dt
ecuación de Bernoulli.
dv
La ecuación m  mg  kv , se puede escribir de la forma
dt
dv k k
 v  g  v(t )  v  g que es una ecuación de la forma (iii)
dt m m
k
12. Resuelva el PVI v(t )  v  g, v (0)  v0
m
13. La ley de Torricelli establece que el agua fluirá desde una abertura en la parte inferior del
tanque con la misma velocidad que alcanzaría al caer desde la superficie del agua a la abertura.
dh
Una de las formas de la ecuación de Torricelli es A( h )   k 2 gh , donde h es la altura del
dt
agua en el tanque, k es el área de la abertura de la parte inferior del tanque, A(h) es el área de la
sección transversal a la altura h, y g es la aceleración debida a la gravedad ( g  32 pies / s 2 ). Un
tanque de agua hemisférico tiene un radio de 6 pies. Cuando el tanque está lleno, una válvula
circular con un radio de 1 pulgada se abre en la parte inferior, como se muestra en la figura.
¿Cuánto tiempo es necesario para que el tanque se vacíe completamente?

14. El tanque cilíndrico de agua mostrado en la figura tiene una altura de 18 pies. Cuando el
tanque está lleno, una válvula circular se abre en la parte inferior del tanque. Después de 30
minutos, la profundidad del agua es de 12 pies.
a) ¿Cuánto tiempo es necesario para que el tanque se vacíe completamente?

b) ¿Cuál es la profundidad del agua en el tanque después de 1 hora?
15. Suponer que el tanque del ejercicio 14 tiene una altura de 20 pies, un radio de 8 pies, y la
válvula circular tiene un radio de 2 pulgadas. El tanque está completamente lleno cuando la
válvula está abierta. ¿Cuánto tiempo es necesario para que el tanque se vacíe completamente?
16. Un cuerpo cuya temperatura es de 30C requiere de dos minutos para descender su
temperatura a 20C si es colocado en un refrigerador con temperatura constante de 10C .
¿Cuánto tiempo tardará el mismo cuerpo en bajar su temperatura de 40C a 35C suponiendo
que el medio se encuentre a la temperatura constante de 15C ?
17. La temperatura de un cuerpo difiere en 20C de la del medio en que se coloca, cuya
temperatura TM se mantiene constante. A los 3 minutos la diferencia es de 10C .¿Cuánto tiempo
deberá transcurrir para que la diferencia sea de 10C ?
18. Un químico desea enfriar desde 80C hasta 60C una solución contenida en un matraz.
Coloca el dispositivo en un recipiente amplio por el que circula agua a 15C . Observa que
después de 2 minutos la temperatura ha descendido a 70C . Estimar el tiempo total de
enfriamiento.

19. Se tienen dos recipientes con agua que se mantiene a una temperatura constante de 20C
el recipiente A y 14C el recipiente B . A las 8:00 a.m un termómetro que marca 25C
es puesto en el recipiente B. A las 8:01 a.m el termómetro marca 20C . Mas tarde el
termómetro es puesto en el recipiente A . Luego a las 8:06 a.m el termómetro marca 23C .
¿Cuándo fue llevado el termómetro del recipiente B al A ?
20. Describir y resolver el sistema cinético de primer orden
A 
k1
C
B 
k2
C
Donde  A0 ,  B0 cantidades iniciales.
21. Aceptando una cinética de primer orden, plantear y resolver las ecuaciones que gobiernan el
sistema cinético A  k1 
k2
 B  C
k3
Donde  A0 ,  B0  C 0  0 son las condiciones iniciales en t  0 .

22. Si el 45% de una sustancia radioactiva se desintegra en 200 años. ¿Cuál es su vida media? En
cuánto tiempo se desintegrará 60% de la cantidad original?
23. Una sustancia radioactiva A se descompone dando lugar a una sustancia radioactiva B , a la
vez se desintegra dando lugar a un producto estable C .
A 
k1
 B 
k2
 C
en el instante t  0 se tienen 10 mg de A mientras que no hay nada de B y C . La vida media

de A es de 2 horas, mientras que de B es 1 hora. ¿Cuál es el valor de B(t ) y C (t ) después de
2 horas?
24. Hallar la solución de PVI recursivo
Pn(t )    [(n  1) Pn (t )  nPn 1 (t )] para n  0,1, 2...
P0 (0)  1 y Pn (0)  0 para n  1, 2...
Sistema que aparece en el estudio de procesos estocásticos

Iedo 1

Cargado por

Copyright:

Formatos disponibles

Iedo 1

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Iedo 1

Cargado por

Copyright:

Formatos disponibles

CAPITULO I

Georg Friedrich Riemann (1826-1866)

Georg Friedrich Riemann

Julio Flores Dionicio Página 1

ECUACION DIFERENCIAL ORDINARIA

ECUACION DIFERENCIAL PARCIAL

DEFINICIÓN 2: ORDEN DE UNA ECUACION DIFERENCIAL

Por ejemplo y  xy es de orden 1 ; y   y   xy  0 es de orden 2.

Julio Flores Dionicio Página 2

DEFINICIÓN 3: DE UNA ECUACION DIFERENCIAL ORDINARIA LINEAL

Donde las funciones ai ( x)  , i 1,...n son llamadas coeficientes de la EDO lineal

DEFINICION 5. Si los coeficientes ai ( x) no dependen de x, se dice que la EDO lineal es a

Se dice que está normalizada.

EJEMPLO 1. xy  sen  x   tg  x  . EDO lineal de orden 1 a coeficientes variables, no homogénea

Julio Flores Dionicio Página 3

SOLUCION DE UNA ECUACION DIFERENCIAL ORDINARIA

Cualquier función   x  , defi nida en un intervalo I y que tiene al menos n derivadas

Julio Flores Dionicio Página 4

EJEMPLO 1. Consideremos la ecuación diferencial y’ = 0 en un intervalo I. recordemos que la

Julio Flores Dionicio Página 5

Julio Flores Dionicio Página 6

Julio Flores Dionicio Página 7

EXISTENCIA Y UNICIDAD DE SOLUCION DE UNA ECUACION DIFERENCIAL

¿Existe una solución al problema? Si la hay, ¿es única?

Para un problema de valor inicial, como

Julio Flores Dionicio Página 8

TEOREMA: Existencia de una solución única

Julio Flores Dionicio Página 9

La función y  Ce x es una familia uniparamétrica de soluciones de la ecuación de primer orden

es una solución del PVI planteado

Julio Flores Dionicio Página 10

INTERVALO DE EXISTENCIA Y UNICIDAD

Julio Flores Dionicio Página 11

El intervalo I depende tanto de F  x, y  como de la condición inicial y  x0   y0 . La condición

Julio Flores Dionicio Página 12

ECUACIONES DIFERENCIALES ORDINARIAS

Julio Flores Dionicio Página 13

ECUACIONES DIFERENCIALES ORDINARIAS

 PRIMER TIPO: INTEGRACION DIRECTA (ID): y  f ( x) (1)

Julio Flores Dionicio Página 14

 TERCER TIPO: LINEALES DE PRIMER ORDEN (LI): y  a t  y  b t  (3)

TEOREMA (TFC). Sea f integrable en  a, b , entonces, dado x0  a, b e y0  , la función y,

Es continua en  a, b y se tiene y ( x0 )  y0 . Si además f es continua en  a, b entonces la función

Nos enfocaremos en la resolución de las ecuaciones y no siempre seremos demasiados rigurosos

PRIMER TIPO: INTEGRACION DIRECTA (ID)

Es una ecuación de la forma

EJEMPLO 1. Las soluciones de la ecuación y  sen( x) son de la forma

y   sen( x)dx  C   cos( x)  C , con C  R

EJEMPLO 2. La ecuación y  x tiene como soluciones a las funciones de la forma

Julio Flores Dionicio Página 15

En el cálculo anterior hemos reemplazado la constante arbitraria C por ln(k) (donde k  e  0 )

PROBLEMA DE CAUCHY (PVI)

 y( x)  f ( x) para todo x  I ,

Integrando la ecuación y  f ( x) entre x0 y x  I obtenemos

Julio Flores Dionicio Página 16

primitiva bien particular de f: aquella que cumple F ( x0 )  0 . Entonces la constante C en (7)

SEGUNDO TIPO: VARIABLES SEPARABLES (vs)

Se dice que es de variables separables