ECUACIONES DIFERENCIALES

2 DE PRIMER ORDEN

2.1 Curvas solución sin una solución

2.1.1 Campos direccionales
2.1.2 ED de primer orden autónomas
2.2 Variables separables
2.3 Ecuaciones lineales
2.4 Ecuaciones exactas
2.5 Soluciones por sustitución
2.6 Un método numérico
REPASO DEL CAPÍTULO 2

La historia de las matemáticas tiene muchos relatos de personas que han dedicado
gran parte de su vida a la solución de ecuaciones, al principio de ecuaciones
algebraicas y después de ecuaciones diferenciales. En las secciones 2.2 a 2.5
estudiaremos algunos de los métodos analíticos más importantes para resolver
ED de primer orden. Sin embargo, antes de que empecemos a resolverlas,
debemos considerar dos hechos: es posible que una ecuación diferencial no tenga
soluciones y que una ecuación diferencial tenga una solución que con los
métodos existentes actuales no se puede determinar. En las secciones 2.1 y 2.6
no resolveremos ninguna ED pero mostraremos cómo obtener información
directamente de la misma ecuación. En la sección 2.1 podemos ver cómo, a partir
de la ED, obtenemos información cualitativa de la misma respecto a sus gráficas,
lo que nos permite interpretar los dibujos de las curvas solución. En la sección 2.6
usamos ecuaciones diferenciales para construir un procedimiento numérico para
soluciones aproximadas.

Página 1 de 48

08367_02_ch02_p034-081-ok.indd 34 6/4/09 12:15:44 PM

 2.1 CURVAS SOLUCIÓN SIN UNA SOLUCIÓN
REPASO DE MATERIAL
● La primera derivada como pendiente de una recta tangente.

● El signo algebraico de la primera derivada indica crecimiento o decrecimiento.

INTRODUCCIÓN Imaginemos por un momento que nos enfrentamos con una ecuación diferencial
de primer orden dy
dx
f (x, y), y que además no podemos encontrar ni inventar un método para
resolverla analíticamente. Esto no es tan malo como se podría pensar, ya que la ecuación diferencial
en sí misma a veces puede “decirnos” concretamente cómo se “comportan” sus soluciones.
Iniciaremos nuestro estudio de las ecuaciones diferenciales de primer orden con dos formas
cualitativas de analizar una ED. Estas dos formas nos permiten determinar, de una manera aproximada,
cómo es una curva solución sin resolver realmente la ecuación.

2.1.1 CAMPOS DIRECCIONALES

ALGUNAS PREGUNTAS FUNDAMENTALES En la sección 1.2 vimos que si f (x,
y) y f
y satisfacen algunas condiciones de continuidad, se pueden responder preguntas
cualitativas acerca de la existencia y unicidad de las soluciones. En esta sección veremos
otras preguntas cualitativas acerca de las propiedades de las soluciones. ¿Cómo se com-
porta una solución cerca de un punto dado? ¿Cómo se comporta una solución cuando x
: ? Con frecuencia, estas preguntas se pueden responder cuando la función f depende
sólo de la variable y. Sin embargo, comenzaremos con un simple concepto de cálculo:
Una derivada dy
dx de una función derivable y
y(x) da las pendientes de las
y rectas tangentes en puntos de su gráfica.
pendiente = 1.2
PENDIENTE Debido a que una solución y
y(x) de una ecuación diferencial de
primer orden
(2, 3)
dy

f (x, y) (1)
dx
es necesariamente una función derivable en su intervalo I de definición, debe también
ser continua en I. Por tanto la curva solución correspondiente en I no tiene cortes y debe
x tener una recta tangente en cada punto (x, y(x)). La función f en la forma normal (1) se
llama función pendiente o función razón. La pendiente de la recta tangente en (x, y(x))
en una curva solución es el valor de la primera derivada dy
dx en este punto y sabemos
a) elemento lineal en un punto. de la ecuación (1) que es el valor de la función pendiente f (x, y(x)). Ahora suponga-
mos que (x, y) representa cualquier punto de una región del plano xy en la que está
y
curva definida la función f. El valor f (x, y) que la función f le asigna al punto representa la
solución pendiente de una recta o que la visualizaremos como un segmento de recta llamado
elemento lineal. Por ejemplo, considere la ecuación dy
dx
0.2xy, donde f (x, y)

(2, 3) 0.2xy. En el punto (2, 3) la pendiente de un elemento lineal es f (2, 3)
0.2(2)(3)

1.2. La figura 2.1.1a muestra un segmento de recta con pendiente 1.2 que pasa por
tangente (2, 3). Como se muestra en la figura 2.1.1b, si una curva solución también pasa por el
punto (2, 3), lo hace de tal forma que el segmento de recta es tangente a la curva; en otras
palabras, el elemento lineal es una recta tangente miniatura en ese punto.
x
CAMPO DIRECCIONAL Si evaluamos sistemáticamente a f en una malla rectan-
gular de puntos en el plano xy y se dibuja un elemento lineal en cada punto (x, y) de la
b) el elemento lineal es tangente malla con pendiente f (x, y), entonces al conjunto de todos estos elementos lineales se
a la curva solución que
le llama campo direccional o campo de pendientes de la ecuación diferencial dy
dx
pasa por el punto.

f (x, y). Visualmente, la dirección del campo indica el aspecto o forma de una familia
FIGURA 2.1.1 El elemento lineal es de curvas solución de la ecuación diferencial dada y, en consecuencia, se pueden ver
tangente a la curva solución en (2, 3). a simple vista aspectos cualitativos de la solución, por ejemplo, regiones en el plano

Página 2 de 48

08367_02_ch02_p034-081-ok.indd 35 6/4/09 12:15:44 PM

 en las que una solución presenta un comportamiento poco común. Una sola curva
solución que pasa por un campo direccional debe seguir el patrón de flujo del campo:
el elemento lineal es tangente a la curva cuando intercepta un punto de la malla. La
figura 2.1.2 muestra un campo direccional generado por computadora de la ecuación
diferencial dy
dx
sen(x  y) en una región del plano xy. Observe cómo las tres cur-
vas solución que se muestran a colores siguen el flujo del campo.

EJEMPLO 1 Campo direccional

El campo direccional para la ecuación diferencial dy
dx
0.2xy que se muestra en la figura
FIGURA 2.1.2 Las curvas solución 2.1.3a se obtuvo usando un paquete computacional en el que se definió una malla 5  5 (mh,
siguen el flujo de un campo direccional. nh) con m y n enteros, haciendo – 5  m  5, 5  n  5, y h
1. Observe en la figura
2.1.3a que en cualquier punto del eje de las x (y
0) y del eje y (x
0), las pendientes son
y
f (x, 0)
0 y f (0, y)
0, respectivamente, por lo que los elementos lineales son horizontales.
Además observe que en el primer cuadrante para un valor fijo de x los valores de f (x, y)

4 0.2xy aumentan conforme crece y; análogamente, para una y los valores de f (x, y)

0.2xy aumentan conforme x aumenta. Esto significa que conforme x y y crecen, los ele-
2 mentos lineales serán casi verticales y tendrán pendiente positiva ( f (x, y)
0.2xy  0 para
x  0, y  0). En el segundo cuadrante,  f (x, y) aumenta conforme crecen x y y crecen,
x por lo que nuevamente los elementos lineales serán casi verticales pero esta vez tendrán
pendiente negativa ( f (x, y)
0.2xy 0 para x 0, y  0). Leyendo de izquierda a dere-
_2 cha, imaginemos una curva solución que inicia en un punto del segundo cuadrante, se
mueve abruptamente hacia abajo, se hace plana conforme pasa por el eje y y después,
_4
conforme entra al primer cuadrante, se mueve abruptamente hacia arriba; en otras palabras,
_4 _2 2 4 su forma sería cóncava hacia arriba y similar a una herradura. A partir de esto se podría
inferir que y :  conforme x :
. Ahora en el tercer y el cuarto cuadrantes, puesto que
a) Campo direccional para
f (x, y)
0.2xy  0 y f (x, y)
0.2xy 0, respectivamente, la situación se invierte: una
dy/dx
0.2xy.
curva solución crece y después decrece conforme nos movamos de izquierda a derecha.
Vimos en la ecuación (1) de la sección 1.1 que y
e0.1x es una solución explícita de
2
y
dy
dx
0.2xy; usted debería comprobar que una familia uniparamétrica de soluciones
de la misma ecuación está dada por: y
ce0.1x . Con objeto de comparar con la figura 2.1.3a,
2
4
c>0
en la figura 2.1.3b se muestran algunos miembros representativos de esta familia.
2

c=0 x EJEMPLO 2 Campo direccional

c<0
_2 Utilice un campo direccional para dibujar una curva solución aproximada para el pro-
blema con valores iniciales dy
dx
sen y, y(0)
 32.
_4

_4 _2 2 4 SOLUCIÓN Antes de proceder, recuerde que a partir de la continuidad de f (x, y)
sen y y

b) Algunas curvas solución
2
f
y
cos y el teorema 1.2.1 garantiza la existencia de una curva solución única que pase
en la familia y
ce 0.1x . por un punto dado (x0, y0) en el plano. Ahora nuevamente seleccionando en nuestro paquete
FIGURA 2.1.3 Campo direccional y computacional la opción para una región rectangular 5  5 y dando puntos (debidos a la
curvas solución. condición inicial) en la región con separación vertical y horizontal de 12 unidad, es decir,
en puntos (mh, nh), h
12 , m y n enteros tales como 10  m  10, 10  n  10. En
la figura 2.1.4 se presenta el resultado. Puesto que el lado derecho de dy
dx
sen y es 0
en y
0, y en y
p, los elementos lineales son horizontales en todos los puntos cuyas
segundas coordenadas son y
0 o y
p. Entonces tiene sentido que una curva solución
que pasa por el punto inicial (0, 32), tenga la forma que se muestra en la figura.

CRECIMIENTO/DECRECIMIENTO La interpretación de la derivada dy
dx como

una función que da la pendiente juega el papel principal en la construcción de un
campo direccional. A continuación se usará otra contundente propiedad de la primera
derivada, es decir, si dy
dx  0 (o dy
dx 0) para toda x en un intervalo I, entonces
una función derivable y
y(x) es creciente (o decreciente) en I.

Página 3 de 48

08367_02_ch02_p034-081-ok.indd 36 6/4/09 12:15:46 PM

 y
COMENTARIOS
4
Dibujar a mano un campo direccional es directo pero tardado; por eso es proba-
2 ble que en la vida solo una o dos veces se realice esta tarea, pero generalmente
es más eficiente realizarlo usando un paquete computacional. Antes de las calcu-
x
ladoras, de las computadoras personales y de los programas se utilizaba el mé-
todo de las isoclinas para facilitar el dibujo a mano de un campo direccional.
_2
Para la ED dy
dx
f (x, y), cualquier miembro de la familia de curvas f (x, y)
_4
c, donde c es una constante, se llama isoclina. Se dibujan elementos lineales
que pasen por los puntos en una isoclina dada, digamos, f (x, y)
c1 todos con la
_4 _2 2 4 misma pendiente c1. En el problema 15 de los ejercicios 2.1 tiene dos oportuni-
dades para dibujar un campo direccional a mano.
FIGURA 2.1.4 Campo direccional
del ejemplo 2.

2.1.2 ED DE PRIMER ORDEN AUTÓNOMAS

ED DE PRIMER ORDEN AUTÓNOMAS En la sección 1.1 dividimos la clase
de las ecuaciones diferenciales ordinarias en dos tipos: lineales y no lineales. Ahora
consideraremos brevemente otra clase de clasificación de las ecuaciones diferenciales
ordinarias, una clasificación que es de particular importancia en la investigación cua-
litativa de las ecuaciones diferenciales. Una ecuación diferencial ordinaria en la que la
variable independiente no aparece explícitamente se llama autónoma. Si el símbolo x
denota a la variable independiente, entonces se puede escribir una ecuación diferencial
autónoma de primer orden como f (y, y)
0 o en la forma normal como
dy

f (y). (2)
dx
Supondremos que la función f en la ecuación (2) y su derivada f  son funciones conti-
nuas de y en algún intervalo I. Las ecuaciones de primer orden
f ( y) f (x, y)
p p
dy dy

1  y2 y
0.2xy
dx dx
son respectivamente autónoma y no autónoma.
Muchas ecuaciones diferenciales que se encuentran en aplicaciones o ecuaciones
que modelan leyes físicas que no cambian en el tiempo son autónomas. Como ya
hemos visto en la sección 1.3, en un contexto aplicado, se usan comúnmente otros
símbolos diferentes de y y de x para representar las variables dependientes e indepen-
dientes. Por ejemplo, si t representa el tiempo entonces al examinar a
dA dx dT dA 1

kA,
kx(n  1  x),
k(T  Tm),
6 A,
dt dt dt dt 100
donde k, n y Tm son constantes, se encuentra que cada ecuación es independiente del
tiempo. Realmente, todas las ecuaciones diferenciales de primer orden introducidas en
la sección 1.3 son independientes del tiempo y por tanto son autónomas.

PUNTOS CRÍTICOS Las raíces de la función f en la ecuación (2) son de especial

importancia. Decimos que un número real c es un punto crítico de la ecuación dife-
rencial autónoma (2) si es una raíz de f, es decir, f (c)
0. Un punto crítico también
se llama punto de equilibrio o punto estacionario. Ahora observe que si sustituimos
la función constante y(x)
c en la ecuación (2), entonces ambos lados de la ecuación
son iguales a cero. Esto significa que:
Si c es un punto crítico de la ecuación (2), entonces y(x)
c es una solución
constante de la ecuación diferencial autónoma.
Una solución constante y(x)
c se llama solución de equilibrio; las soluciones de
equilibrio son las únicas soluciones constantes de la ecuación (2).

Página 4 de 48

08367_02_ch02_p034-081-ok.indd 37 6/4/09 12:15:47 PM

 Como ya lo hemos mencionado, podemos decir cuándo una solución no constante
y
y(x) de la ecuación (2) está creciendo o decreciendo determinando el signo alge-
braico de la derivada dy
dx; en el caso de la ecuación (2) hacemos esto identificando
los intervalos del eje y en los que la función f (y) es positiva o negativa.

EJEMPLO 3 Una ED autónoma

La ecuación diferencial

dP

P(a  bP),
dt
eje P donde a y b son constantes positivas, tiene la forma normal dP
dt
f (P), la de la ecua-
ción (2) con t y P jugando los papeles de x y y respectivamente y por tanto es autónoma.
a De f (P)
P(a – bP)
0 vemos que 0 y a
b son puntos críticos de la ecuación, así que
b
las soluciones de equilibrio son P(t)
0 y P(t)
a
b. Poniendo los puntos críticos en
una recta vertical, dividimos esta recta en tres intervalos definidos por  P 0, 0
0 P a
b, a
b P . Las flechas en la recta que se presenta en la figura 2.1.5 indican
el signo algebraico de f (P)
P(a – bP) en estos intervalos y si una solución constante
P(t) está creciendo o decreciendo en un intervalo. La tabla siguiente explica la figura:
FIGURA 2.1.5 Esquema de fase de
dP
dt
P(a  bP). Intervalo Signo de f (P) P(t) Flecha

(, 0) menos decreciente apunta hacia abajo

(0, a
b) más creciente apunta hacia arriba
(a
b, ) menos decreciente apunta hacia abajo

La figura 2.1.5 se llama un esquema de fase unidimensional, o simplemente

esquema de fase, de la ecuación diferencial dP
dt
P(a  bP). La recta vertical se
y llama recta de fase.
R

CURVAS SOLUCIÓN Sin resolver una ecuación diferencial autónoma, normalmen-

(x0, y0) te podemos decir gran cantidad de detalles respecto a su curva solución. Puesto que
I la función f en la ecuación (2) es independiente de la variable x, podemos suponer
que f está definida para  x  o para 0  x . También, puesto que f y su
derivada f  son funciones continuas de y en algún intervalo I del eje y, los resultados
x principales del teorema 1.2.1 valen en alguna franja o región R en el plano xy corres-
pondiente a I, y así pasa por algún punto (x0, y0) en R por el que pasa una curva solución
de la ecuación (2). Véase la figura 2.1.6a. Para realizar nuestro análisis, supongamos
a) región R. que la ecuación (2) tiene exactamente dos puntos críticos c1 y c2 y que c1 c2. Las gráfi-
cas de las soluciones y(x)
c1 y y(x)
c2 son rectas horizontales y estas rectas dividen
y
y(x) = c2 R3 la región R en tres subregiones R1, R2 y R3, como se muestra en la figura 2.1.6b. Aquí se
presentan sin prueba alguna de nuestras conclusiones respecto a una solución no cons-
R2 tante y(x) de la ecuación (2):
(x0, y0)
• Si (x0, y0) es una subregión Ri , i
1, 2, 3, y y(x) es una solución cuya gráfica
y(x) = c1
pasa a través de este punto, entonces y(x) permanece en la subregión Ri para
toda x. Como se muestra en la figura 2.1.6b, la solución y(x) en R2 está acotada
R1
x por debajo con c1 y por arriba con c2, es decir, c1 y(x) c2 para toda x. La
curva solución está dentro de R2 para toda x porque la gráfica de una solución no
constante de la ecuación (2) no puede cruzar la gráfica de cualquier solución de
b) subregiones R1, R2, y R3 de R. equilibrio y(x)
c1 o y(x)
c2. Véase el problema 33 de los ejercicios 2.1.
FIGURA 2.1.6 Las rectas y(x)
c1 y • Por continuidad de f debe ser f (y)  0 o f (y) 0 para toda x en una
y(x)
c2 dividen a R en tres subregiones subregión Ri , i
1, 2, 3. En otras palabras, f (y) no puede cambiar de signo
horizontales. en una subregión. Véase el problema 33 de los ejercicios 2.1.

Página 5 de 48

08367_02_ch02_p034-081-ok.indd 38 6/4/09 12:15:48 PM

 • Puesto que dy
dx
f (y(x)) es ya sea positiva o negativa en una subregión Ri ,
i
1, 2, 3, una solución y(x) es estrictamente monótona, es decir, y(x) está
creciendo o decreciendo en la subregión Ri. Por tanto y(x) no puede oscilar, ni
puede tener un extremo relativo (máximo o mínimo). Véase el problema 33
de los ejercicios 2.1.
• Si y(x) está acotada por arriba con un punto crítico c1 (como en la subregión
R1 donde y(x) c1 para toda x), entonces la gráfica de y(x) debe tender a la
gráfica de la solución de equilibrio y(x)
c1 conforme x :  o x : . Si
y(x) está acotada, es decir, acotada por arriba y por debajo por dos puntos
críticos consecutivos (como en la subregión R2 donde c1 y(x) c2 para
toda x), entonces la gráfica de y(x) debe tender a las gráficas de las soluciones
de equilibrio y(x)
c1 y y(x)
c2, conforme x :  en una y x :  en
la otra. Si y(x) está acotada por debajo por un punto crítico (como en la
subregión R3 donde c2 y(x) para toda x), entonces la gráfica de y(x) debe
tender a la gráfica de la solución de equilibrio y(x)
c2 conforme ya sea
x :  o x : . Véase el problema 34 de los ejercicios 2.1.
Considerando estos hechos, analicemos la ecuación diferencial del ejemplo 3.

EJEMPLO 4 Volver a tratar el ejemplo

Los tres intervalos determinados en el eje P o recta de fase con los puntos críticos P

0 y P
a
b ahora corresponden en el plano tP a tres subregiones definidas por:
R1:  P 0, R 2: 0 P a
b, y R 3: a
b P ,
donde  t . El esquema de fase de la figura 2.1.7 nos dice que P(t) está de-
P P
creciendo en R1, creciendo en R2 y decreciendo en R3. Si P(0)
P0 es un valor inicial,
R3
entonces en R1, R2 y R3 tenemos, respectivamente, que:
decreciente P
0 i) Para P0 0, P(t) está acotada por arriba. Puesto que P(t) está decreciendo
a sin límite conforme aumenta t, y así P(t) : 0 conforme t : . Lo que
b
P0 significa que en el eje t negativo, la gráfica de la solución de equilibrio P(t)
creciente R2
0, es una asíntota horizontal para una curva solución.
0 ii) Para 0 P0 a
b, P(t) está acotada. Puesto que P(t) está creciendo,
t P(t) : a
b conforme t :  y P(t) : 0 conforme t : . Las gráficas
decreciente P0 R1 de las dos soluciones de equilibrio, P(t)
0 y P(t)
a
b, son rectas
horizontales que son asíntotas horizontales para cualquier curva solución
recta de fase Plano tP que comienza en esta subregión.
iii) Para P0  a
b, P(t) está acotada por debajo. Puesto que P(t) está
FIGURA 2.1.7 Esquema de fase y
decreciendo, P(t) : a
b conforme t : . La gráfica de la solución de
curvas solución en cada una de las tres
subregiones.
equilibrio P(t)
a
b es una asíntota horizontal para una curva solución.
En la figura 2.1.7 la recta de fase es el eje P en el plano tP. Por claridad la recta de
fase original de la figura 2.1.5 se ha reproducido a la izquierda del plano en el cual se
han sombreado las regiones R1, R2 y R3. En la figura se muestran las gráficas de las
soluciones de equilibrio P(t)
a
b y P(t)
0 (el eje t) como las rectas punteadas
azules; las gráficas sólidas representan las gráficas típicas de P(t) mostrando los tres
casos que acabamos de analizar.

En una subregión tal como R1 en el ejemplo 4, donde P(t) está decreciendo y no

está acotada por debajo, no se debe tener necesariamente que P(t) : . No inter-
prete que este último enunciado significa que P(t) :  conforme t : ; podríamos
tener que P(t) :  conforme t : T, donde T  0 es un número finito que depende
de la condición inicial P(t0)
P0. Considerando términos dinámicos, P(t) “explota”
en un tiempo finito; considerando la gráfica, P(t) podría tener una asíntota vertical en
t
T  0. Para la subregión R3 vale una observación similar.
La ecuación diferencial dy
dx
sen y en el ejemplo 2 es autónoma y tiene un nú-
mero infinito de puntos críticos, ya que sen y
0 en y
np, con n entero. Además, sabe-

Página 6 de 48

08367_02_ch02_p034-081-ok.indd 39 6/4/09 12:15:48 PM

 mos que debido a que la solución y(x) pasa por (0, 32) está acotada por arriba y por debajo
por dos puntos críticos consecutivos (p y(x) 0) y decrece (sen y 0 para p
y 0), la gráfica de y(x) debe tender a las gráficas de las soluciones de equilibrio como
asíntotas horizontales: y(x) : p conforme x :  y y(x) : 0 conforme x : .

EJEMPLO 5 Curvas solución de una ED autónoma

La ecuación autónoma dy
dx
(y  1)2 tiene un solo punto crítico 1. Del esquema
de fase de la figura 2.1.8a concluimos que una solución y(x) es una función creciente
en las subregiones definidas por  y 1 y 1 y , donde  x . Para
una condición inicial y(0)
y0 1, una solución y(x) está creciendo y está acotada
por arriba por 1 y así y(x) : 1 conforme x : ; para y(0)
y0  1, una solución y(x)
está creciendo y está acotada.
Ahora y(x)
1 1
(x  c) es una familia uniparamétrica de soluciones de la
ecuación diferencial (vea el problema 4 de los ejercicios 2.2). Una condición ini-
cial dada determina un valor para c. Para las condiciones iníciales, y(0)
1 1
y y(0)
2  1, encontramos, respectivamente, que y(x)
1 − 1/(x  12), y(x)
1 − 1/(x
− 1). Como se muestra en las figuras 2.1.8b y 2.1.8c, la gráfica de cada una de estas
y y y
x =1
creciente (0, 2)
y=1 y =1
1

x x
(0, −1)
creciente

1
x= −
2

a) recta de fase b) plano xy c) plano xy

y(0) 1 y(0)  1
FIGURA 2.1.8 Comportamiento de las soluciones cerca de y
1.

funciones racionales tienen una asíntota vertical. Pero tenga en mente que las solucio-
nes de los problemas con valores iniciales
dy dy

( y  1) 2, y(0)
1 y
( y  1) 2, y(0)
2 .
dx dx
están definidas en intervalos especiales. Éstos son, respectivamente,
1 1
y(x) 1 1
, 21 x y y(x) 1 , x 1.
x 2 x 1
Las curvas solución son las partes de las gráficas de las figuras 2.1.8b y 2.1.8c que
se muestran en azul. Como lo indica el esquema de fase, para la curva solución de la
y0 y0 figura 2.1.8b, y(x) : 1 conforme x :  para la curva solución de la figura 2.1.8c, y(x)
c c c c :  conforme x : 1 por la izquierda.
y0 y0
ATRACTORES Y REPULSORES Suponga que y(x) es una solución no constante de
la ecuación diferencial autónoma dada en (1) y que c es un punto crítico de la ED.
Básicamente hay tres tipos de comportamiento que y(x) puede presentar cerca de c. En
a) b) c) d) la figura 2.1.9 hemos puesto a c en las cuatro rectas verticales. Cuando ambas puntas
FIGURA 2.1.9 El punto crítico c es un de flecha en cualquier lado del punto c apuntan hacia c, como se muestra en la figura
atractor en a) y un repulsor en b) y semi- 2.1.9a, todas las soluciones y(x) de la ecuación (1) que comienzan en el punto inicial
estable en c) y d). (x0, y0) suficientemente cerca de c presentan comportamiento asintótico límx→ y(x)
c.

Página 7 de 48

08367_02_ch02_p034-081-ok.indd 40 6/4/09 12:15:49 PM

 Por esta razón se dice que el punto crítico c es asintóticamente estable. Utilizando una
analogía física, una solución que comienza en c se parece a una partícula cargada que,
con el tiempo, se transforma en una partícula de carga contraria y así c también se conoce
como un atractor. Cuando ambas puntas de flecha a los lados de la flecha del punto c
apuntan alejándose de c, como se muestra en la figura 2.1.9b, todas las soluciones y(x) de
la ecuación (1) que comienzan en un punto inicial (x0, y0) se alejan de c conforme crece x.
las pendientes de los varían las pendientes En este caso se dice que el punto crítico c es inestable. Un punto crítico inestable se co-
elementos lineales de los elementos sobre
noce como un repulsor, por razones obvias. En las figuras 2.1.9c y 2.1.9d se muestra el
sobre una recta una recta vertical.
punto crítico c que no es ni un atractor ni un repulsor. Pero puesto que c presenta carac-
horizontal son
todas iguales. terísticas tanto de atractor como de repulsor, es decir, una solución que comienza desde
y un punto inicial (x0, y0) que está suficientemente cerca de c es atraída hacia c por un lado
y repelida por el otro, este punto crítico se conoce como semiestable. En el ejemplo 3 el
punto crítico a
b es asintóticamente estable (un atractor) y el punto crítico 0 es inestable
(un repulsor). El punto crítico 1 del ejemplo 5 es semiestable.

x ED AUTÓNOMAS Y CAMPOS DIRECCIONALES Si una ecuación diferencial de

primer orden es autónoma, entonces vemos del miembro derecho de su forma normal
dy
dx
f (y) que las pendientes de los elementos lineales que pasan por los puntos en
la malla rectangular que se usa para construir un campo direccional para la ED que sólo
depende de la coordenada y de los puntos. Expresado de otra manera, los elementos li-
neales que pasan por puntos de cualquier recta horizontal deben tener todos la misma
pendiente; por supuesto, pendientes de elementos lineales a lo largo de cualquier recta
FIGURA 2.1.10 Campo direccional
vertical, variarán. Estos hechos se muestran examinando la banda horizontal amarilla y
para una ED autónoma.
la banda vertical azul de la figura 2.1.10. La figura presenta un campo direccional para la
ecuación autónoma dy
dx
2y – 2. Recordando estos hechos, examine nuevamente
la figura 2.1.4.

EJERCICIOS 2.1
2.1.1 CAMPOS DIRECCIONALES dy

e0.01x y
2
2.
dx
En los problemas 1 a 4 reproduzca el campo direccional dado ge- a) y(6)
0 b) y(0)
1
nerado por computadora. Después dibuje a mano, una curva solu-
c) y(0)
4 d) y(8)
4
ción aproximada que pase por cada uno de los puntos indicados.
Utilice lápices de colores diferentes para cada curva solución.
y
dy
1.
x2  y2 8
dx
a) y(2)
1 b) y(3)
0 4

c) y(0)
2 d) y(0)
0 x
y
_4
3

2 _8

1 _8 _4 4 8
x FIGURA 2.1.12 Campo direccional del problema 2.
_1

_2 dy
3.
1  xy
_3 dx
_3 _2 _1 1 2 3 a) y(0)
0 b) y(1)
0
FIGURA 2.1.11 Campo direccional del problema 1. c) y(2)
2 d) y(0)
4

Página 8 de 48

08367_02_ch02_p034-081-ok.indd 41 6/4/09 12:15:49 PM

 42 ● CAPÍTULO 2 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

y En los problemas 13 y 14 la figura dada representa la gráfica

4
de f (y) y de f (x), respectivamente. Dibuje a mano un campo
direccional sobre una malla adecuada para dy
dx
f (y) (pro-
2 blema 13) y después para dy
dx
f (x) (problema 14).

x 13. f

_2 1

_4 y
_4 _2 2 4 1

FIGURA 2.1.13 Campo direccional del problema 3.

dy FIGURA 2.1.15 Gráfica del problema 13.

4.
(sen x) cos y
dx
14.
a) y(0)
1 b) y(1)
0 f

c) y(3)
3 d) y(0)
52

y

4
1
2
x
1
x

_2
FIGURA 2.1.16 Gráfica del problema 14.
_4
15. En los incisos a) y b) dibuje isoclinas f (x, y)
c (vea los
_4 _2 2 4 omentarios de la página 37) para la ecuación diferencial
FIGURA 2.1.14 Campo direccional del problema 4. dada usando los valores de c indicados. Construya un campo
direccional sobre una malla dibujando con cuidado elemen-
En los problemas 5 a 12 use un paquete computacional para tos lineales con la pendiente adecuada en los puntos elegi-
obtener un campo direccional para la ecuación diferencial dos de cada isoclina. En cada caso, utilice esta dirección
dada. Dibuje a mano una curva solución aproximada que pase para dibujar una curva solución aproximada para el PVI que
por los puntos dados. consiste en la ED y en la condición inicial y (0)
1.
5. y
x 6. y
x  y a) dy
dx
x  y; c un entero que satisface 5  c  5
a) y(0)
0 a) y(2)
2 b) dy
dx
x 2  y 2; c
14, c
1, c
94, c
4
b) y(0)
3 b) y(1)
3
Problemas para analizar
dy dy 1
7. y
x 8.
16. a) Considere el campo direccional de la ecuación dife-
dx dx y
rencial dy
dx
x(y – 4)2 – 2, pero no use tecnología
a) y(1)
1 a) y(0)
1 para obtenerlo. Describa las pendientes de los elemen-
b) y(0)
4 b) y(2)
1 tos lineales en las rectas x
0, y
3, y
4 y y
5.
dy dy b) Considere el PVI dy
dx
x(y – 4)2 – 2, y(0)
y0,
9.
0.2x 2  y 10.
xey
dx dx donde y0 4. Analice, basándose en la información
a) y(0)
12 a) y(0)
2 del inciso a), ¿sí puede una solución y(x) :  con-
forme x : ?
b) y(2)
1 b) y(1)
2.5
 dy y 17. Para la ED de primer orden dy
dx
f (x, y) una curva en
11. y
y  cos x 12.
1 el plano definido por f (x, y)
0 se llama ceroclina de
2 dx x
la ecuación, ya que un elemento lineal en un punto de la
a) y(2)
2 ( )
a) y 12
2 curva tiene pendiente cero. Use un paquete computacional
b) y(1)
0 b) y (32)
0 para obtener un campo direccional en una malla rectangu-

Página 9 de 48

08367_02_ch02_p034-081-ok.indd 42 6/4/09 12:15:50 PM

 lar de puntos dy
dx
x2  2y y después superponga la 29. f
gráfica de la ceroclina y
12 x 2 sobre el campo direccional.
Analice el campo direccional. Analice el comportamiento
de las curvas solución en regiones del plano definidas por
y 12 x 2 y por y  12 x 2. Dibuje algunas curvas solución c y
aproximadas. Trate de generalizar sus observaciones.
18. a) Identifique las ceroclinas (vea el problema 17) en los
problemas 1, 3 y 4. Con un lápiz de color, circule FIGURA 2.1.17 Gráfica del problema 29.
todos los elementos lineales de las figuras 2.1.11,
2.1.13 y 2.1.14, que usted crea que pueden ser un ele-
mento lineal en un punto de la ceroclina. 30. f

b) ¿Qué son las ceroclinas de una ED autónoma de pri-

mer orden?
1

2.1.2 ED DE PRIMER ORDEN AUTÓNOMAS 1 y

19. Considere la ecuación diferencial de primer orden dy
dx

y – y3 y la condición inicial y(0)
y0. A mano, dibuje
la gráfica de una solución típica y(x) cuando y0 tiene los
valores dados.
FIGURA 2.1.18 Gráfica del problema 30.
a) y 0  1 b) 0 y 0 1
c) 1 y 0 0 d) y 0 1
Problemas para analizar
20. Considere la ecuación diferencial autónoma de primer
orden dy
dx
y2 – y4 y la condición inicial y(0)
y0. A 31. Considere la ED autónoma dy
dx
(2
p)y – sen y.
mano, dibuje la gráfica de una solución típica y(x) cuando Determine los puntos críticos de la ecuación. Proponga
y0 tiene los valores dados. un procedimiento para obtener un esquema de fase de la
ecuación. Clasifique los puntos críticos como asintótica-
a) y 0  1 b) 0 y 0 1
mente estable, inestable o semiestable.
c) 1 y 0 0 d) y 0 1
32. Un punto crítico c de una ED de primer orden autónoma
En los problemas 21 a 28 determine los puntos críticos y el es- se dice que está aislada si existe algún intervalo abierto
quema de fase de la ecuación diferencial autónoma de primer que contenga a c pero no otro punto crítico. ¿Puede exis-
orden dada. Clasifique cada punto crítico como asintótica- tir una ED autónoma de la forma dada en la ecuación (1)
mente estable, inestable o semiestable. Dibuje a mano curvas para la cual todo punto crítico no esté aislado? Analice:
solución típicas en las regiones del plano xy determinadas por no considere ideas complicadas.
las gráficas de las soluciones de equilibrio. 33. Suponga que y(x) es una solución no constante de la ecua-
ción diferencial autónoma dy
dx
f (y) y que c es un punto
dy
21.
y 2  3y 22. dy
y 2  y 3 crítico de la ED. Analice. ¿Por qué no puede la gráfica de
dx dx y(x) cruzar la gráfica de la solución de equilibrio y
c?
dy ¿Por qué no puede f (y) cambiar de signo en una de las re-
23.
(y  2)4 24. dy
10  3y  y 2 giones analizadas de la página 38? ¿Por qué no puede y(x)
dx dx
oscilar o tener un extremo relativo (máximo o mínimo)?
dy
25.
y 2(4  y 2) 26. dy
y(2  y)(4  y) 34. Suponga que y(x) es una solución de la ecuación autó-
dx dx
noma dy
dx
f (y) y está acotada por arriba y por debajo
28. dy
ye  9y
dy y
27.
y ln( y  2) por dos puntos críticos consecutivos c1 c2, como una
dx dx ey subregión R2 de la figura 2.1.6b. Si f (y)  0 en la región,
entonces límx: y(x)
c 2. Analice por qué no puede exis-
En los problemas 29 y 30 considere la ecuación diferencial tir un número L c2 tal que límx: y(x)
L. Como parte
autónoma dy
dx
f (y), donde se presenta la gráfica de f.
de su análisis, considere qué pasa con y (x) conforme
Utilice la gráfica para ubicar los puntos críticos de cada una
x : .
de las ecuaciones diferenciales. Dibuje un esquema de fase de
cada ecuación diferencial. Dibuje a mano curvas solución 35. Utilizando la ecuación autónoma (1), analice cómo se
típicas en las subregiones del plano xy determinadas por las puede obtener información respecto a la ubicación de
gráficas de las soluciones de equilibrio. puntos de inflexión de una curva solución.

Página 10 de 48

08367_02_ch02_p034-081-ok.indd 43 6/4/09 12:15:53 PM

 36. Considere la ED dy
dx
y2 – y – 6. Use sus ideas del pro- cuerpo de masa m que está cayendo bajo la influencia de
blema 35 para encontrar los intervalos en el eje y para los la gravedad. Debido a que el término –kv representa la
que las curvas solución son cóncavas hacia arriba y en los que resistencia del aire, la velocidad de un cuerpo que cae de
las curvas solución son cóncavas hacia abajo. Analice por una gran altura no aumenta sin límite conforme pasa el
qué cada curva solución de un problema con valores ini- tiempo t. Utilice un esquema de fase de la ecuación dife-
ciales dy
dx
y2  y – 6, y(0)
y0, donde 2 y0 rencial para encontrar la velocidad límite o terminal del
3, tiene un punto de inflexión con la misma coordenada cuerpo. Explique su razonamiento.
y. ¿Cuál es la coordenada y? Con cuidado dibuje la curva
41. Suponga que el modelo del problema 40 se modifica de tal
solución para la que y(0)
1. Repita para y(2)
2.
manera que la resistencia del aire es proporcional a v2, es
37. Suponga que la ED autónoma en la ecuación (1) no tiene decir
puntos críticos. Analice el comportamiento de las solu- dv
ciones. m mg kv2 .
dt
Vea el problema 17 de los ejercicios 1.3. Utilice un es-
Modelos matemáticos quema de fase para determinar la velocidad terminal del
cuerpo. Explique su razonamiento.
38. Modelo de población La ecuación diferencial en el
ejemplo 3 es un muy conocido modelo de población. 42. Reacciones químicas Cuando se combinan ciertas cla-
Suponga que la ED se cambia por ses de reacciones químicas, la razón con la que se forman
los nuevos componentes se modela por la ecuación dife-
dP rencial autónoma

P(aP  b),
dt dX

k(  X)(   X),
dt
donde a y b son constantes positivas. Analice qué le pasa
a la población P conforme pasa el tiempo. donde k  0 es una constante de proporcionalidad y b 
a  0. Aquí X(t) denota el número de gramos del nuevo
39. Modelo de población Otro modelo de población está componente al tiempo t.
dado por
a) Utilice un esquema de fase de la ecuación diferencial
dP para predecir el comportamiento de X(t) conforme

kP  h,
dt t : .
b) Considere el caso en que a
b. Utilice un esquema
donde h y k son constantes positivas. ¿Para qué valor ini-
de fase de la ecuación diferencial para predecir el
cial P(0)
P0 este modelo predice que la población des-
comportamiento de X(t) conforme t :  cuando X(0)
aparecerá?
a. Cuando X(0)  a.
40. Velocidad terminal En la sección 1.3 vimos que la c) Compruebe que una solución explícita de la ED en
ecuación diferencial autónoma el caso en que k
1 y a
b es X(t)
a  1
(t 
dv c). Determine una solución que satisfaga que X(0)

m mg kv. a
2. Después determine una solución que satisfaga
dt que X(0)
2a. Trace la gráfica de estas dos solucio-
donde k es una constante positiva y g es la aceleración nes. ¿El comportamiento de las soluciones conforme
de la gravedad, es un modelo para la velocidad v de un t :  concuerdan con sus respuestas del inciso b)?

2.2 VARIABLES SEPARABLES

REPASO DE MATERIAL
● Fórmulas básicas de integración (véase al final del libro).

● Técnicas de integración: integración por partes y por descomposición en fracciones parciales.

INTRODUCCIÓN Comenzaremos nuestro estudio de cómo resolver las ecuaciones diferenciales

con las más simple de todas las ecuaciones diferenciales: ecuaciones diferenciales de primer orden con
variables separables. Debido a que el método que se presenta en esta sección y que muchas de las técni-
cas para la solución de ecuaciones diferenciales implican integración, consulte su libro de cálculo para
recordar las fórmulas importantes (como  du
u) y las técnicas (como la integración por partes).

Página 11 de 48

08367_02_ch02_p034-081-ok.indd 44 6/4/09 12:15:54 PM

 SOLUCIÓN POR INTEGRACIÓN Considere la ecuación diferencial de primer
orden dy
dx
f (x, y). Cuando f no depende de la variable y, es decir, f (x, y)
g(x),
la ecuación diferencial
dy

g(x) (1)
dx
se puede resolver por integración. Si g(x) es una función continua, al integrar ambos
lados de la ecuación (1) se obtiene y
g(x) dx = G(x)  c, donde G(x) es una anti-
derivada (integral indefinida) de g(x). Por ejemplo, si dy
dx
1  e2x, entonces su
solución es y (1 e 2x ) dx o y x 12 e2x c.

UNA DEFINICIÓN La ecuación (l) así como su método de solución, no son más
que un caso especial en el que f, en la forma normal dy
dx
f (x, y) se puede factori-
zar como el producto de una función de x por una función de y.

DEFINICIÓN 2.2.1 Ecuación separable

Una ecuación diferencial de primer orden de la forma

dy

g(x)h(y)
dx
Se dice que es separable o que tiene variables separables.

Por ejemplo, las ecuaciones

dy dy

y 2xe3x4y y
y  sen x
dx dx
son respectivamente, separable y no separable. En la primera ecuación podemos fac-
torizar f (x, y)
y 2xe 3x4y como
g(x) h(y)
p p
f (x, y)
y2xe3x4y
(xe3x )( y2e4y ),

pero en la segunda ecuación no hay forma de expresar a y  sen x como un producto

de una función de x por una función de y.
Observe que al dividir entre la función h(y), podemos escribir una ecuación sepa-
rable dy
dx
g(x)h(y) como

dy
p( y)
g(x), (2)
dx

donde, por conveniencia p(y) representa a l
h(y). Podemos ver inmediatamente que la
ecuación (2) se reduce a la ecuación (1) cuando h(y)
1.
Ahora si y
f(x) representa una solución de la ecuación (2), se tiene que
p(f(x))f (x)
g(x), y por tanto

 p( (x))(x) dx
 g(x) dx. (3)

Pero dy
f (x)dx, por lo que la ecuación (3) es la misma que

p( y) dy g(x) dx o H(y) G(x) c, (4)

donde H(y) y G(x) son antiderivadas de p(y)
1
h(y) y g(x), respectivamente.

Página 12 de 48

08367_02_ch02_p034-081-ok.indd 45 6/4/09 12:15:55 PM

 MÉTODO DE SOLUCIÓN La ecuación (4) indica el procedimiento para resolver
ecuaciones separables. Al integrar ambos lados de p(y) dy
g(x) dx, se obtiene una fa-
milia uniparamétrica de soluciones, que usualmente se expresa de manera implícita.

NOTA No hay necesidad de emplear dos constantes cuando se integra una ecuación
separable, porque si escribimos H(y)  c1
G(x)  c2, entonces la diferencia c2 – c1 se
puede reemplazar con una sola constante c, como en la ecuación (4). En muchos casos
de los capítulos siguientes, sustituiremos las constantes en la forma más conveniente
para una ecuación dada. Por ejemplo, a veces se pueden reemplazar los múltiplos o las
combinaciones de constantes con una sola constante.

EJEMPLO 1 Solución de una ED separable

Resuelva (1  x) dy  y dx
0.

SOLUCIÓN Dividiendo entre (1  x)y, podemos escribir dy
y
dx
(1  x), de

donde tenemos que

  dy
y


dx
1x
ln y 
ln 1  x   c1
y
eln1xc1
eln1x ⴢ ec1 ; leyes de exponentes

 1  x  ec1


ec1(1  x).
;  11  xx 

1(1 x, x), x 1
x <1

Haciendo c igual a ⫾ec1 se obtiene y
c(1  x).

SOLUCIÓN ALTERNATIVA Como cada integral da como resultado un logaritmo, la

elección más prudente para la constante de integración es lnc, en lugar de c. Rees-
cribiendo el segundo renglón de la solución como lny
ln1  x  lnc nos permi-
te combinar los términos del lado derecho usando las propiedades de los logaritmos.
De lny
lnc(1  x) obtenemos inmediatamente que y
c(1  x). Aun cuando no
todas las integrales indefinidas sean logaritmos, podría seguir siendo más conveniente
usar lnc. Sin embargo, no se puede establecer una regla firme.

En la sección 1.1 vimos que una curva solución puede ser sólo un tramo o un arco
de la gráfica de una solución implícita G(x, y)
0.

EJEMPLO 2 Curva solución

dy x
Resuelva el problema con valores iniciales
 , y(4)
3.
dx y
SOLUCIÓN Si reescribe la ecuación como y dy
x dx, obtiene

 
y dy
 x dx y
y2
2
x2

  c1.
2
Podemos escribir el resultado de la integración como x 2  y 2
c 2, sustituyendo a la
constante 2c1 por c2. Esta solución de la ecuación diferencial representa una familia de
circunferencias concéntricas centradas en el origen.
Ahora cuando x
4, y
3, se tiene 16  9
25
c2. Así, el problema con valo-
res iniciales determina la circunferencia x 2  y 2
25 de radio 5. Debido a su sencillez
podemos escribir de esta solución implícita como una solución explícita que satisfaga la

Página 13 de 48

08367_02_ch02_p034-081-ok.indd 46 6/4/09 12:15:56 PM

 y condición inicial. Vimos en el ejemplo 3 de la sección 1.1, esta solución como y

f2(x) o y
 125  x2, 5 x 5. Una curva solución es la gráfica de una fun-
ción derivable. En este caso la curva solución es la semicircunferencia inferior, que se
muestra en azul oscuro en la figura 2.2.1 que contiene al punto (4, 3).

x 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,
(4, −3) si r es una raíz de la función h(y), entonces sustituyendo y
r en dy
dx
g(x)h(y)
se encuentra que ambos lados son iguales a cero; es decir, y
r es una solución
constante de la ecuación diferencial. Pero después de que las variables se separan,
dy
FIGURA 2.2.1 Curvas solución para el lado izquierdo de
g (x) dx está indefinido en r. Por tanto, y
r podría no
el PVI del ejemplo 2. h( y)
representar a 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 3 Pérdida de una solución

dy
Resuelva
y 2  4.
dx
SOLUCIÓN Poniendo la ecuación en la forma

y  2  dy
dx.
1 1
dy

dx o 4
 4
(5)
y2  4 y2

La segunda ecuación en la ecuación (5) 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
ln y 2 ln y 2 x c1
4 4

y 2 y 2
o ln 4x c2 o e4x c2
.
y 2 y 2

Aquí hemos sustituido 4c1 por c2. Por último, después de sustituir
ec2 por c y despe-
jando y de la última ecuación, obtenemos una familia uniparamétrica de soluciones

1  ce4x
y
2 . (6)
1  ce4x

Ahora, si factorizamos el lado derecho de la ecuación diferencial como dy
dx

(y  2)(y  2), sabemos del análisis de puntos críticos de la sección 2.1 que y
2 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 (6) correspondiendo al valor
c
0. Sin embargo, y
2 es una solución singular; ésta no se puede obtener de la
ecuación (6) 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 (5) indica claramente que
debemos excluir y

2 en estos pasos.

EJEMPLO 4 Un problema con valores iniciales

dy
Resuelva (e2y  y) cos x
ey sen 2x, y(0)
0.
dx

Página 14 de 48

08367_02_ch02_p034-081-ok.indd 47 6/4/09 12:15:57 PM

 SOLUCIÓN Dividiendo la ecuación entre ey cos x se obtiene

e2y  y sen 2x
dy
dx.
ey cos x
Antes de integrar se realiza la división del lado izquierdo y utilizamos la identidad
trigonométrica sen 2x
2 sen x cos x en el lado derecho. Entonces tenemos que

integración de partes  (ey ye y) dy 2 sen x dx

se obtiene e y  yey  ey
2 cos x  c. (7)

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  yey  ey
4  2 cos x. (8)

USO DE COMPUTADORA Los Comentarios al final de la sección 1.1 mencionan

que puede ser difícil utilizar una solución implícita G(x, y)
0 para encontrar una solu-
ción explícita y
f(x). La ecuación (8) muestra que la tarea de despejar y en términos
de x puede presentar más problemas que solamente el aburrido trabajo de presionar
símbolos; ¡en algunos casos simplemente no se puede hacer! Las soluciones implícitas
tales como la ecuación (8) son un poco frustrantes; ya que no se aprecia ni en la gráfica
y
de la ecuación ni en el intervalo una solución definida que satisfaga que y(0)
0. El
2
problema de “percibir” cuál es la solución implícita en algunos casos se puede resol-
ver mediante la tecnología. Una manera* de proceder es utilizar la aplicación contour
1 plot de un sistema algebraico de computación (SAC). Recuerde del cálculo de varias
variables que para una función de dos variables z
G(x, y) las curvas bi-dimensionales
x definidas por G(x, y)
c, donde c es una constante, se llaman las curvas de nivel de la
función. En la figura 2.2.2 se presentan algunas de las curvas de nivel de la función G(x,
_1 y)
ey  yey  ey  2 cos x que se han reproducido con la ayuda de un SAC. La fa-
milia de soluciones definidas por la ecuación (7) son las curvas de nivel G(x, y)
c. En
la figura 2.2.3 se muestra en color azul la curva de nivel G(x, y)
4, que es la solución
_2
_2 _1 1 2
particular de la ecuación (8). La otra curva de la figura 2.2.3 es la curva de nivel G(x, y)

2, que es miembro de la familia G(x, y)
c que satisface que y(p
2)
0.
FIGURA 2.2.2 Curvas de nivel Si al determinar un valor específico del parámetro c en una familia de soluciones
G(x, y)
c, donde de una ecuación diferencial de primer orden llegamos a una solución particular, hay una
G(x, y)
ey  yey  ey  2 cos x.
inclinación natural de la mayoría de los estudiantes (y de los profesores) a relajarse y estar
satisfechos. Sin embargo, una solución de un problema con valores iniciales podría no ser
y única. Vimos en el ejemplo 4 de la sección 1.2 que el problema con valores iniciales
2
dy

x y1/2, y(0)
0 (9)
c=4 dx
1
1 4
tiene al menos dos soluciones, y
0 y y
16 x . Ahora ya podemos resolver esa ecua-
x ción. Separando las variables e integrando y1
2 dy
x dx obtenemos
(0, 0) (π /2,0)
x2 x2 2
_1 c=2 2y1/2 c1 o y c .
2 4
_2 1 4
_2 _1 1 2 Cuando x
0, entonces y
0, así que necesariamente, c
0. Por tanto y
16 x . Se
perdió la solución trivial y
0 al dividir entre y1
2. Además, el problema con valores
FIGURA 2.2.3 Curvas de nivel iniciales, ecuación (9), tiene una cantidad infinitamente mayor de soluciones porque
c
2 y c
4.
para cualquier elección del parámetro a  0 la función definida en tramos

*
En la sección 2.6 analizaremos algunas otras maneras de proceder que están basadas en el concepto de una
solución numérica.

Página 15 de 48

08367_02_ch02_p034-081-ok.indd 48 6/4/09 12:15:58 PM

 y
0, (x  a ) , x a
a=0 a>0
y
1
16
2 2 2
xa

satisface tanto a la ecuación diferencial como a la condición inicial. Véase la fi-

(0, 0) x gura 2.2.4.

SOLUCIONES DEFINIDAS POR INTEGRALES Si g es una función continua en

FIGURA 2.2.4 Soluciones de la un intervalo abierto I que contiene a a, entonces para toda x en I,
ecuación (9) definida en tramos.
d
dx
 x

a
g(t) dt
g(x).

Usted podría recordar que el resultado anterior es una de las dos formas del teorema
fundamental del cálculo. Es decir, ax g(t) dt es una antiderivada de la función g. En
ocasiones esta forma es conveniente en la solución de ED. Por ejemplo, si g es continua
en un intervalo I que contiene a x0 y a x, entonces una solución del sencillo problema
con valores iniciales dy
dx
g(x), y(x0)
y0, que está definido en I está dado por

y(x)
y0   x

x0
g(t) dt

Usted debería comprobar que y(x) definida de esta forma satisface la condición inicial.
Puesto que una antiderivada de una función continua g no siempre puede expresarse
en términos de las funciones elementales, esto podría ser lo mejor que podemos hacer
para obtener una solución explícita de un PVI. El ejemplo siguiente ilustra esta idea.

EJEMPLO 5 Un problema con valores iniciales

dy

ex , y(3)
5.
2
Resuelva
dx

SOLUCIÓN La función g(x)
e−x2 es continua en (, ), pero su antiderivada

no es una función elemental. Utilizando a t como una variable muda de integración,
podemos escribir

3
x
dy
dt
dt
x

3
et dt
2

y(t) ]x
3
x

3
et dt
2

y(x)  y(3)
x

3
et dt
2

y(x)
y(3)  
3
x
et dt.
2

Utilizando la condición inicial y(3)
5, obtenemos la solución

y(x)
5  3
x
et dt.
2

El procedimiento que se mostró en el ejemplo 5 también funciona bien en las ecua-

ciones separables dy
dx
g(x) f (y) donde, f (y) tiene una antiderivada elemental pero g(x)
no tiene una antiderivada elemental. Véanse los problemas 29 y 30 de los ejercicios 2.2.

Página 16 de 48

08367_02_ch02_p034-081-ok.indd 49 6/4/09 12:15:59 PM

 COMENTARIOS

i) Como acabamos de ver en el ejemplo 5, algunas funciones simples no tienen

una antiderivada que es una función elemental. Las integrales de estas clases de
funciones se llaman no elementales. Por ejemplo 3x e−t dt y sen x2 dx son integra-
2

les no elementales. Retomaremos nuevamente este concepto en la sección 2.3.

ii) En algunos de los ejemplos anteriores vimos que la constante de la familia
uniparamétrica de soluciones de una ecuación diferencial de primer orden se
puede redefinir cuando sea conveniente. También se puede presentar con faci-
lidad el caso de que dos personas obtengan distintas expresiones de las mismas
respuestas resolviendo correctamente la misma ecuación. Por ejemplo, sepa-
rando variables se puede demostrar que familias uniparamétricas de soluciones
de la ED (l  y2) dx  (1 x2) dy
0 son
xy
arctan x  arctan y
c o
c.
1  xy
Conforme avance en las siguientes secciones, considere que las familias de so-
luciones pueden ser equivalentes, en el sentido de que una se puede obtener de
otra, ya sea por redefinición de la constante o utilizando álgebra o trigonometría.
Vea los problemas 27 y 29 de los ejercicios 2.2.

EJERCICIOS 2.2
En los problemas 1 a 22 resuelva la ecuación diferencial dada dy dy
21.
x 1  y2 22. (ex  ex )
y2
por separación de variables. dx dx
En los problemas 23 a 28 encuentre una solución explícita del
dy dy problema con valores iniciales dado.
1. sen 5x 2. (x 1)2
dx dx dx
23.
4(x2  1), x(> 4)
1
3. dx  e 3xdy
0 4. dy  (y  1) 2dx
0 dt
dy dy dy y2  1
5. x 4y 6. 2xy 2 0 24.
, y(2)
2
dx dx dx x2  1
dy dy dy
7. e3x 2y
8. e x y e y
e 2x y
25. x2
y  xy, y(1)
1
dx dx dx
2 2
dx y 1 dy 2y 3 dy
9. y ln x 10. 26.  2y
1, y(0)
52
dy x dx 4x 5 dt

11. csc y dx  sec 2x dy
0 13

27. 11  y2 dx  11  x2 dy
0, y(0)

2
12. sen 3x dx  2y cos 33x dy
0 28. (1  x 4) dy  x(1  4y 2) dx
0, y(1)
0
13. (e y  1) 2ey dx  (e x  1) 3ex dy
0 En los problemas 29 y 30 proceda como en el ejemplo 5 y de-
termine una solución explícita del problema con valores ini-
14. x(1  y 2) 1
2 dx
y(1  x 2) 1
2 dy ciales dado.
dy

yex , y(4)
1
2

dQ 29.
dS dx
15. kS 16. k(Q 70)
dr dt
dy
30.
y 2 sen x 2, y(2)
13
dP dN dx
17. P P2 18. N Ntet 2
dt dt 31. a) Encuentre una solución al problema con valores inicia-
les que consiste en la ecuación diferencial del ejemplo
19.
dy xy 3x y 3
20.
dy xy 2y x 2 3 y de las condiciones iniciales y(0)
2, y(0)
2,
dx xy 2x 4y 8 dx xy 3y x 3 ()
y y 14
1.

Página 17 de 48

08367_02_ch02_p034-081-ok.indd 50 6/4/09 12:16:00 PM

 b) Encuentre la solución de la ecuación diferencial en el y1(0)
4, y2(0)
2, y3(1)
2 y y4(1)
4. Trace
ejemplo 4 cuando se utiliza In c1 como la constante de la gráfica de cada solución y compare con sus dibu-
integración del lado izquierdo en la solución y 4 In c1 jos del inciso a). Indique el intervalo de definición
se sustituye por In c. Después resuelva los mismos pro- exacto de cada solución.
blemas con valores iniciales que en el inicio a).
41. a) Determine una solución explícita del problema con
dy valores iniciales
32. Encuentre una solución de x
y2  y que pase por
dx dy 2x 1
los puntos indicados. , y( 2) 1.
a) (0, 1) b) (0, 0) ( )
c) 12, 12 ( )
d) 2, 14 dx 2y
b) Utilice un programa de graficación para dibujar la
33. Encuentre una solución singular del problema 21 y del gráfica de la solución del inciso a). Use la gráfica para
problema 22. estimar el intervalo I de definición de la solución.
34. Demuestre que una solución implícita de c) Determine el intervalo I de definición exacto me-
2x sen 2 y dx  (x2  10) cos y dy
0 diante métodos analíticos.

está dada por ln(x2  10)  csc y
c. Determine las so- 42. Repita los incisos a) al c) del problema 41 para el PVI que
luciones constantes si se perdieron cuando se resolvió la consiste en la ecuación diferencial del problema 7 y de la
ecuación diferencial. condición inicial y(0)
0.

Con frecuencia, un cambio radical en la forma de la solución Problemas para analizar

de una ecuación diferencial corresponde a un cambio muy 43. a) Explique por qué el intervalo de definición de la solu-
pequeño en la condición inicial o en la ecuación misma. En ción explícita y
f2(x) del problema con valores ini-
los problemas 35 a 38 determine una solución explícita del ciales en el ejemplo 2 es el intervalo abierto (5, 5).
problema con valores iniciales dado. Utilice un programa de b) ¿Alguna solución de la ecuación diferencial puede
graficación para dibujar la gráfica de cada solución. Compare cruzar el eje x? ¿Usted cree que x2  y2
1 es una
cada curva solución en una vecindad de (0, 1). solución implícita del problema con valores iniciales
dy dy
dx
x
y, y(1)
0?
35.
(y  1)2, y(0)
1
dx 44. a) Si a  0 analice las diferencias, si existen, entre las
dy soluciones de los problemas con valores iniciales
36.
(y  1)2, y(0)
1.01 que consisten en la ecuación diferencial dy
dx
x
y
dx
y de cada una de las condiciones iniciales y(a)
a,
dy
37.
(y  1)2  0.01, y(0)
1 y(a)
a, y(a)
a y y(a)
a.
dx
b) ¿Tiene una solución el problema con valores iniciales
dy dy
dx
x
y, y(0)
0?
38.
(y  1)2  0.01, y(0)
1
dx c) Resuelva dy
dx
x
y, y(1)
2 e indique el inter-
39. Toda ecuación autónoma de primer orden dy
dx
f (y) es valo de definición exacto de esta solución.
separable. Encuentre las soluciones explícitas y1(x), y2(x), 45. En los problemas 39 y 40 vimos que toda ecuación di-
y3(x) y y4(x) de la ecuación diferencial dy
dx
y – y3, que ferencial autónoma de primer orden dy
dx
f(y) es se-
satisfagan, respectivamente las condiciones iniciales y1(0)

parable. ¿Ayuda este hecho en la solución del problema
2, y2(0)
12 , y3(0)
 12 y y4(0)
2. Utilice un programa
dy
de graficación para cada solución. Compare estas gráficas con valores iniciales
11  y2 sen2 y, y(0)
12?
con las bosquejadas en el problema 19 de los ejercicios 2.1. dx
Analice. A mano, dibuje una posible curva solución del
Dé el intervalo de definición exacto para cada solución. problema.
40. a) La ecuación diferencial autónoma de primer orden
46. Sin usar tecnología. ¿Cómo podría resolver
dy
dx
1
(y 3) no tiene puntos críticos. No obs-
tante, coloque 3 en la recta de fase y obtenga un es- dy
quema de fase de la ecuación. Calcule d 2 y
dx2 para ( 1x  x) dx
1y  y?
determinar dónde las curvas solución son cóncavas Lleve a cabo sus ideas.
hacia arriba y dónde son cóncavas hacia abajo (vea
los problemas 35 y 36 de los ejercicios 2.1). Utilice 47. Determine una función cuyo cuadrado más el cuadrado
el esquema de fase y la concavidad para que, a mano, de su derivada es igual a 1.
dibuje algunas curvas solución típicas. 48. a) La ecuación diferencial del problema 27 es equiva-
b) Encuentre las soluciones explícitas y1(x), y2(x), y3(x) lente a la forma normal
y y4(x) de la ecuación diferencial del inciso a) que dy 1  y2
satisfagan, respectivamente las condiciones iniciales

dx B1  x 2

Página 18 de 48

08367_02_ch02_p034-081-ok.indd 51 6/4/09 12:16:05 PM

 en la región cuadrada del plano xy definida por x de la familia de soluciones de la ecuación diferencial
1, y 1. Pero la cantidad dentro del radical es no ne- dy 8x 5
. Experimente con diferentes números
gativa también en las regiones definidas por x  1, dx 3y 2 1
y  1. Dibuje todas las regiones del plano xy para las de las curvas de nivel así como con diferentes regiones
que esta ecuación diferencial tiene soluciones reales. rectangulares definidas por a  x  b, c  y  d.
b) Resuelva la ED del inciso a) en las regiones definidas b) En diferentes ejes coordenados dibuje las gráficas
por x  1, y  1. Después determine una solución de las soluciones particulares correspondientes a las
implícita y una explícita de la ecuación diferencial su- condiciones iniciales: y(0)
1; y(0)
2; y(1)

jeta a y(2)
2. 4; y(1)
3.

51. a) Determine una solución implícita del PVI

Modelo matemático
(2y  2) dy  (4x3  6x) dx
0, y(0)
3.
49. Puente suspendido En la ecuación (16) de la sección
1.3 vimos que un modelo matemático para la forma de un b) Utilice el inciso a) para encontrar una solución explí-
cable flexible colgado de dos postes es cita y
f(x) del PVI.
dy W c) Considere su respuesta del inciso b) como una sola

, (10) función. Use un programa de graficación o un SAC
dx T1
para trazar la gráfica de esta función y después utilice
donde W denota la porción de la carga vertical total entre la gráfica para estimar su dominio.
los puntos P1 y P2 que se muestran en la figura 1.3.7. La
ED, ecuación (10) es separable bajo las siguientes condi- d) Con la ayuda de un programa para determinar raíces
ciones que describen un puente suspendido. de un SAC, determine la longitud aproximada del in-
Supongamos que los ejes x y y están como se mues- tervalo de definición más grande posible de la solu-
tra en la figura 2.2.5, es decir, el eje x va a lo largo de la ción y
f(x) del inciso b). Utilice un programa de
superficie de la carretera y el eje y pasa por (0, a), que graficación o un SAC para trazar la gráfica de la curva
es el punto más bajo de un cable en la región que abarca solución para el PVI en este intervalo.
el puente, que coincide con el intervalo [L
2, L
2]. En el
caso de un puente suspendido, la suposición usual es que la 52. a) Utilice un SAC y el concepto de curvas de nivel para
carga vertical en (10) es sólo una distribución uniforme de dibujar las gráficas representativas de los miembros
la superficie de la carretera a lo largo del eje horizontal. En de la familia de soluciones de la ecuación diferencial
otras palabras, se supone que el peso de todos los cables es dy x(1  x)

. Experimente con diferentes núme-
despreciable en comparación con el peso de la superficie de dx y(2  y)
la carretera y que el peso por unidad de longitud de la super- ros de curvas de nivel así como en diferentes regiones
ficie de la carretera (digamos, libras por pie horizontal) es rectangulares del plano xy hasta que su resultado se
una constante r. Utilice esta información para establecer y parezca a la figura 2.2.6.
resolver un adecuado problema con valores iniciales a par- b) En diferentes ejes coordenados, dibuje la gráfica de
tir del cual se determine la forma (una curva con ecuación la solución implícita correspondiente a la condición
y
f(x)) de cada uno de los dos cables en un puente sus- inicial y(0)
23. Utilice un lápiz de color para indicar
pendido. Exprese su solución del PVI en términos del pan- el segmento de la gráfica que corresponde a la curva
deo h y de la longitud L. Véase la figura 2.2.5. solución de una solución f que satisface la condición
inicial. Con ayuda de un programa para determinar raí-
y ces de un SAC, determine el intervalo I de definición
cable aproximado más largo de la solución f. [Sugerencia:
h (pandeo) Primero encuentre los puntos en la curva del inciso a)
(0, a) donde la recta tangente es vertical.]
c) Repita el inciso b) para la condición inicial y(0)
2.
L/2 L/2 x y
L longitud
superficie de la carretera (carga)

FIGURA 2.2.5 Forma de un cable del problema 49.

x

Tarea para el laboratorio de computación

50. a) Utilice un SAC y el concepto de curvas de nivel para
dibujar las gráficas representativas de los miembros FIGURA 2.2.6 Curvas de nivel del problema 52.

Página 19 de 48

08367_02_ch02_p034-081-ok.indd 52 6/4/09 12:16:07 PM

 2.3 ECUACIONES LINEALES
REPASO DE MATERIAL
● Repase la definición de las ED en las ecuaciones (6 ) y (7) de la sección 1.1

INTRODUCCIÓN Continuamos con nuestra búsqueda de las soluciones de las ED de primer orden
examinando ecuaciones lineales. Las ecuaciones diferenciales lineales son una familia especialmente
“amigable” de ecuaciones diferenciales en las que, dada una ecuación lineal, ya sea de primer orden
o de un miembro de orden superior, siempre hay una buena posibilidad de que podamos encontrar
alguna clase de solución de la ecuación que podamos examinar.

UNA DEFINICIÓN En la ecuación (7) de la sección 1.1, se presenta la forma de

una ED lineal de primer orden. Aquí, por conveniencia, se reproduce esta forma en la
ecuación (6) de la sección 1.1, para el caso cuando n
1.

DEFINICIÓN 2.3.1 Ecuación lineal

Una ecuación diferencial de primer orden de la forma

dy
a1(x)  a0(x)y
g(x) (1)
dx
se dice que es una ecuación lineal en la variable dependiente y.

Se dice que la ecuación lineal (1) es homogénea cuando g(x)
0; si no es no

homogénea.

FORMA ESTÁNDAR Al dividir ambos lados de la ecuación (1) entre el primer coefi-
ciente, a1(x), se obtiene una forma más útil, la forma estándar de una ecuación lineal:

dy
 P(x)y
f (x). (2)
dx
Buscamos una solución de la ecuación (2) en un intervalo I, en el cual las dos funcio-
nes P y f sean continuas.
En el análisis que se presenta a continuación ilustraremos una propiedad y un proce-
dimiento y terminaremos con una fórmula que representa la forma de cada solución de la
ecuación (2). Pero más importantes que la fórmula son la propiedad y el procedimiento,
porque ambos conceptos también se aplican a ecuaciones lineales de orden superior.

LA PROPIEDAD La ecuación diferencial (2) tiene la propiedad de que su solución

es la suma de las dos soluciones, y
yc  yp, donde yc es una solución de la ecuación
homogénea asociada
dy
 P(x)y
0 (3)
dx
y yp es una solución particular de ecuación no homogénea (2). Para ver esto, observe que

d
dx
dy
dx [ ] [
dy
]
––– [yc  yp]  P(x)[ yc  yp]
–––c  P(x)yc  –––p  P(x)yp
f (x).
dx

0 f (x)

Página 20 de 48

08367_02_ch02_p034-081-ok.indd 53 6/4/09 12:16:07 PM

 Ahora la ecuación (3) es también separable. Por lo que podemos determinar yc al es-
cribir la ecuación (3) en la forma
dy
 P(x) dx
0
y
e integramos. Despejando y, se obtiene yc
ceP(x)dx. Por conveniencia escribimos
yc
cy1(x), donde y1
eP(x)dx. A continuación se utiliza el hecho de que dy1
dx 
P(x)y1
0, para determinar yp.

EL PROCEDIMIENTO Ahora podemos definir una solución particular de la ecua-

ción (2), siguiendo un procedimiento llamado variación de parámetros. Aquí, la idea
básica es encontrar una función, u tal que yp
u(x)y1(x)
u(x)e−P(x)dx sea una solución
de la ecuación (2). En otras palabras, nuestra suposición para yp es la misma que yc

cy1(x) excepto que c se ha sustituido por el “parámetro variable” u. Sustituyendo yp

uy1 en la ecuación (2) se obtiene
Regla del producto cero

dy
u –––1
dx
du
y1–––
dx
P(x)uy1 f(x) o [dy
u –––1
dx
P(x)y1 ] du
y1 –––
dx
f (x)

du
por tanto y1 f (x).
dx
Entonces separando las variables e integrando se obtiene

du

f (x)
y1(x)
dx y u
 f (x)
y1(x)
dx.

Puesto que y1(x)
eP(x)dx, vemos que 1
y1(x)
eP(x)dx. Por tanto

yp
uy1
 f (x)
y1(x)
dx eP(x)d x
eP(x)d x  eP(x)d x f (x) dx,

y y
ceP (x)dx  eP (x)dx eP (x) dxf(x) dx.  (4)

yc yp
Por tanto, si la ecuación (2) tiene una solución, debe ser de la forma de la ecuación (4).
Recíprocamente, es un ejercicio de derivación directa comprobar que la ecuación (4)
es una familia uniparamétrica de soluciones de la ecuación (2).
No memorice la fórmula que se presenta en la ecuación (4). Sin embargo recuerde
el término especial
e ∫P(x)dx (5)
ya que se utiliza para resolver la ecuación (2) de una manera equivalente pero más
fácil. Si la ecuación (4) se multiplica por (5),

eP(x)d xy
c   e P(x)d x f (x) dx, (6)

y después se deriva la ecuación (6),

d P(x)d x
dx
e [ ]
y
e P(x)d x f (x), (7)

dy
se obtiene eP(x)dx  P(x)eP(x)dx y
eP(x)dx f (x). (8)
dx

Dividiendo el último resultado entre e P(x)dx se obtiene la ecuación (2).

Página 21 de 48

08367_02_ch02_p034-081-ok.indd 54 6/4/09 12:16:08 PM

 MÉTODO DE SOLUCIÓN El método que se recomienda para resolver la ecuación
(2) consiste en realidad en trabajar con las ecuaciones (6) a (8) en orden inverso. En otras
palabras, si la ecuación (2) se multiplica por la ecuación (5), obtenemos la ecuación (8). Se
reconoce que el lado izquierdo de la ecuación (8) es la derivada del producto de e P(x)dx por
y. Esto nos conduce a la ecuación (7). Entonces, integrando ambos lados de la ecuación
(7) se obtiene la solución (6). Como podemos resolver la ecuación (2) por integración,
después de multiplicar por e P(x)dx, esta función se llama factor integrante de la ecuación
diferencial. Por conveniencia resumiremos estos resultados. Nuevamente le indicamos
que no debe memorizar la fórmula (4) sino seguir cada vez el siguiente procedimiento.

SOLUCIÓN DE UNA ECUACIÓN LINEAL DE PRIMER ORDEN

i) Ponga la ecuación lineal de la forma (1) en la forma estándar (2).

ii) Identifique de la identidad de la forma estándar P(x) y después
determine el factor integrante e P(x)dx.
iii) Multiplique la forma estándar de la ecuación por el factor integrante. El
lado izquierdo de la ecuación resultante es automáticamente la derivada
del factor integrante y y:
d P(x)dx
dx
e [
y
e P(x)dx f (x).]
iv) Integre ambos lados de esta última ecuación.

EJEMPLO 1 Solución de una ED lineal homogénea

dy
Resuelva  3y
0.
dx

SOLUCIÓN Esta ecuación lineal se puede resolver por separación de variables. En

otro caso, puesto que la ecuación ya está en la forma estándar (2), vemos que P(x)

3 y por tanto el factor integrante es e (3)dx
e3x. Multiplicando la ecuación por este
factor y reconociendo que

3x
dy 3x
d 3x
e 3e y 0 es la misma que
[e y] 0.
dx dx
Integrando ambos lados de la última ecuación se obtiene e3xy
c. Despejando y se
obtiene la solución explícita y
ce 3x,  x .

EJEMPLO 2 Solución de una ED lineal no homogénea

dy
Resuelva  3y
6.
dx

SOLUCIÓN La ecuación homogénea asociada a esta ED se resolvió en el ejemplo 1.

Nuevamente la ecuación está ya en la forma estándar (2) y el factor integrante aún es
e (3)dx
e3x. Ahora al multiplicar la ecuación dada por este factor se obtiene

3x
dy 3x 3x
d 3x 3x
e 3e y 6e , que es la misma que [e y] 6e .
dx dx
Integrando ambos lados de la última ecuación se obtiene e3xy
2e3x  c o
y
2  ce 3x,  x .

Página 22 de 48

08367_02_ch02_p034-081-ok.indd 55 6/4/09 12:16:09 PM

 La solución final del ejemplo 2 es la suma de dos soluciones: y
yc  yp, donde yc

ce3x es la solución de la ecuación homogénea del ejemplo 1 y yp
2 es una solu-
ción particular de la ecuación no homogénea y  – 3y
6. No necesita preocuparse de
si una ecuación lineal de primer orden es homogénea o no homogénea; cuando sigue
y el procedimiento de solución que se acaba de describir, la solución de una ecuación
no homogénea necesariamente produce y
yc  yp. Sin embargo, la diferencia entre
1 resolver una ED homogénea y una no homogénea será más importante en el capítulo 4,
x donde se resolverán ecuaciones lineales de orden superior.
_1 Cuando a1, a0 y g en la ecuación (1) son constantes, la ecuación diferencial es
_2
y =_2 autónoma. En el ejemplo 2 podemos comprobar de la forma normal dy
dx
3(y  2)
que 2 es un punto crítico y que es inestable (un repulsor). Así, una curva solución
_3
con un punto inicial ya sea arriba o debajo de la gráfica de la solución de equilibrio
_1 1 2 3 4 y
2 se aleja de esta recta horizontal conforme x aumenta. La figura 2.3.1, obtenida
con la ayuda de una aplicación para trazo de gráficas, muestra la gráfica de y
2
FIGURA 2.3.1 Algunas soluciones junto con otras curvas solución.
y  3y
6.

CONSTANTE DE INTEGRACIÓN Observe que en el análisis general y en los

ejemplos 1 y 2 no se ha considerado una constante de integración en la evaluación de
la integral indefinida en el exponente e P(x)dx. Si consideramos las leyes de los expo-
nentes y el hecho de que el factor integrante multiplica ambos lados de la ecuación
diferencial, usted podría explicar por qué es innecesario escribir P(x)dx  c. Vea el
problema 44 de los ejercicios 2.3.

SOLUCIÓN GENERAL Suponga que las funciones P y f en la ecuación (2) son

continuas en un intervalo I. En los pasos que conducen a la ecuación (4) mostramos
que si la ecuación (2) tiene una solución en I, entonces debe estar en la forma dada en
la ecuación (4). Recíprocamente, es un ejercicio directo de derivación comprobar que
cualquier función de la forma dada en (4) es una solución de la ecuación diferencial (2)
en I. En otras palabras (4) es una familia uniparamétrica de soluciones de la ecuación
(2) y toda solución de la ecuación (2) definida en I es un miembro de esta familia. Por
tanto llamamos a la ecuación (4) la solución general de la ecuación diferencial en
el intervalo I. (Véase los Comentarios al final de la sección 1.1.) Ahora escribiendo la
ecuación (2) en la forma normal y 
F(x, y), podemos identificar F(x, y)
P(x)y
 f (x) y F
y
P(x). De la continuidad de P y f en el intervalo I vemos que F y
F
y son también continuas en I. Con el teorema 1.2.1 como nuestra justificación,
concluimos que existe una y sólo una solución del problema con valores iniciales
dy
 P(x)y
f(x), y(x0)
y0 (9)
dx
definida en algún intervalo I0 que contiene a x0. Pero cuando x0 está en I, encontrar una
solución de (9) es exactamente lo mismo que encontrar un valor adecuado de c en la
ecuación (4), es decir, a toda x0 en I le corresponde un distinto c. En otras palabras,
el intervalo de existencia y unicidad I0 del teorema 1.2.1 para el problema con valores
iniciales (9) es el intervalo completo I.

EJEMPLO 3 Solución general

dy
Resuelva x  4y
x 6e x.
dx

SOLUCIÓN Dividiendo entre x, obtenemos la forma estándar

dy 4
 y
x5e x. (10)
dx x

Página 23 de 48

08367_02_ch02_p034-081-ok.indd 56 6/4/09 12:16:10 PM

 En esta forma identificamos a P(x)
4
x y f (x)
x5ex y además vemos que P y f son
continuas en (0, ). Por tanto el factor integrante es
podemos utilizar ln x en lugar de ln
x
ya que x  0

e4dx/x
e4ln x
eln x4
x4.

Aquí hemos utilizado la identidad básica blogbN
N, N  0. Ahora multiplicamos la

ecuación (10) por x4 y reescribimos

4
dy d
x 4x 5y xex como
[x 4y] xex.
dx dx

De la integración por partes se tiene que la solución general definida en el intervalo (0,
) es x4y
xe x  e x  c o y
x 5e x  x 4e x  cx 4.

Excepto en el caso en el que el coeficiente principal es 1, la reformulación de la

ecuación (1) en la forma estándar (2) requiere que se divida entre a1(x). Los valores
de x para los que a1(x)
0 se llaman puntos singulares de la ecuación. Los puntos
singulares son potencialmente problemáticos. En concreto, en la ecuación (2), si P(x)
(que se forma al dividir a0(x) entre a1(x)) es discontinua en un punto, la discontinuidad
puede conducir a soluciones de la ecuación diferencial.

EJEMPLO 4 Solución general

dy
Determine la solución general de (x 2  9)  xy
0.
dx

SOLUCIÓN Escribimos la ecuación diferencial en la forma estándar

dy x
 y
0 (11)
dx x 2  9

e identificando P(x)
x
(x2 – 9). Aunque P es continua en (, 3), (3, 3) y (3,
), resolveremos la ecuación en el primer y tercer intervalos. En estos intervalos el
factor integrante es
ex d x/(x 9)
e2 2x d x/(x 9)
e2 lnx 9
1x2  9 .
2 1 2 1 2

Después multiplicando la forma estándar (11) por este factor, obtenemos

d
dx 
1x2  9 y
0. 
Integrando ambos lados de la última ecuación se obtiene 1x  9 y
c. Por
2

tanto para cualquiera x  3 o x 3 la solución general de la ecuación es

c
y
.
1x 2  9
Observe en el ejemplo 4 que x
3 y x
3 son puntos singulares de la ecuación
y que toda función en la solución general y
c
1x 2  9 es discontinua en estos pun-
tos. Por otra parte, x
0 es un punto singular de la ecuación diferencial en el ejemplo
3, pero en la solución general y
x5ex – x4ex  cx4 es notable que cada función de esta
familia uniparamétrica es continua en x
0 y está definida en el intervalo (, ) y no
sólo en (0, ), como se indica en la solución. Sin embargo, la familia y
x5ex – x4ex  cx4
definida en (, ) no se puede considerar la solución general de la ED, ya que el punto
singular x
0 aún causa un problema. Véase el problema 39 en los ejercicios 2.3.

Página 24 de 48

08367_02_ch02_p034-081-ok.indd 57 6/4/09 12:16:11 PM

 EJEMPLO 5 Un problema con valores iniciales

dy
Resuelva  y
x, y(0)
4.
dx

SOLUCIÓN La ecuación está en forma estándar, y P(x)
1 y f(x)
x son continuas

en (, ). El factor integrante es e dx
e x, entonces integrando
d x
[e y]
xex
dx
se tiene que exy
xex – ex  c. Despejando y de esta última ecuación se obtiene la
solución general y
x  1  ce x. Pero de la condición general sabemos que y
4
cuando x
0. El sustituir estos valores en la solución general implica que c
5. Por
tanto la solución del problema es
y
x  1  5ex,  x . (12)
La figura 2.3.2, que se obtuvo con la ayuda de un programa de graficación, mues-
y tra la gráfica de (12) en azul oscuro, junto con las gráficas, de las otras soluciones re-
presentativas de la familia uniparamétrica y
x – 1 cex. En esta solución general
4
identificamos yc
cex y yp
x – 1. Es interesante observar que conforme x aumenta,
c>0 las gráficas de todos los miembros de la familia son cercanas a la gráfica de la solución
2
particular yp
x – 1 que se muestra con una línea sólida de la figura 2.3.2. Esto es de-
x bido a que la contribución de yc
cex a los valores de una solución es despreciable al
aumentar los valores de x. Decimos que yc
cex es un término transitorio, ya que yc
_2
c<0
: 0 conforme x : . Mientras que este comportamiento no es característico de todas
las soluciones generales de las ecuaciones lineales (véase el ejemplo 2), el concepto de
_4 c=0 un transitorio es frecuentemente importante en problemas aplicados.
_4 _2 2 4

FIGURA 2.3.2 Algunas soluciones COEFICIENTES DISCONTINUOS En aplicaciones, los coeficientes P(x) y f(x)
y  y
x. en (2) pueden ser continuos por tramos. En el siguiente ejemplo f(x) es continua por
tramos en [0, ) con una sola discontinuidad, en particular un salto (finito) discontinuo
en x
1. Resolvemos el problema en dos partes correspondientes a los dos interva-
los en los que f está definida. Es entonces posible juntar las partes de las dos soluciones
en x
1 así que y(x) es continua en [0, ).

EJEMPLO 6 Un problema con valores iniciales

Resuelvaa
dy
dx
y
f (x), y(0)
0 donde f (x)
1,
0,
0
x
x 1,
1.

y SOLUCIÓN En la figura 2.3.3 se muestra la gráfica de la función discontinua f.

Resolvemos la ED para y(x) primero en el intervalo [0, 1] y después en el intervalo
(1, ). Para 0  x  1 se tiene que

x
dy
dx
y 1 o, el equivalente, d x
dx
[e y] ex.

FIGURA 2.3.3 f(x) discontinua. Integrando esta última ecuación y despejando y se obtiene y
1  c1ex. Puesto que
y(0)
0, debemos tener que c1
1 y por tanto y
1  ex, 0  x  1. Entonces
para x  1 la ecuación
dy
y
0
dx

Página 25 de 48

08367_02_ch02_p034-081-ok.indd 58 6/4/09 12:16:12 PM

 conduce a y
c2ex. Por tanto podemos escribir

1c e e,
x
, 0  x  1,
y y
x
2 x  1.
Invocando a la definición de continuidad en un punto, es posible determinar c2 así la
última función es continua en x
1. El requisito de límx→1 y(x)
y(1) implica que
c2e1
1 – e1 o c2
e1. Como se muestra en la figura 2.3.4, la función
x
 1  ex, 0  x  1,
1
y
(13)
FIGURA 2.3.4 Gráfica de la función (e  1)ex, x1
de (13). es continua en (0, ).
Es importante considerar la ecuación (13) y la figura 2.3.4 como un bloque pe-
queño; le pedimos que lea y conteste el problema 42 de los ejercicios 2.3.

FUNCIONES DEFINIDAS POR INTEGRALES Al final de la sección 2.2 analiza-

mos 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 y sen x 2 dx
2

no son integrales elementales. En matemáticas aplicadas algunas funciones importan-

tes están definidas en términos de las integrales no elementales. Dos de esas funciones
especiales son la función error y la función error complementario:

erf(x)

2
1
0
x
et dt
2
y erfc(x)

2
1
 x

et dt.
2
(14)

Del conocido resultado 0 et dt
1
2* podemos escribir (2
1 ) 0 et dt
1.
2 2

Entonces de la forma 0
0x  x se ve de la ecuación (14) que la función error
complementario, erfc(x), se relaciona con erf(x) por erf(x)  erfc(x)
1. Debido a su
importancia en probabilidad, estadística y en ecuaciones diferenciales parciales apli-
cadas se cuenta con extensas tablas de la función error. Observe que erf(0)
0 es un
valor obvio de la función. Los valores de erf(x) se pueden determinar con un sistema
algebraico de computación (SAC).

EJEMPLO 7 La función error

dy
Resuelva el problema con valores iniciales  2xy
2, y(0)
1.
dx

SOLUCIÓN Puesto que la ecuación ya se encuentra en la forma normal, el factor

integrante es e−x dx, y así de
2

x
y d x2 x2 2
t2 2
[e y] 2e obtenemos
y 2ex e dt cex . (15)
dx 0

Aplicando y(0)
1 en la última expresión obtenemos c
1. Por tanto, la solución del

problema es


x
x
y
2ex et dt  ex o y
ex [1  1 erf(x)].
2 2 2 2

En la figura 2.3.5 se muestra en azul oscuro, la gráfica de esta solución en el intervalo

FIGURA 2.3.5 Algunas soluciones (, ) junto con otros miembros de la familia definida en la ecuación (15), obtenida
de y  2xy
2. con la ayuda de un sistema algebraico de computación.

*
Este resultado normalmente se presenta en el tercer semestre de cálculo.

Página 26 de 48

08367_02_ch02_p034-081-ok.indd 59 6/4/09 12:16:13 PM

 USO DE COMPUTADORAS Algunos sistemas algebraicos de computación como
Mathematica y Maple permiten obtener soluciones implícitas o explícitas para algunos
tipos de ecuaciones diferenciales, usando la instrucción dsolve.*

COMENTARIOS

i) En general, una ED lineal de cualquier orden se dice que es homogénea cuando

g(x)
0 en la ecuación (6) de la sección 1.1. Por ejemplo, la ED lineal de se-
gundo orden y – 2y   6y
0 es homogénea. Como se puede ver en este ejem-
plo y en el caso especial de la ecuación (3) de esta sección, la solución trivial
y
0 es siempre una solución de una ED lineal homogénea.
ii) A veces, una ecuación diferencial de primer orden es no lineal en una variable
pero es lineal en la otra variable. Por ejemplo, la ecuación diferencial
dy 1


dx x  y 2
es no lineal en la variable y. Pero su recíproca
dx dx

x  y2 o  x
y2
dy dy
se reconoce como lineal en la variable x. Usted debería comprobar que el factor
integrante es e (1)dy
ey e integrando por partes se obtiene la solución ex-
plícita x
y2  2y  2  cey para la segunda ecuación. Esta expresión es,
entonces, una solución implícita de la primera ecuación.
iii) Los matemáticos han adoptado como propias algunas palabras de ingeniería
que consideran adecuadas para describir. La palabra transitorio, que ya hemos
usado, es uno de estos términos. En futuros análisis ocasionalmente se presenta-
rán las palabras entrada y salida. La función f en la ecuación (2) es la función
de entrada o de conducción; una solución y(x) de la ecuación diferencial para
una entrada dada se llama salida o respuesta.
iv) El término funciones especiales mencionado en relación con la función error
también se aplica a la función seno integral y a la integral seno de Fresnel
introducidas en los problemas 49 y 50 de los ejercicios 2.3. “Funciones especia-
les” es una rama de las matemáticas realmente bien definidas. En la sección 6.3
se estudian funciones más especiales.

*Ciertas instrucciones se deletrean igual, pero las instrucciones en Mathematica inician con una letra
mayúscula (Dsolve) mientras que en Maple la misma instrucción comienza con una letra minúscula
(dsolve). Cuando analizamos la sintaxis de las instrucciones, nos comprometimos y escribimos, por
ejemplo dsolve.

EJERCICIOS 2.3
En los problemas 1 a 24 determine la solución general de la 5. y  3x 2y
x 2 6. y  2xy
x 3
ecuación diferencial dada. Indique el intervalo I más largo en
7. x 2y  xy
1 8. y
2y  x 2  5
el que está definida la solución general. Determine si hay al-
gunos términos transitorios en la solución general. dy dy
9. x y x 2 senx 10. x 2y 3
dx dx
dy dy
1. 5y 2. 2y 0 dy dy
dx dx x3 x2
11. x 4y x 12. (1 x) xy x
dx dx
dy dy
3. y e3x 4. 3 12y 4 13. x 2y  x(x  2)y
e x
dx dx

Página 27 de 48

08367_02_ch02_p034-081-ok.indd 60 6/4/09 12:16:15 PM

 14. xy  (1  x)y
ex sen 2x
dy
15. y dx  4(x  y 6) dy
0 33.
dx
2xy f (x), y(0) 2, donde
16. y dx
(ye y  2x) dy
x, 0 x 1
17. cos x
dy
 (sen x)y
1
f (x)
0, x 1
dx
dy dy
34. (1 x 2) 2xy f (x), y(0) 0, donde
18. cos2x sen x  (cos3x)y
1 dx
dx
dy x, 0 x 1
19. (x  1)
dx
 (x  2)y
2xex f (x)
x, x 1
dy 35. Proceda en una forma similar al ejemplo 6 para resolver el pro-
20. (x  2)2
5  8y  4xy
dx blema con valores iniciales y  P(x)y
4x, y(0)
3, donde

2>x, 0  x  1,
dr 2,
21.  r sec 
cos  P(x)

d x  1.
dP Utilice un programa de graficación para para trazar la grá-
22.  2tP
P  4t  2
dt fica de la función continua y(x).
dy
23. x  (3x  1)y
e3x 36. Considere el problema con valores iniciales y   exy

dx
f (x), y(0)
1. Exprese la solución del PVI para x  0
24. (x 2  1)
dy
 2y
(x  1)2 como una integral no elemental cuando f (x)
1. ¿Cuál
dx es la solución cuando f (x)
0? ¿Y cuándo f (x)
ex?
En los problemas 25 a 30 resuelva el problema con valores ini- 37. Exprese la solución del problema con valores iniciales
ciales. Indique el intervalo I más largo en el que está definida y  – 2xy
1, y(1)
1, en términos de erf(x).
la solución.
25. xy  y
ex, y(1)
2 Problemas para analizar
38. Lea nuevamente el análisis siguiente al ejemplo 2. Cons-
dx
26. y  x
2y2, y(1)
5 truya una ecuación diferencial lineal de primer orden
dy para la que todas las soluciones no constantes tienden a la
di asíntota horizontal y
4 conforme x : .
27. L  Ri
E, i(0)
i0,
dt 39. Lea nuevamente el ejemplo 3 y después analice, usando
L, R, E e i 0 constantes
el teorema 1.2.1, la existencia y unicidad de una solución
dT del problema con valores iniciales que consiste en xy 
28.
k(T  Tm ); T(0)
T0,
dt – 4y
x6ex y de la condición inicial dada.
k, T m y T 0 constantes a) y(0)
0 b) y(0)
y 0, y 0  0
dy c) y(x 0)
y 0, x 0  0, y 0  0
29. (x  1)  y
ln x, y(1)
10
dx
40. Lea nuevamente el ejemplo 4 y después determine la solu-
30. y  (tan x)y
cos 2x, y(0)
1 ción general de la ecuación diferencial en el intervalo (3, 3).
En los problemas 31 a 34 proceda como en el ejemplo 6 para 41. Lea nuevamente el análisis siguiente al ejemplo 5.
resolver el problema con valores iniciales dado. Utilice un pro- Construya una ecuación diferencial lineal de primer orden
grama de graficación para trazar la función continua y(x). para la que todas las soluciones son asintóticas a la recta
y
3x  5 conforme x : .
dy
31. 2y f (x), y(0) 0, donde
dx 42. Lea nuevamente el ejemplo 6 y después analice por qué
es técnicamente incorrecto decir que la función en (13) es
1, 0 x 3 una “solución” del PVI en el intervalo [0, ).
f (x)
0, x 3
43. a) Construya una ecuación diferencial lineal de primer
dy orden de la forma xy   a0(x)y
g(x) para la cual yc
32. y f (x), y(0) 1, donde
c
x3 y yp
x3. Dé un intervalo en el que y
x3 
dx
c
x3 es la solución general de la ED.
1, 0 x 1 b) Dé una condición inicial y(x0)
y0 para la ED que
f (x)

1, x 1 se determinó en el inciso a) de modo que la solución

Página 28 de 48

08367_02_ch02_p034-081-ok.indd 61 6/4/09 12:16:16 PM

 del PVI sea y
x3  1
x3. Repita si la solución es y
el corazón se está estimulando, el voltaje E a través del
x3  2
x3. Dé un intervalo de definición I de cada una corazón satisface la ecuación diferencial lineal
de estas soluciones. Trace la gráfica de las curvas so- dE 1
lución. ¿Hay un problema con valores iniciales cuya
 E.
dt RC
solución esté definida en (, )?
c) ¿Es único cada PVI encontrado en el inciso b)? Es decir, Resuelva la ED sujeta a E(4)
E0.
puede haber más de un solo PVI para el cual, digamos,
y
x3  1
x3, x en algún intervalo I, es la solución?
Tarea para el laboratorio de computación
44. Al determinar el factor integrante (5), no usamos una
constante de integración en la evaluación de P(x) dx. 48. a) Exprese la solución del problema con valores inicia-
Explique por qué usar P(x) dx  c no tiene efecto en la les y  2xy
1, y(0) 1 2, en términos de
solución de (2). erfc(x).
45. Suponga que P(x) es continua en algún intervalo I y a es un b) Utilice las tablas de un SAC para determinar el valor
número en I. ¿Qué se puede decir acerca de la solución del de y(2). Use un SAC para trazar la gráfica de la curva
problema con valores iniciales y   P(x)y
0, y(a)
0? solución para el PVI en (, ).
49. a) La función seno integral está definida por
Modelos matemáticos Si(x) x
, donde el integrando está defi-
0 (sent>t) dt
46. Series de decaimiento radiactivo El siguiente siste- nido igual a 1 en t
0. Exprese la solución y(x) del
ma de ecuaciones diferenciales se encuentra en el estudio problema con valores iniciales x3y   2x2y
10 sen
del decaimiento de un tipo especial de series de elemen- x, y(1)
0 en términos de Si(x).
tos radiactivos: b) Use un SAC para trazar la gráfica de la curva solu-
ción para el PVI para x  0.
dx

1x c) Use un SAC para encontrar el valor del máximo ab-
dt soluto de la solución y(x) para x  0.
dy

1x  2 y, 50. a) La integral seno de Fresnel está definida por
dt x
0 sen(pt >2) dt.
2
S(x) . Exprese la solución y(x) del
donde l1 y l2 son constantes. Analice cómo resolver este sis- problema con valores iniciales y  – (sen x2)y
0,
tema sujeto a x(0)
x0, y(0)
y0. Lleve a cabo sus ideas. y(0)
5, en términos de S(x).
b) Use un SAC para trazar la gráfica de la curva solu-
47. Marcapasos de corazón Un marcapasos de corazón
ción para el PVI en (, ).
consiste en un interruptor, una batería de voltaje cons-
tante E0, un capacitor con capacitancia constante C y c) Se sabe que S(x) : 12 conforme x :  y S(x) :  12
un corazón como un resistor con resistencia constante conforme x :  . ¿A dónde tiende la solución y(x)
R. Cuando se cierra el interruptor, el capacitor se carga; cuando x : ? ¿Y cuando x : ?
cuando el interruptor se abre, el capacitor se descarga en- d) Use un SAC para encontrar los valores del máximo
viando estímulos eléctricos al corazón. Todo el tiempo absoluto y del mínimo absoluto de la solución y(x).

2.4 ECUACIONES EXACTAS

REPASO DE MATERIAL
● Cálculo de varias variables.

● Derivación parcial e integración parcial.

● Diferencial de una función de dos variables.

INTRODUCCIÓN Aunque la sencilla ecuación diferencial de primer orden

y dx  x dy
0
es separable, podemos resolver la ecuación en una forma alterna al reconocer que la expresión del
lado izquierdo de la ecuación es la diferencial de la función f (x, y)
xy, es decir
d(xy)
y dx  x dy.
En esta sección analizamos ecuaciones de primer orden en la forma diferencial M(x, y) dx  N(x, y) dy

0. Aplicando una prueba simple a M y a N, podemos determinar si M(x, y) dx  N(x, y) dy es una
diferencial de una función f (x, y). Si la respuesta es sí, construimos f integrando parcialmente.

Página 29 de 48

08367_02_ch02_p034-081-ok.indd 62 6/4/09 12:16:18 PM

 DIFERENCIAL DE UNA FUNCIÓN DE DOS VARIABLES Si z
f (x, y) es una
función de dos variables con primeras derivadas parciales continuas en una región R
del plano xy, entonces su diferencial es
f f
dz
dx  dy. (1)
x y
En el caso especial cuando f (x, y)
c, donde c es una constante, entonces la ecuación
(1) implica que
f f
dx  dy
0. (2)
x y
En otras palabras, dada una familia de curvas f (x, y)
c, podemos generar una ecua-
ción diferencial de primer orden si calculamos la diferencial de ambos lados de la
igualdad. Por ejemplo, si x2  5xy  y3
c, entonces la ecuación (2) da la ED de
primer orden
(2x  5y) dx  (5x  3y 2 ) dy
0. (3)

UNA DEFINICIÓN Por supuesto, que no todas las ED de primer orden escritas en
la forma M(x, y) dx  N(x, y) dy
0 corresponden a una diferencial de f (x, y)
c. Por
tanto para nuestros objetivos es muy importante regresar al problema anterior; en par-
ticular, si nos dan una ED de primer orden tal como la ecuación (3), ¿hay alguna forma
de reconocer que la expresión diferencial (2x  5y) dx  (5x  3y 2) dy es la diferen-
cial d(x 2  5xy  y 3)? Si la hay, entonces una solución implícita de la ecuación (3) es
x 2  5xy  y 3
c. Podemos contestar esta pregunta después de la siguiente definición.

DEFINICIÓN 2.4.1 Ecuación exacta

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 ésta corresponde a la diferencial de alguna función
f (x, y) definida en R. Una ecuación diferencial de primer orden de la forma
M(x, y) dx  N(x, y) dy
0
se dice que es una ecuación exacta si la expresión del lado izquierdo es una
diferencial exacta.

Por ejemplo x 2y 3 dx  x 3y 2 dy
0 es una ecuación exacta, ya que su lado iz-

quierdo es una diferencial exacta:
d 13 x3 y3
x2 y3 dx  x3y2 dy.

Observe que si hacemos las identificaciones M(x, y)
x 2y 3 y N(x, y)
x 3y 2, entonces

M
y
3x 2y 2
N
x. El teorema 2.4.1, que se presenta a continuación, muestra
que la igualdad de las derivadas parciales M
y y N
x no es una coincidencia.

TEOREMA 2.4.1 Criterio para una diferencial exacta

Sean M(x, y) y N(x, y) continuas y que tienen primeras derivadas parciales con-
tinuas en una región rectangular R definida por a x b, c y d. Entonces
una condición necesaria y suficiente para que M(x, y) dx  N(x, y) dy sea una
diferencial exacta es
M N

. (4)
y x

Página 30 de 48

08367_02_ch02_p034-081-ok.indd 63 6/4/09 12:16:19 PM

 PRUEBA DE LA NECESIDAD Por simplicidad suponemos que M(x, y) y N(x, y) tie-
nen primeras derivadas parciales continuas para todo (x, y). Ahora si la expresión
M(x, y) dx  N(x, y) dy es exacta, existe alguna función f tal que para toda x en R,

f f
M(x, y) dx  N(x, y) dy
dx  dy.
x y

f f
Por tanto M(x, y)
, N(x, y)
,
x y

y
M
y


 f
y x


2 f


 f
y x x y
N

.
x 
La igualdad de las parciales mixtas es una consecuencia de la continuidad de las pri-
meras derivadas parciales de M(x, y) y N(x, y).
La parte de suficiencia del teorema 2.4.1 consiste en mostrar que existe una fun-
ción f para la que f
x
M(x, y) y f
y
N(x, y) siempre que la ecuación (4) sea
válida. La construcción de la función f en realidad muestra un procedimiento básico
para resolver ecuaciones exactas.

MÉTODO DE SOLUCIÓN Dada una ecuación en la forma diferencial M(x, y) dx 

N(x, y) dy
0, determine si la igualdad de la ecuación (4) es válida. Si es así, entonces
existe una función f para la que
f

M(x, y).
x
Podemos determinar f integrando M(x, y) respecto a x mientras y se conserva cons-
tante:

f (x, y)
 M(x, y) dx  g(y), (5)

donde la función arbitraria g(y) es la “constante” de integración. Ahora derivando

(5) respecto a y y suponiendo que f
y
N(x, y):

f


y y

 M(x, y) dx  g(y)
N(x, y).

Se obtiene g( y)
N(x, y) 

y
 M(x, y) dx. (6)

Por último, se integra la ecuación (6) respecto a y y se sustituye el resultado en la

ecuación (5). La solución implícita de la ecuación es f (x, y)
c.
Haremos algunas observaciones en orden. Primero, es importante darse cuenta de
que la expresión N(x, y)  (
y)  M(x, y) dx en (6) es independiente de x, ya que


x 
N(x, y) 

y
 M(x, y) dx 
Nx  y x  M(x, y) dx
Nx  My
0.
Segunda, pudimos iniciar bien el procedimiento anterior con la suposición de que f
y

N(x, y). Después, integrando N respecto a y y derivando este resultado, encontraría-
mos las ecuaciones que, respectivamente, son análogas a las ecuaciones (5) y (6),

f (x, y)
 N(x, y) dy  h(x) y h(x)
M(x, y) 


x
 N(x, y) dy.

En ninguno de ambos casos se deben memorizar estas fórmulas.

Página 31 de 48

08367_02_ch02_p034-081-ok.indd 64 6/4/09 12:16:20 PM

 EJEMPLO 1 Resolviendo una ED exacta

Resuelva 2xy dx  (x 2  1) dy
0.

SOLUCIÓN Con M(x, y)
2xy y N(x, y)
x 2  1 tenemos que

M N

2x
.
y x

Así la ecuación es exacta y por el teorema 2.4.1 existe una función f (x, y) tal que

f f

2xy y
x2  1.
x y

Al integrar la primera de estas ecuaciones, se obtiene:

f (x, y)
x 2y  g (y).
Tomando la derivada parcial de la última expresión con respecto a y y haciendo el
resultado igual a N(x, y) se obtiene
f

x2  g(y)
x2  1. ; N(x, y)
y

Se tiene que g(y)
1 y g(y)
y. Por tanto f (x, y)
x 2y  y, así la solución de
la ecuación en la forma implícita es x 2y  y
c. La forma explícita de la solución se
ve fácilmente como y
c
(1  x 2) y está definida en cualquier intervalo que no con-
tenga ni a x
1 ni a x
1.

NOTA La solución de la ED en el ejemplo 1 no es f (x, y)
x 2y  y. Sino que es

f (x, y)
c; si se usa una constante en la integración de g (y), podemos escribir la
solución como f (x, y)
0. Observe que la ecuación también se podría haber resuelto
por separación de variables.

EJEMPLO 2 Solución de una ED exacta

Resuelva (e 2y  y cos xy) dx  (2xe 2y  x cos xy  2y) dy
0.

SOLUCIÓN La ecuación es exacta ya que

M N

2e 2y  xy sen xy  cos xy
.
y x

Por tanto existe una función f (x, y) para la cual

f f
M(x, y)
y N(x, y)
.
x y
Ahora, para variar, comenzaremos con la suposición de que f
y
N(x, y); es decir

f

2xe2y  x cos xy  2y
y

f (x, y)
2x  e2y dy  x  cos xy dy  2  y dy.

Página 32 de 48

08367_02_ch02_p034-081-ok.indd 65 6/4/09 12:16:21 PM

 Recuerde que la razón por la que x sale del símbolo  es que en la integración respecto
a y se considera que x es una constante ordinaria. Entonces se tiene que

f (x, y)
xe 2y  sen xy  y 2  h(x)

f

e2y  y cos xy  h(x)
e 2y  y cos xy, ; M(x, y)
x
y así h (x)
0 o h(x)
c. Por tanto una familia de soluciones es
xe 2y  sen xy  y 2  c
0.

EJEMPLO 3 Problema con valores iniciales

dy xy2  cos x sen x

Resuelva
, y(0)
2.
dx y(1  x2)

SOLUCIÓN Al escribir la ecuación diferencial en la forma

(cos x sen x  xy 2) dx  y(1  x 2) dy
0,
reconocemos que la ecuación es exacta porque

M N

2xy
.
y x
f
Ahora
y(1  x2)
y
y2
f(x, y)
(1  x 2 )  h(x)
2
f

xy2  h(x)
cos x sen x  x y 2.
x

La última ecuación implica que h (x)
cos x sen x. Integrando se obtiene

1
h(x) (cos x)( sen x dx) cos 2 x.
2

y2 1
Por tanto (1 x2) cos2 x c1 o y2 (1 x2) cos2 x c, (7)
2 2
y
donde se sustituye 2c1 por c. La condición inicial y
2 cuando x
0 exige que
4(1)  cos 2 (0)
c, y por tanto c
3. Una solución implícita del problema es enton-
ces y 2(1  x 2)  cos 2 x
3.
x En la figura 2.4.1, la curva solución del PVI es la curva dibujada en azul oscuro, y
forma parte de una interesante familia de curvas. Las gráficas de los miembros de la fa-
milia uniparamétrica de soluciones dadas en la ecuación (7) se puede obtener de diferen-
tes maneras, dos de las cuales son utilizando un paquete de computación para trazar grá-
ficas de curvas de nivel (como se analizó en la sección 2.2) y usando un programa de
FIGURA 2.4.1 Algunas gráficas graficación para dibujar cuidadosamente la gráfica de las funciones explícitas obtenidas
de los miembros de la familia para diferentes valores de c despejando a y de y 2
(c  cos 2 x)
(1  x 2) para y.
y 2(1  x 2)  cos 2x
c.
FACTORES INTEGRANTES Recuerde de la sección 2.3 que el lado izquierdo de la
ecuación lineal y   P(x)y
f (x) se puede transformar en una derivada cuando mul-
tiplicamos la ecuación por el factor integrante. Esta misma idea básica algunas veces
funciona bien para una ecuación diferencial no exacta M(x, y) dx  N(x, y) dy
0.

Página 33 de 48

08367_02_ch02_p034-081-ok.indd 66 6/4/09 12:16:22 PM

 Es decir, algunas veces es posible encontrar un factor integrante m(x, y) de manera
que, después de multiplicar el lado izquierdo de
m(x, y)M(x, y) dx  m(x, y)N(x, y) dy
0 (8)
es una diferencial exacta. En un intento por encontrar m, regresamos al criterio (4) de
la exactitud. La ecuación (8) es exacta si y sólo si (mM)y
(mN)x, donde los subíndi-
ces denotan derivadas parciales. Por la regla del producto de la derivación la última
ecuación es la misma que mMy  m y M
mNx  m x N o
m x N  m y M
(My  Nx)m. (9)
Aunque M, N, My y Nx son funciones conocidas de x y y, la dificultad aquí al determinar
la incógnita m(x, y) de la ecuación (9) es que debemos resolver una ecuación diferencial
parcial. Como no estamos preparados para hacerlo, haremos una hipótesis para simpli-
ficar. Suponga que m es una función de una variable; por ejemplo, m depende sólo de x.
En este caso, m x
dm
dx y m y
0, así la ecuación (9) se puede escribir como
d My  Nx

. (10)
dx N
Estamos aún en un callejón sin salida si el cociente (My  Nx )
N depende tanto de x
como de y. Sin embargo, si después de que se hacen todas las simplificaciones alge-
braicas el cociente (My  Nx )
N resulta que depende sólo de la variable x, entonces la
ecuación (10) es separable así como lineal. Entonces de la sección 2.2 o de la sección
2.3 tenemos que m(x)
e ((MyNx)/N)dx. Análogamente, de la ecuación (9) tenemos que
si m depende sólo de la variable y, entonces
d Nx  My

. (11)
dy M
En este caso, si (N x  My)
M es una función sólo de y, podemos despejar m de la
ecuación (11).
Resumiendo estos resultados para la ecuación diferencial.
M(x, y) dx  N(x, y) dy
0. (12)
• Si (My  Nx)
N es una función sólo de x, entonces un factor integrante para
la ecuación (12) es

(x)
e

MyNx
N
dx
. (13)
• Si (Nx  My)
M es una función sólo de y, entonces un factor integrante de (12) es

(y)
e
 NxMy
M
dy
. (14)

EJEMPLO 4 Una ED no exacta hecha exacta

La ecuación diferencial no lineal de primer orden

xy dx  (2x 2  3y 2  20) dy
0
es no exacta. Identificando M
xy, N
2x 2  3y 2  20, encontramos que las deriva-
das parciales My
x y Nx
4x. El primer cociente de la ecuación (13) no nos conduce
a nada, ya que
My  Nx x  4x 3x

2
2
N 2x  3y  20 2x  3y 2  20
2

depende de x y de y. Sin embargo, la ecuación (14) produce un cociente que depende

sólo de y:
Nx  My 4x  x 3x 3



.
M xy xy y

Página 34 de 48

08367_02_ch02_p034-081-ok.indd 67 6/4/09 12:16:23 PM

 El factor integrante es entonces e 3dy
y
e 3lny
e lny
y 3. Después de multiplicar la
3

ED dada por m(y)
y3, la ecuación resultante es

xy 4 dx  (2x 2y 3  3y 5  20y 3) dy
0.
Usted debería comprobar que la última ecuación es ahora exacta así como mostrar,
usando el método que se presentó en esta sección, que una familia de soluciones es
1 2 4
2x y  12 y 6  5y 4
c.

COMENTARIOS

i) Cuando pruebe la exactitud de una ecuación, se debe asegurar que tiene exac-
tamente la forma M(x, y) dx  N(x, y) dy
0. Algunas veces una ecuación dife-
rencial se escribe como G(x, y) dx
H(x, y) dy . En este caso, primero rescriba
como G(x, y) dx  H(x, y) dy
0 y después identifique M(x, y)
G(x, y) y N(x,
y)
H(x, y) antes de utilizar la ecuación (4).
ii) En algunos libros de ecuaciones diferenciales el estudio de las ecuaciones
exactas precede al de las ED lineales. Entonces el método que acabamos de des-
cribir para encontrar los factores integrantes se puede utilizar para deducir un
factor integrante para y   P(x) y
f (x). Reescribiendo la última ecuación en la
forma diferencial (P(x)y  f (x)) dx  dy
0, vemos que
M y  Nx

P(x).
N

A partir de la ecuación (13) hemos obtenido el conocido factor integrante e P(x)dx,

utilizado en la sección 2.3.

EJERCICIOS 2.4
En los problemas 1 a 20 determine si la ecuación diferencial 12. (3x 2y  e y ) dx  (x 3  xe y  2y) dy
0
dada es exacta. Si lo es, resuélvala.
dy
1. (2x  1) dx  (3y  7) dy
0 13. x
2xe x  y  6x 2
dx
2. (2x  y) dx  (x  6y) dy
0
3. (5x  4y) dx  (4x  8y 3) dy
0
14. 1  3y  x dydx  y
3x  1
4. (sen y  y sen x) dx  (cos x  x cos y  y) dy
0
5. (2xy 2  3) dx  (2x 2y  4) dy
0
15. x y
2 3

1

dx
1  9x 2 dy
 x 3y 2
0

6. 2y  1x  cos 3x dxdy  xy  4x  3y sen 3x
0

2
3 16. (5y  2x)y  2y
0

17. (tan x  sen x sen y) dx  cos x cos y dy
0

7. (x 2  y 2) dx  (x 2  2xy) dy
0
18. (2y sen x cos x  y  2y 2e xy ) dx
2

8. 1  ln x 
y
x
dx
(1  ln x) dy

(x  sen2 x  4xye xy ) dy
2

9. (x  y 3  y 2 sen x) dx
(3xy 2  2y cos x) dy

10. (x 3  y 3) dx  3xy 2 dy
0 19. (4t 3y  15t 2  y) dt  (t 4  3y 2  t) dy
0

11. (y ln y  e xy) dx  1y  x ln y dy
0 20. 1t  t1  t

2 2
y
y 2 
dt  ye y  2
t
t  y2
dy
0

Página 35 de 48

08367_02_ch02_p034-081-ok.indd 68 6/4/09 12:16:24 PM

 En los problemas 21 a 26 resuelva el problema con valores b) Demuestre que las condiciones iniciales y(0)
2 y
iniciales. y(1)
1 determinan la misma solución implícita.
21. (x  y)2 dx  (2xy  x 2  1) dy
0, y(1)
1 c) Encuentre las soluciones explícitas y1(x) y y2(x) de la
ecuación diferencial del inciso a) tal que y1(0)
2
22. (e  y) dx  (2  x  ye ) dy
0,
x y
y(0)
1 y y2(1)
1. Utilice un programa de graficación para
23. (4y  2t  5) dt  (6y  4t  1) dy
0, y(1)
2 trazar la gráfica de y1(x) y y2(x).

24.  3y 2  t 2 dy
y 5 t
 4
0, y(1)
1
dt 2y
25. ( y 2 cos x  3x 2y  2x) dx Problemas para analizar
 (2y sen x  x 3  ln y) dy
0, y(0)
e 40. Considere el concepto de factor integrante utilizado en
los problemas 29 a 38. ¿Son las dos ecuaciones Mdx  N
26. 1 1 y  cos x  2xy dxdy
y(y  sen x), y(0)
1
2
dy
0 y mM dx  mN dy
0 necesariamente equivalen-
tes en el sentido de que la solución de una es también una
En los problemas 27 y 28 determine el valor de k para el que solución de la otra? Analice.
la ecuación diferencial es exacta.
41. Lea nuevamente el ejemplo 3 y después analice por qué
27. ( y 3  kxy 4  2x) dx  (3xy 2  20x 2y 3) dy
0 podemos concluir que el intervalo de definición de la so-
28. (6xy 3  cos y) dx  (2kx 2y 2  x sen y) dy
0 lución explícita del PVI (curva azul de la figura 2.4.1) es
(1, 1).
En los problemas 29 y 30 compruebe que la ecuación dife-
rencial dada es no exacta. Multiplique la ecuación diferencial 42. Analice cómo se pueden encontrar las funciones M(x, y) y
dada por el factor integrante indicado m(x, y) y compruebe que N(x, y) tal que cada ecuación diferencial sea exacta. Lleve
la nueva ecuación es exacta. Resuelva. a cabo sus ideas.

29. (xy sen x  2y cos x) dx  2x cos x dy
0;

m(x, y)
xy

a) M(x, y) dx  xe x y  2xy 
1
x dy
0

30. (x 2  2xy  y 2) dx  (y 2  2xy  x 2) dy
0;

m(x, y)
(x  y)2
b) x 1/2 1/2
y 
x
x y
2
dx  N(x, y) dy
0

43. Algunas veces las ecuaciones diferenciales se re-

En los problemas 31 a 36 resuelva la ecuación diferencial dada
suelven con una idea brillante. Este es un pe-
determinando, como en el ejemplo 4, un factor integrante ade-
queño ejercicio de inteligencia: aunque la ecuación
cuado. (x  1x2  y2) dx  y dy
0 no es exacta, demuestre
31. (2y 2  3x) dx  2xy dy
0 cómo el reacomodo (x dx  y dy)
1x2  y2
dx y la
observación 12 d(x 2  y 2)
x dx  y dy puede conducir a
32. y(x  y  1) dx  (x  2y) dy
0 una solución.
33. 6xy dx  (4y  9x 2) dy
0 44. Verdadero o falso: toda ecuación de primer orden separa-
ble dy
dx
g(x)h( y) es exacta.
34. cos x dx  1   2
y
sen x dy
0

35. (10  6y  e3x ) dx  2 dy
0

Modelos matemáticos
36. (y 2  xy 3) dx  (5y 2  xy  y 3 sen y) dy
0
45. Cadena cayendo Una parte de una cadena de 8 pies de
En los problemas 37 y 38 resuelva el problema con valores longitud está enrollada sin apretar alrededor de una cla-
iniciales determinando, como en el ejemplo 5, un factor inte- vija en el borde de una plataforma horizontal y la parte
grante adecuado. restante de la cadena cuelga descansando sobre el borde
37. x dx  (x 2y  4y) dy
0, y(4)
0 de la plataforma. Vea la figura 2.4.2. Suponga que la lon-
gitud de la cadena que cuelga es de 3 pies, que la cadena
38. (x 2  y 2  5) dx
(y  xy) dy, y(0)
1 pesa 2 lb
pie y que la dirección positiva es hacia abajo.
39. a) Demuestre que una familia de soluciones uniparamé- Comenzando en t
0 segundos, el peso de la cadena que
trica de soluciones de la ecuación cuelga causa que la cadena sobre la plataforma se desenro-
lle suavemente y caiga al piso. Si x(t) denota la longitud de
(4xy  3x 2) dx  (2y  2x 2) dy
0 la cadena que cuelga de la mesa al tiempo t  0, entonces
es x 3  2x 2y  y 2
c. v
dx
dt es su velocidad. Cuando se desprecian todas las

Página 36 de 48

08367_02_ch02_p034-081-ok.indd 69 6/4/09 12:16:26 PM

 fuerzas de resistencia se puede demostrar que un modelo Tarea para el laboratorio de computación
matemático que relaciona a v con x está dado por
46. Líneas de flujo
dv a) La solución de la ecuación diferencial
xv  v2
32x.
dx
a) Rescriba este modelo en forma diferencial. Proceda
como en los problemas 31 a 36 y resuelva la ED para
2xy
(x2  y2 ) 2
dx  1
y2  x2
(x2  y2) 2 
dy
0

v en términos de x determinando un factor integrante es una familia de curvas que se pueden interpretar
adecuado. Determine una solución explícita v(x). como líneas de flujo de un fluido que discurre alrede-
b) Determine la velocidad con que la cadena abandona dor de un objeto circular cuya frontera está descrita
la plataforma. por la ecuación x2  y2
1. Resuelva esta ED y ob-
clavija serve que la solución f (x, y)
c para c
0.
b) Use un SAC para dibujar las líneas de flujo para c
0,

0.2,
0.4,
0.6 y
0.8 de tres maneras diferentes.
borde de la
plataforma Primero, utilice el contourplot de un SAC. Segundo,
despeje x en términos de la variable y. Dibuje las dos
x(t)
funciones resultantes de y para los valores dados de
c, y después combine las gráficas. Tercero, utilice el
SAC para despejar y de una ecuación cúbica en térmi-
FIGURA 2.4.2 Cadena desenrollada del problema 45. nos de x.

2.5 SOLUCIONES POR SUSTITUCIÓN

REPASO DE MATERIAL
● Técnicas de integración.

● Separación de variables.

● Solución de ED.

INTRODUCCIÓN Normalmente resolvemos una ecuación diferencial reconociéndola dentro de

una cierta clase de ecuaciones (digamos separables, lineales o exactas) y después aplicamos un proce-
dimiento, que consiste en pasos matemáticos específicos para el tipo de ecuación que nos conducen
a la solución de la misma. Pero no es inusual que nos sorprenda el tener una ecuación diferencial que
no pertenece a alguna de las clases de ecuaciones que sabemos cómo resolver. Los procedimientos
que se analizan en esta sección pueden ser útiles en este caso.

SUSTITUCIONES Con frecuencia el primer paso para resolver una ecuación diferen-
cial es transformarla en otra ecuación diferencial mediante una sustitución. Por ejemplo,
suponga que se quiere transformar la ecuación diferencial de primer orden dy
dx
f (x,
y) sustituyendo y
g(x, u), donde u se considera una función de la variable x. Si g tiene
primeras derivadas parciales, entonces, usando la regla de la cadena
dy g dx g du dy du
obtenemos
gx (x, u) gu(x, u) .
dx x dx u dx dx dx
Al sustituir dy
dx por la derivada anterior y sustituyendo y en f(x, y) por g (x, u), obte-
du
nemos la ED dy
dx
f (x, y) que se convierten en g x (x, u)  gu(x, u)
f (x, g (x, u)), la
du du dx
cual, resuelta para , tiene la forma
F(x, u). Si podemos determinar una solu-
dx dx
ción u
f(x) de esta última ecuación, entonces una solución de la ecuación diferen-
cial original es y(x)
g(x, f(x)).
En el análisis siguiente examinaremos tres clases diferentes de ecuaciones dife-
renciales de primer orden que se pueden resolver mediante una sustitución.

Página 37 de 48

08367_02_ch02_p034-081-ok.indd 70 6/4/09 12:16:28 PM

 ECUACIONES HOMÓGENEAS Si una función f tiene la propiedad f (tx, ty)

t a f (x, y) para algún número real a, entonces se dice que es una función homogénea de
grado a. Por ejemplo f (x, y)
x 3  y 3 es una función homogénea de grado 3, ya que
f (tx, ty)
(tx) 3  (ty) 3
t 3(x 3  y 3)
t 3f (x, y),

mientras que f (x, y)
x 3  y 3  1 es no homogénea. Una ED de primer orden en

forma diferencial

M(x, y) dx  N(x, y) dy
0 (1)

se dice que es homogénea* si ambas funciones coeficientes M y N son ecuaciones ho-

mogéneas del mismo grado. En otras palabras, la ecuación (1) es homogénea si

M(tx, ty)
t␣M(x, y) y N(tx, ty) = t␣N(x, y).

Además, si M y N son funciones homogéneas de grado a, podemos escribir

M(x, y)
x␣M(1, u) y N(x, y)
x␣N(1, u) donde u
y/x, (2)

M(x, y)
y␣M(v, 1) y N(x, y)
y␣N(v, 1) donde v
x/y. (3)

Vea el problema 31 de los ejercicios 2.5. Las propiedades (2) y (3) sugieren las sus-
tituciones que se pueden usar para resolver una ecuación diferencial homogénea. En
concreto, cualquiera de las sustituciones y
ux o x
vy, donde u y v son las nuevas
variables dependientes, reducirán una ecuación homogénea a una ecuación diferencial
de primer orden separable. Para mostrar esto, observe que como consecuencia de (2)
una ecuación homogénea M(x, y)dx  N(x, y)dy
0 se puede reescribir como

x␣M(1, u) dx  x␣N(1, u) dy
0 o bien M(1, u) dx  N(1, u) dy
0,

donde u
y
x o y
ux. Sustituyendo la diferencial dy
u dx  x du en la última

ecuación y agrupando términos, obtenemos una ED separable en las variables u y x:

M(1, u) dx  N(1, u)[u dx  x du]
0

[M(1, u)  uN(1, u)] dx  xN(1, u) du
0
dx N(1, u) du
o 
0.
x M(1, u)  uN(1, u)

En este momento le damos el mismo consejo que en las secciones anteriores. No memo-
rice nada de aquí (en particular la última fórmula); más bien, cada vez siga el procedi-
miento. Pruebe a partir de la ecuación (3) que las sustituciones x
vy y dx
v dy  y dv
también conducen a una ecuación separable siguiendo un procedimiento similar.

EJEMPLO 1 Solución de una ED homogénea

Resuelva (x 2  y 2) dx  (x 2  xy) dy
0.

SOLUCIÓN Examinando a M(x, y)
x 2  y 2 y a N(x, y)
x 2  xy se muestra que

estas funciones coeficientes son homogéneas de grado 2. Si hacemos y
ux, entonces

*
Aquí la palabra homogénea no significa lo mismo que en la sección 2.3. Recuerde que una ecuación lineal
de primer orden a1(x)y a 0 (x)y g(x) es homogénea cuando g(x)
0.

Página 38 de 48

08367_02_ch02_p034-081-ok.indd 71 6/4/09 12:16:29 PM

 dy
u dx  x du, de modo que después de sustituir, la ecuación dada se convierte en
(x2 u2x2) dx (x2 ux2)[u dx x du] 0
x2 (1 u) dx x3(1 u) du 0
1 u dx
du 0
1 u x
2 dx
du 1 0.  división larga
1 u x
Después de integrar la última ecuación se obtiene
u  2 ln 1  u   ln x 
ln c 
y y
 
  2 ln 1   ln x 
ln c.
x x
; sustituyendo de nuevo u
y
x

Utilizando las propiedades de los logaritmos, podemos escribir la solución anterior como

(x y) 2 y
ln o (x y) 2 cxey/x.
cx x

Aunque cualquiera de las soluciones indicadas se puede usar en toda ecuación

diferencial homogénea, en la práctica se intenta con x
vy cuando la función M(x, y)
sea más fácil que N(x, y). También podría ocurrir que después de utilizar una sustitu-
ción, podemos encontrar integrales que son difíciles o imposibles de evaluar en forma
cerrada; y el cambiar las sustituciones puede facilitar el problema.

ECUACIÓN DE BERNOULLI La ecuación diferencial

dy
 P(x)y
f (x)y n, (4)
dx

donde n es cualquier número real, se llama ecuación de Bernoulli. Observe que para
n
0 y n
1, la ecuación (4) es lineal. Para n ã 0 y n ã 1 la sustitución u
y 1n
reduce cualquier ecuación de la forma (4) a una ecuación lineal.

EJEMPLO 2 Solución de una ED de Bernoulli

dy
Resuelva x  y
x 2 y 2.
dx

SOLUCIÓN Primero reescribimos la ecuación como

dy 1
 y
xy 2
dx x

al dividir entre x. Con n
2 tenemos u
y1 o y
u1. Entonces sustituimos

dy dy du du


u2 ; Regla de la cadena
dx du dx dx

en la ecuación dada y simplificando. El resultado es

du 1
 u
x.
dx x

Página 39 de 48

08367_02_ch02_p034-081-ok.indd 72 6/4/09 12:16:30 PM

 El factor integrante para esta ecuación lineal en, digamos, (0, ) es
1
ed x/x
eln x
eln x
x1.
d 1
Integrando [x u]
1
dx

se obtiene x1u
x  c o u
x 2  cx. Puesto que u
y1, tenemos que y
1
u,

así una solución de la ecuación dada es y
1
(x 2  cx).

Observe que no hemos obtenido una solución general de la ecuación diferencial

no lineal original del ejemplo 2 ya que y
0 es una solución singular de la ecuación.

REDUCCIÓN A SEPARACIÓN DE VARIABLES Una ecuación diferencial de la

forma
dy

f (Ax  By  C) (5)
dx
Se puede siempre reducir a una ecuación con variables separables por medio de la
sustitución u
Ax  By  C, B ã 0. El ejemplo 9 muestra la técnica.

EJEMPLO 3 Un problema con valores iniciales

dy
Resuelva
(2x  y) 2  7, y(0)
0.
dx

SOLUCIÓN Si hacemos u
2x  y, entonces du
dx
2  dy
dx, por lo que la

ecuación diferencial se expresa como
du du
 2
u2  7 o
u 2  9.
dx dx

La última ecuación es separable. Utilizando fracciones parciales

du
(u  3)(u  3)

dx o
1 1


1
6 u3 u3
du
dx 
y después de integrar se obtiene

1 u 3 u 3  sustitu yendo e6c1 por c

ln x c1 o e6x 6c1
ce6x.
y 6 u 3 u 3

Despejando u de la última ecuación y resustituyendo a u en términos de x y y, se ob-

tiene la solución
x 3(1  ce6x ) 3(1  ce6x)
u
o y
2x  . (6)
1  ce6x 1  ce6x

Por último, aplicando la condición inicial y(0)
0 a la última ecuación en (6) se ob-
tiene c
1. La figura 2.5.1, obtenida con la ayuda de un programa de graficación,
3(1  e6x)
muestra en azul oscuro la gráfica de la solución particular y
2x  junto
FIGURA 2.5.1 Algunas soluciones de 1  e6x
y
(2x  y) 2  7. con las gráficas de algunos otros miembros de la familia de soluciones (6).

Página 40 de 48

08367_02_ch02_p034-081-ok.indd 73 6/4/09 12:16:31 PM

 EJERCICIOS 2.5
Cada una de las ED de los problemas 1-14 es homogénea. Cada una de las ED de los problemas 23 a 30 es de la forma
En los problemas 1 a 10 resuelva la ecuación diferencial dada dada en la ecuación (5).
usando las sustituciones adecuadas. En los problemas 23 a 28 resuelva la ecuación diferencial
dada usando una sustitución adecuada.
1. (x  y) dx  x dy
0 2. (x  y) dx  x dy
0
dy dy 1  x  y
3. x dx  ( y  2x) dy
0 4. y dx
2(x  y) dy 23.
(x  y  1) 2 24.

dx dx xy
5. (y 2  yx) dx  x 2 dy
0
dy dy
6. (y 2  yx) dx  x 2 dy
0 25.
tan2 (x  y) 26.
sen(x  y)
dx dx
dy y  x
7.
dy dy
dx y  x 27.
2  1y  2x  3 28.
1  eyx5
dx dx
dy x  3y
8.
En los problemas 29 y 30 resuelva el problema con valores
dx 3x  y
iniciales dado.
(
9. y dx  x  1xy dy
0 ) dy

cos(x  y), y(0)
>4
29.
dx
dy
10. x
y  1x2  y2, x0 dy 3x  2y
dx 30.
, y (1)
1
dx 3x  2y  2
En los problemas 11 a 14 resuelva el problema con valores
iniciales dado.
dy Problemas para analizar
11. xy
2

y3  x3, y(1)
2
dx 31. Explique por qué es posible expresar cualquier ecuación di-
dx ferencial homogénea M(x, y) dx  N(x, y) dy
0 en la forma
12. (x  2y )
xy, y(1)
1
2 2


dy dy y

F .
13. (x  ye y
x) dx  xe y
x dy
0, y(1)
0 dx x

14. y dx  x(ln x  ln y  1) dy
0, y(1)
e Podría comenzar por demostrar que

M(x, y)
x␣M(1, y/x) y N(x, y)
x␣N(1, y/x).
Cada una de las ED de los problemas 15 a 22 es una ecuación
32. Ponga la ecuación diferencial homogénea
de Bernoulli.
(5x 2  2y 2) dx  xy dy
0
En los problemas 15 a 20 resuelva la ecuación diferencial
dada usando una sustitución adecuada. en la forma dada en el problema 31.

dy 1 dy 33. a) Determine dos soluciones singulares de la ED en el

15. x y
2 16.  y
ex y2
dx y dx problema 10.
b) Si la condición inicial y(5)
0 es como se indicó para
dy dy el problema 10, entonces ¿cuál es el intervalo I de de-
17.
y (xy 3  1) 18. x  (1  x)y
xy2
dx dx finición más grande en el cual está definida la solu-
ción? Utilice un programa de graficación para obtener
dy dy la gráfica de la curva solución para el PVI.
19. t2  y2
ty 20. 3(1  t2)
2ty( y3  1)
dt dt
34. En el ejemplo 3 la solución y(x) es no acotada conforme
En los problemas 21 y 22 resuelva el problema con valores x :
. Sin embargo, y(x) es asintótica a una curva con-
iniciales dado. forme x :  y a una diferente curva conforme x : .
¿Cuáles son las ecuaciones de estas curvas?
dy
21. x2  2xy
3y4, y(1)
12 35. La ecuación diferencial dy
dx
P(x)  Q(x)y  R(x)y2
dx
se conoce como la ecuación de Riccati.
dy a) Una ecuación de Riccati se puede resolver por dos
22. y1/2  y3/2
1, y(0)
4 sustituciones consecutivas, siempre y cuando conoz-
dx

Página 41 de 48

08367_02_ch02_p034-081-ok.indd 74 6/4/09 12:16:32 PM

 camos una solución particular, y1, de la ecuación. de una cadena que se desliza por el borde de una plata-
Muestre que la sustitución y
y1  u reduce la ecua- forma horizontal es
ción de Riccati a una ecuación de Bernoulli (4) con
dv
n
2. La ecuación de Bernoulli se puede entonces xv  v 2
32x.
reducir a una ecuación lineal sustituyendo w
u1. dx
b) Determine una familia uniparamétrica de soluciones En ese problema se le pidió que resolviera la ED convir-
de la ecuación diferencial tiéndola en una ecuación exacta usando un factor inte-
grante. Esta vez resuelva la ED usando el hecho de que es
dy 4 1

 2  y  y2 una ecuación de Bernoulli.
dx x x
38. Crecimiento de la población En el estudio de la pobla-
donde y1
2
x es una solución conocida de la ecua-
ción dinámica uno de los más famosos modelos para un
ción.
crecimiento poblacional limitado es la ecuación logística
36. Determine una sustitución adecuada para resolver
dP
xy
y ln(xy).
P(a  bP),
dt

donde a y b son constantes positivas. Aunque retomaremos

Modelos matemáticos
esta ecuación y la resolveremos utilizando un método al-
37. Cadena cayendo En el problema 45 de los ejercicios ternativo en la sección 3.2, resuelva la ED por esta primera
2.4 vimos que un modelo matemático para la velocidad v vez usando el hecho de que es una ecuación de Bernoulli.

2.6 UN MÉTODO NUMÉRICO

INTRODUCCIÓN Una ecuación diferencial dy
dx
f (x, y) es una fuente de información. Comen-
zaremos este capítulo observando que podríamos recolectar información cualitativa de una ED de
primer orden respecto a sus soluciones aun antes de intentar resolver la ecuación. Entonces en las sec-
ciones 2.2 a 2.5 examinamos a las ED de primer orden analíticamente, es decir, desarrollamos algunos
procedimientos para obtener soluciones explícitas e implícitas. Pero una ecuación diferencial puede
tener una solución aun cuando no podamos obtenerla analíticamente. Así que para redondear el esquema
de los diferentes tipos de análisis de las ecuaciones diferenciales, concluimos este capítulo con un mé-
todo con el cual podemos “resolver” la ecuación diferencial numéricamente; esto significa que la ED se
utiliza como el principio básico de un algoritmo para aproximar a la solución desconocida.
En esta sección vamos a desarrollar únicamente el más sencillo de los métodos numéricos, un
método que utiliza la idea de que se puede usar una recta tangente para aproximar los valores de una
función en una pequeña vecindad del punto de tangencia. En el capítulo 9 se presenta un tratamiento
más extenso de los métodos numéricos.

USANDO LA RECTA TANGENTE Suponemos que el problema con valores iniciales

y′
f (x, y), y(x0)
y0 (1)
tiene una solución. Una manera de aproximar esta solución es usar rectas tangentes. Por
ejemplo, sea que y(x) denote la solución incógnita para el problema con valores inicia-
les y 0.1 1y 0.4x2, y(2) 4. La ecuación diferencial no lineal en este PVI no
se puede resolver directamente por cualquiera de los métodos considerados en las sec-
ciones 2.2, 2.4 y 2.5; no obstante, aún podemos encontrar valores numéricos aproxi-
mados de la incógnita y(x). En concreto, supongamos que deseamos conocer el valor
de y(2, 5). El PVI tiene una solución y como el flujo del campo direccional de la ED
en la figura 2.6.1a sugiere, una curva solución debe tener una forma similar a la curva
que se muestra en azul.
El campo direccional de la figura 2.6.1a se generó con elementos lineales que pasan
por puntos de una malla de coordenadas enteras. Puesto que la curva solución pasa por el

Página 42 de 48

08367_02_ch02_p034-081-ok.indd 75 6/4/09 12:16:34 PM

 punto inicial (2, 4), el elemento lineal en este punto es una recta tangente con pendiente
dada por f (2, 4) 0.1 4 0.4(2) 2 1.8. Como se muestra en la figura 2.6.1a y el
“zoom in” (acercamiento) de la figura 2.6.1b, cuando x está cerca de 2, los puntos en la
curva solución están cerca de los puntos de la recta tangente (el elemento lineal). Utilizando
el punto (2, 4), la pendiente f (2, 4)
1.8 y la forma punto pendiente de una recta, encon-
tramos que una ecuación de la recta tangente es y
L(x), donde L(x)
1.8x  0.4. Esta
última ecuación se llama linealización de y(x) en x
2 que se puede utilizar para aproxi-
mar los valores dentro de una pequeña vecindad de x
2. Si y1
L(x1) denota la coorde-
nada y en la recta tangente y y(x1) es la coordenada y de la curva solución correspondiente
a una coordenada x, x1 que está cerca de x
2, entonces y(x1)  y1. Si elegimos, x1
2.1,
entonces y1
L(2.1)
1.8(2.1)  0.4
4.18, entonces y(2.1)  4.18.

curva
4 solución

(2, 4) pendiente
2 m = 1.8

x
_2 2
a) campo direccional para y  0. b) elemento lineal
en (2, 4).
FIGURA 2.6.1 Amplificación de una vecindad del punto (2, 4).

y MÉTODO DE EULER Para generalizar el procedimiento que acabamos de ilustrar,

curva solución
usamos la linealización de una solución incógnita y(x) de (1) en x
x0:
L(x)
y0  f (x0 , y0)(x  x0). (2)
(x1, y(x1))
La gráfica de esta linealización es una recta tangente a la gráfica de y
y (x) en el punto
error
(x0, y0). Ahora hacemos que h sea un incremento positivo del eje x, como se muestra en
(x1, y1)
(x0, y0) la figura 2.6.2. Entonces sustituyendo x por x1
x0  h en la ecuación (2), obtenemos
pendiente = f(x0, y0)
L(x1)
y0  f (x0, y0)(x0  h  x0) o y 1
y0  hf(x1, y1),
donde y1
L(x1). El punto (x1, y1) en la recta tangente es una aproximación del
L(x) h
punto (x1, y(x1)) sobre la curva solución. Por supuesto, la precisión de la aproxima-
x0 x1 = x 0 + h x
ción L(x1)  y(x1) o y1  y(x1) depende fuertemente del tamaño del incremento h.
Normalmente debemos elegir este tamaño de paso para que sea “razonablemente
FIGURA 2.6.2 Aproximación de y(x1)
usando una recta tangente. pequeño”. Ahora repetimos el proceso usando una segunda “recta tangente” en (x1,
y1).* Identificando el nuevo punto inicial como (x1, y1) en lugar de (x0, y0) del análisis
anterior, obtenemos una aproximación y2  y(x 2) correspondiendo a dos pasos de lon-
gitud h a partir de x0, es decir, x 2
x1  h
x 0  2h, y
y(x2)
y(x0  2h)
y(x1  h)  y2
y1  hf (x1, y1).
Continuando de esta manera, vemos que y1, y2, y3, . . . , se puede definir recursivamente
mediante la fórmula general
yn1
yn  hf (xn, yn), (3)
donde x n
x 0  nh, n
0, 1, 2, . . . Este procedimiento de uso sucesivo de las “rectas
tangentes” se llama método de Euler.

*
Esta no es una recta tangente real, ya que (x1, y1) está sobre la primera tangente y no sobre la curva solución.

Página 43 de 48

08367_02_ch02_p034-081-ok.indd 76 6/4/09 12:16:35 PM

 EJEMPLO 1 Método de Euler

Considere el problema con valores iniciales y
0.1 1y  0.4x2, y(2)
4 Utilice

el método de Euler para obtener una aproximación de y(2.5) usando primero h
0.1
y después h
0.05.
TABLA 2.1 h
0.1
xn yn SOLUCIÓN Con la identificación f (x, y)
0.11y  0.4x2 la ecuación (3) se con-
vierte en
2.00 4.0000
2.10 4.1800
(
yn1
yn  h 0.11yn  0.4x2n . )
2.20 4.3768 Entonces para h
0.1, x0
2, y0
4 y n
0 encontramos
2.30 4.5914
2.40 4.8244 ( ) (
y1
y0  h 0.11y0  0.4x20
4  0.1 0.114  0.4(2) 2
4.18, )
2.50 5.0768
que, como ya hemos visto, es una estimación del valor y(2.1). Sin embargo, si usamos el
paso de tamaño más pequeño h
0.05, le toma dos pasos alcanzar x
2.1. A partir de
(
y1
4  0.05 0.114  0.4(2)2
4.09 )
TABLA 2.2 h
0.05
xn yn
(
y2
4.09  0.05 0.114.09  0.4(2.05)2
4.18416187 )
tenemos y1  y(2.05) y y 2  y(2.1). El resto de los cálculos fueron realizados usando
2.00 4.0000
un paquete computacional. En las tablas 2.1 y 2.2 se resumen los resultados, donde
2.05 4.0900
cada entrada se ha redondeado a cuatro lugares decimales. Vemos en las tablas 2.1 y
2.10 4.1842
2.15 4.2826
2.2 que le toma cinco pasos con h
0.1 y 10 pasos con h
0.05, respectivamente,
2.20 4.3854 para llegar a x
2.5. Intuitivamente, esperaríamos que y10
5.0997 correspondiente
2.25 4.4927 a h
0.05 sea la mejor aproximación de y(2.5) que el valor y5
5.0768 correspon-
2.30 4.6045 diente a h
0.1.
2.35 4.7210
2.40 4.8423 En el ejemplo 2 aplicamos el método de Euler para una ecuación diferencial para
2.45 4.9686 la que ya hemos encontrado una solución. Hacemos esto para comparar los valores de
2.50 5.0997 las aproximaciones yn en cada caso con los valores verdaderos o reales de la solución
y(xn) del problema con valores iniciales.

EJEMPLO 2 Comparación de los valores aproximados y reales

Considere el problema con valores iniciales y 
0.2xy, y(1)
1. Utilice el método de Euler
para obtener una aproximación de y (1.5) usando primero h
0.1 y después h
0.05.

SOLUCIÓN Con la identificación f (x, y)
0.2xy, la ecuación (3) se convierte en

yn1
yn  h(0.2xnyn)
donde x 0
1 y y 0
1. De nuevo con la ayuda de un paquete computacional obtenga
los valores de las tablas 2.3 y 2.4.
TABLA 2.4 h
0.05
xn yn Valor real Error absoluto % Error relativo

1.00 1.0000 1.0000 0.0000 0.00

1.05 1.0100 1.0103 0.0003 0.03
TABLA 2.3 h
0.1 1.10 1.0206 1.0212 0.0006 0.06
1.15 1.0318 1.0328 0.0009 0.09
xn yn Valor real Error absoluto % Error relativo
1.20 1.0437 1.0450 0.0013 0.12
1.00 1.0000 1.0000 0.0000 0.00 1.25 1.0562 1.0579 0.0016 0.16
1.10 1.0200 1.0212 0.0012 0.12 1.30 1.0694 1.0714 0.0020 0.19
1.20 1.0424 1.0450 0.0025 0.24 1.35 1.0833 1.0857 0.0024 0.22
1.30 1.0675 1.0714 0.0040 0.37 1.40 1.0980 1.1008 0.0028 0.25
1.40 1.0952 1.1008 0.0055 0.50 1.45 1.1133 1.1166 0.0032 0.29
1.50 1.1259 1.1331 0.0073 0.64 1.50 1.1295 1.1331 0.0037 0.32

Página 44 de 48

08367_02_ch02_p034-081-ok.indd 77 6/4/09 12:16:36 PM

 En el ejemplo 1 se calcularon los valores verdaderos o reales de la solución cono-
cida y
e0.1(x −1) . (Compruebe.) El error absoluto se define como
2

valor real – aproximado .

El error relativo y el error relativo porcentual son, respectivamente,

error absoluto error absoluto × 100.
y
valor real valor real
Es evidente de las tablas 2.3 y 2.4 que la precisión de las aproximaciones mejora
conforme disminuye el tamaño del paso h. También nosotros vemos esto aun cuando
el error relativo porcentual esté creciendo en cada paso, no parece estar mal. Pero no
debe engañarse por un ejemplo. Si simplemente cambiamos el coeficiente del lado de-
recho de la ED del ejemplo 2 de 0.2 a 2 entonces en xn
1.5 los errores relativos por-
centuales crecen dramáticamente. Véase el problema 4 del ejercicio 2.6.

UNA ADVERTENCIA El método de Euler es sólo uno de los diferentes métodos en

los que se puede aproximar una solución de una ecuación diferencial. Aunque por su
sencillez es atractivo, el método de Euler rara vez se usa en cálculos serios. Aquí se ha
presentado sólo para dar un primer esbozo de los métodos numéricos. En el capítulo 9
trataremos en detalle el análisis de los métodos numéricos que tienen mucha precisión,
en especial el método de Runge-Kutta conocido como el método RK4.

SOLUCIONADORES NUMÉRICOS Independientemente de si se puede realmente

encontrar una solución explícita o implícita, si existe una solución de una ecuación
diferencial, ésta se representa por una curva suave en el plano cartesiano. La idea bá-
sica detrás de cualquier método numérico para las ecuaciones diferenciales ordinarias
de primer orden es de alguna manera aproximar los valores de y de una solución para
valores de x preseleccionados. Comenzamos con un punto inicial dado (x0, y0) de una
curva solución y procedemos a calcular en un modelo paso por paso una secuencia
de puntos (x1, y1), (x2, y2),…, (xn, yn) cuyas coordenadas y, yi se aproximan a las coor-
denadas y, y(xi) de los puntos (x1, y(x1)), (x2, y(x2)), …, (xn, y(xn)) que yacen sobre la
gráfica de la solución normalmente desconocida y(x). Tomando las coordenadas x más
cercanas (es decir, para valores pequeños de h) y uniendo los puntos (x1, y1), (x2, y2),…,
(xn, yn) con segmentos de recta cortos, obtenemos una curva poligonal cuyas caracte-
rísticas cualitativas esperamos sean cercanas a las de una curva solución real. El dibujo
de curvas es muy adecuado en una computadora. A un programa de cómputo escrito
para implementar un método numérico o para mostrar una representación visual de
una solución aproximada que ajusta los datos numéricos producidos por este segundo
y método se le conoce como un solucionador numérico. Comercialmente hay disponi-
5
bles muchos solucionadores numéricos ya sea que estén integrados en un gran paquete
método computacional, tal como en un sistema algebraico computacional o que sean un pa-
4
RK4
solución quete autónomo. Algunos paquetes computacionales simplemente dibujan las aproxi-
3 exacta maciones numéricas generadas, mientras que otros generan pesados datos numéricos
así como la correspondiente aproximación o curvas solución numérica. En la figura
2
2.6.3 se presenta a manera de ilustración la conexión natural entre los puntos de las
(0,1)
1 método gráficas producidas por un solucionador numérico, las gráficas poligonales pintadas
Euler con dos colores son las curvas solución numérica para el problema con valores inicia-
x les y 
0.2xy, y(0)
1 en el intervalo [0, 4] obtenidas de los métodos de Euler y RK4
_1 usando el tamaño de paso h
1. La curva suave en azul es la gráfica de la solución
_1 1 2 3 4 5 exacta y
e0.1x del PVI. Observe en la figura 2.6.3 que, aun con el ridículo tamaño
2

FIGURA 2.6.3 Comparación de los de paso de h
1, el método RK4 produce la “curva solución” más creíble. La curva
métodos de Runge-Kutta (RK4) y de solución numérica obtenida del método RK4 es indistinguible de la curva solución real
Euler. en el intervalo [0, 4] cuando se usa el tamaño de paso usual de h
0.1.

Página 45 de 48

08367_02_ch02_p034-081-ok.indd 78 6/4/09 12:16:37 PM

 y USANDO UN SOLUCIONADOR NUMÉRICO No es necesario conocer los di-
6 ferentes métodos numéricos para utilizar un solucionador numérico. Un solucionador
5 usualmente requiere que la ecuación diferencial se pueda expresar en la forma normal
dy
dx
f (x, y). Los solucionadores numéricos que sólo generan curvas requieren que se
4
les proporcione f (x, y) y los datos iniciales x0 y y0 y que se indique el método numérico
3 deseado. Si la idea es aproximarse al valor numérico de y(a), entonces un solucionador
2 numérico podría requerir además expresar un valor de h o, del mismo modo, dar el nú-
1
mero de pasos que quiere tomar para llegar de x
x0 a x
a. Por ejemplo, si queremos
aproximar y(4) para el PVI que se muestra en la figura 2.6.3, entonces, comenzando en
x x
0 le tomaría cuatro pasos llegar a x
4 con un tamaño de paso de h
1; 40 pasos
_1 son equivalentes a un tamaño de paso de h
0.1. Aunque aquí no investigaremos todos
_2 _1 1 2 3 4 5 los problemas que se pueden encontrar cuando se intenta aproximar cantidades matemá-
FIGURA 2.6.4 Una curva solución ticas, al menos debe estar consciente del hecho de que el solucionador numérico puede
que no ayuda mucho. dejar de funcionar cerca de ciertos puntos o dar una incompleta o engañosa imagen
cuando se aplica a ciertas ecuaciones diferenciales en la forma normal. La figura 2.6.4
muestra la gráfica que se obtuvo al aplicar el método de Euler a un problema con valores
iniciales de primer orden dy
dx
f (x, y), y(0)
1. Se obtuvieron resultados equiva-
lentes utilizando tres diferentes solucionadores numéricos, sin embargo la gráfica di-
fícilmente es una posible curva solución. (¿Por qué?) Hay diferentes caminos de solución
cuando un solucionador numérico tiene dificultades; las tres más obvias son disminuir el
tamaño del paso, usar otro método numérico e intentar con un solucionador diferente.

EJERCICIOS 2.6
En los problemas 1 y 2 use el método de Euler para obtener En los problemas 11 y 12 utilice un solucionador para obtener
una aproximación a cuatro decimales del valor indicado, una curva solución numérica para el problema con valores iniciales
ejecute a mano la ecuación de recursión (3), usando primero dado. Primero utilice el método de Euler y después el método RK4.
h
0.1 y después usando h
0.05. Utilice h
0.25 en cada caso. Superponga ambas curvas solución
en los mismos ejes coordenados. Si es posible, utilice un color
1. y
2x  3y  1, y(1)
5; y(1.2)
diferente para cada curva. Repita, usando h
0.1 y h
0.05.
2. y
x  y 2, y(0)
0; y(0.2)
11. y
2(cos x)y, y(0)
1
En los problemas 3 y 4 use el método de Euler para obte-
ner una aproximación a cuatro decimales del valor indicado. 12. y
y(10  2y), y(0)
1
Primero utilice h
0.1 y después utilice h
0.05. Determine
una solución explícita para cada problema con valores inicia-
les y después construya tablas similares a las tablas 2.3 y 2.4. Problemas para analizar
3. y
y, y(0)
1; y(1.0) 13. Use un solucionador numérico y el método de Euler para
aproximar y(0.1), donde y(x) es la solución de y
2xy 2,
4. y
2xy, y(1)
1; y(1.5)
y(0)
1. Primero use h
0.1 y después use h
0.05.
En los problemas 5 a 10 use un solucionador numérico y el Repita, usando el método RK4. Analice qué podría cau-
método de Euler para obtener una aproximación a cuatro de- sar que las aproximaciones a y(1.0) difieran mucho.
cimales del valor indicado. Primero utilice h
0.1 y después
utilice h
0.05.
Tarea para el laboratorio de computación
5. y
e , y(0)
0; y(0.5)
y
14. a) Utilice un solucionador numérico y el método RK4
6. y
x  y , y(0)
1;
2 2
y(0.5) para trazar la gráfica de la solución del problema con
valores iniciales y
2xy  1, y(0)
0.
7. y
(x  y) 2, y(0)
0.5; y(0.5) b) Resuelva el problema con valores iniciales por uno de
los procedimientos analíticos desarrollados en las
8. y
xy  1y, y(0)
1; y(0.5)
secciones anteriores en este capítulo.
y c) Use la solución analítica y(x) que encontró en el in-
9. y
xy 2  , y(1)
1; y(1.5)
x ciso b) y un SAC para determinar las coordenadas de
10. y
y  y 2, y(0)
0.5; y(0.5) todos los extremos relativos.

Página 46 de 48

08367_02_ch02_p034-081-ok.indd 79 6/4/09 12:16:37 PM

 REPASO DEL CAPÍTULO 2
Responda los problemas 1 a 4 sin consultar las respuestas del libro. f
Llene los espacios en blanco o responda si es verdadero o falso.

1. La ED lineal, y   ky
A, donde k y A son constantes,

es autónomo. El punto crítico de la ecuación 1
es un (atractor o repulsor) para k  0 y un 1 P
(atractor o repulsor) para k 0.
dy
2. El problema x  4y
0, y(0)
k , tiene un número
dx
infinito de soluciones para k
y no tiene so-
lución para k
. FIGURA 2.R.3 Gráfica del problema 8.
3. La ED lineal, y   k1y
k2, donde k1 y k2 son constantes
distintas de cero, siempre tiene una solución constante. 9. La figura 2.R.4 es una parte de un campo direccional de
una ecuación diferencial dy
dx
f (x, y). Dibuje a mano
dos diferentes curvas solución, una que es tangente al ele-
4. La ED lineal, a1(x)y  a2(x)y
0 es también separable. mento lineal que se muestra en negro y el otro que es tan-
gente al elemento lineal que se muestra de color (rojo).

En los problemas 5 y 6 construya una ecuación diferencial de

primer orden dy
dx
f (y) cuyo esquema de fase es consis-
tente con la figura dada.

5. y

FIGURA 2.R.4 Parte de un campo direccional del problema 9.

FIGURA 2.R.1 Gráfica del problema 5. 10. Clasifique cada ecuación diferencial como separable,
6. y exacta, lineal, homogénea o Bernoulli. Algunas ecuacio-
nes pueden ser de más de una clase. No las resuelva.
4
dy x  y dy 1
a)
b)

2 dx x dx y  x

0 dy dy 1
c) (x  1)
y  10 d)

dx dx x(x  y)

FIGURA 2.R.2 Gráfica del problema 6. dy y 2  y dy

e)
f)
5y  y 2
dx x 2  x dx
7. El número 0 es un punto crítico de la ecuación diferen- dy
cial autónoma dx
dt
xn, donde n es un entero positivo. g) y dx
( y  xy 2) dy h) x
ye x/y  x
dx
¿Para qué valores de n es 0 asintóticamente estable?
¿Semiestable? ¿Inestable? Repita para la ecuación dife- i) xy y  y 2
2x j) 2xy y  y 2
2x 2
rencial dx
dt
xn.
k) y dx  x dy
0
8. Considere la ecuación diferencial dP / dt
f (P), donde
f (P)
0.5P3  1.7P  3.4. l) x 2

2y
x dx
(3  ln x ) dy
2

La función f (P) tiene una raíz real, como se muestra en la

figura 2.R.3. Sin intentar resolver la ecuación diferencial, dy x y y dy
m)
 1 n)  e 2x y
0
3 2

estime el valor de límt→ P(t). dx y x x 2 dx

Página 47 de 48

08367_02_ch02_p034-081-ok.indd 80 6/4/09 12:16:38 PM

 En los problemas resuelva la ecuación diferencial dada. y

11. ( y 2  1) dx
y sec2 x dy
12. y(ln x  ln y) dx
(x ln x  x ln y  y) dy
dy
13. (6x  1)y2  3x2  2y3
0 x
dx
dx 4y2  6xy
14.
 2
dy 3y  2x
dQ
15. t  Q
t 4 ln t
dt
FIGURA 2.R.5 Gráfica para el problema 23.
16. (2x  y  1)y
1
17. (x 2  4) dy
(2x  8xy) dx 24. Utilice el método de Euler con tamaño de paso h
0.1
para aproximar y(1.2), donde y(x) es una solución del pro-
18. (2r 2 cos u sen u  r cos u) du blema con valores iniciales y 1 x 1y , y(1) 9.
 (4r  sen u  2r cos2 u) dr
0
En los problemas 25 y 26 cada figura representa una parte de
En los problemas 19 y 20 resuelva el problema con valores un campo direccional de una ecuación diferencial de primer
iniciales dado e indique el intervalo I más largo en el que la orden dy
dx
f (y). Reproduzca esta figura en una hoja y des-
solución está definida. pués termine el campo direccional sobre la malla. Los puntos
de la malla son (mh, nh) donde h
21, m y n son enteros, 7
dy 7  m  7, 7  n  7. En cada campo direccional dibuje a
19. senx (cos x)y 0, y 2
dx 6 mano una curva solución aproximada que pase por cada uno
de los puntos sólidos mostrados en rojo. Analice: ¿parece que
dy
20. 2(t 1)y 2 0, y(0) 1
8
la ED tiene puntos críticos en el intervalo 3.5  m  3.5?
dt Si es así, clasifique los puntos críticos como asintóticamente
21. a) Sin resolver, explique por qué el problema con valores estables, inestables o semiestables.
iniciales
dy 25. y

1y, y(x0)
y0
dx
3
no tiene solución para y0 0. 2
b) Resuelva el problema con valores iniciales del inciso 1
a) para y0  0 y determine el intervalo I más largo en x
el que la solución está definida.
_1
22. a) Determine una solución implícita del problema con _2
valores iniciales
_3
dy y 2  x 2

, y(1)
 12. _3 _2 _1 1 2 3
dx xy

b) Determine una solución explícita del problema del FIGURA 2.R.6 Parte de un campo direccional del problema 25.
inciso a) e indique el intervalo de solución más largo y
de I en el que la solución está definida. Aquí puede 26.
ser útil un programa de graficación. 3

23. En la figura 2.R.5 se presentan las gráficas de algunos miem- 2

bros de una familia de soluciones para una ecuación dife- 1
rencial de primer orden dy
dx
f (x, y). Las gráficas de dos x
soluciones implícitas, una que pasa por el punto (1, 1) y la _1
otra que pasa por (1, 3) se muestran en rojo. Reproduzca
_2
la figura en una hoja. Con lápices de colores trace las curvas
solución para las soluciones y
y1(x) y y
y2(x) definidas _3
por las soluciones implícitas tales como y1(1)
1 y y2(1)
_3 _2 _1 1 2 3

3, respectivamente. Estime los intervalos en los que las
soluciones y
y1(x) y y
y2(x) están definidas. FIGURA 2.R.7 Parte de un campo direccional del problema 26.

Página 48 de 48

08367_02_ch02_p034-081-ok.indd 81 6/4/09 12:16:41 PM