Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

TEMA 1.
Ecuaciones Diferenciales de Primer Orden y sus
aplicaciones
1.1 Conceptos básicos de ecuaciones diferenciales.

1.1.1 Definición de ecuación diferencial.

Definición 1.1 Si una ecuación contiene las derivadas o diferenciales de una o más variables
dependientes con respecto a una o más variables independientes, se dice que es
una ecuación diferencial.

Ejemplo. Las siguientes ecuaciones son ejemplos de ecuaciones diferenciales

dy
a)  cos x
dx
d2y
b) 2
 k2y  0
dx
c) (x 2  y 2)dx  2xydy  0
u 2  u
2
 2u 
d)  h  2  2 

t  x y 
d 2i di 1
e) L 2
 R  i  E cos ωt
dt dt C
3
 d2w  dw
f)  2   xy w0
 dx  dx
3
d2y  dy 
g) 2
 7   8y  0
dx  dx 
d y d2x
2
h)  x
dt 2 dt 2

Una ecuación que contiene sólo derivadas ordinarias de una o más variables dependientes con
respecto a una sola variable independiente es llamada ecuación diferencial ordinaria

Ejemplo. Las siguientes ecuaciones son ejemplos de ecuaciones diferenciales ordinarias

dy
a)  5y  1
dx
b) (x  y)dx  4 ydy  0
du dv
c)  x
dx dx
d) y''  2 y'  6 y  0

Una ecuación que contiene las derivadas parciales de una o más variables dependientes con

1
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

respecto a dos o más variables independientes se llama ecuación diferencial parcial.

Ejemplo. Las siguientes ecuaciones son ejemplos de ecuaciones diferenciales parciales

u u
a) 
x y
u u
b) x y u
x y

1.1.2 Orden de una ecuación diferencial.

Orden de una ecuación diferencial

Definición 1.2 El orden de una ecuación diferencial es igual al orden de la más alta derivada de
las variables dependientes que aparecen en la ecuación.

Grado de una ecuación diferencial

El grado de una ecuación diferencial es el grado algebraico o exponente de la derivada de más
alto orden de la ecuación.

Ejemplo. Determina el orden y grado de las ecuaciones diferenciales

orden grado
x 2
y' ' '2e y' ' x y' y  x 4 3 1
( y' ' )3  xy( y' ) 4  y  x 2 3
yyIV  y'  20 4 1

1.1.3 Solución de una ecuación diferencial.

Definición 1.3 Una solución de una ecuación diferencial ordinaria de orden n, sobre el
intervalo de definición de la ecuación, es una función que posee por lo menos n
derivadas en dicho intervalo y satisface la ecuación

Ejemplo.
a) Verificar que y  e 2 x  3 es solución en todos los números reales de
y'2 y  6
b) Verificar que y  x ln x es solución en (0 , ) de
2

x 2 y' '3xy'4 y  0
c) Verificar que y  x 2 ln x es solución en (- , 0) de
x 2 y' '3xy'4 y  0
x

 e dt  Ce es solución de
x 2 2
x2
d) Verificar que y  e t

y'2xy  1

2
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

e) Verificar que x 2 y  y 2  C es solución de

2xydx  ( x 2  2 y)dy  0
Solución
a) La función y  e 2 x  3 está definida en todos los números reales y tiene derivada en todos los
reales. Derivamos a y con respecto a x, se tiene
y'  2e 2 x
sustituimos en la ecuación diferencial ( y'2 y  6 ), obtenemos

2e 2 x  2(e 2 x  3)  2e 2 x  2e2 x  6  6
como se tiene la identidad, entonces si es solución.

b) La función y  x 2 ln x está definida en el intervalo (0 , ) y tiene primer y segunda derivada

en este intervalo. Obtenemos la primera y segunda deriva con respecto a x
y'  2x ln x  x
y' '  2 ln x  2  1  2 ln x  3

sustituimos en la ecuación ( x 2 y' '3xy'4 y  0 ), se tiene

x 2 (2 ln x  3)  3x (2 x ln x  x )  4x 2 ln x
 2x 2 ln x  3x 2  6x 2 ln x  3x 2  4x 2 ln x  0
como se tiene la identidad, entonces si es solución.

c) Como la función y  x 2 ln x no está definida en el intervalo ( , 0) , entonces no es solución

de la ecuación

e
x2 t2 2
d) La función y  e dt  Ce  x está definida en todos los números reales. Para obtener la
0

primera derivada utilizamos el teorema fundamental de cálculo

x

e
x 2 t2 2 2 2
y'  2xe dt  e  x e x  2Cxe  x
0
x

e
2
t2 2
y'  2xe  x dt  1  2Cxe  x
0

sustituimos en la ecuación y'2xy  1

 x 2 x t 2 2 
x

0  e  e dt  Ce  x 
2
t2 x 2
 2xe  x e dt  1  2Cxe  2 x  
 0 
x x

e e
x 2 t2 x 2 x 2 t2 2
 2xe dt  1  2Cxe  2 xe dt  2Cxe  x
0 0

11
como se tiene la identidad, entonces si es solución.

3
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

e) Derivamos implícitamente la expresión x 2 y  y 2  C

dy dy
2 xy  x 2  2y 0
dx dx
dy
2 xy  ( x 2  2 y) 0
dx
2 xydx  ( x 2  2 y)dy  0
como se obtiene la ecuación diferencial, entonces si es solución
Tarea 1
1.1.4 Problema de valores iniciales.

Condiciones iniciales
Definición 1.4 Para una ecuación diferencial ordinaria de orden n de la forma

y ( n )  f x , y, y' , ..., y ( n 1)  (1.1)
son comunes condiciones del tipo
y(xo) = k0, y'(xo) = k1, . . . ,y(n - 1)(xo) = kn - 1 (1.2)
donde los ki son números reales y xo es un punto del intervalo donde está
definida la ecuación. Al conjunto de condiciones de la forma (1.2), para la
ecuación (1.1) se llama condiciones iniciales. La ecuación (1.1) y las
condiciones iniciales constituyen un problema con valores iniciales.

Condiciones a la frontera
Cuando las condiciones no son dadas en un mismo punto xo, las condiciones se llaman
condiciones a la frontera.

1.2 Solución de ecuaciones diferenciales de primer orden.

Veremos la como encontrar la solución para seis tipos de ecuaciones diferenciales ordinarias de
primer orden. Estoa seis tipos de ecuaciones diferenciales no son mutuamente excluyentes. Esto
quiere decir que, una ecuación diferencial puede ser de dos o más tipos por lo que puede ser
resuelta por varias formas.

1.2.1 Ecuaciones de variables separables.

Una ecuación diferencial ordinaria de primer orden se puede escribir de la forma

M ( x , y)dx  N( x , y)dy  0 (1.3)

Ecuaciones de variables separables

Definición 1.5 Si M de la ecuación (1.3) es una función de x únicamente y N de y, entonces la
ecuación (1.3) toma la forma
M ( x )dx  N( y)dy  0

4
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

la cual se llama ecuación de variables separables.

Método de solución
Dada la ecuación diferencial
M ( x , y)dx  N( x , y)dy  0
mediante pasos algebraicos debemos llegar a una ecuación diferencial de la forma
P( x )dx  Q( y)dy  0
y la solución es
 P(x)dx   Q( y)dy  C
Al aplicar el método de variables separables en una ecuación diferencial ordinaria de primer
orden debemos contestar en forma afirmativa la pregunta ¿somos capaces de juntar las x con dx y
juntar las y con dy?, en caso contrario debemos intentar resolver la ecuación diferencial por otro
método.

Ejemplo. Resolver las siguientes ecuaciones diferenciales

dy
a) y  x(y 2  1)
dx
b) x y'  1  x 2  y 2  x 2 y 2
2

dy 1 2x  1
c) 2  
dx y y
5
d) y'2 y  1 , y(0) 
2
Solución
dy
a) y  x(y 2  1)
dx

Hagámonos la pregunta ¿podemos separar las variables?, es decir, somos capaces de escribir la
ecuación de modo que el coeficiente de dx sea una expresión que contenga únicamente la
variable x y el coeficiente de dy sea una expresión que contenga únicamente la variable y
mediante pasos algebraicos. Si la respuesta es si, podemos encontrar la solución de la ecuación
diferencial por el método de variables separables, si la respuesta es no resolveremos la ecuación
diferencial por otro método. Nuestra respuesta es si, separamos las variables
ydy  x(y 2  1)dx
y
dy  xdx
y 1
2

integramos
y
 y2  1 dy   xdx
1 1
ln(y2  1)  x 2  C1
2 2

5
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

multiplicamos por 2 toda la expresión

1 1 
2 ln(y2  1)  x 2  C1 
2 2 
ln(y2  1)  x 2  2C1
la constante C1 tiene un valor desconocido, por lo que la constante 2C1 también tiene un valor
desconocido que lo podemos escribir como C. La solución es
ln(y 2  1)  x 2  C

b) x 2 y'  1  x 2  y 2  x 2 y 2

Separamos las variables

dy
x2  1  x 2  y2  x 2 y2
dx
factorizamos y2 en los dos últimos términos
dy
x2  (1  x 2)  (1  x 2)y 2
dx
2
factorizamos 1  x
dy
x2  (1  x 2 )(1  y 2 )
dx
x 2 dy  (1  x 2 )(1  y 2 )dx
x2
dy  (1  x 2 )dx
1 y 2

1 1  x2 
dy   2 dx
1  y2  x 
dy  1 
  2  1dx
1  y 2
x 
integramos
dy  1 
1  y 2
   2  1dx  C
x 
1
arctan y    x  C
x
la solución es
1
arctan y  C  x
x

dy 1 2x  1
c) 2  
dx y y

Separamos las variables

dy 2 x  1 1 2x  2
2   
dx y y y
2 ydy  (2 x  2)dx

6
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

integramos
 2ydy   (2x  2)dx  C
y 2  x 2  2x  C
la solución es
y2  x 2  2x  C

5
d) y'2 y  1 , y(0) 
2

Separamos las variables

y'  1  2 y
dy
 1  2y
dx
dy
 dx
(1  2 y)
integramos
dy
 (1  2 y)   dx  C
1
 ln 1  2 y  x  C
2
sustituimos la condición inicial
1 5
 ln 1  2   0  C
2 2
despejamos C
1 1
C   ln  4   ln 4
2 2
sustituimos el valor de C en la solución de la ecuación
1 1
 ln 1  2 y  x  ln 4
2 2
simplificamos
ln 1  2 y  2x  ln 4
ln 1  2 y  ln 4  2x
a
aplicamos la propiedad de los logaritmos ln a  ln b  ln
b
 1  2y 
ln   2 x
 4 
1  2y
 e 2x
4
1  2 y  4e  2 x
despejamos y
1  4e 2 x  2 y

7
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

1  4e 2 x  2 y
1  4e  2 x
y
2
la solución es
1  4e 2 x
y
2
Tarea 2
Ecuación homogénea
Antes de ver la definición de una ecuación diferencial homogénea necesitamos definir función
homogénea de grado k en x y y.

Función homogénea de grado k en x y y

Definición 1.6 Se dice que la función f(x , y) es homogénea de grado k en x y y si y sólo si
f ( tx , ty)  t k f ( x , y)

Ejemplo. Determine si las siguientes funciones son homogéneas y cual es su grado

y3
a) f ( x , y)  x 2  2xy 
x
x
x3
b) f(x , y)  e y  xy
y
x3 x
c) f(x , y)   xy cos
y y
d) f(x , y)  x 2  y 2
x3
e) f(x , y)   xy
y2
Solución
y3
a) f ( x , y)  x 2  2xy 
x

Valuamos la función en ( tx , ty)

( ty)3
f ( tx , ty)  ( tx )  2( tx )( ty) 
2

( tx )
y3  y3 
 t 2 x 2  2t 2 xy  t 2  t 2  x 2  2 xy  
x  x
f ( tx , ty)  t 2f ( x , y)
entonces f es homogénea de grado 2.

Otra forma de ver si es homogénea. Es hacer la suma algebraica de los exponentes en cada uno de
los términos y si todas las sumas de todos los términos son iguales, entonces la función es
homogénea del grado igual a la suma, es decir,

8
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

f ( x , y)  y3
2 xy
x2 + –
x
exponentes 2 11 3 1
2 2 2
las sumas de los exponentes de cada término son iguales, entonces f es homogénea de grado 2,
como ya sabíamos.

x
x3
b) f(x , y)  e y  xy
y

Valuamos la función en ( tx , ty)

tx
(tx)3 ty
x
t 3x 3 y 2 2 x
3 x  3 xy
2 x

f(tx , ty)  e  (tx)(ty)  e  t xy  t e  t xy  t
y 2
e  xy 
ty ty y  y 
 
la función f es homogénea de grado 2

Al usar el método del inciso a), es necesario que los argumentos de las funciones trascendentales
sean de grado 0, como en la exponencial se observa que su argumento es de grado cero
x3 x / y
f(x , y) = e + xy
y
exponentes 3 – 1 1+1
2 2

x3 x
c) f(x , y)   xy cos
y y

Valuamos la función en ( tx , ty)

(tx)3 tx
f(tx,ty)   (tx)(ty) cos
ty ty
t 3x 3 2 x x3 2 x
  t xy cos  t 2  t xy cos
ty y y y
 x3 x
f(tx,ty)  t 2   xy cos 
 y y
f(tx,ty)  t 2f(x,y)
entonces f s homogénea de grado 2

Al usar el método del inciso a), es necesario que los argumentos de las funciones trigonométricas
sean de grado 0, como en el coseno se observa que su argumento es de grado cero, se tiene
x3 x
f ( x , y)  + xy cos
y y
exponentes 3  1 11
2 2
como las sumas de los términos son iguales la función es homogénea de grado 2.

9
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

d) f(x , y)  x 2  y 2

Valuamos la función en ( tx , ty)

f ( tx , ty)  ( tx ) 2  ( ty) 2  t 2 x 2  t 2 y 2  t 2 ( x 2  y 2 )  t x 2  y 2
la función es homogénea de grado 1.

Usamos el método del inciso a), se tiene

 
1
f ( x , y)  x 2  y 2  x 2  y 2 2
el grado de la expresión dentro de la raíz es 2. Empleamos la propiedad de los exponentes,
multiplicamos
1
2   1
2
dando como grado de la función igual a 1.

x3
e) f(x , y)   xy
y2

Usamos el método del inciso a), se tiene

x3
f ( x , y) 
+ xy
y2
exponentes 32 11
1 2
como las sumas no son iguales la función no es homogénea.

Ecuación diferencial homogénea

Definición 1.7 La ecuación
M( x , y)dx  N( x , y)  0
es homogénea si y sólo si M(x , y) y N(x , y) son funciones homogéneas del
mismo grado.

Método de solución
Una ecuación de la forma
M( x , y)dx  N( x , y)  0
homogénea, se reduce a una ecuación de variables separables usando un cambio de variable de la
forma y  ux o x  vy , donde u y v son nuevas variables dependientes, ysa sólo un cambio de
variable.

y x
Al resolver la ecuación resultante y reemplazando u por o v  obtenemos la solución de la
x y
ecuación diferencial.

10
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

Ejemplo. Resolver las siguientes ecuaciones

a) ( x  y)dx  xdy  0
 2  xy 
b)  x e  y 2 dx  xydy
 
 
c)  
ydx  2 xy  x dy  0
 x 
d) ydx   y cos  x dy  0 , y(0) = 2
 y 
Solución
a) ( x  y)dx  xdy  0

Se verá si las funciones M y N son homogéneas del mismo grado

M ( x , y)  x  y N ( x , y)  x
exponentes 1 1 1
como todos los términos tienen la misma suma de exponentes, la ecuación diferencial es
homogénea.

Al escoger el cambio de variable observaremos que diferencial tiene menos términos (esto nos
puede facilitar las integrales). Como dy tiene menos términos, se escogerá el cambio y = ux,
calculamos el diferencial de y
dy  xdu  udx
sustituimos en la ecuación, se tiene
( x  ux )dx  x ( xdu  udx )  0
separamos las variables
( x  ux )dx  x 2du  xudx  0
( x  2ux )dx  x 2du  0
x (1  2u )dx  x 2du  0
1 1
dx  du  0
x 1  2u
integramos
dx du
 x

1  2u
 C1

1
ln x  ln 1  2u  C1
2
2 ln x  ln 1  2u  2C1  C 2
ln x 2  ln 1  2u  C2
ln x 2 1  2u  C2
x 2 (1  2u )  eC 2  C
regresamos a las variables originales
 y
x 2 1  2   C
 x

11
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

x 2  2xy  C

 2  xy 
b)  x e  y 2 dx  xydy
 
 

Se verá si las funciones M y N son homogéneas del mismo grado

M(x , y)  x 2e - y / x  y 2 N(x , y)  xy
exponentes 2 2 11
2 2 2
como todos los términos tienen el mismo exponente, entonces la ecuación es homogénea.

Como la diferencial dy tiene menos términos se escogerá el cambio de variable y = ux.

Calculamos el diferencial
dy  xdu  udx
sustituimos en la ecuación
 2  uxx 
 x e  (ux) 2 dx  x(ux)(xdu  udx)
 
 
separamos las variables
 2  uxx 
 x e  u 2 x 2 dx  x 3udu  x 2 u 2dx
 
 
2 u
x e dx  x udu
3

dx
 ue u du
x
integramos
dx
 x   ue du  C
u

integramos por partes

Uu dv  e u du
dU  du v  e u
ln x  ue u   e u du  C
ln x  ue u  e u  C
regresamos a las variables originales, se tiene
y y
y
ln x  e x  e x  C
x

c)  
ydx  2 xy  x dy  0

Se verá si las funciones M y N son homogéneas del mismo grado

M ( x, y)  y , N( x, y)  2 xy  x
1 1 1

12
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

como todos los términos tiene el mismo exponente, entonces la ecuación es homogénea. Como la
diferencial dx tiene menos términos se escogerá el cambio de variable x  vy . Calculamos el
diferencial
x  vy , dx  vdy  ydv
sustituimos en la ecuación diferencial

y( vdy  ydv)  2 ( vy) y  vy dy  0 
separamos variables

yvdy  y 2 dv  2 vy 2  vy dy  0 

y 2 dv  2 vy 2  vy  yv dy  0 
y 2 dv  2 vy 2 dy  0
y 2 dv  2 y vdy  0
dv dy
 0
2 v y
integramos
dv dy
 0
2 v y
dv dy
2 v  y B
v1 / 2  ln y  B
regresamos a las variables originales
x
 ln y  B
y
x
ln y  B 
y
x x
B 
ye y
 Ce y

x

y  Ce y

la solución es
x
y
ye C

 x 
d) ydx   y cos  x dy  0
 y 

Se verá si las funciones M y N son homogéneas del mismo grado

x
M(x , y)  y N(x , y)  y cos  x
y
exponentes 1 1 1
como todos los términos tiene la misma suma, entonces la ecuación es homogénea.

13
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

Como la diferencial dx tiene menos términos se escogerá el cambio de variable x  vy .

Calculamos el diferencial
dx  ydv  vdy
sustituimos en la ecuación diferencial
 vy 
y(ydv  vdy)   y cos  vy dy  0
 y 
separamos variables
y 2dv  yvdy  y cos v  vydy  0
y 2dv  y cos vdy  0
dv dy
 0
cos v y
dy
sec vdv  0
y
integramos
dy
 sec vdv   y
C

ln sec v  tan v  ln y  C
regresamos a las variables originales
x x
ln sec  tan  ln y  C
y y
aplicamos la condición inicial
0 0
ln sec  tan  ln 2  C
2 2
C  ln 2
sustituimos el valor de C
x x
ln sec  tan  ln y  ln 2
y y
 x x
ln y sec  tan   ln 2
 y y
 x x
y sec  tan   2
 y y
la solución es
 x x
y sec  tan   2
 y y
Tarea 3
1.2.2 Ecuaciones diferenciales exactas.

Antes de ver la definición de una ecuación diferencial exacta, debemos ver la definición de una
diferencial exacta.

14
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

Diferencial exacta
Definición 1.8 Una expresión diferencial
M( x , y)dx  N( x , y)dy
es una diferencial exacta en una región R del plano xy si corresponde a la
diferencial total de alguna función F(x , y), es decir,
F(x , y) F(x , y)
dF  dx  dy  M(x , y)dx  N(x , y)dy
x y

Ecuación diferencial exacta

Definición 1.9 Una ecuación
M( x , y)dx  N( x , y)dy  0
es una ecuación diferencial exacta si la expresión del miembro izquierdo de la
igualdad es una diferencial exacta.
Si
M ( x , y)dx  N( x , y)dy  0 (1.4)
es exacta, entonces por definición existe F tal que
dF  Mdx  Ndy
y F( x , y)  C es la solución general de (1.4)

M N
Teorema 1.1 Si M ( x , y) , N( x , y) y y son funciones continuas de x y y. Entonces la
y x
ecuación M ( x , y)dx  N( x , y)dy  0 es una ecuación diferencial exacta si y sólo
si
M(x , y) N(x , y)

y x

Método de solución
Es importante tener la ecuación diferencial en la forma
M ( x, y)dx  N( x, y)dy  0
para poder identificar claramente las funciones M y N.

Se integra la función M con respecto a x, es decir,

F( x , y)   M ( x , y)dx  g( y)

para encontrar g(y). Derivemos a F(x , y) con respecto a y y la igualamos con N(x , y), es decir,
F(x , y)   M(x , y)dx
  g' (y)  N(x , y)
y y

15
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

despejamos g'(y) (observe que al despejar g'(y), se debe obtener una función que depende
solamente de y) e integrando obtenemos g(y). Y la solución es
F( x , y)   M( x , y)dx  g( y)  C

También se puede integrar la función N con respecto a y, es decir,

F( x , y)   N( x , y)dy  h ( x )
para encontrar h(x). Derivemos a F(x , y) con respecto a x y la igualamos con M(x , y), es decir,
F(x , y)   N(x , y)dx
  h' (x)  M(x , y)
x x
despejamos h'(x) (observe que al despejar h'(x), se debe obtener una función que depende
solamente de x) e integrando obtenemos h(x). Y la solución es
F( x , y)   N( x , y)dy  h ( x )  C

Ejemplo. Resuelva las siguientes ecuaciones diferenciales

a) (5x  4 y)dx  (4x  8 y3 )dy  0
b) (2 x  y)dx  ( x  6 y)dy  0
 s
c) 1  ln r  dr  (1  ln r)ds
 r
d) 2ysenx cos x  y  2y e dx  x  sen x  4xye  dy
2 xy 2 2 xy 2

2x y 2  3x 2
e) dx  dy  0 , y(1) = 1
y3 y4
Solución
a) (5x  4 y)dx  (4x  8 y3 )dy  0

La ecuación está en la forma requerida, por lo que

M ( x , y)  5 x  4 y y N ( x , y)  4 x  8 y 3

Al derivar parcialmente con respecto a una variable, las demás variables se consideran constantes
(Si se deriva con respecto a x, se considera que y es constante. Si se deriva con respecto a y, se
considera que x es constante). Calculemos las derivadas parciales
M  (5x  4y) N  (4x - 8y3 )
  4,  4
y y x x
como las dos derivadas parciales son iguales, entonces la ecuación diferencial es exacta.
Integramos M con respecto a x
F( x , y)   (5x  4 y)dx  g ( y)

Al integrar con respecto a x, se considera que y es constante (si se integra con respecto a y, se
considera que x es constante
5
F(x , y)  x 2  4xy  g(y)
2

16
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

derivamos a F con respecto a y e igualamos con N

F
 4x  g'(y)  4x  8y3
y
despejamos g', se tiene
g' ( y)  8y3
observe que g' depende solamente de y.

Integramos
g ( y)    8y3dy  2 y 4
sustituimos g(y) en F(x,y) y la igualamos a C
5 2
x  4xy  2 y 4  C
2
entonces la solución es
5x 2  8xy  4 y 4  C

b) (2x  y)dx  ( x  6 y)dy  0

La ecuación no está en la forma requerida, por lo que la escribiremos en la forma requerida

(2x  y)dx  ( x  6 y)dy  0
por lo que
M ( x , y)  2 x  y N ( x , y)   x  6 y
calculamos las derivadas parciales
M N
 1,  1
y x
como las derivadas parciales no son iguales, la ecuación diferencial no es exacta. Por lo que, la
ecuación no se puede resolver por el método de ecuaciones diferenciales exactas, pero se puede
resolver por el método de ecuaciones diferenciales homogéneas.

 s
c) 1  ln r  dr  (1  ln r)ds
 r

Para resolver esta ecuación diferencial podemos reemplazar a r por x y a s por y. O considerar
mentalmente a r como x y a s como y, y seguir en método de solución, esto es lo que haremos.

La ecuación no está en la forma requerida, por lo que la escribiremos en la forma requerida

 s
1  ln r  dr  ( ln r  1)ds  0
 r
por lo que
s
M(r, s)  1  ln r  , N(r, s)  ln r  1
r
calculamos las derivadas parciales
M 1 N 1
 
s r r r

17
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

como las derivadas parciales son iguales, entonces la ecuación diferencial es exacta. Integramos a
N con respecto a r, s se considera una constante
F(r, s)   (ln r  1)ds  h (r )  s ln r  s  h (r )
derivamos a F con respecto a r e igualándolo con M
F s s
  h' (r)  1  lnr 
r r r
despejamos h ' (r )
h ' (r )  1  ln r
observe que h ' (r ) depende solamente de r. Integramos
h (r )  r  r ln r  r  r ln r
sustituimos h (r ) en F(r, s)
s ln r  s  r ln r  C
(r  s) ln r  C  s
entonces la solución es
(r  s) ln r  C  s

d) 2ysenx cos x  y  2y e dx  x  sen x  4xye  dy

2 xy 2 2 xy 2

La ecuación no está en la forma requerida, por lo que la escribiremos en la forma requerida

 2
  2

2 ysenx cos x  y  2 y 2e xy dx  sen 2 x  4 xyexy  x dy  0
por lo que
2 2
M(x , y)  2 ysenx cos x  y  2 y 2e xy , N(x , y)  sen 2 x  4 xyexy  x
calculamos las derivadas parciales
M 2 2 N 2 2
 2senx cos x  1  4 yexy  4xy3e xy ,  2senx cos x  4 yexy  4 xy3e xy  1
y x
como las derivadas parciales son iguales, entonces la ecuación diferencial es exacta. Integramos a
M con respecto a x
 
F(x , y)   2 ysenx cos x  y  2 y 2e xy dx  g(y)  ysen 2 x  xy  2e xy  g(y)
2 2

derivamos a F con respecto a y e igualamos con N

F 2 2
 sen 2 x  x  4xyexy  g'(y)  sen 2 x  4xyexy  x
x
despejamos g'(x)
g'(x) = 0
observe que g'(x) depende solamente de x. Integramos
g(x) = 0
entonces la solución es
2
ysen 2 x - xy  2e xy  C

2x y 2  3x 2
e) dx  dy  0
y3 y4

18
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

Calculamos las derivadas parciales

 2x   y 2  3x 2 
 3   
M  y  6x N  y 4  6x
  4 ,   4
y y y x x y
como las derivadas parciales son iguales, entonces la ecuación diferencial es exacta. Integramos a
M con respecto a x
2x x2
F( x , y)   3 dx  g ( y)  3  g( y)
y y
derivamos a F con respecto a y e igualamos con N
F 3x 2 y 2  3x 2
  4  g ' ( y) 
y y y4
despejamos g' ( y)
y 2  3x 2 3x 2 y 2  3x 2  3x 2 y 2
g ' ( y)   4   4
y4 y y4 y
1
g ' ( y)  2
y
integramos a g' ( y) con respecto a y
1
g ( y)  
y
sustituimos g(y) en F(x , y)
x2 1
F( x , y)  3   C
y y

Aplicamos la condición inicial

(1) 2 1
 C
(1) 3 1
obtenemos
C=0
sustituimos el valor de C en F(x , y) y simplificamos
x2 1
 0
y3 y
x 2  y2  0
la solución es
x 2  y2
Tarea 4

1.2.3 Ecuaciones diferenciales reducibles a exactas mediante un factor integrante

que dependa de una sola de las variables.

Factores integrantes

19
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

Definición 1.10 Si la ecuación

M ( x , y)dx  N( x , y)dy  0
no es exacta, escogeremos una función u = u(x , y) tal que
uM( x , y)dx  uN( x , y)dy  0
sea exacta. La función u se llama factor integrante.

1
Ejemplo. Demostrar que la función u ( x , y)  es un factor integrante de la ecuación
x2y
 y 2dx  ( x 2  xy)dy  0
y encontrar la solución
Solución
Primero veamos si la ecuación diferencial no es exacta. Como la ecuación está en la forma
requerida, entonces
M ( x , y)   y 2 N ( x , y)  x 2  xy
calculamos la derivadas parciales
M N
 2 y  2x  y
y x
las derivadas parciales no son iguales multiplicamos la ecuación por el factor integrante
1 1
2
(  y 2 )dx  2 (x 2  xy)dy  0
x y x y
y 1 1
 2
dx    dy  0
x y x
calculamos la derivadas parciales
M 1 N 1
 2  2
y x x x

Las derivadas parciales son iguales. Integramos a M con respecto a x

y y
F(x , y)    2 dx  g(y)   g(y)
x x
derivamos a F con respecto a y e igualamos con N
F 1 1 1
  g'(y)  
y x y x
despejamos g'(y) e integrando
1
g' (y) 
y
1
g(y)   dy  ln y
y
la solución es
y
 ln y  C
x

20
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

Factores integrantes que dependen únicamente de x ó de y

Si el factor integrante u depende solamente de x. Multiplicamos la ecuación por el factor
integrante u ( x )
u ( x )M ( x , y)dx  N ( x , y)dy   0
u ( x )M( x , y)dx  u ( x ) N( x , y)dy  0
calculamos las derivadas parciales
[u(x)M(x , y)] M(x , y)
 u(x)
y y
[u(x)N(x , y)] N(x, y) du(x)
 u(x)  N(x , y)
x x dx
igualamos las dos derivadas parciales, esto es para que la ecuación resultante sea exacta, se tiene
M(x , y) N(x, y) du(x)
u(x)  u(x)  N(x , y)
y x dx
M(x , y) N(x, y) du(x)
u(x)  u(x)  N(x , y)
y x dx
 M N  du
u     N
 y x  dx
M N

du y x
 u
dx N
o
M N

u' y x

u N
el término del lado derecho de la igualdad debe depender solamente de x. Entonces el factor de
integración es
M N

y x
 dx
ue N

Si u depende solamente de y. De una manera similar, se tiene que

N M

x y
M
depende solamente de y, entonces
N M

x y
 dy
ue M

es un factor integrante

21
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

Método de solución
Dada la ecuación diferencial
M ( x , y)dx  N( x.y)dy  0
no exacta (se debe verificar). Se busca un factor integrante comprobando que
M N

y x
N
depende únicamente de x (no aparece la variable y). O
N M

x y
M
depende únicamente de y (no aparece la variable x).

Se determina el factor integrante, u, se multiplica la ecuación

M( x , y)dx  N( x.y)dy  0
por el factor integrante, es decir,
uM ( x , y)dx  uN( x.y)dy  0
se comprueba que la nueva ecuación diferencial es exacta y se resuelve por el método de
ecuaciones exactas.

Ejemplo. Encontrar el factor integrante y la solución de las siguientes ecuaciones:

a) 6 xydx  (4 y  9 x 2 )dy  0
b) ( x  y)dx  x ln x dy  0
c) ( x 4 ln x  2 xy 3 )dx  3x 2 y 2 dy  0 , y(1) = 2
Solución
a) 6xydx  (4 y  9 x 2 )dy  0

Primero veamos si es exacta. Calculamos las derivadas parciales

M N
 6x  18x
y x
no son iguales. A priori, no sabemos si existe un factor integrante que dependa únicamente de x ó
dependa únicamente de y. Por lo que, veamos primero si existe un factor integrante que dependa
únicamente de x, se tiene
M N

y x 6 x  18x  12 x
 
N 4 y  9x 2
4 y  9x 2
vemos que el miembro del lado derecho de la igualdad no es una función que depende
únicamente de x. Por lo que, concluimos que no existe un factor integrante que depende
únicamente de x. Ahora veamos si existe un factor que depende solamente de y, se tiene
N M

x y 18x  6x 12x 2
  
M 6xy 6 xy y
vemos que el miembro del lado derecho de la igualdad si es una función que depende únicamente
de y. Por lo que, concluimos que si existe un factor integrante que depende únicamente de y, éste
22
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

es
2
 y dy 2
u(y)  e  e 2lny  e lny  y 2
multiplicamos la ecuación por el factor integrante encontrado y comprobamos que es exacta la
ecuación resultante
y 2 (6 xy)dx  y 2 (4 y  9 x 2 )dy  0
6 xy3dx  (4 y3  9 x 2 y 2 )dy  0
comprobamos que es exacta
M N
 18xy 2  18xy 2
y x
la ecuación es exacta. Integramos N con respecto a y
F(x , y)   (4y 3  9x 2 y 2 )dy  h(x)  y 4  3x 2 y 3  h(x)
derivamos F con respecto a x e igualándola con M
F
 6xy 3  h' (x)  6xy 3
x
despejamos h ' ( x ) e integrándola
h' (x)  0, h(x)  0
y la solución es
y 4  3x 2 y3  C

b) ( x  y)dx  x ln x dy  0

Primero veamos si es exacta. Calculamos las derivadas parciales

M N
1  lnx  1
y x
no son iguales.

Veamos primero si existe un factor integrante que dependa únicamente de x, se tiene

M N

y x 1  ( ln x  1) ln x 1
  
N x ln x x ln x x
vemos que el miembro del lado derecho de la igualdad es una función que depende únicamente
de x. Por lo que, concluimos que si existe un factor integrante que depende únicamente de x, éste
es
1
  dx 1 1
u(x)  e x  e ln x  e ln x  x 1 
x
multiplicamos la ecuación por el factor integrante y comprobamos que es exacta la ecuación
resultante
1 1
(x  y)dx  (xlnx)dy  0
x x
 y
1  dx  lnx dy  0
 x

23
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

comprobamos que es exacta

M 1 N 1
 , 
y x x x
la ecuación es exacta. Integramos M con respecto a x
 y
F(x , y)   1  dx  g(y)  x  ylnx  g(y)
 x
derivamos F con respecto a y e igualamos con N
F
 lnx  g' (y)  lnx
y
g' (y)  0, g(y)  0
y la solución es
x  y ln x  C

c) ( x 4 ln x  2xy 3 )dx  3x 2 y 2 dy  0

Primero veamos si es exacta. Calculamos las derivadas parciales

M  ( x 4 ln x  2 xy 3 ) N  (3x 2 y 2 )
  6 xy 2 ,   6xy 2
y y x x
no son iguales. Veamos primero si existe un factor integrante que dependa únicamente de x, se
tiene
M N

y x  6 xy 2  6 xy 2  12 xy 2 4
 2 2
 2 2

N 3x y 3x y x
vemos que la expresión del lado derecho de la igualdad si es una función que depende
únicamente de x. Por lo que, concluimos que si existe un factor integrante que depende
únicamente de x, éste es
4
  dx 4
u ( x )  e x  e 4 ln x  e ln x  x 4
multiplicamos la ecuación por el factor integrante y comprobamos que es exacta la ecuación
resultante
x 4 ( x 4 ln x  2xy 3 )dx  x 4 (3x 2 y 2 )dy  0
(ln x  2x 3 y3 )dx  3x 2 y 2dy  0
comprobamos que es exacta
M  (ln x  2 x 3 y 3 ) N  (3x 2 y 2 )
  6 x 3 y 2 ,   6 x  3 y 2
y y x x
la ecuación es exacta. Integramos M con respecto a x
F( x, y)   (ln x  2x 3 y 3 )dx  g ( y)
F( x, y)  x ln x  x  x 2 y 3  g( y)  C
derivamos F con respecto a y e igualamos con N
F
 3x 2 y 2  g' ( y)  3x 2 y 2
y
despejamos g' ( y) e integramos
g ' ( y)  0 , g ( y)  0

24
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

sustituimos g(y) en F(x , y)

x ln x  x  x 2 y 3  C

Aplicamos la condición inicial

(1) ln(1)  (1)  (1) 2 (2) 3  C
C7
la solución es
x ln x  x  x 2 y 3  7
Tarea 6
1.2.4 Ecuaciones diferenciales lineales.

Ecuaciones lineales
La forma general de una ecuación diferencial ordinaria lineal de primer orden es:
a0 ( x ) y' a1 ( x ) y  g ( x )
o
y' p( x ) y  f ( x ) (1.5)

Método de solución
Es necesario tener la ecuación diferencial en la forma
y' p( x ) y  f ( x )
para poder identificar claramente las funciones p( x ) y f ( x ) . Observe que el coeficiente de y'
debe ser 1.

La solución de (1.5) en un intervalo en el cual p(x) y f(x) son continuas es

y(x) 
1
u(x) 

u(x)f(x)dx  C 
donde
u(x)  e 
p(x)dx

Ejemplo. Resolver las siguientes ecuaciones:

dy
a) x  2y  3
dx
b) (1  x 2 )dy  ( xy  x 3  x )dx  0
c) ydx  ( xy  2x  yey )dy  0
d) y' y cos x  senx cos x , y(0) = 1
Solución
dy
a) x  2y  3
dx

Siempre tenemos que tener la ecuación diferencial de la forma requerida y' p( x ) y  f ( x ) (o

x ' p( y) x  f ( y) ), es decir el coeficiente de y' (o x') es uno. Observamos que el coeficiente de y'

25
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

es x, dividimos la ecuación entre x

2 3
y'  y
x x
la función u es
2
 x dx 2
ue  e 2lnx  e lnx  x 2
entonces
y(x) 
1  23
2 
x 
x
x
 1
dx  C  2
 x
 3xdx  C x1  32 x
2
2 
 C

la solución es
3
y(x)   Cx -2
2

b) (1  x 2 )dy  ( xy  x 3  x )dx  0

Reescribimos la ecuación en la forma deseada, dividimos la ecuación entre dx

dy
(1  x 2)  xy  x 3  x  0
dx
(1  x )y'  xy   x 3  x
2

dividimos la ecuación entre x 2  1

x x(x 2  1)
y'  y  
1 x2 x2 1
x
y'  y  x
1 x2
la función u es
x 1 1 1
 1  x 2 dx ln(1  x 2 )
ln(1  x 2 ) 2
ue e 2
e  (1  x 2 ) 2
entonces
1 
1
 1  1 2
3

1 
y(x)   (1  x 2 2
) xdx  C   1 
 (1  x 2 2
)  C
(1  x) 2   (1  x) 2  2 3 
la solución es
1
1 -
y(x)   (1  x )  C(1  x )
2 2 2

c) ydx  ( xy  2 x  yey )dy  0

Reescribimos la ecuación en la forma requerida, dividimos por dy

dx
y  (y  2)x  yey
dy
dividimos por y
y2
x'  x  ey
y
note que en esta ecuación se toma como variable dependiente a x, es decir, x es la función

26
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

desconocida. La función u es
y2
 dy 2
ue y
 e y  2lny  e y  lny  y 2 e y
la solución es

ye
1
x(y)  2 y
 y e e dy  C y 1e  y e
2 y y
2 y
2 2y
dy  C 
integramos por partes, tomamos
u dv
y2 e2 y
2y 1 2y
e
2
2 1 2y
e
4
0 1 2y
e
8
la solución es
1 ey ey C
x(y)  e y   2 2 y
2 2y 4y ye

d) y' y cos x  senx cos x

La función u es
u  e
cos xdx
 e senx
la solución es
y( x ) 
e
1
senx
 e senx
senx cos xdx  C 
integramos por partes
u  senx dv  e senx cos xdx
du  cos xdx v  e senx

y( x )  e senx e senx senx   e senx cos xdx  C 

 e senx e senx senx  e senx  C 
y( x )  senx  1  Ce senx

Aplicamos la condición inicial

y(0)  sen (0)  1  Ce sen ( 0)  1
1  C  1
C2
la solución es
y( x )  senx  1  2e senx
Tarea 27

27
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

1.2.5 Ecuaciones diferenciales de Bernoulli.

Ecuación de Bernoulli
Definición 1.11 Una ecuación diferencial ordinaria de la forma
y' p( x ) y  f ( x ) y n
donde n es cualquier número real, se llama ecuación de Bernoulli.

Si n  0 o n  1 , entonces la ecuación de Bernoulli es una ecuación diferencial ordinaria lineal.

Método de solución

Siempre debemos tener la ecuación diferencial en l forma y' p( x ) y  f ( x ) y n (o

x ' p( y) x  f ( y) x n ), es decir, el coeficiente de y' (o x ' ) es uno.

Si n  0 y n  1 , entonces el cambio de variable w  y1 n reduce la ecuación de Bernoulli en

una ecuación diferencial ordinaria lineal de la forma
w '(1  n )p( x ) w  (1  n )f ( x )
se resuelve la ecuación diferencial para w y se regresa a la variable y, sustituimos w por y1 n .

Ejemplo. Resolver las siguientes ecuaciones

dy
a)  y  ex y2
dx
dy
b) 3(1  x 2 )  2xy( y 3  1)
dx
1 3
dy
c) y2  y 2  1 , y(0) = 4
dx
Solución
dy
a)  y  ex y2
dx

Siempre tenemos que tener la ecuación diferencial de la forma requerida y' p( x ) y  f ( x ) y n (o

x ' p( y) x  f ( y) x n ), es decir el coeficiente de y' (o x') es uno. El signo negativo de (  y ) lo
podemos considerar mentalmente, pero escribiremos la ecuación en la forma
y'(1) y  e x y 2

De la ecuación diferencial vemos que n  2 , entonces 1  n  1  2  1 . Hacemos el cambio de

variable w  y1 2  y 1 , da

w '(1) w  (1)e x
w ' w  e x

28
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

la función u es
u  e
dx
 ex
la solución de w es

e
1
    1
 2


 1
w  x  e x (  e x )dx  C  e-x   e 2 x dx  C  e -x  e 2 x  C    Ce-x
2
regresamos a la variable y, se tiene
1
y 1   e x  Ce  x
2
finalmente, la solución es
1
y
1
 e x  Ce  x
2

dy
b) 3(1  x 2 )  2xy( y 3  1)
dx

Reescribimos la ecuación para tener la forma requerida

dy
3(1  x 2)  2xy 4  2xy
dx
3(1  x )y'  2xy  2xy 4
2

2x 2x
y'  y y4
3(1  x )
2
3(1  x )
2

De la ecuación diferencial vemos que n  4 , entonces 1  n  1  4  3 . Hacemos el cambio de

variable w  y1 4  y 3
2x 2x
w'  (  3) w  (  3)
3(1  x )
2
3(1  x 2)
2x 2x
w'  w
1 x 2
1  x2
la función u es
2 x
 dx 2 2 1
ue 1 x2
 e ln(1  x )  e ln(1  x )
 (1  x 2 ) 1
la solución en w es
1  2 1  2x   2  2x 
(1  x 2 )-1   
w  (1  x )   2  dx  C   (1  x )   2 2
dx  C 
 1 x    (1  x ) 
 1 
w  (1  x 2 )  2
 C  1  (1  x 2 )C
1  x 
regresamos a la variable y
y 3  1  (1  x 2 )C
finalmente, la solución es


y  1  (1  x 2 )C 
1 / 3

29
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

1 3
dy
c) y2  y 2 1
dx

Reescribimos la ecuación
1 3
y y' y  1
2 2

1

y' y  y 2

1  1 3
De la ecuación diferencial vemos que n   , entonces 1  n  1      . Hacemos el
2  2 2
 1 3
1   
cambio de variable w  y  2
 y 2 , da
3 3
w ' w 
2 2
la función u es
3 3
 2 dx x
ue  e2
la solución en w es
1  3 2x   x 
3 3 3 3
x  x
3 
w e dx  C   e 2
 e 2
 C   1  Ce 2
2
e2    
x

regresamos a la variable y
3 3
 x
y 2  1  Ce 2

aplicamos la condición inicial

3 3
 ( 0)
(4) 2  1  Ce 2
 1 C
C7
sustituimos la constante
3 3
 x
y  1  7e
2 2

finalmente, la solución es
2
  x 3
3
y  1  7e 2 
 
Tarea 8
1.3 Aplicaciones de ecuaciones diferenciales de primer orden.

1.3.1 Trayectorias ortogonales.

Definición 1.12 Dos curvas se dicen ortogonales en un punto de intersección si y solo si sus
rectas tangentes son perpendiculares en el punto de intersección.

Esto significa que la pendiente de la recta tangente de una curva es la recíproca negativa de la
otra, a excepción que las rectas tangentes sean paralelas a los ejes coordenados.

30
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

Definición 1.13 Cuando todas las curvas de una familia de curvas G ( x , y , C1 )  0 intersecan
ortogonalmente a todas las curvas de otra familia H( x , y , C 2 )  0 , se dice que
las familias son, cada una, trayectorias ortogonales a la otra.

dy
Supóngase que  f ( x, y) es la ecuación diferencial que describe a una familia de curvas dada.
dx
Entonces la ecuación diferencial de la familia ortogonal es
dy 1

dx f ( x , y)
Ejemplo. Halle las trayectorias ortogonales para la familia
a) cos y  Ce  x
b) 2x 2  y 2  4C1 x
Solución
a) cos y  Ce  x

Hallamos la ecuación diferencial de la familia, para hacer esto se deriva implícitamente con
respecto a x la expresión dada
dy
 seny  Ce  x
dx
despejamos C de la expresión de la familia
cos y
C   x  e x cos y
e
sustituimos C en la derivada
dy
 seny  (e x cos y)e  x   cos y
dx
dy
seny  cos y
dx
dy cos y

dx seny
la ecuación diferencial de las trayectorias ortogonales es
dy 1 seny
 
dx cos y cos y
seny
dy seny

dx cos y
resolvemos la ecuación diferencial por el método de variables separables
cos y
dy  dx
seny
cos y
 seny dy   dx  C1
ln seny   x  C1

31
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

La familia ortogonal es
ln seny   x  C1

b) 2x 2  y 2  4C1x

Hallamos la ecuación diferencial de la familia, para hacer esto se deriva implícitamente con
respecto a x la expresión dada
dy
4x  2 y  4C1
dx
sustituimos 4C1 en la ecuación de la familia dada ( 2x 2  y 2  4C1 x )
 dy 
2x 2  y 2   4x  2 y  x
 dx 
despejamos la derivada de y con respecto a x
2x 2  y 2 dy
 4x  2 y
x dx
2 2
dy 2x  y 2x 2  y 2  4x 2 y 2  2x 2
2y   4x  
dx x x x
2 2
dy y  2x

dx 2xy
la ecuación diferencial de la familia de trayectorias ortogonales es
dy 2xy
 2
dx y  2x 2
esta ecuación diferencial es homogénea
( y 2  2 x 2 )dy  2xydx
realizamos el cambio de variable
x  uy, dx  udy  ydu
entonces
( y 2  2(uy) 2 )dy  2(uy) y(udy  ydu)
( y 2  2u 2 y 2 )dy  2u 2 y 2dy  2uy3du
y 2dy  2uy3du
1
dy  2udu
y
integramos
1
 y dy    2udu  C
ln y   u 2  C
x
sustituimos u 
y
2
x
ln y     C
 y

32
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

la familia ortogonal es
2
x
ln y     C
 y
Tarea 9
1.3.2 Decaimiento radiactivo.

Si una molécula tiene la tendencia a descomponerse espontáneamente en moléculas menores, a

un ritmo no afectado por la presencia de otras sustancias. Entonces es natural que el número de
 dA 
moléculas a descomponerse en una unidad de tiempo   será proporcional al número total
 dt 
presente (A(t)). Una reacción química de este tipo se denomina reacción de primer orden, cuya
ecuación diferencial es:
dA
 kA
dt
La vida media o semivida es una medida de la estabilidad de una sustancia radiactiva. La vida
media es simplemente el tiempo necesario para que se desintegre la mitad de los átomos de una
cantidad inicial A0.

Una aplicación de este modelo de decaimiento es la que encontró Willard Libby para la
determinación de edades de los fósiles usando el carbono radiactivo. La teoría se basa en que el
isótopo carbono 14 se produce en la atmósfera por la acción de la radiación cósmica sobre el
nitrógeno. El cociente de la cantidad de C-14 y la cantidad de carbono ordinario presentes en la
atmósfera es constante, y en consecuencia, la proporción de isótopo presente en todos los
organismos vivos es la misma que en la atmósfera. Cuando un organismo muere, la absorción de
C-14 cesa. Así, comparando la proporción de C-14 que hay en un fósil con la proporción
constante encontrada en la atmósfera es posible obtener una estimación razonable de su edad. El
método se basa en que la semivida del C-14 radiactivo es de aproximadamente 5600 años.

Ejemplo. La vida media del carbono C-14 es alrededor de 5600 años. Se encuentra que un hueso
fosilizado contiene una milésima de la concentración de C-14 que se encuentra en la materia
viva. a) Determine el valor de la constante de proporcionalidad de la ecuación diferencial. b)
Estime la edad del fósil.
Solución
La ecuación diferencial es
dA
 kA
dt
resolvemos la ecuación diferencial
dA
 kdt
A
dA
 A   kdt
ln A  kt  C
A( t )  e kt  C  Ce kt
suponemos que inicialmente la cantidad de carbono C-14 es A(0) = A0, tenemos
A(0)  Ce k ( 0 )  A 0

33
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

C  A0
A( t )  A 0e kt

a) Consideramos la vida media del carbono C-14

1
A(5600)  A 0e k ( 5600)  A0
2
1
e k ( 5600) 
2
1
k (5600)  ln
2
1
ln
k 2
5600
k  0.000123776

b) Para estimar la edad del fósil

A( t )  A 0e 0.000123776 t
A( t )  A 0e  0.000123776 t  0.001A 0
e  0.000123776 t  0.001
 0.000123776t  ln 0.001
ln 0.001
t
 0.000123776
t = 55808.52 años
Tarea 10

1.3.3 Modelos de población.

Consideremos una población homogénea aislada. La rapidez con que la población en cada
instante se incrementa depende del número de habitantes en el instante, es decir,
dN
 f ( N( t ))
dt
dos de estos modelos son:
dN
 kN
dt
donde k es una constante
dN
 N(a  bN) ecuación logística
dt
donde a y b son constantes.

Ejemplo. Un cultivo tiene inicialmente una cantidad N 0 de bacterias. Para t  1 hora, el número
de bacterias es 2 N 0 . Si la rapidez de crecimiento es proporcional al número de bacterias
presente, determine el tiempo necesario para que el número de bacterias sea el triple de su
cantidad inicial.

34
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

Solución. El modelo de crecimiento para este ejemplo es

dN
 kN
dt
sujeta a las condiciones N(0) = N0 y N(1) = 2N0

Resolvemos la ecuación diferencial (variables separables), se tiene

dN
 kdt
N
integramos
ln( N)  kt  c
N( t )  e kt  c  ec e kt  Ce kt
donde C  ec . Aplicamos las condiciones, se tiene
N(0)  Ce k ( 0 )  C  N 0
N( t )  N 0e kt
N(1)  N 0e k (1)  2 N 0
e k (1)  2
k  ln 2  0.693
entonces la solución es
N( t )  N 0e0.693t

Para determinar el valor de t para el que las bacterias se triplican, despejamos t de

N( t )  N 0e0.693t  3N 0
e0.693t  3
obtenemos el logaritmo, se tiene
0.693t  ln 3  1.0986
1.086
t
0.693
t = 1.58 horas

Ejemplo. El número N(t) de personas en una comunidad que están expuestas a un anuncio
particular se rige mediante la ecuación logística. Al inicio, N(0)  500 , y se observa que
N(1)  1000 . Determine N(t) si se predice que el número límite de personas en la comunidad que
verán el anuncio es 50000.
Solución
La ecuación logística es
dN
 N(a  bN)
dt
donde a y b son constantes a determinar. Y esta ecuación es de Bernoulli, la escribiremos de
forma diferente para visualizar la ecuación de Bernoulli
N'  Na  bN 2
N'aN  bN 2

35
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

Observamos que: n  2 , 1  n  1  2  1 . Realizando el cambio de variable w  N 1 ,

obtenemos la ecuación
w '(1)aw  (1)b
w ' aw  b
la función u es
u  e   e at
adt

w es
w
e
1

at 
 b 
e at bdt  C  e -at  e at  C
a 
b
w   Ce -at
a

La solución de la ecuación logística es

1 b
  Ce -at
N( t ) a

Aplicamos las condiciones a la frontera

1 b 1
  Ce -a ( 0)   0.002
N(0) a 500
b
 C  0.002 (A)
a
1 b 1
  Ce -a (1)   0.001
N(1) a 1000
b
 Ce -a  0.001 (B)
a
1 b  1
lím  lím   Ce -at    0.00002
t  N ( t )

t  a
 50000
b
 0.00002 (C)
a

Con las ecuaciones (A), (B) y (C) formamos el sistema

b
 C  0.002 (A)
a
b
 Ce -a  0.001 (B)
a
b
 0.00002 (C)
a
b
sustituimos  0.00002 en la ecuación (A)
a
0.00002  C  0.002
C  0.002  0.00002
C  0.00198

36
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

b
sustituimos  0.00002 y C = 0.00098 en la ecuación (B)
a
0.00002  0.00198e -a  0.001
0.00198e -a  0.001  0.00002  0.00098
0.00098
e -a 
0.00198
 0.00098 
 a  ln 
 0.00198 
 0.00098 
a   ln 
 0.00198 
a  0.703
b
sustituimos  0.00002 , C = 0.00098 y a = 0.703 en el número de personas N(t)
a
1
 0.00002  0.00198e -0.703 t
N( t )
entonces
1
N( t ) 
0.00002  0.00198e -0.703t
Tarea 11
1.3.4 Ley de enfriamiento de Newton.

De acuerdo con la ley de empírica de enfriamiento de Newton –o de calentamiento–, la rapidez

con que la temperatura de un cuerpo cambia es proporcional a la diferencia que hay entre la
temperatura del cuerpo y la del medio que la rodea, la denominada temperatura ambiental. Si
T( t ) representa la temperatura de un cuerpo en el instante t, Tm la temperatura del medio que lo
dT
rodea y la rapidez a la que cambia la temperatura del cuerpo, la ley de Newton de
dt
enfriamiento y calentamiento se traduce en el enunciado matemático
dT
 (T  Tm )
dt
o
dT
 k (T  Tm )
dt
donde k es una contante de proporcionalidad. En cualquier caso, calentamiento o enfriamiento, si
Tm es una constante, se sostiene que k  0 .

Ejemplo. La temperatura de una taza de café es de 95C cuando se acaba de servir. Un minuto
después se enfrió a 88C en un cuarto que está a 20C . ¿Cuánto demorará para que el café
alcance una temperatura de 65C ?
Solución. Se debe resolver el problema de valores iniciales
dT
 k (T  20)
dt
T(0)  95

37
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

y determinar el valor de k de modo que T(1)  88 .

La ecuación diferencial es separable, separamos las variables

dT
 kdt
T  20
integramos
dT
 T  20   kdt  C
ln(T  20)  kt  C
T  20  e kt  C  Ce kt
T( t )  Ce kt  20
consideremos la condición inicial
T(0)  Ce k ( 0)  20  95
C  75
T( t )  75e kt  20

Consideramos la condición T(1)  88

T(1)  75e k (1)  20  88
75e k  68
68
ek 
75
 68 
k  ln   0.098
 75 

La solución es
T( t )  75e 0.1098 t  20

Ahora determinaremos el tiempo necesario para que T( t )  65

T( t )  75e 0.1098 t  20  65
e 0.1098 t  45
45 3
e 0.1098 t 
75 5
3
 0.1098t  ln 
5
1 3
t ln   5.2
 0.1098  5 

El tiempo necesario es
t = 5.2 min
Tarea 12

38
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

1.3.5 Mezclas.

Considere un tanque que inicialmente

contiene una cantidad A 0 de sustancia
soluble disuelta homogéneamente en un
volumen V0 de sustancia solvente.
Suponga que entra al tanque sustancia
solvente que contiene Ce unidades de
sustancia soluble por unidad de sustancia Figura 1.1
solvente, a una rapidez re de unidades de sustancia solvente por unidad de tiempo. Y está
saliendo del tanque la solución bien homogeneizada a una rapidez rs unidades de volumen por
unidad de tiempo.

Consideremos que A( t ) es la cantidad de sustancia soluble y V( t ) es el volumen en el tanque, en

el instante t.

La rapidez de cambio de sustancia soluble en el tanque, en el instante t, A' ( t ) ; debe ser igual a la
rapidez con que la sustancia entra, R e ( t ) , menos la rapidez con la que sale, R s ( t ) , es decir,
A' ( t )  R e ( t )  R s ( t ) (1.6)
por otro lado
A
R e ( t )  C e re y R s (t)  rs
V
sustituimos en (1.6)
A
A'  C e re  rs
V
o
rs
A' A  C e re (1.7)
V
y el volumen en el instante t está dado por
V( t )  V0  (re  rs ) t (1.8)
sustituimos en (1.7), se tiene
rs
A' A  Ce re (1.9)
V0  (re  rs ) t

existen dos casos: a) volumen constante y b) volumen variable.

a) Para mantener el volumen constante, se requiere que rs  re . Entonces la ecuación (1.9) queda
como
r
A' s A  Ce re
V0

39
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

b) Como el volumen es variable, entonces rs  re . Y en este caso la ecuación es

rs
A' A  Ce re
V0  (re  rs ) t
la ecuación anterior es valida hasta que el nivel del liquido no supere la capacidad del tanque ó
que el tanque quede vacío.

Ejemplo. Un depósito grande se llena al máximo con 1500 litros de agua pura. Se bombea al
depósito salmuera que contiene 0.3 kg de sal por cada 1 litro a razón de 15 litros/min. La solución
bien mezclada se bombea hacia fuera a una rapidez de 30 litros/min. a) ¿Cuándo se vacía el
depósito? b) Calcule el número A(t) de kg de sal en el depósito en el tiempo t. c) Calcule la
concentración de sal en el depósito en el tiempo t. d) ¿En qué momento la concentración de la sal
en el depósito es igual a 0.15 kg por litro?
Solución
a) Como la rapidez de salmuera que entra, re, no es igual a la rapidez con que sale, rs. El volumen
de salmuera en el depósito es
V( t )  V0  (re  rs ) t (A)
igualamos a cero el volumen y despejamos el tiempo
V( t )  V0  (re  rs ) t  0
V0  (re  rs ) t  0
V
t 0
re  rs
sustituimos valores
1500
t
15  30
t = 100 min.

b) Como la rapidez de salmuera que entra, re, no es igual a la rapidez con que sale, rs. La ecuación
diferencial es
rs
A' A  Ce re
V0  (re  rs ) t
sustituimos valores
30
A' A  (0.3)(15)
1500  (15  30) t
30
A' A  4.5
1500  15t
2
A' A  4.5
100  t
resolvemos la ecuación diferencial
2
 100  t dt 2
ue  e  2 ln(100  t )  eln(100  t )  (100  t )  2

A( t ) 
1
(100  t )  2
 (100  t) 2
4.5dt  C 
 4.5 
A( t )  (100  t ) 2   C
100  t 

40
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

A( t )  4.5(100  t )  C(100  t ) 2

Como inicialmente el depósito contiene agua pura, es decir, A(0) = 0, tenemos

A(0)  4.5(100  0)  C(100  0) 2  0
450  10000C  0
450 9
C 
10000 200
9
A( t )  4.5(100  t )  (100  t ) 2
200

c) La concentración de sal está dada por

A( t )
C( t ) 
V( t )
sustituimos A(t) del inciso (b) y V(t) dado por la expresión (A) y sustituimos valores
9 9
4.5(100  t )  (100  t ) 2 4.5(100  t )  (100  t ) 2
C( t )  200  200
Vo  (re  rs ) t 1500  (15  30) t
9 9
4.5(100  t )  (100  t ) 2 4.5(100  t )  (100  t ) 2
 200  200
1500  15t 15(100  t )
4.5 9
C( t )   (100  t )
15 200(15)
3
C( t )  0.3  (100  t )
1000

d) El tiempo en que la concentración es igual a 0.15 kg es

3
C( t )  0.3  (100  t )  0.15
1000
3
(100  t )  0.15
1000
1000
100  t  (0.15)  50
3
t  100  50  50
t = 50 min.
Tarea 13
1.3.6 Mecánica elemental.

Supongamos que un cuerpo de masa m cae libremente, tan solo bajo la

acción de la gravedad, como se ilustra en la figura 1.2. Consideremos un
sistema de referencia con el origen en una posición arbitraria en un punto
arriba del cuerpo y el eje y positivo hacia abajo, como se ilustra en la
figura 1.3, la distancia recorrida por el cuerpo es positiva, su velocidad
hacia abajo es positiva, así como la gravedad. Supongamos que en el
instante t0 el cuerpo está en una posición y0 con una velocidad v0. En este
Figura 1.2

41
Víctor Gerardo Serrano Domínguez Apuntes de Ecuaciones Diferenciales

caso, la única fuerza que actúa sobre el cuerpo es mg. Su aceleración será
d2y
y por la segunda ley de Newton (F = ma), se tiene
dt 2
d2y d2y
mg  m 2 o g  2
dt dt
ésta no es una ecuación de primer orden. Pero realizamos el cambio de
variable Figura 1.3
dy dv d 2 y
v , 
dt dt dt 2
entonces
dv
g
dt
que es un ecuación de primer orden de variables separables. Separamos las variables
gdt  dv
t v

 gdt   dv
t0 v0

( t  t 0 )g  v  v 0
v  ( t  t 0 )g  v 0
como
dy
v
dt
dy
 ( t  t 0 )g  v 0
dt
dy  [( t  t 0 )g  v 0 ]dt
y t

 dy   [(t  t 0 )g  v0 ]dt
y0 t0

1
y  y0  ( t  t 0 )2 g  v0 ( t  t 0 )
2
1
y  (t  t 0 )2 g  v0 (t  t 0 )  y0
2

42