Tema 2 (Ecuaciones Diferenciales Ordinarias de Primer Orden)

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 94

2

Ecuaciones Diferenciales
Ordinarias de Primer Orden

ESQUEMA

2.1. ECUACIÓN DIFERENCIAL ORDINARIA DE PRIMER ORDEN. EXISTENCIA Y


UNICIDAD DE SOLUCIONES
2.1.1. Ecuación diferencial ordinaria de primer orden. Formas de expresión
2.1.2. Teorema básico de existencia y unicidad de soluciones

2.2. ECUACIONES DE VARIABLES SEPARADAS

2.3. ECUACIONES HOMOGÉNEAS


2.3.1. Definición de ecuación homogénea
2.3.2. Caracterización de las ecuaciones homogéneas
2.3.3. Resolución de ecuaciones homogéneas
2.3.4. Ecuaciones reducibles a homogéneas

2.4. ECUACIONES DIFERENCIALES EXACTAS


2.4.1. Definición de ecuación diferencial exacta
2.4.2. Caracterización de las ecuaciones diferenciales exactas
2.4.3. Resolución de ecuaciones diferenciales exactas

2.5. FACTORES INTEGRANTES


2.5.1. Definición de factor integrante
2.5.2. Cálculo de factores integrantes
2.5.3. Resolución de ecuaciones diferenciales mediante factores integrantes
2.5.4. Multiplicidad del factor integrante

2.6. ECUACIÓN LINEAL


2.6.1. Definición de ecuación lineal
2.6.2. Resolución de ecuaciones lineales homogéneas

- 39 -
40 Ecuaciones Diferenciales

2.6.3. Resolución de ecuaciones lineales completas


• Solución general de la ecuación lineal completa
• Resolución por el método del factor integrante
• Resolución por el método de Lagrange o de variación de
constantes
2.6.4. Ecuación de Bernoulli
2.6.5. Ecuación de Riccati

2.7. ECUACIÓN DE LAGRANGE


2.7.1. Definición de la ecuación de Lagrange
2.7.2. Resolución de la ecuación de Lagrange

2.8. ECUACIÓN DE CLAIRAUT


2.8.1. Definición de la ecuación de Clairaut
2.8.2. Resolución de la ecuación de Clairaut

2.9. APLICACIONES GEOMÉTRICAS DE LAS ECUACIONES DIFERENCIALES


ORDINARIAS DE PRIMER ORDEN
2.9.1. Trayectorias ortogonales a una familia de curvas
2.9.2. Trayectorias oblicuas a una familia de curvas

EJERCICIOS
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 41

2.1. ECUACIÓN DIFERENCIAL ORDINARIA DE PRIMER ORDEN.


EXISTENCIA Y UNICIDAD DE SOLUCIONES

2.1.1. Ecuación diferencial ordinaria de primer orden. Formas de


expresión

El estudio cuantitativo de las ecuaciones diferenciales ordinarias debe comenzar


lógicamente por las más sencillas, es decir, las de primer orden. Así, en este tema se
describen varios tipos fundamentales de ecuaciones diferenciales de primer orden, junto
con los métodos analíticos específicos que existen para calcular sus soluciones, y
diferentes ejemplos de aplicaciones. Por tanto, el objetivo del tema contempla tres
aspectos importantes del estudio de las EDO de primer orden: la clasificación en tipos
fundamentales, el desarrollo de técnicas analíticas de resolución apropiadas para cada
tipo, y la descripción de ejemplos prácticos de aplicación.

Definición (Ecuación diferencial ordinaria de primer orden). Se denomina ecuación


diferencial ordinaria de primer orden a una ecuación cuya incógnita es una función
f : D ⊆ o → o de una variable independiente de la forma y = f (x) , y en cuya expresión
aparece la derivada de y de primer orden.

Ejemplo: Ecuación diferencial ordinaria lineal de primer orden:


y ′ + y = −2x + 1

Ejemplo: Ecuación diferencial ordinaria polinómica de primer orden y grado 2:


dy
y − y2 + x +1 = 0
dx

Ejemplo: Ecuación diferencial ordinaria algebraica irracional de primer orden:

y′ = 4 y − x

Definición (Forma general o implícita de una EDO de primer orden). Dada una ecuación
diferencial ordinaria de primer orden de incógnita y (x ) , se llama forma general o implícita
de la ecuación a la forma de expresión en la que la variable independiente x, la función y,
y su derivada y' aparecen en un solo miembro de la ecuación, siendo el otro miembro
igual a cero:
F (x, y, y ′) = 0

o en notación diferencial:

 dy 
F  x, y,  = 0
 dx 

donde F es una expresión analítica elemental.

Ejemplo: La siguiente EDO de primer orden de incógnita y (x ) está expresada en


forma general o implícita:

y′ + y + 3x2 = 0
42 Ecuaciones Diferenciales

Ejemplo: La siguiente EDO de primer orden de incógnita x (t ) está expresada en


forma general o implícita:
2
 dx  dx
x  + x2 =0
 dt  dt

Definición (Forma normal o explícita de una EDO de primer orden). Dada una ecuación
diferencial ordinaria de primer orden de incógnita y (x) , se llama forma normal o explícita
de la ecuación a la forma de expresión en la que la derivada aparece despejada en un
miembro de la ecuación:
y ′ = F (x, y )

o en notación diferencial:
dy
= F (x, y )
dx
donde F es una expresión analítica elemental.

Ejemplo: La siguiente EDO de primer orden de incógnita y (x ) está expresada en


forma normal o explícita:

y ′ = 2xy 2 + 2x 2 y + x + y − 2

Ejemplo: La siguiente EDO de primer orden de incógnita x (t ) está expresada en


forma normal o explícita:
dx
= cos 2 (t )
dt

Definición (Forma diferencial de una EDO de primer orden). Dada una ecuación
diferencial ordinaria de primer orden de incógnita y (x ) , se llama forma diferencial de la
dy
ecuación a la forma de expresión en la que la derivada y ′ = aparece en forma de
dx
suma de diferenciales:
G (x, y )dx + H ( x, y )dy = 0

donde G y H son expresiones analíticas elementales.

Ejemplo: La siguiente EDO de primer orden de incógnita y (x ) está expresada en


forma diferencial:

(y + 3x )dx + 2ydy = 0
2

Ejemplo: La siguiente EDO de primer orden de incógnita x (t ) está expresada en


forma diferencial:

 sen2 (x) 

 1− t 
( )
dt + x 2 + ln(t ) dx = 0
 
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 43

Observación. Una EDO de primer orden en forma diferencial siempre puede reescribirse
en forma normal, ya que:
dy −G (x, y )
G (x, y )dx + H (x, y )dy = 0 ⇔ = = F (x, y )
dx H (x, y )
Asimismo, en el caso de que la EDO de primer orden sea lineal respecto a la derivada y',
es decir, cuando y' aparezca solamente en forma de potencia de primer grado, también
es posible pasar fácilmente la ecuación a forma normal, tal como se muestra en los
siguientes ejemplos. Precisamente, las EDO de primer orden lineales en y' constituyen
los tipos más básicos a considerar en el estudio de las ecuaciones diferenciales de
primer orden.

Ejemplo: La EDO de primer orden lineal en y':

x ⋅ y ′ − xy = arcsen( x)

resulta en forma normal:


xy + arcsen(x )
y′ =
x

Ejemplo: La EDO de primer orden lineal en y' expresada en forma implícita:


cotg (x) y ′ + sen(x ) y − x = 0

resulta en forma normal:


x − sen(x ) y
y′ =
cotg (x)

2.1.2. Teorema básico de existencia y unicidad de soluciones

En el tema anterior se ha introducido el problema de valor inicial o problema de Cauchy


como uno de los tipos fundamentales de problemas relativos a la resolución de las
ecuaciones diferenciales ordinarias de orden arbitrario p. Este problema consiste en
encontrar una solución particular de la EDO considerando p condiciones iniciales o
valores que deben tomar la función solución y sus derivadas hasta el orden p−1 en un
punto a concreto. Además, a través de los ejemplos mostrados de integración directa ha
podido observarse que la función solución del problema existe y es única, al menos bajo
ciertas condiciones. Con respecto a esto, existen varios teoremas de existencia y
unicidad de la solución del problema de valor inicial para la EDO de primer orden, los
cuales especifican las condiciones que debe cumplir la EDO para que la función solución
del problema exista y sea única. Así, el teorema de Peano plantea una condición muy
general que garantiza la existencia de solución del problema de Cauchy, pero no su
unicidad, mientras que los teoremas de Picard y Cauchy contemplan condiciones
generales de existencia y unicidad. Sin embargo, el estudio de estos teoremas queda
fuera del ámbito de este texto, de modo que en su lugar se enuncia a continuación el
teorema básico de existencia y unicidad de la solución del problema de valor inicial para
la EDO de primer orden. Se trata de un teorema más restrictivo que los anteriores, pero
de gran utilidad práctica, ya que las condiciones que plantea son más fáciles de verificar
que las de Picard o Cauchy.
44 Ecuaciones Diferenciales

Teorema (Teorema básico de existencia y unicidad de la solución del problema de


valor inicial para la EDO de primer orden). Dada una ecuación diferencial ordinaria de
primer orden en forma normal y ′ = F (x, y ) , siendo F : D ⊆ o 2 → o una función derivable
∂F
respecto a y en D, si F y son funciones continuas en D, entonces para todo punto
∂y
(a, b) ∈ D existe una única función f : E(a, δ) → o solución de la EDO que verifica
y = f (a) = b , es decir, una única función solución del problema de valor inicial:

y ′ = F (x, y ) ; y (a) = b

Ejemplo: Sea la EDO de primer orden:

y2
y′ =
x −1
Aplicando el teorema básico de existencia y unicidad, es posible
determinar el dominio D del plano en cada uno de cuyos puntos (x, y ) ∈ D
existe y es única la solución de la EDO. En efecto, se considera la función
y2
F (x, y ) = , con dom(F ) = {(x, y ) ∈ o 2 | x ≠ 1} , así como la derivada
x −1
∂F (x, y ) 2y
parcial = . Se tiene que:
∂y x −1

• F es una función elemental, continua en todo su dominio


dom(F ) = {(x, y ) ∈ o 2 | x ≠ 1} .

∂F
• es una función elemental, continua en todo el dominio de
∂y
derivabilidad parcial dom(F y′ ) = {(x, y ) ∈ o 2 | x ≠ 1} .

Por tanto, el dominio D en que se cumple el teorema de existencia y


unicidad es D = dom(F ) ∩ dom(F y′ ) = {(x, y ) ∈ o 2 | x ≠ 1} .

Ejemplo: Sea la EDO de primer orden:


3 2
y′ = x + y2

Aplicando el teorema básico de existencia y unicidad, es posible


determinar el dominio D del plano en cada uno de cuyos puntos (x, y ) ∈ D
existe y es única la solución de la EDO. En efecto, se considera la función
3 2
F (x, y ) = x + y 2 , con dom(F ) = o 2 , así como la derivada parcial
∂F (x, y ) 2y
= . Ahora bien, dado que esta expresión no está
∂y
3⋅
3
(x 2
+y )
2 2

definida en el punto (x, y ) = (0,0) del dominio de F, es necesario


∂F
comprobar si en él también existe la derivada parcial mediante la
∂y
definición como límite:
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 45

3
∂F (0,0) F (0,0 + h) − F (0,0) h2 − 0
= lim = lim =
∂y h→ 0 h h→ 0 h

 1
3 2  lim+ 3
= 1+ = +∞
h 1 h 0 ∂F (0,0)
= lim = lim 3 = 0 = h→0
1 ⇒ ∃/
h→ 0 h h→0 h 1 ∂y
 lim = 1− = −∞
h→0 −
3
h 0

Se tiene que:
• F es una función elemental, continua en todo su dominio
dom(F ) = o 2 .

∂F
• es una función elemental, continua en todo el dominio de
∂y
derivabilidad parcial dom(F y′ ) = o 2 − {(0,0)} .

Por tanto, el dominio D en que se cumple el teorema de existencia y


unicidad es D = dom(F ) ∩ dom(F y′ ) = o 2 − {(0,0)} .

Ejemplo: Sea el problema de valor inicial:

y ′ = 1 − cos(xy ) ; y (0) = 0

Se trata de verificar si se cumplen las condiciones de existencia y


unicidad de la solución. En efecto, se considera la función
F (x, y ) = 1 − cos(xy ) , así como la derivada parcial
∂F (x, y ) xsen(xy )
= . Ahora bien, dado que esta expresión no está
∂y 2 1 − cos(xy )
definida en el punto (x, y ) = (0,0) especificado en la condición inicial, es
necesario comprobar si en él existe esta derivada parcial:

∂F (0,0) F (0,0 + h) − F (0,0) 1 − cos(0) − 0 0


= lim = lim = lim = 0
∂y h→ 0 h h→ 0 h h→ 0 h

 xsen( xy )
∂F (x, y )  si (x, y ) ≠ (0,0)
Por tanto, =  2 1 − cos(xy ) .
∂y 
 0 si (x, y ) = (0,0)

Continuidad de F en (0,0):

• lim F ( x, y ) = lim 1 − cos(xy ) = 1 − cos(0) = 0


x →0 x →0
y →0 y →0

• F (0,0) = 1 − cos(0) = 0

∴ F continua en (0,0)
46 Ecuaciones Diferenciales

∂F
Continuidad de en (0,0):
∂y

∂F (x, y ) xsen(xy ) x 1 − cos2 (xy )


• lim = lim = lim =
x →0 ∂y x →0 2 1 − cos(xy ) x →0 2 1 − cos(xy )
y →0 y →0 y →0

x 1 + cos(xy ) 1 − cos(xy ) x 1 + cos(xy ) 0⋅ 2


= lim = lim = 2 =0
x →0 2 1 − cos(xy ) x →0 2
y →0 y →0

∂F (0,0)
• =0
∂y

∂F
∴ continua en (0,0)
∂y

Por tanto, se cumple el teorema de existencia y unicidad, y la solución del


problema de valor inicial y ′ = 1 − cos(xy ) ; y (0) = 0 es única.

Observación. La interpretación geométrica del teorema de existencia y unicidad para la


∂F
EDO de primer orden y ′ = F (x, y ) es la siguiente. Si F y son funciones continuas en D,
∂y
entonces para todo punto (a, b) ∈ D existe una única curva integral que pasa por el punto
(a,b) del plano, y por tanto, es solución del problema de valor inicial y ′ = F (x, y ) ;
y (a) = b .

2.2. ECUACIONES DE VARIABLES SEPARADAS

Se trata de uno de los tipos más sencillos de EDO de primer orden, ya que tras una
sencilla manipulación algebraica se convierten en ecuaciones de integración directa.

Definición (Ecuación de variables separadas). Una ecuación diferencial ordinaria de


primer orden es una ecuación de variables separadas si y sólo si puede expresarse en la
forma diferencial:
G1 (x )G2 ( y )dx + H1 (x)H2 ( y )dy = 0

donde G1, G2, H1 y H2 son expresiones analíticas elementales.

Ejemplo: La siguiente EDO de primer orden de incógnita y (x ) es una ecuación de


variables separadas:
2y
3x 2 ydx + dy = 0
cos(x )

1
donde G1 (x ) = 3x 2 , G2 ( y ) = y , H1 (x ) = , H2 ( y ) = 2 y .
cos( x )
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 47

Ejemplo: La siguiente EDO de primer orden de incógnita x (t ) :

x x ′ = et + t + 2
puede ponerse en forma de ecuación de variables separadas:
dx
x = et + t + 2
dt

( )
xdx = et + t + 2 dt

(− e − t − 2)dt + xdx = 0
t

donde G1 (t ) = −e t − t − 2 , G2 (x) = 1 , H1 (t ) = 1 , H2 (x ) = x.

La siguiente proposición ofrece una caracterización alternativa de las ecuaciones de


variables separadas cuando están expresadas en forma normal sin necesidad de
pasarlas a forma diferencial.

Proposición (Caracterización de las ecuaciones de variables separadas). Una ecuación


diferencial ordinaria de primer orden es una ecuación de variables separadas si y sólo si
puede expresarse en la forma normal:
y ′ = G (x )H ( y )

donde G y H son expresiones analíticas elementales.

Demostración:

dy
En efecto, llamando y ′ = , resulta:
dx

dy 1 1
y ′ = G (x )H ( y ) ⇔ = G (x)H ( y ) ⇔ dy = G ( x)dx ⇔ −G (x )dx + dy = 0
dx H( y) H( y)

Ejemplo: La siguiente EDO de primer orden de incógnita y (x ) es una ecuación de


variables separadas:
y ′ = (x + 3)tg ( y )

donde F (x) = x + 3 , G ( y ) = tg ( y ) .

Procedimiento (Resolución de la ecuación de variables separadas). Para resolver una


ecuación de variables separadas G1 (x )G2 ( y )dx + H1 (x)H2 ( y )dy = 0 se procede de la
siguiente manera:
1. Se separan los diferenciales en la igualdad, poniendo dy a la izquierda y dx a la
derecha:
H1 (x)H2 ( y )dy = −G1 (x )G2 ( y )dx

2. Se divide la ecuación por G2 ( y )H1 (x ) :

H2 ( y ) G (x)
dy = − 1 dx
G2 ( y ) H1 (x)
48 Ecuaciones Diferenciales

3. Se integra en cada miembro respecto a su variable:


H2 ( y ) G (x)
∫G dy = ∫ − 1 dx
2 ( y ) H1 (x)

Ejemplo: Integrar la EDO de primer orden de variables separadas:


1
dx + 5xdy = 0
y2

Siguiendo el procedimiento descrito, se tiene:


−1
5xdy = dx
y2

−1
y 2dy = dx
5x

2 −1
∫y dy = ∫ 5x dx
1
3
y3 = −1 ln
5
(x )+ C
y 3 = −53 ln( x ) + C

y= 3 −3 ln
5
(x )+ C
Ejemplo: Integrar la EDO de primer orden de variables separadas:

(
− 2xdx + 1 + x 4 e y dy = 0 )
Siguiendo el procedimiento descrito, se tiene:

(1 + x )e dy = 2xdx
4 y

2x
e ydy = dx
1 + x4

y 2x
∫e dy = ∫ 1 + x4 dx
( )
e y = arctg x2 + C

y = ln(arctg(x ) + C ) 2

Ejemplo: Resolver el problema de valor inicial:

( )
− e2t dt + 1 + e2t xdx = 0 ; x (0) = 1

Se trata de una ecuación de variables separadas, y su solución general es:

(1 + e )xdx = e
2t 2t
dt

e2t
xdx = dt
1 + e2t
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 49

e2t
∫ xdx = ∫ 1 + e2t dt
1 x2
2
= 2
(
1 ln 1 + e2t )+ C
(
x 2 = ln 1 + e 2t + 2C )
Sustituyendo ahora la condición x(0) = 1 en la solución general, se
obtiene una ecuación que permite determinar el valor de la constante de
integración C:
x( 0 ) = 1

(1)2 = ln(1 + 1) + 2C

1 = ln(2) + 2C

C= 1
2
− 12 ln(2)

Por tanto, la solución del problema de valor inicial es:

( )
x2 = ln 1 + e2t + 1 − ln(2)

Ejemplo: Problema de mezclas. Un depósito contiene inicialmente 5 m3 de una


solución compuesta por un 80% de agua y un 20% de alcohol. En el
depósito se vierte otra solución que contiene un 60% de agua y un 40%
de alcohol a razón de 0'001 m3 por segundo, y al mismo tiempo del
depósito se evacua líquido a igual velocidad. Suponiendo que la solución
en el depósito es siempre homogénea (se mezcla constantemente), hallar
la función que determina la cantidad de alcohol presente en el depósito
en cada instante, así como la cantidad de alcohol en el depósito
transcurridos 60 segundos.

Figura 2.1: Problema de mezclas.


Sea q = f (t ) la cantidad de alcohol en la solución del depósito en un
instante t. Dado que la velocidad de entrada de líquido (con un porcentaje
de alcohol del 40%) es de 0'001 m3/s, y ésta es igual a la velocidad de
salida del líquido mezcla (con un porcentaje de alcohol del
100 q = 20q %), puede ponerse que la variación de alcohol en el depósito
5
40 − (0'001) 20 q = 0'0004 − 0'0002q m3/s. Además, por
es v = (0'001) 100 100
dq
definición de derivada se cumple que v = , de donde se obtiene la EDO
dt
50 Ecuaciones Diferenciales

de primer orden de variables separadas, dada en forma normal, que


determina la cantidad de alcohol en un instante t:
dq
= 0'0004 − 0'0002q
dt
Además, como la cantidad inicial de alcohol es conocida, q(0) = 1 m3, el
problema a resolver es un problema de valor inicial:
dq
= 0'0004 − 0'0002q ; q(0) = 1
dt
Poniendo las variables y los diferenciales por separado en cada miembro
de la igualdad, con dq a la izquierda y dt a la derecha, se puede calcular
la solución general por integración directa:
1
dq = dt
0'0004 − 0'0002q

1
∫ 0'0004 − 0'0002q dq = ∫ dt
−1 ln
0'0002
( 0'0004 − 0'0002q ) = t + C
ln( 0'0004 − 0'0002q ) = (−0'0002)t − 0'0002C

0'0004 − 0'0002q = e ( −0'0002)t − 0'0002C

0'0004 − 0'0002q = ±e −0'0002C e ( −0'0002)t

0'0004 − 0'0002q = Ce ( −0'0002)t

q = 2 − 5000Ce ( −0'0002)t

Sustituyendo ahora la condición q(0) = 1 en la solución general, se


obtiene una ecuación que permite determinar el valor de la constante de
integración C:
q(0) = 1

2 − 5000C = 1

C= 1
5000

Por tanto, la solución del problema de valor inicial es:

q = 2 − e ( −0'0002)t

Finalmente, la cantidad de alcohol en el depósito transcurridos 60


segundos es:

q(60) = 2 − e ( −0'0002)(60) = 2 − e −0'012 ≅ 1'011928 m3


2. Ecuaciones Diferenciales Ordinarias de Primer Orden 51

2.3. ECUACIONES HOMOGÉNEAS

2.3.1. Definición de ecuación homogénea

Se trata de un tipo de EDO de primer orden que se puede reducir a una ecuación de
variables separadas mediante un cambio de variable específico. Previamente a su
definición es necesario recordar el concepto de función homogénea de varias variables.

Definición (Función homogénea). Una función escalar de varias variables f : D ⊆ o n → o


se llama función homogénea de grado α ∈ o si y sólo si para cualquier x = (x1,..., xn ) ∈ D
y cualquier t ∈ o tal que tx ∈ D , se verifica que:

f (tx ) = f (tx1,..., txn ) = t αf (x1,..., xn ) = t αf (x )

Ejemplo: La función de dos variables:

f (x, y ) = x2 y + xy 2

es una función homogénea de grado α = 3 , ya que:

( )
f (tx, ty ) = (tx)2 (ty ) + (tx)(ty )2 = t 3x2 y + t 3xy 2 = t 3 x2 y + xy 2 = t 3f (x, y )

Ejemplo: La función de tres variables:

xy + xz + yz
f (x, y, z ) =
x

es una función homogénea de grado α = 12 , ya que:

(tx)(ty ) + (tx )(tz ) + (ty )(tz ) t 2xy + t 2xz + t 2 yz


f (tx, ty, tz ) = = =
tx tx

t 2 (xy + xz + yz ) xy + xz + yz 1
= = t = t 2 f (x, y, z )
tx x

Ejemplo: La función de dos variables:


x
f (x, y ) =
y

es una función homogénea de grado α = 0 , ya que:


tx x
f (tx, ty ) = = = f (x, y ) = t 0f (x, y )
ty y

Ejemplo: La función de dos variables:

f (x, y ) = x2 + xy − y

no es una función homogénea, ya que para todo α ∈ o :

f (tx, ty ) = (tx )2 + (tx )(ty ) − ty = t 2x2 + t 2xy − ty


52 Ecuaciones Diferenciales

t αf (x, y ) = t α x2 + xy − y = t 2α x2 + t 2α xy − t 2α y

luego f (tx, ty ) ≠ t αf (x, y ) .

Definición (Ecuación diferencial homogénea). Una ecuación diferencial ordinaria de


primer orden en forma diferencial G (x, y )dx + H (x, y )dy = 0 es una ecuación diferencial
homogénea si y sólo si G(x,y) y H(x,y) son funciones homogéneas de igual grado α ∈ o :

G (tx, ty ) = t α G (x, y ) ; H (tx, ty ) = t α H (x, y )

Ejemplo: La siguiente EDO de primer orden de incógnita y (x ) es una ecuación


diferencial homogénea:

(x 2
)
+ xy dx + 4 y 2dy = 0

ya que G (x, y ) = x2 + xy y H (x, y ) = 4 y 2 son funciones homogéneas de


grado α = 2 :

( )
G (tx, ty ) = (tx ) 2 + (tx )(ty ) = t 2 x 2 + t 2 xy = t 2 x 2 + xy = t 2G (x, y )

H (tx, ty ) = 4(ty ) 2 = 4t 2 y 2 = t 2 (4 y ) = t H(x, y)


2 2

Ejemplo: La siguiente EDO de primer orden de incógnita y (x ) es una ecuación


diferencial homogénea:
3x − y dx + 3 x dy = 0

ya que G (x, y ) = 3 x − y y H (x, y ) = 3 x son funciones homogéneas de grado


α = 13 :

1
G (tx, ty ) = 3 tx − ty = 3 t ( x − y ) = 3 t 3 x − y = t 3 G (x, y )
1
H (tx, ty ) = 3 tx = 3 t 3 x = t 3 H (x, y )

2.3.2. Caracterización de las ecuaciones homogéneas

La siguiente proposición ofrece una caracterización alternativa de las ecuaciones


diferenciales homogéneas cuando están expresadas en forma normal sin necesidad de
pasarlas a forma diferencial.

Proposición (Caracterización de las ecuaciones diferenciales homogéneas). Una


ecuación diferencial ordinaria de primer orden en forma normal y ′ = F (x, y ) es una
ecuación diferencial homogénea si y sólo si F es una función homogénea de grado cero, es
y
decir, si y sólo si F puede expresarse en la forma F   :
x

y
F (tx, ty ) = F (x, y ) = F  
x
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 53

Demostración:

Considérese en primer lugar que la EDO y ′ = F (x, y ) es homogénea, es decir, puede expresarse en forma
diferencial:

G (x, y )dx + H (x, y )dy = 0

donde G(x,y) y H(x,y) son funciones homogéneas de igual grado α ∈ o . Si se pasa esta ecuación a forma
normal, se tiene que:

−G (x, y )
y′ = = F ( x, y )
H (x, y )

Por tanto,

− G (tx, ty ) − t αG (x, y ) − G (x, y )


F (tx, ty ) = = = = F (x, y )
H (tx, ty ) t αH (x, y ) H (x, y )

1
luego F es una función homogénea de grado cero. Además, haciendo t = en la igualdad anterior, se tiene:
x

x y  y y
F (x, y ) = F (tx, ty ) = F  ,  = F 1,  = F  
x x  x x

y
Recíprocamente, si la EDO y′ = F (x, y ) verifica que F puede expresarse en la forma F   , es decir, F es una
x
función homogénea de grado cero, entonces pasando la ecuación a forma diferencial, resulta:

dy
y ′ = F (x, y ) ⇔ = F (x, y ) ⇔ dy = F (x, y )dx ⇔ −F (x, y )dx + dy = 0
dx

donde G (x, y ) = −F (x, y ) y H (x, y ) = 1 son funciones homogéneas de igual grado α = 0 , luego la EDO y′ = F (x, y )
es homogénea.

Ejemplo: La siguiente EDO de primer orden de incógnita y (x ) es una ecuación


diferencial homogénea:
3xy
y′ =
x − y2
2

3xy
ya que F (x, y ) = 2 es una función homogénea de grado cero:
x − y2

3(tx )(ty ) 3t 2 xy 3t 2 xy 3xy


F (tx, ty ) = = = = = F (x, y )
2
(tx ) − (ty ) 2 2 2
t x −t y 2 2 2
(
t x −y 2 2
) x − y2
2

y
Además, F puede expresarse en la forma F   :
x

3y 3xy y
3 
3xy x x 2 x  = F y 
F (x, y ) =
2 2
=
2 2
=
2
=
2
 
x −y x −y y  y  x
1− 1−  
2 2
x x x
54 Ecuaciones Diferenciales

Ejemplo: La siguiente EDO de primer orden de incógnita x (t ) es una ecuación


diferencial homogénea:

dx t2
=
dt tx

t2
ya que F (t, x ) = es una función homogénea de grado cero:
tx

(rt ) 2 r 2t 2 r2 t2 t2
F (rt, rx) = = = = = F (t, x )
(rt )(rx) r 2tx r 2 tx tx

x
Además, H puede expresarse en la forma F   :
t 

t2
t2 t2 1 x
F (t, x ) = = = = F 
tx tx x t 
t2 t

2.3.3. Resolución de ecuaciones homogéneas

Proposición (Cambio de variable en la ecuación diferencial homogénea). Dada una


ecuación diferencial homogénea G (x, y )dx + H ( x, y )dy = 0 , con el cambio de variable
y = tx , se transforma en una ecuación de variables separadas de la forma
G1 (t )G2 (x )dt + H1 (t )H2 (x)dx = 0 :

 y = tx 
G (x, y )dx + H (x, y )dy = 0 ⇒   ⇒ G1 (t )G2 (x)dt + H1 (t )H2 (x)dx = 0
dy = xdt + tdx 

Demostración:

En efecto, sea la ecuación homogénea:

G (x, y )dx + H (x, y )dy = 0

 y = tx 
donde G(x,y) y H(x,y) son funciones homogéneas de igual grado α ∈ o . Si se aplica el cambio  ,
dy = xdt + tdx 
resulta:

G (x, tx )dx + H (x, tx )(xdt + tdx ) = 0

x αG (1, t )dx + x αH (1, t )(xdt + tdx) = 0

H (t )xdt + (G (t ) + tH (t ))dx = 0

que es una ecuación de variables separadas de la forma G1(t )G2 (x )dt + H1(t )H2 (x )dx = 0 con G1(t ) = H (t ) ,
G2 (x ) = x , H1(t ) = G (t ) + tH (t ) , H2 (x ) = 1 .
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 55

Procedimiento (Resolución de la ecuación diferencial homogénea). Para resolver una


ecuación diferencial homogénea G (x, y )dx + H ( x, y )dy = 0 se procede de la siguiente
manera:
1. Se aplica el cambio de variable y = tx , obteniéndose una ecuación de variables
separadas:
 y = tx 
G (x, y )dx + H (x, y )dy = 0 ⇒   ⇒ G1 (t )G2 (x)dt + H1 (t )H2 (x)dx = 0
dy = xdt + tdx 
2. Se separan las variables y los diferenciales en la igualdad, poniendo dx a la
izquierda y dt a la derecha:
H2 (x ) G (t )
dx = − 1 dt
G2 (x) H1 (t )

3. Se integra en cada miembro respecto a su variable:


H2 (x ) G (t )
∫G dx = ∫ − 1 dt
2 (x) H1 (t )

y
4. Se deshace el cambio en la solución sustituyendo t = .
x

Ejemplo: Integrar la EDO homogénea:

(x 2
)
− 3 y 2 dx + 2xydy = 0

 y = tx 
En primer lugar, se aplica el cambio   , obteniéndose:
dy = xdt + tdx 

(x 2
)
− 3t 2x2 dx + 2tx2 (xdt + tdx ) = 0

(
2tx 3dt + 1 − t 2 x2dx = 0 )
A continuación, se integra la ecuación de variables separadas obtenida:

(1 − t )x dx = −2tx dt
2 2 3

1 −2t
dx = dt
x 1 − t2
1 −2t
∫ x dx = ∫ 1 − t 2 dt

ln x = ln 1 − t 2  + C
( )
 

x = eC 1 − t 2

(
x = ±eC 1 − t 2 )
(
x = C 1 − t2 )
Finalmente, se deshace el cambio de variable:
56 Ecuaciones Diferenciales

 y2 
x = C1 − 2 
 x 

Ejemplo: Integrar la EDO homogénea:


dy y y
= + cos2  
dx x x

 y = tx 
En primer lugar, se aplica el cambio   , obteniéndose:
dy = xdt + tdx 

y  y 
dy =  + cos2   dx
x  x 

 tx  tx  
xdt + tdx =  + cos2   dx
x  x 

(
xdt + tdx = t + cos2 (t ) dx )
cos2 (t )dx = xdt
A continuación, se integra la ecuación de variables separadas obtenida:
1 1
dx = dt
x cos2 (t )

1 1
∫ x dx = ∫ cos2 (t ) dt
ln( x ) = tg(t ) + C

x = etg(t ) +C

x = ±eCetg(t )

x = Cetg(t )
Finalmente, se deshace el cambio de variable:
y
tg  
x= Ce  x 

Ejemplo: Resolver el problema de valor inicial:


(x + y )dx − xdy = 0 ; y (1) = 1

Se trata de una ecuación diferencial homogénea. Para determinar su


 y = tx 
solución general se aplica el cambio   , obteniéndose:
dy = xdt + tdx 
(x + tx )dx − x (xdt + tdx ) = 0

xdx = x2dt
A continuación, se integra la ecuación de variables separadas obtenida:
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 57

xdx = x2dt
1
dx = dt
x
1
∫ x dx = ∫ dt
( )
ln x = t + C

Se deshace el cambio de variable:


y
ln( x ) = +C
x
y = x ln( x ) − Cx

Sustituyendo ahora la condición y (1) = 1 en la solución general, se obtiene


una ecuación que permite determinar el valor de la constante de
integración C:
y (1) = 1

1 ⋅ ln(1 ) − C = 1

C = −1
Por tanto, la solución del problema de valor inicial es:
y = x ln( x ) + x

2.3.4. Ecuaciones reducibles a homogéneas

A continuación se estudian ciertas ecuaciones diferenciales de primer orden que no son


homogéneas, pero que pueden reducirse a éstas mediante cambios de variable
específicos.

Definición (Ecuación diferencial reducible a homogénea). Una ecuación diferencial


ordinaria de primer orden es una ecuación reducible a homogénea si y sólo si puede
expresarse en la forma normal:

 ax + by + c 
y ′ = F  
 a′x + b′y + c′ 
donde F es una expresión analítica elemental.

Ejemplo: La siguiente EDO de primer orden de incógnita y (x ) es una ecuación


diferencial reducible a homogénea:
2x + 3 y − 1
y′ =
3x − y + 4
58 Ecuaciones Diferenciales

Ejemplo: La siguiente EDO de primer orden de incógnita x (t ) es una ecuación


diferencial reducible a homogénea:
dx t + x + 2
= ln 
dt t − x + 2

 ax + by + c 
Observación. Nótese que la ecuación y ′ = F   no es homogénea, porque la
 a′x + b′y + c′ 
función F no es homogénea de grado cero, debido a la presencia de los términos c y c'
(siempre que al menos uno de ellos no sea nulo). Además, en una ecuación reducible a
homogénea las expresiones ax + by + c , a′x + b′y + c′ son expresiones afines de dos
variables x e y que corresponden a sendas rectas del plano OXY. Precisamente, la
dependencia o independencia lineal entre los vectores directores de estas dos rectas
permite definir los cambios de variable apropiados para su resolución.

Proposición (Cambio de variable en la ecuación diferencial reducible a homogénea).


Dada una ecuación diferencial reducible a homogénea de la forma:

 ax + by + c 
y ′ = F  
 a′x + b′y + c′ 
se verifica que:
• Si los vectores directores de las rectas ax + by + c = 0 , a ′x + b ′y + c ′ = 0 son
libres, es decir, si ab′ ≠ a′b , las rectas se cortan en un punto (x0,y0), por lo que
 u = x − x0
mediante el cambio de variable  la EDO se transforma en una
v = y − y0
ecuación homogénea:
dv  au + bv 
= F 
du  a ′u + b ′v 
• Si los vectores directores de las rectas ax + by + c = 0 , a ′x + b ′y + c ′ = 0 son
ligados, es decir, si ab′ = a′b , las rectas son paralelas, por lo que mediante el
cambio de variable y = −ba x + u la EDO se transforma en una ecuación de
variables separadas:
du a  bu + c 
= b
+ F 
dx  b ′u + c ′ 

Demostración:

En efecto, sea la ecuación reducible a homogénea:

 ax + by + c 
y ′ = F  
 a′x + b′y + c′ 

y considérese en primer lugar el caso ab′ ≠ a′b , tal que las rectas ax + by + c = 0 , a ′x + b ′y + c ′ = 0 se cortan en
 u = x − x0
un punto (x0,y0). Si se aplica el cambio  , se produce una traslación del origen de coordenadas al
v = y − y0
punto (x0,y0) que elimina los términos c y c' y transforma la EDO en una ecuación homogénea:
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 59

dy  ax + by + c  u = x − x0 ; du = dx  dv  a(u + x0 ) + b(v + y0 ) + c 
= F   ⇒  ⇒ = F  ⇒

dx  a′x + b′y + c′  v = y − y 0 ; dv = dy  du  a′(u + x0 ) + b′(v + y0 ) + c′ 

dv  au + bv + ax0 + ay0 + c  dv  au + bv 
⇒ = F  ⇒
 = F  
du  a′u + b′v + a′x0 + b′y0 + c′  du  a′u + b′v 

ya que, por ser (x0,y0) punto de corte de ambas rectas ax + by + c = 0 , a ′x + b ′y + c ′ = 0 , se cumple que
ax0 + by0 + c = 0 , a′x0 + b′y0 + c′ = 0 .

Por otra parte, si ab′ = a′b , donde k = a′ = b' , se aplica el cambio y = −a x + u , con lo que se obtiene:
a b b

 ax + by + c   ax + by + c   y = −a x + u 
dy dy  b 
= F   ⇒ = F   ⇒  ⇒
′ ′ ′ −a
dx  a x + b y + c  dx  kax + kby + c′  dy = b dx + du

⇒ b
−a dx + du
= F 
( b
)
 ax + b −a x + u + c 
 ⇒ −a + du = F  bu + c  ⇒ du = a + F  bu + c 
   
dx
 b
( )
 kax + kb −a x + u + c′ 

b dx  kbu + c′  dx b  b′u + c′ 

que es una EDO de variables separadas.

Observación. En la práctica, para resolver las ecuaciones reducibles a homogéneas


 ax + by + c 
y ′ = F   del primer tipo, es decir, tales que las rectas ax + by + c = 0 ,
 a' x + b′y + c′ 
a ′x + b ′y + c ′ = 0 se cortan en un punto (x0,y0), se emplea el cambio de variable mixto
 u = x − x0
y − y 0 = t (x − x 0 ) que combina los cambios  (transformación a homogénea), y
v = y − y0
v = tu (transformación a variables separadas), de manera que la ecuación se transforma
directamente de una EDO reducible a homogénea a una EDO de variables separadas.

Procedimiento (Resolución de la ecuación diferencial reducible a homogénea tipo I).


 ax + by + c 
Para resolver una ecuación diferencial reducible a homogénea y ′ = F   , tal
 a′x + b′y + c′ 
que ab′ ≠ a′b (las rectas ax + by + c = 0 , a ′x + b ′y + c ′ = 0 se cortan en un punto (x0,y0)),
se procede de la siguiente manera:
1. Se aplica el cambio de variable y − y 0 = t (x − x0 ) , obteniéndose una ecuación de
variables separadas:

dy  ax + by + c   y − y 0 = t (x − x 0 )  (x − x 0 )dt + tdx  a + bt 
= F   ⇒  ⇒ = F 
dx  a ′x + b ′y + c ′  dy = ( x − x 0 )dt + tdx  dx  a ′ + b ′t 

2. Se separan las variables y los diferenciales en la igualdad, poniendo dx a la


izquierda y dt a la derecha:
1 1
dx = dt
x − x0  a + bt 
F −t
 a ′ + b ′t 
3. Se integra en cada miembro respecto a su variable:
1 1
∫x−x dx = ∫  a + bt 
dt
0
F −t
 a ′ + b ′t 
60 Ecuaciones Diferenciales

y − y0
4. Se deshace el cambio en la solución sustituyendo t = .
x − x0

Ejemplo: Integrar la EDO reducible a homogénea:


x − y +1
y′ =
x + y +1

Se trata de una EDO reducible a homogénea del tipo I, ya que


ab ′ = 1 ≠ a ′b = −1 . El punto de corte (x0,y0) entre las dos rectas
x − y + 1 = 0 , x + y + 1 = 0 es:

x − y + 1 = 0
 ⇒ 2x + 2 = 0 ⇒ x = −1, y = 0
x + y + 1 = 0
 y = t (x + 1) 
Se aplica el cambio   , obteniéndose:
dy = ( x + 1)dt + tdx 
x − y +1
dy = dx
x + y +1

x − t (x + 1) + 1
(x + 1)dt + tdx = dx
x + t (x + 1) + 1

(x + 1) − t (x + 1)
(x + 1)dt + tdx = dx
(x + 1) + t (x + 1)

1− t
(x + 1)dt + tdx = dx
1+ t

1 − 2t − t 2
(x + 1)dt = dx
1+ t
A continuación, se integra la ecuación de variables separadas obtenida:
1 1+ t
dx = dt
x +1 1 − 2t − t 2
1 −2 − 2t
∫ x + 1 dx = −1
2 ∫ 1 − 2t − t 2 dt

ln x + 1 = −21 ln 1 − 2t − t 2  + C
( )  

eC
x +1 =
1 − 2t − t 2

± eC
x +1 =
1 − 2t − t 2

C
x +1 =
1 − 2t − t 2
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 61

Finalmente, se deshace el cambio de variable:


C
x +1 =
2y y2
1− −
x + 1 (x + 1) 2

Procedimiento (Resolución de la ecuación diferencial reducible a homogénea tipo II).


 ax + by + c 
Para resolver una ecuación diferencial reducible a homogénea y′ = F   , tal
 a′x + b′y + c′ 
que ab′ = a′b (las rectas ax + by + c = 0 , a ′x + b ′y + c ′ = 0 son paralelas), se procede de la
siguiente manera:

1. Se aplica el cambio de variable y = −ba x + u , obteniéndose una ecuación de


variables separadas:

dy  ax + by + c   y = −a x + u  du  bu + c 
= F   ⇒ b ⇒ a + F
 − a dx + du
= b

dx  a ′x + b ′y + c ′  dy = b  dx  b ′u + c ′ 

2. Se separan las variables y los diferenciales en la igualdad, poniendo du a la


izquierda y dx a la derecha:
1
du = dx
a + F  bu + c 
b
 
 b ′u + c ′ 
3. Se integra en cada miembro respecto a su variable:
1
∫ du = ∫ dx
a + F  bu + c 
b
 
 b ′u + c ′ 
a
4. Se deshace el cambio en la solución sustituyendo u = y + x.
b

Ejemplo: Integrar la EDO reducible a homogénea:


x + y +1
y′ =
x + y −1

Se trata de una EDO reducible a homogénea del tipo II, ya que


 y = −x + u 
ab ′ = a ′b = 1 . Se aplica el cambio   , obteniéndose:
dy = −dx + du
x + y +1
dy = dx
x + y −1

x − x + u +1
− dx + du = dx
x − x + u −1
u +1
− dx + du = dx
u −1
62 Ecuaciones Diferenciales

2u
du = dx
u −1
A continuación, se integra la ecuación de variables separadas obtenida:
u −1
du = 2dx
u

 1
∫ 1 − u du = ∫ 2dx
( )
u − ln u = 2x + C

Finalmente, se deshace el cambio de variable:


(
x + y − ln x + y = 2x + C)
( )
ln x + y = y − x − C

x + y = e −C e y − x

x + y = ±e −C e y − x

x + y = Ce y − x

2.4. ECUACIONES DIFERENCIALES EXACTAS

2.4.1. Definición de ecuación diferencial exacta

Las ecuaciones diferenciales exactas son EDO de primer orden cuya expresión en forma
diferencial coincide con la diferencial (total) de una función de dos variables. Su
resolución puede llevarse a cabo mediante un proceso de integración parcial respecto a
la primera variable, seguido de una derivación parcial respecto a la segunda variable, y
de una nueva integración simple respecto a la segunda variable.

Definición (Ecuación diferencial exacta). Una ecuación diferencial ordinaria de primer


orden en forma diferencial G (x, y )dx + H (x, y )dy = 0 es una ecuación diferencial exacta en
un conjunto D ⊆ o 2 si y sólo si existe una función F : D ⊆ o 2 → o , tal que para todo
(x, y ) ∈ D verifica:

dF ( x, y ) = G (x, y )dx + H (x, y )dy

es decir,
∂F (x, y ) ∂F (x, y )
= G ( x, y ) ; = H (x, y )
∂x ∂y

Ejemplo: La siguiente EDO de primer orden de incógnita y (x ) es una ecuación


diferencial exacta:

(3x 2
) ( )
+ 4xy 2 dx + 4x 2 y − 3 y 2 dy = 0
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 63

ya que considerando la función F (x, y ) = x 3 + 2x 2 y 2 − y 3 se cumple que:

∂F (x, y ) ∂F (x, y )
= 3x 2 + 4xy 2 ; = 4x 2 y − 3 y 2
∂x ∂y

es decir,

( ) (
dF ( x, y ) = 3x 2 + 4xy 2 dx + 4x 2 y − 3 y 2 dy )
Ejemplo: La siguiente EDO de primer orden de incógnita x (t ) es una ecuación
diferencial exacta:

 2t   8x 3 
 dt +  dx = 0
 1 + t + 2x 4
2
  1 + t 2 + 2x 4 
 

(
ya que considerando la función F (t, x ) = ln 1 + t 2 + 2x 4 se cumple que: )
∂F (t, x ) 2t ∂F (t, x) 8x3
= ; =
∂t 1 + t 2 + 2x 4 ∂x 1 + t 2 + 2x 4
es decir,

 2t   8x 3 
dF (t, x) =  dt +  dx
 1 + t + 2x 4
2
  1 + t 2 + 2x 4 
 

2.4.2. Caracterización de las ecuaciones diferenciales exactas

La siguiente proposición proporciona una condición necesaria y suficiente que


caracteriza a las ecuaciones diferenciales exactas, y que además permite determinar de
forma sencilla si una EDO de primer orden es exacta o no.

Proposición (Caracterización de las ecuaciones diferenciales exactas). Sea una


ecuación diferencial ordinaria de primer orden en forma diferencial
G (x, y )dx + H ( x, y )dy = 0 , tal que las funciones G y H poseen derivadas parciales
continuas en un conjunto D ⊆ o 2 . La ecuación G (x, y )dx + H ( x, y )dy = 0 es una ecuación
diferencial exacta en D ⊆ o 2 si y sólo si para todo (x, y ) ∈ D se verifica:

∂G (x, y ) ∂H (x, y )
=
∂y ∂x

Demostración:

1. Si la ecuación G (x, y )dx + H(x, y )dy = 0 es una ecuación diferencial exacta en D ⊆ o2 , entonces para todo
∂G (x, y ) ∂H (x, y )
(x, y ) ∈ D se verifica = .
∂y ∂x

En efecto, por definición de ecuación diferencial exacta existe una función F : D ⊆ o2 → o , tal que para todo
(x, y ) ∈ D verifica:

∂F ( x, y ) ∂F (x, y )
= G ( x, y ) ; = H (x, y )
∂x ∂y
64 Ecuaciones Diferenciales

y de aquí:

∂ 2F (x, y ) ∂G (x, y ) ∂ 2F (x, y ) ∂H (x, y )


= ; =
∂x∂y ∂y ∂y∂x ∂x

Ahora bien, como las derivadas parciales de G y H son continuas en (x, y ) ∈ D y coinciden con las derivadas
parciales segundas de F, podemos aplicar el teorema de Schwarz a la función F, de donde resulta para todo
(x, y ) ∈ D :

∂ 2F (x, y ) ∂ 2F (x, y )
=
∂x∂y ∂y∂x

y en consecuencia:

∂G (x, y ) ∂H (x, y )
=
∂y ∂x

∂G (x, y ) ∂H (x, y )
2. Si para todo (x, y ) ∈ D se verifica = , entonces la ecuación G (x, y )dx + H (x, y )dy = 0 es una
∂y ∂x

ecuación diferencial exacta en D ⊆ o2 .

Se trata de probar que existe una función F : D ⊆ o2 → o , tal que para todo (x, y ) ∈ D verifica
dF (x, y ) = G (x, y )dx + H (x, y )dy . Dado que las funciones G y H son continuas en (x, y ) ∈ D , ambas son
parcialmente continuas en (x, y ) ∈ D respecto a cada variable, luego existe la función F : D ⊆ o2 → o integral
parcial de G respecto a la variable x:

F (x, y ) = ∫ G (x, y )dx + ϕ( y )

siendo ϕ una función arbitraria que depende sólo de la variable y, equivalente a la constante de integración de
las integrales simples, pues su derivada respecto a x es nula. Derivando respecto a y la función F obtenida,
resulta:

∂F (x, y )
∂y
=

∂y
(∫ G(x, y )dx ) + ddyϕ
Por tanto, para que la función F verifique dF (x, y ) = G (x, y )dx + H (x, y )dy y la EDO sea exacta, debe cumplirse
para (x, y ) ∈ D que:


∂y
(∫ G(x, y)dx) + ddyϕ = H(x, y)
es decir,


dy
= H ( x, y ) −

∂y
(∫ G(x, y)dx )

y como ϕ es una función arbitraria que depende sólo de la variable y, su derivada también deberá ser
dy

independiente de x. En consecuencia, se trata de verificar que H (x, y ) −



∂y
(∫ G(x, y )dx) es una función arbitraria
que no depende de x, es decir, una función cuya derivada respecto a x es nula. Se tiene que:

2
∂ 
 H (x, y ) −
∂x 

∂y
(∫ G(x, y)dx) = ∂H∂(xx, y ) − ∂∂y∂x (∫ G(x, y )dx) =

∂G (x, y ) ∂H (x, y )
de donde, por el teorema de Schwarz y la hipótesis = , se obtiene:
∂y ∂x

∂2
=
∂H (x, y )
∂x

∂x∂y
(∫ G(x, y)dx) = ∂H∂(xx, y ) − ∂∂y  ∂∂x (∫ G(x, y )dx ) = ∂H∂(xx, y ) − ∂G∂(xy, y ) = 0
 
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 65

En conclusión, la función F verifica dF (x, y ) = G (x, y )dx + H (x, y )dy y la EDO es exacta.

Ejemplo: La siguiente EDO de primer orden de incógnita y (x ) es una ecuación


diferencial exacta:
y cos( x)dx + (sen(x ) + 2y )dy = 0

ya que, llamando G (x, y ) = y cos( x ) , H (x, y ) = sen(x ) + 2y , se cumple que:

∂G (x, y ) ∂H ( x, y )
= cos(x ) ; = cos(x)
∂y ∂x

Ejemplo: La siguiente EDO de primer orden de incógnita x (t ) es una ecuación


diferencial exacta:

(2x 2
)
− t dt + 4txdx = 0

ya que, llamando G (t, x) = 2x 2 , H (t, x) = 4tx , se cumple que:

∂G (t, x ) ∂H (t, x)
= 4x ; = 4x
∂x ∂t

2.4.3. Resolución de ecuaciones diferenciales exactas

La solución de la ecuación diferencial exacta G (x, y )dx + H ( x, y )dy = 0 es la familia


uniparamétrica de funciones F (x, y ) = C , donde F : D ⊆ o 2 → o es una función que para
todo (x, y ) ∈ D verifica dF ( x, y ) = G (x, y )dx + H (x, y )dy , lo que puede comprobarse de
forma inmediata sin más que diferenciar la igualdad F (x, y ) = C . El procedimiento que se
describe a continuación permite obtener la solución de la EDO exacta mediante un
proceso de tres etapas integración-derivación-integración. Precisamente, este
procedimiento se basa en la demostración de la condición de suficiencia del teorema
anterior de caracterización de las EDO exactas.

Procedimiento (Resolución de la ecuación diferencial exacta). Para resolver una


ecuación diferencial exacta G (x, y )dx + H (x, y )dy = 0 se debe encontrar una función
∂F (x, y ) ∂F (x, y )
F : D ⊆ o 2 → o que verifique = G (x, y ) , = H (x, y ) para todo (x, y ) ∈ D ,
∂x ∂y
para lo cual se procede de la siguiente forma:
∂F (x, y )
1. Se integra = G (x, y ) parcialmente respecto a x:
∂x

F (x, y ) = ∫ G (x, y )dx + ϕ( y )

donde ϕ es una función desconocida dependiente de la variable y.


2. Para determinar ϕ, la expresión F(x,y) obtenida se deriva parcialmente respecto a
∂F (x, y )
y, y se iguala a = H (x, y ) :
∂y

∂F (x, y )
∂y
=

∂y
(∫ G(x, y)dx ) + ddyϕ = H(x, y)
66 Ecuaciones Diferenciales


3. Se despeja en la ecuación formada y se integra respecto a y. La función ϕ
dy
calculada se lleva a la expresión de F(x,y), y con ello se obtiene la solución
general de la EDO exacta.

Ejemplo: Integrar la ecuación diferencial exacta:

(x 2
) (
+ 2xy dx + x 2 + 2y 2 dy = 0 )
Se trata efectivamente de una EDO exacta, ya que llamando
G (x, y ) = x 2 + 2xy , H (x, y ) = x 2 + 2y 2 , se cumple que:

∂G (x, y ) ∂H (x, y )
= 2x ; = 2x
∂y ∂x

Primero, se integra G (x, y ) = x 2 + 2xy parcialmente respecto a x:

∫ (x )
2
F (x, y ) = + 2xy dx = 13 x 3 + x2 y + ϕ( y )

donde ϕ es una función que depende de y. A continuación, se deriva F(x,y)


parcialmente respecto a y, y se iguala a H (x, y ) = x 2 + 2y 2 :

∂F (x, y )
∂y
=
∂ 1 3
∂y 3
(
x + x2y +

dy
= x2 +

dy
)

x2 + = x 2 + 2y 2
dy


Por último, se despeja en la ecuación formada y se integra respecto a
dy
la variable y:

= 2y 2
dy

ϕ′( y ) = 2y 2
2
∫ ϕ′( y)dy = ∫ 2y dy

ϕ( y ) = 2
3
y3 + C

En conclusión, la solución general de la ecuación diferencial exacta es:


1
3
x3 + x2y + 2
3
y3 + C = 0

Observación. El procedimiento de integración de la ecuación diferencial exacta también


puede realizarse intercambiando las variables en el proceso de integración-derivación-
∂F (x, y )
integración descrito, es decir, primero se trataría de integrar = H (x, y )
∂y
parcialmente respecto a la variable y, después derivar la expresión obtenida
F (x, y ) = ∫ H (x, y )dy + ϕ(x ) parcialmente respecto a la variable x igualando el resultado a
∂F (x, y ) dϕ
= G (x, y ) , y finalmente despejar en la ecuación formada e integrarla respecto
∂x dx
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 67

a x. Normalmente suele comenzarse por la integración parcial de la función G o H que


sea más fácil de integrar.

Ejemplo: Integrar la ecuación diferencial exacta:

(e x
)
+ 3 y dx + (3x + cos( y ))dy = 0

Se trata efectivamente de una EDO exacta, ya que llamando


x
G (x, y ) = e + 3 y , H (x, y ) = 3x + cos( y ) , se cumple que:

∂G (x, y ) ∂H (x, y )
= 3; =3
∂y ∂x

Primero, se integra H (x, y ) = 3x + cos( y ) parcialmente respecto a y:

F (x, y ) = ∫ (3x + cos( y ))dy = 3xy + sen( y ) + ϕ(x)


donde ϕ es una función que depende de x. A continuación, se deriva F(x,y)
parcialmente respecto a x, y se iguala a G (x, y ) = e x + 3 y :

∂F (x, y ) ∂ dϕ dϕ
= (3xy + sen( y )) + = 3y +
∂x ∂x dx dx

3y + = e x + 3y
dx

Por último, se despeja en la ecuación formada y se integra respecto a
dx
la variable x:

= ex
dx

ϕ′(x) = e x
x
∫ ϕ′(x)dx = ∫ e dx

ϕ(x) = e x + C

En conclusión, la solución general de la ecuación diferencial exacta es:

3xy + sen( y ) + e x + C = 0

Ejemplo: Resolver el problema de valor inicial:

(y cos(xy ) − xy sen(xy ))dx + (x cos(xy ) − x ysen(xy ))dy = 0 ; y ( ) = 1


2 2 π
3

Se trata de una ecuación diferencial exacta, ya que llamando


2 2
G (x, y ) = y cos(xy ) − xy sen( xy ) , H ( x, y ) = x cos(xy ) − x ysen(xy ) , se cumple
que:
∂G ( x, y )
= cos(xy ) − 3xysen( xy ) − x2 y 2 cos(xy )
∂y

∂H (x, y )
= cos(xy ) − 3xysen(xy ) − x2 y 2 cos(xy )
∂x
68 Ecuaciones Diferenciales

Primero, se integra G (x, y ) = y cos(xy ) − xy 2sen(xy ) parcialmente respecto


a x:

F (x, y ) = ∫ (y cos(xy ) − xy
2
)
sen(xy ) dx =

= ∫ y cos(xy )dx − ∫ xy 2sen(xy )dx

• ∫ y cos(xy )dx = sen(xy ) + ϕ( y )


 u = xy dv = ysen(xy )dx 
• − ∫ xy 2sen(xy )dx =  =
du = y v = − cos (xy) 

= xy cos(xy ) − ∫ y cos(xy )dx = xy cos(xy ) − sen(xy ) + ϕ( y )

∴ F (x, y ) = xy cos(xy ) + ϕ( y )

donde ϕ es una función que depende de y. A continuación, se deriva F(x,y)


parcialmente respecto a y, y se iguala a H (x, y ) = x cos(xy ) − x2ysen(xy ) :

∂F (x, y ) ∂ dϕ dϕ
= (xy cos(xy )) + = x cos(xy ) − x2 ysen(xy ) +
∂y ∂y dy dy


x cos(xy ) − x2 ysen( xy ) + = x cos(xy ) − x2 ysen( xy )
dy


Por último, se despeja en la ecuación formada y se integra respecto a
dy
la variable y:

=0
dy

ϕ′( y ) = 0

∫ ϕ′( y )dy = ∫ 0 ⋅ dy
ϕ( y ) = C

La solución general de la ecuación diferencial exacta es:


xy cos( xy ) + C = 0

()
Sustituyendo ahora la condición y 3π = 1 en la solución general, se obtiene
una ecuación que permite determinar el valor de la constante de
integración C:

()
y 3π = 1

(3π )(1) cos(3π ) + C = 0


π +C = 0
6

C = −6π

Por tanto, la solución del problema de valor inicial es:


2. Ecuaciones Diferenciales Ordinarias de Primer Orden 69

xy cos(xy ) − π =0
6

2.5. FACTORES INTEGRANTES

2.5.1. Definición de factor integrante

A continuación se estudian ciertas ecuaciones diferenciales no exactas, pero que se


pueden transformar en exactas de una forma relativamente sencilla multiplicándolas
adecuadamente por una función auxiliar µ, llamada factor integrante, es decir,
ecuaciones diferenciales de la forma G (x, y )dx + H ( x, y )dy = 0 , para las que se cumple:

∂G (x, y ) ∂H ( x, y )

∂y ∂x

tales que al multiplicarlas por µ(x, y ) se transforman en ecuaciones diferenciales de la


forma µ(x, y )G (x, y )dx + µ(x, y )H ( x, y )dy = 0 , que verifican:

∂ ∂
(µ(x, y )G (x, y )) = (µ(x, y )H (x, y ))
∂y ∂x

esto es, EDO exactas que pueden resolverse por el método visto de integración-
derivación-integración. La utilidad de esta técnica está en que la solución general de la
ecuación original y de la transformada en exacta es esencialmente la misma, ya que al
multiplicar la EDO por el factor integrante sólo se añadirán algunas soluciones a su
solución general, concretamente el conjunto de soluciones particulares que hacen
µ( x, y ) = 0 (siempre que existan, y no coincidan con las de la EDO original), y
recíprocamente sólo se eliminarán algunas soluciones de la EDO original, aquellas para
las cuales el factor integrante no esté definido (si existen).

Definición (Factor integrante). Dada una EDO de primer orden no exacta


G (x, y )dx + H ( x, y )dy = 0 , es decir, una ecuación diferencial tal que:

∂G (x, y ) ∂H ( x, y )

∂y ∂x

se llama factor integrante de la ecuación diferencial a cualquier función µ : D ⊆ o 2 → o


tal que al multiplicar la ecuación por µ( x, y ) se transforme en una ecuación diferencial
exacta µ(x, y )G (x, y )dx + µ(x, y )H ( x, y )dy = 0 , para la cual se verifica:

∂ ∂
(µ(x, y )G (x, y )) = (µ(x, y )H (x, y ))
∂y ∂x

Ejemplo: La siguiente EDO de primer orden de incógnita y (x ) es una ecuación


diferencial no exacta:

(2x 2
) (
y + y 2 dx + x 3 + 3
2
)
xy dy = 0

ya que, llamando G (x, y ) = 2x 2 y + y 2 , H (x, y ) = x 3 + 32 xy , se cumple que:


70 Ecuaciones Diferenciales

∂G ( x, y ) ∂H (x, y )
= 2x2 + 2 y ≠ = 3x2 + 32 y
∂y ∂x

La función µ(x, y ) = xy es un factor integrante de la ecuación diferencial,


ya que al multiplicarla por él, se obtiene la EDO exacta:

( ) (
xy 2x 2 y + y 2 dx + xy x 3 + 3
2
)
xy dy = 0

(2x 3 2
) ( )
y + xy 3 dx + x 4 y + 32 x 2 y 2 dy = 0

puesto que se cumple:


∂ ∂
(µ(x, y )G (x, y )) = (µ(x, y )H (x, y )) = 4x 3 y + 3xy 2
∂y ∂x

Ejemplo: La siguiente EDO de primer orden de incógnita x (t ) es una ecuación


diferencial no exacta:

 tsen(tx ) 2 
(sen(tx ) − t cos(x ))dt +  + t sen(x ) dx = 0
 x 

tsen(tx )
ya que, llamando G (t, x ) = sen(tx ) − t cos( x) , H (t, x) = + t 2sen(x) ,
x
se cumple que:
∂G (t, x ) ∂H (t, x ) sen(tx )
= t cos(tx ) + tsen(x ) ≠ = + t cos(tx ) + 2tsen( x )
∂x ∂t x
1
La función µ(t, x ) = es un factor integrante de la ecuación diferencial, ya
t
que al multiplicarla por él, se obtiene la EDO exacta:
1 1  tsen(tx ) 2 
(sen(tx ) − t cos(x ))dt +  + t sen( x ) dx = 0
t t x 

 sen(tx )   sen(tx ) 
 − cos( x) dt +  + tsen(x) dx = 0
 t   x 
puesto que se cumple:
∂ ∂
(µ(t, x )G (t, x)) = (µ(t, x)H (t, x)) = cos(tx) + sen(x)
∂x ∂t

Observación. En la práctica, la búsqueda de factores integrantes se realiza para


ecuaciones diferenciales no exactas de la forma G (x, y )dx + H ( x, y )dy = 0 , en las que las
∂G (x, y ) ∂H (x, y )
derivadas parciales , sean suficientemente parecidas, como ocurre en
∂y ∂x
los ejemplos precedentes.

2.5.2. Cálculo de factores integrantes

Sea una EDO de primer orden no exacta G (x, y )dx + H (x, y )dy = 0 , es decir, una ecuación
diferencial tal que:
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 71

∂G (x, y ) ∂H ( x, y )

∂y ∂x

y sea µ : D ⊆ o 2 → o un factor integrante que convierte a la ecuación en una ecuación


diferencial exacta de la forma µ(x, y )G (x, y )dx + µ(x, y )H ( x, y )dy = 0 . Por definición de
ecuación diferencial exacta se debe cumplir:
∂ ∂
(µ(x, y )G (x, y )) = (µ(x, y )H (x, y ))
∂y ∂x
de donde, aplicando la regla de derivación del producto, resulta:
∂µ(x, y ) ∂G ( x, y ) ∂µ( x, y ) ∂H (x, y )
G ( x, y ) + µ(x, y ) = H (x, y ) + µ( x, y )
∂y ∂y ∂x ∂x

∂G ( x, y ) ∂H ( x, y ) ∂µ(x, y ) ∂µ(x, y )
µ(x, y ) − µ( x, y ) = H ( x, y ) − G ( x, y )
∂y ∂x ∂x ∂y

 ∂G ( x, y ) ∂H ( x, y )  ∂µ(x, y ) ∂µ( x, y )
µ( x, y ) −  = H ( x, y ) − G ( x, y )
 ∂ y ∂ x  ∂x ∂y

que es la ecuación diferencial general de cálculo del factor integrante. Esta ecuación es
una ecuación en derivadas parciales, normalmente mucho más difícil de resolver que la
EDO original. De ahí que resulte inviable considerar el cálculo general de factores
integrantes para ecuaciones diferenciales ordinarias. Sin embargo, la EDD de cálculo del
factor integrante es resoluble en algunos casos particulares, por ejemplo, si se posee
alguna información adicional respecto al factor integrante µ, tal como que µ dependa
sólo de la variable x, o sólo de la variable y, o de una expresión concreta de dos variables
como la suma x+y o el producto xy. Precisamente, estos casos son los que se van a
considerar en el cálculo práctico de factores integrantes.

Proposición (Cálculo del factor integrante dependiente de la variable x). La condición


necesaria y suficiente para que una EDO de primer orden no exacta
G (x, y )dx + H ( x, y )dy = 0 admita un factor integrante µ dependiente sólo de la variable x
es que se verifique:
∂G (x, y ) ∂H (x, y )

∂y ∂x
= f (x )
H (x, y )

siendo el factor integrante igual a:

µ(x) = Ce ∫ f ( x )dx

Demostración:

En efecto, sea µ un factor integrante dependiente sólo de la variable x que convierte a la ecuación en una
ecuación diferencial exacta de la forma µ(x )G (x, y )dx + µ(x )H (x, y )dy = 0 . Por definición de ecuación diferencial
exacta se debe cumplir:


(µ(x)G (x, y )) = ∂ (µ(x)H(x, y ))
∂y ∂x

y de aquí:
72 Ecuaciones Diferenciales

∂G (x, y ) dµ ∂H (x, y )
µ( x ) = H (x, y ) + µ( x )
∂y dx ∂x

∂G (x, y ) ∂H ( x, y ) dµ
µ( x ) − µ( x ) = H (x, y )
∂y ∂x dx

 ∂G (x, y ) ∂H (x, y )  dµ
µ(x ) −  = H (x, y )
 ∂y ∂x  dx

∂G (x, y ) ∂H (x, y )

1 dµ ∂y ∂x
⋅ =
µ(x ) dx H (x, y )

1 dµ
y como ⋅ es una función que depende sólo de la variable x, de la igualdad anterior resulta la condición
µ( x ) dx
referida:

∂G (x, y ) ∂H (x, y )

∂y ∂x
= f (x )
H (x, y )

Para hallar el factor integrante, basta con integrar la condición:

1 dµ
⋅ = f (x )
µ(x ) dx

µ′(x )
= f (x )
µ( x )

µ′(x )
∫ µ(x) dx = ∫ f (x)dx

( )
ln µ(x ) = ∫ f (x )dx + C

µ(x ) = ece ∫ f ( x )dx

µ(x ) = ±ece ∫ f ( x )dx

µ(x ) = Ce ∫ f ( x )dx

Ejemplo: Hallar un factor integrante dependiente sólo de la variable x para la


ecuación diferencial no exacta:

 2 y
 x
2
 xy − dx + 2x y + 1 dy = 0 ( )
En primer lugar, se comprueba que la EDO admite un factor integrante de
y
la forma µ(x ) . Llamando G (x, y ) = xy 2 − , H (x, y ) = 2x2 y + 1 , se cumple
x
que:
∂G (x, y ) ∂H (x, y ) 1 1
− 2xy − − 4 xy − − 2xy
∂y ∂x x x − 1 − 2x2 y −1
= = = = = f (x )
H ( x, y ) 2
2x y + 1 2
2x y + 1 2
x 2x y + 1 x ( )
El factor integrante es:
−1
∫ dx −ln( x ) C 1
µ( x ) = Ce ∫ f ( x )dx = Ce x = Ce = = {C = 1} =
x x
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 73

Proposición (Cálculo del factor integrante dependiente de la variable y). La condición


necesaria y suficiente para que una EDO de primer orden no exacta
G (x, y )dx + H ( x, y )dy = 0 admita un factor integrante µ dependiente sólo de la variable y
es que se verifique:
∂H (x, y ) ∂G (x, y )

∂x ∂y
= f (y)
G (x, y )

siendo el factor integrante igual a:

µ( y ) = Ce ∫ f ( y )dy

Demostración:

En efecto, sea µ un factor integrante dependiente sólo de la variable y que convierte a la ecuación en una
ecuación diferencial exacta de la forma µ( y )G ( x, y )dx + µ( y )H (x, y )dy = 0 . Por definición de ecuación diferencial
exacta se debe cumplir:


(µ( y )G (x, y )) = ∂ (µ( y )H(x, y ))
∂y ∂x

y de aquí:

dµ ∂G (x, y ) ∂H ( x, y )
G (x, y ) + µ( y ) = µ( y )
dy ∂y ∂x

dµ ∂H (x, y ) ∂G (x, y )
G (x, y ) = µ( y ) − µ( y )
dy ∂x ∂y

dµ  ∂H (x, y ) ∂G (x, y ) 
G (x, y ) = µ( y ) − 
dy  ∂x ∂y 

∂H (x, y ) ∂G (x, y )

1 dµ ∂x ∂y
⋅ =
µ( y ) dy G (x, y )

1 dµ
y como ⋅ es una función que depende sólo de la variable y, de la igualdad anterior resulta la condición
µ( y ) dy
referida:

∂H (x, y ) ∂G (x, y )

∂x ∂y
= f ( y)
G (x, y )

Para hallar el factor integrante, basta con integrar la condición:

1 dµ
⋅ = f (y)
µ( y ) dy

µ′( y )
= f ( y)
µ( y )

µ′( y )
∫ µ( y ) dy = ∫ f ( y )dy

( )
ln µ( y ) = ∫ f ( y )dy + C

µ( y ) = ece ∫ f ( y )dy
74 Ecuaciones Diferenciales

µ( y ) = ±ece ∫ f ( y )dy

µ( y ) = Ce ∫ f ( y )dy

Ejemplo: Hallar un factor integrante dependiente sólo de la variable y para la


ecuación diferencial no exacta:

(xy 2
) (
− y 3 dx + 1 − xy 2 dy = 0 )
En primer lugar, se comprueba que la EDO admite un factor integrante de
la forma µ(y ) . Llamando G (x, y ) = xy 2 − y 3 , H (x, y ) = 1 − xy 2 , se cumple
que:
∂H (x, y ) ∂G (x, y )

∂x ∂y − y 2 − 2xy + 3 y 2 − 2xy + 2 y 2
= = =
G (x, y ) xy 2 − y 3 xy 2 − y 3

2y ( y − x) −2
= 2
= = f ( y)
− y ( y − x) y

El factor integrante es:


−2
∫ dy
−2 ln( y ) C 1
µ( y ) = Ce ∫ f ( y )dy = Ce y
= Ce = = {C = 1} =
y2 y2

Proposición (Cálculo del factor integrante dependiente de una expresión w(x,y)). La


condición necesaria y suficiente para que una EDO de primer orden no exacta
G (x, y )dx + H (x, y )dy = 0 admita un factor integrante µ dependiente de una expresión
sencilla w(x,y) es que se verifique:
∂G ( x, y ) ∂H (x, y )

∂y ∂x
= f (w (x, y ))
∂w (x, y ) ∂w (x, y )
H (x, y ) − G (x, y )
∂x ∂y

siendo el factor integrante igual a:

µ(w ) = Ce ∫ f (w )dw

Demostración:

En efecto, sea µ un factor integrante dependiente de una expresión w(x,y) que convierte a la ecuación en una
ecuación diferencial exacta de la forma µ(w (x, y ))G (x, y )dx + µ(w (x, y ))H (x, y )dy = 0 . Por definición de ecuación
diferencial exacta se debe cumplir:


(µ(w (x, y ))G (x, y )) = ∂ (µ(w (x, y ))H(x, y ))
∂y ∂x

y de aquí:

∂µ(w (x, y )) ∂G (x, y ) ∂µ(w (x, y )) ∂H (x, y )


G (x, y ) + µ(w (x, y )) = H (x, y ) + µ(w (x, y ))
∂y ∂y ∂x ∂x

dµ ∂w (x, y ) ∂G (x, y ) dµ ∂w (x, y ) ∂H (x, y )


G (x, y ) ⋅ + µ(w (x, y )) = H (x, y ) ⋅ + µ(w (x, y ))
dw ∂y ∂y dw ∂x ∂x
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 75

∂G (x, y ) ∂H (x, y ) dµ ∂w (x, y ) dµ ∂w (x, y )


µ(w (x, y )) − µ(w (x, y )) = H (x, y ) ⋅ − G (x, y ) ⋅
∂y ∂x dw ∂x dw ∂y

 ∂G (x, y ) ∂H (x, y )  dµ  ∂w (x, y ) ∂w (x, y ) 


µ(w (x, y )) −  =  H (x, y ) − G ( x, y ) 
 ∂y ∂ x  dw  ∂ x ∂y 

∂G (x, y ) ∂H (x, y )

1 dµ ∂y ∂x
⋅ =
µ(w (x, y )) dw H (x, y ) ∂w (x, y ) − G (x, y ) ∂w (x, y )
∂x ∂y

1 dµ
y como ⋅ es una función que depende de la expresión w(x,y), de la igualdad anterior resulta la
µ(w (x, y )) dw
condición referida:

∂G ( x, y ) ∂H (x, y )

∂y ∂x
= f (w (x, y ))
∂w (x, y ) ∂w (x, y )
H (x, y ) − G (x, y )
∂x ∂y

Para hallar el factor integrante, basta con integrar la condición:

1 dµ
⋅ = f (w )
µ(w ) dw

µ′(w )
= f (w )
µ(w )

µ′(w )
∫ µ(w ) dw = ∫ f (w )dw

( )
ln µ(w ) = ∫ f (w )dw + C

µ(w ) = ece ∫ f (w )dw

µ(w ) = ±ece ∫ f (w )dw

µ(w ) = Ce ∫ f (w )dw

Ejemplo: Hallar un factor integrante dependiente del producto xy para la ecuación


diferencial no exacta:

(y + xy )dx + (x − x y )dy = 0
2 2

En primer lugar, se comprueba que la EDO admite un factor integrante de


la forma µ(xy ) . Llamando G (x, y ) = y + xy 2 , H (x, y ) = x − x2 y , se cumple
que:
∂G (x, y ) ∂H (x, y )

∂y ∂x 1 + 2xy − 1 + 2xy
= =
H (x, y )
∂w ( x, y )
− G (x, y )
∂w (x, y ) x − x 2
y y − y (
+ xy 2
x ) ( )
∂x ∂y

4xy 4xy −2 −2
= 2 2 2 2
= 2 2
= = = f (w (x, y ))
xy − x y − xy − x y − 2x y xy w

−2
∫ dw −2ln(w ) C C 1
µ(w ) = Ce ∫ f (w )dw = Ce w = Ce = 2
= 2 2
= {C = 1} = 2 2
w x y x y
76 Ecuaciones Diferenciales

2.5.3. Resolución de ecuaciones diferenciales mediante factores


integrantes

La resolución de ecuaciones diferenciales de primer orden por transformación en


ecuaciones exactas mediante factores integrantes conlleva en primer lugar la búsqueda
de un factor integrante apropiado, lo que es un problema no trivial, tal como se ha visto
en el subapartado anterior. El siguiente procedimiento intuitivo describe el método
práctico de cálculo de factores integrantes y resolución de la EDO dada.

Procedimiento (Resolución de una EDO de primer orden mediante un factor


integrante). Para resolver una ecuación diferencial de primer orden no exacta
G (x, y )dx + H (x, y )dy = 0 mediante un factor integrante se procede de la siguiente forma:

∂G (x, y ) ∂H (x, y )
1. Se calculan las derivadas parciales , y se estudia la existencia
∂y ∂x
del factor integrante considerando tres casos:
∂G (x, y ) ∂H (x, y )

∂y ∂x
• Si = f (x ) , entonces el factor integrante es de la forma
H (x, y )
µ(x) = Ce ∫ f ( x )dx .

∂H (x, y ) ∂G (x, y )

∂x ∂y
• Si = f ( y ) , entonces el factor integrante es de la forma
G (x, y )
µ( y ) = Ce ∫ f ( y )dy .

• Si es posible encontrar una expresión sencilla w(x,y), tal como x+y, x−y,
∂G ( x, y ) ∂H (x, y )

∂y ∂x
xy, x/y, para la cual se verifica = f (w (x, y )) ,
∂w (x, y ) ∂w (x, y )
H (x, y ) − G (x, y )
∂x ∂y
entonces el factor integrante es de la forma µ(w ) = Ce ∫ f (w )dw .

2. Se multiplica la EDO original G (x, y )dx + H (x, y )dy = 0 por el factor integrante
µ( x, y ) , obteniéndose la EDO exacta:

µ(x, y )G (x, y )dx + µ(x, y )H ( x, y )dy = 0

3. Se resuelve la EDO exacta, cuya solución general es equivalente a la de la


original, mediante el procedimiento conocido de integración-derivación-
integración:

• F (x, y ) = ∫ µ(x, y )G (x, y )dx + ϕ( y )

∂F (x, y )

∂y
=

∂y
(∫ µ(x, y)G(x, y)dx ) + ddyϕ = µ(x, y)H(x, y)

• ϕ( y ) = ∫  µ(x, y )H (x, y ) −

∂y
(∫ µ(x, y)G(x, y)dx )dy
 
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 77

Ejemplo: Integrar la siguiente EDO de primer orden mediante un factor integrante:


 
(
 y + xy 2 − 2x dx + 3x + 2x 2 y dy = 0
 y 2 
)

Se trata efectivamente de una EDO no exacta, ya que llamando
2x
G (x, y ) = y + xy 2 − , H (x, y ) = 3x + 2x 2 y , se cumple que:
y2

∂G (x, y ) 4x ∂H (x, y )
= 1 + 2xy + ≠ = 3 + 4xy
∂y y3 ∂x

Cálculo del factor integrante:


∂G (x, y ) ∂H (x, y )

∂y ∂x
=
∂w ( x, y ) ∂w (x, y )
H (x, y ) − G (x, y )
∂x ∂y

4x
1 + 2xy +
− 3 − 4xy
y3
= =
 2x  ∂w (x, y )
( ∂x
)
2 ∂w ( x, y ) 
3x + 2x y 2
− y + xy −
 y 2  ∂y

y 3 + 2xy 4 + 4 x − 3 y 3 − 4 xy 4
= =
 2x  ∂w (x, y )
3
( 2 ∂w ( x, y )
y 3x + 2x y
∂x
)
3
− y y + xy −

2
y 2  ∂y

− 2 y 3 − 2xy 4 + 4 x
=
(3xy 3
+ 2x 2 y 4 ) ∂w (x, y )
∂x
(
− y 4 + xy 5 − 2xy
∂w (x, y )
∂y
)
∂w ∂w
• Caso factor integrante µ(x ) ( w = x , = 1, = 0 ):
∂x ∂y

− 2 y 3 − 2xy 4 + 4 x
=
(3xy 3
) (
+ 2x 2 y 4 (1) − y 4 + xy 5 − 2xy (0) )
− 2y 3 − 2xy 4 + 4x
= ≠ f (x)
3xy 3 + 2x 2 y 4

∂w ∂w
• Caso factor integrante µ(y ) ( w = y , = 0, = 1 ):
∂x ∂y

− 2 y 3 − 2xy 4 + 4 x
=
(3xy 3
) (
+ 2x 2 y 4 (0) − y 4 + xy 5 − 2xy (1) )
2y 3 + 2xy 4 − 4x 2
= = = f ( y)
(y 3
+ xy 4
− 2x y ) y

Por tanto, la EDO admite un factor integrante dependiente de la variable


y, de la forma:
78 Ecuaciones Diferenciales

2
∫ y dy 2 ln( y )
µ( y ) = Ce ∫ f ( y )dy = Ce = Ce = Cy 2 = {C = 1} = y 2

Se multiplica la EDO original por el factor integrante µ( y ) = y 2 ,


obteniéndose la EDO exacta:

(y 3
) (
+ xy 4 − 2x dx + 3xy 2 + 2x 2 y 3 dy = 0 )
Se integra µ( y )G (x, y ) = y 3 + xy 4 − 2x parcialmente respecto a x:

∫ (y )
3
F (x, y ) = + xy 4 − 2x dx = xy 3 + 1
2
x 2 y 4 − x 2 + ϕ( y )

donde ϕ es una función que depende de y. A continuación, se deriva F(x,y)


parcialmente respecto a y, y se iguala a µ( y )H (x, y ) = 3xy 2 + 2x 2 y 3 :

∂F (x, y )
∂y
=

∂y
xy 3 + ( 1
2
x2y 4 − x2 + ) dϕ
dy
= 3xy 2 + 2x 2 y 3 +

dy


3xy 2 + 2x 2 y 3 + = 3xy 2 + 2x 2 y 3
dy


Por último, se despeja en la ecuación formada y se integra respecto a
dy
la variable y:

=0
dy

ϕ′( y ) = 0

∫ ϕ′( y )dy = ∫ 0dy


ϕ( y ) = C

En conclusión, la solución general de la ecuación diferencial es:

xy 3 + 12 x 2 y 4 − x 2 + C = 0

Ejemplo: Integrar la siguiente EDO de primer orden mediante un factor integrante:

 2
(x ) 
(2xsen(x) cos(x) + xysen(x) )dx +  xsen − 2x cos( x) dy = 0

 y 
Se trata efectivamente de una EDO no exacta, ya que llamando
xsen2 (x )
G (x, y ) = 2xsen(x ) cos( x) + xysen(x) , H (x, y ) = − 2x cos( x ) , se
y
cumple que:
∂G (x, y )
= xsen(x) ≠
∂y

∂H (x, y ) sen2 (x) + 2xsen(x) cos( x)


≠ = − 2 cos( x ) + 2xsen(x )
∂x y

Cálculo del factor integrante:


2. Ecuaciones Diferenciales Ordinarias de Primer Orden 79

∂G (x, y ) ∂H (x, y )

∂y ∂x
=
∂w ( x, y ) ∂w (x, y )
H (x, y ) − G (x, y )
∂x ∂y

sen2 (x) + 2xsen(x) cos( x)


xsen(x) − + 2 cos( x ) − 2xsen( x)
y
= =
 xsen2 (x)  ∂w (x, y ) ∂w ( x, y )
 − 2x cos( x)  − (2xsen(x) cos( x ) + xysen( x) )
 y  ∂x ∂y
 

sen2 (x) + 2xsen(x) cos(x )


2 cos(x ) − xsen(x) −
y
= =
 xsen2 (x )  ∂w (x, y ) ∂w (x, y )
 − 2x cos(x )  − (2xsen(x) cos( x) + xysen(x) )
 y  ∂x ∂y
 

2 y cos( x ) − xysen(x) − sen2 (x) − 2xsen(x) cos( x )


= =
(
xsen2 ( x) − 2xy cos( x )
∂x
)
∂w (x, y )
(
− 2xysen(x ) cos( x) + xy 2sen(x)
∂w (x, y )
∂y
)
− sen2 (x) + 2y cos( x) − 2xsen(x ) cos( x) − xysen(x )
= =
(
x sen2 ( x) − 2y cos( x) )
∂w (x, y )
∂x
− y (2xsen( x) cos( x ) + xysen(x ) )
∂w ( x, y )
∂y

∂w ∂w
• Caso factor integrante µ(x ) ( w = x , = 1, = 0 ):
∂x ∂y

− sen2 (x) + 2 y cos( x) − 2xsen(x) cos( x ) − xysen(x)


=
( )
x sen2 (x) − 2y cos( x) (1) − y (2xsen(x ) cos( x) + xysen(x ) )(0)

2y cos( x) − xysen(x ) − sen2 (x ) − 2xsen(x ) cos(x)


≠ f (x)
(
x sen2 (x) − 2y cos(x ) )
∂w ∂w
• Caso factor integrante µ(y ) ( w = y , = 0, = 1 ):
∂x ∂y

− sen2 (x) + 2 y cos( x) − 2xsen(x) cos(x ) − xysen(x)


=
( )
x sen2 (x) − 2y cos(x) (0) − y (2xsen(x) cos(x ) + xysen( x) )(1)

sen2 (x) + 2xsen( x) cos( x) − 2 y cos( x ) + xysen(x )


= ≠ f ( y)
y (2xsen(x) cos( x ) + xysen(x ) )

y y ∂w − y ∂w 1
• Caso factor integrante µ  ( w = , = , = ):
 
x x ∂x x 2 ∂y x

− sen2 (x ) + 2 y cos( x) − 2xsen(x) cos( x ) − xysen(x)


=
− y
(
x sen2 (x ) − 2 y cos( x )  ) 1
 − y (2xsen(x) cos( x ) + xysen(x ) ) 
 x2  x
80 Ecuaciones Diferenciales

− sen2 (x ) + 2 y cos( x ) − 2xsen(x ) cos( x ) − xysen(x )


= =
( y
x
) y
− sen2 (x ) + 2 y cos( x )   − (2xsen(x ) cos( x ) + xysen(x ) ) 
x

− sen2 (x) + 2y cos( x) − 2xsen(x) cos( x) − xysen(x ) 1 1


= = = = f (w (x, y ) )
( 2  y
− sen (x ) + 2y cos( x ) − 2xsen(x ) cos( x) − xysen(x )  
x
 y
x
w
)
Por tanto, la EDO admite un factor integrante dependiente de y/x, de la
forma:
1
∫ w dw ln(w ) y y
µ(w ) = Ce ∫ f (w )dw = Ce = Ce = Cw = C = {C = 1} =
x x

y
Se multiplica la EDO original por el factor integrante µ(x, y ) = ,
x
obteniéndose la EDO exacta:

(2ysen(x) cos(x) + y sen(x))dx + (sen (x) − 2y cos(x))dy = 0


2 2

Se integra µ(x, y )G (x, y ) = 2 ysen(x) cos( x) + y 2sen( x) parcialmente respecto


a x:

∫ (2ysen(x) cos(x) + y )
2
F (x, y ) = sen(x ) dx = ysen2 (x ) − y 2 cos( x) + ϕ( y )

donde ϕ es una función que depende de y. A continuación, se deriva F(x,y)


parcialmente respecto a y, y se iguala a
µ( x, y )H (x, y ) = sen2 (x ) − 2y cos( x) :

∂F (x, y )
∂y
=

∂y
(
ysen2 (x) − y 2 cos( x) +

dy
)
= sen2 (x) − 2y cos( x) +

dy


sen2 (x ) − 2y cos( x) + = sen2 (x) − 2y cos(x )
dy


=0
dy


Por último, se despeja en la ecuación formada y se integra respecto a
dy
la variable y:

=0
dy

ϕ′( y ) = 0

∫ ϕ′( y )dy = ∫ 0dy


ϕ( y ) = C

En conclusión, la solución general de la ecuación diferencial es:

ysen2 (x ) − y 2 cos( x ) + C = 0
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 81

Ejemplo: Resolver el siguiente problema de valor inicial mediante un factor


integrante dependiente del producto xy:

   
 2 + 1 dx +  x − 1 dy = 0 ; y (1) = 1
 xy 2   y y3 
  
Se trata efectivamente de una EDO no exacta, ya que llamando
1 x 1
G (x, y ) = 2 + 2 , H (x, y ) = − 3 , se cumple que:
xy y y

∂G (x, y ) − 2 ∂H (x, y ) 1
= 3 ≠ =
∂y xy ∂x y

Cálculo del factor integrante:


∂G (x, y ) ∂H (x, y ) −2 1
− −
∂y ∂x xy 3 y
= =
∂w (x, y ) ∂w (x, y )  x 1   1 
H (x, y ) − G (x, y )  − y − 2 + x
∂x ∂y  y y3   xy 2 
  

−2 1 −2 1 − 2 − xy 2
− −
xy 3 y xy 3 xy 3 y
1 1
= = = 2
= = = f (w (x, y ))
1 1 2 − 2 − xy xy w
x − 2 − 2x − 2 − 2 −x
y y y y2

Efectivamente, la EDO admite un factor integrante dependiente del


producto xy, de la forma:
1
∫ dw ln(w )
µ(w ) = Ce ∫ f (w )dw = Ce w = Ce = C w = {C = 1} = w = xy

Se multiplica la EDO original por el factor integrante µ( x, y ) = xy ,


obteniéndose la EDO exacta:

 1  x 
 2xy + dx +  x2 − dy = 0
 y  y 2 

1
Se integra µ(x, y )G (x, y ) = 2xy + parcialmente respecto a x:
y

 1 x
F (x, y ) = ∫  2xy + dx = x2 y + + ϕ( y )
 y  y

donde ϕ es una función que depende de y. A continuación, se deriva F(x,y)


x
parcialmente respecto a y, y se iguala a µ( x, y )H ( x, y ) = x2 − 2 :
y

∂F (x, y ) ∂  2 x  dϕ x dϕ
=  x y +  + = x2 − 2 +
∂y ∂y  y  dy y dy

x dϕ x
x2 − 2
+ = x2 − 2
y dy y
82 Ecuaciones Diferenciales


Por último, se despeja en la ecuación formada y se integra respecto a
dy
la variable y:

=0
dy

ϕ′( y ) = 0

∫ ϕ′( y )dy = ∫ 0dy


ϕ( y ) = C

luego, la solución general de la ecuación diferencial es:


x
x2 y + +C = 0
y

Sustituyendo ahora la condición y (1) = 1 en la solución general, se obtiene


una ecuación que permite determinar el valor de la constante de
integración C:
y (1) = 1

1
(1)2 (1) + +C = 0
1
2+C = 0
C = −2
Por tanto, la solución del problema de valor inicial es:
x
x2 y + −2 = 0
y

2.5.4. Multiplicidad del factor integrante

En este subapartado se demuestra un importante teorema sobre factores integrantes


que proporciona una forma de resolver una EDO de primer orden no exacta sin
integrarla, siempre que se conozcan dos factores integrantes distintos de dicha
ecuación.

Teorema (Multiplicidad del factor integrante). Dada una EDO de primer orden no exacta
G (x, y )dx + H (x, y )dy = 0 , si µ1(x, y ) y µ2 ( x, y ) son dos factores integrantes distintos de la
misma, entonces la solución general de la ecuación es:
µ1(x, y )
=C
µ2 ( x, y )

Demostración:

Diferenciando la expresión dada como solución general de la EDO se obtiene:

∂  µ1(x, y )  ∂  µ1(x, y ) 
 dx +  dy = 0
 
∂x  µ2 (x, y )  ∂y  µ2 (x, y ) 
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 83

∂µ1(x, y ) ∂µ (x, y ) ∂µ (x, y ) ∂µ (x, y )


µ2 (x, y ) − µ1(x, y ) 2 µ2 (x, y ) 1 − µ1(x, y ) 2
∂x ∂x ∂y ∂y
dx + dy = 0
µ22
( x, y ) µ 2
2
( x, y )

 ∂µ (x, y ) ∂µ (x, y )   ∂µ (x, y ) ∂µ (x, y ) 


 µ2 (x, y ) 1 − µ1(x, y ) 2 dx +  µ2 (x, y ) 1 − µ1(x, y ) 2 dy = 0
 ∂x ∂x   ∂y ∂y 

Por tanto, se trata de probar que la ecuación diferencial anterior es equivalente a la original
G (x, y )dx + H (x, y )dy = 0 , es decir, poniendo ambas ecuaciones en forma normal debe verificarse que:

∂µ ( x, y ) ∂µ (x, y )
µ1(x, y ) 2 − µ2 (x, y ) 1
dy − G (x, y ) ∂x ∂x
= =
dx H (x, y ) ∂µ ( x, y ) ∂µ (x, y )
µ2 (x, y ) 1 − µ1(x, y ) 2
∂y ∂y

En efecto, como µ1(x, y ) y µ2 (x, y ) son dos factores integrantes de la ecuación original, al multiplicarla por
cada uno de ellos se obtienen sendas ecuaciones diferenciales exactas:

µ1(x, y )G (x, y )dx + µ1(x, y )H ( x, y )dy = 0

µ2 (x, y )G (x, y )dx + µ2 ( x, y )H ( x, y )dy = 0

para las cuales se verifica por definición de EDO exacta que:

∂ ∂
(µ1(x, y )G (x, y )) = (µ1(x, y )H (x, y ))
∂y ∂x

∂ ∂
(µ2 (x, y )G (x, y )) = (µ2 (x, y )H (x, y ))
∂y ∂x

de donde resulta:

 ∂G (x, y ) ∂H (x, y )  ∂µ (x, y ) ∂µ (x, y )


µ1(x, y ) −  = H (x, y ) 1 − G (x, y ) 1
 ∂y ∂x  ∂x ∂y

 ∂G (x, y ) ∂H (x, y )  ∂µ (x, y ) ∂µ (x, y )


µ2 (x, y ) −  = H (x, y ) 2 − G (x, y ) 2
 ∂y ∂ x  ∂x ∂y

multiplicando la primera de estas igualdades por µ2 (x, y ) , la segunda por µ1(x, y ) , y restándolas miembro a
miembro, se obtiene:

 ∂µ (x, y ) ∂µ (x, y )   ∂µ (x, y ) ∂µ (x, y ) 


H (x, y ) µ2 (x, y ) 1 − µ1(x, y ) 2  − G (x, y ) µ2 (x, y ) 1 − µ1(x, y ) 2  = 0
 ∂x ∂x   ∂y ∂y 

−G
y despejando resulta:
H

 ∂µ (x, y ) ∂µ (x, y )   ∂µ (x, y ) ∂µ (x, y ) 


− G (x, y ) µ2 (x, y ) 1 − µ1(x, y ) 2  = H(x, y ) µ1(x, y ) 2 − µ2 (x, y ) 1 
 ∂y ∂y   ∂x ∂x 

∂µ (x, y ) ∂µ (x, y )
µ1(x, y ) 2 − µ2 (x, y ) 1
− G (x, y ) ∂x ∂x
=
H (x, y ) ∂µ (x, y ) ∂µ (x, y )
µ2 ( x, y ) 1 − µ1(x, y ) 2
∂y ∂y

como se quería demostrar.


84 Ecuaciones Diferenciales

 y    y 
Ejemplo: Dada la EDO  y cos  − x dx +  − x cos  dy = 0 , deducir la solución
  
x    x 
general a partir de un factor integrante dependiente sólo de x, y un factor
y
integrante dependiente de .
x
Se trata efectivamente de una EDO no exacta, ya que llamando
y y
G (x, y ) = y cos  − x , H (x, y ) = − x cos  , se cumple que:
x x

∂G (x, y ) y y  y  ∂H (x, y ) y y y


= cos  − sen  ≠ = − cos  − sen 
∂y x x x ∂x x x x
Cálculo del factor integrante dependiente sólo de x:
∂G (x, y ) ∂H (x, y ) y y y y y y
− cos  − sen  + cos  + sen 
∂y ∂x
= x x x x x x =
H ( x, y ) y
− x cos 
x

y
2 cos 
=  x  = − 2 = f (x )
y x
− x cos 
 
x

Efectivamente, la EDO admite un factor integrante dependiente de x, de la


forma:
−2
∫ dx −2 ln( x ) C 1
µ1(x ) = Ce ∫ f ( x )dx = Ce x = Ce = = {C = 1} =
x2 x2
Cálculo del factor integrante dependiente de y/x:
∂G (x, y ) ∂H (x, y ) y y y y y y
− cos  − sen  + cos  + sen 
∂y ∂x
=  
x x  
x  
x x x =
H (x, y )
∂w (x, y )
− G (x, y )
∂w (x, y ) y  y 1 y 
cos  −  y cos  − x 
∂x ∂y x x x x 

y
2 cos 
x y
= = 2 cos  = 2 cos(w ) = f (w ( x, y ))
y y y y x
cos  − cos  + 1
x x x x
Efectivamente, la EDO admite un factor integrante dependiente del
cociente y/x, de la forma:
y
2sen 
= Ce ∫
2 cos(w )dw
µ2 (w ) = Ce ∫ f (w )dw = Ce2sen(w ) = {C = 1} = e2sen(w ) = e x

Por tanto, la solución general de la EDO es:


µ 1 (x, y )
=C
µ 2 (x, y )
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 85

1
x2 =C
y
2sen 
e x

y
−2sen 
1 x
e =C
x2

2.6. ECUACIÓN LINEAL

2.6.1. Definición de ecuación lineal

Definición (Ecuación lineal de primer orden). Una ecuación diferencial ordinaria de


primer orden es una ecuación lineal si y sólo si se expresa en la forma:
a1(x ) ⋅ y ′ + a0 (x ) ⋅ y = b(x )

donde a0, a1 y b son expresiones analíticas elementales, con a1(x ) ≠ 0 . Asimismo,


dividiendo por a1(x) se obtiene otra forma equivalente de la ecuación lineal:
y ′ + G (x) ⋅ y = H (x )

a0 (x) b( x )
donde G (x ) = y H (x ) = . Si b(x ) = 0 (respectivamente H ( x ) = 0 ), la ecuación
a1(x) a1(x )
lineal se llama homogénea, y en caso contrario, b( x) ≠ 0 (respectivamente H ( x) ≠ 0 ), la
ecuación lineal se llama completa. Además, si a0 y a1 son constantes (respectivamente G
es constante) la ecuación lineal se llama ecuación lineal de coeficientes constantes, y en
caso contrario, si a0 o a1 dependen de x (respectivamente G depende de x), la ecuación
lineal se llama ecuación lineal de coeficientes variables.

Ejemplo: La ecuación diferencial ordinaria de incógnita y (x ) :

y′ + 5y = 0

es una ecuación lineal homogénea de coeficientes constantes.

Ejemplo: La ecuación diferencial ordinaria de incógnita y (x ) :

dy
tg (x ) + xy = 0
dx
es una ecuación lineal homogénea de coeficientes variables.

Ejemplo: La ecuación diferencial ordinaria de incógnita x (t ) :

dx
+ x = 2t
dt
es una ecuación lineal completa de coeficientes constantes.
86 Ecuaciones Diferenciales

Ejemplo: La ecuación diferencial ordinaria de incógnita y (x ) :

y
y′ + = ln(1 + x )
x
es una ecuación lineal completa de coeficientes variables.

2.6.2. Resolución de ecuaciones lineales homogéneas

La resolución de la ecuación lineal homogénea y ′ + G (x ) ⋅ y = 0 es directa, ya que se trata


de una ecuación de variables separadas.

Procedimiento (Resolución de una EDO lineal homogénea). Para resolver una ecuación
lineal homogénea y ′ + G (x) ⋅ y = 0 (EDO de variables separadas) se procede de la
siguiente forma:
1. Se separan los diferenciales en la igualdad, poniendo dy a la izquierda y dx a la
derecha:
dy
+ G (x) ⋅ y = 0
dx
dy
= −G (x ) ⋅ y
dx
dy = −G (x) ⋅ ydx

2. Se divide la ecuación por y :

1
dy = −G (x )dx
y

3. Se integra en cada miembro respecto a su variable:


1
∫ y dy = ∫ − G (x)dx
ln( y ) = ∫ − G (x)dx + C

y = eCe ∫ −G ( x )dx

y = Ce ∫ −G ( x )dx

Ejemplo: Integrar la EDO lineal homogénea:


dy 2xy
− 2 =0
dx x + 4x + 5
Siguiendo el procedimiento descrito, se tiene:
dy 2xy
= 2
dx x + 4x + 5

2xy
dy = 2
dx
x + 4x + 5
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 87

1 2x
dy = 2 dx
y x + 4x + 5

1 2x
∫ y dy = ∫ x2 + 4x + 5 dx
2x + 4 − 4 2x + 4 4
ln( y ) = ∫ x2 + 4x + 5 dx = ∫ x2 + 4x + 5 dx − ∫ x2 + 4x + 5 dx =

(
= ln x2 + 4x + 5 − 4 ∫) 1
(x + 2)2 + 1
(
dx = ln x2 + 4x + 5 − 4arctg(x + 2) + C )
( )
∴ ln( y ) = ln x2 + 4x + 5 − 4arctg(x + 2) + C

( )
y = eC x2 + 4x + 5 e−4arctg( x +2)

( )
y = C x2 + 4 x + 5 e−4arctg( x +2)

Ejemplo: Integrar la EDO lineal homogénea:

y ′ − x2e2x y = 0

Siguiendo el procedimiento descrito, se tiene:


dy
= x2e2x y
dx

dy = x2e2x ydx

1
dy = x2e2xdx
y

1 2 2x
∫ y dy = ∫ x e dx

 u = x2 dv = e2xdx 
ln( y ) = ∫ x2e2xdx =  1 2x 
= 1 x2e2x
2
− ∫ xe2xdx =
du = 2x v = 2 e 

 u = x dv = e2xdx 
= = 1 x2e2x − 1 xe2x + ∫ 12 e2xdx =
1 2x  2 2
du = 1 v = 2 e 

= 1 x2e2x
2
− 1 xe2x
2
+ 1 e2x
4
+C = (
1 x2
2
− 1x
2
+ 1
4
)e 2x
+C

∴ ln( y ) = ( 1 x2
2
− 1x
2
+ 1
4
)e
2x
+C

y = eCe 2
(1 x2 − 12 x + 14 )e2x

y = Ce 2
(1 x2 − 12 x + 14 )e2x
Ejemplo: Resolver el problema de valor inicial:
dx
+
tx
dt sen2 (t )
= 0; x (2π ) = 2
88 Ecuaciones Diferenciales

Se trata de una ecuación lineal homogénea, y su solución general es:


dx tx
+ =0
dt sen2 (t )

dx −tx
=
dt sen2 (t )

−tx
dx = dt
sen2 (t )

1 −t
dx = dt
x sen2 (t )

1 −t
∫ x dx = ∫ sen2 (t ) dt

 −1 
−t  u=t dv = dt 
ln( x ) = ∫ dt =  sen 2
(t )  = tcotg(t ) − ∫ cotg(t )dt =
sen2 (t ) du = dt v = cotg(t ) 

cos(t )
= tcotg(t ) − ∫ dt = tcotg(t ) − ln(sen(t )) + C
sen(t )

∴ ln( x ) = tcotg(t ) − ln(sen(t )) + C

eCetcotg(t )
x =
sen(t )

Cetcotg(t )
x=
sen(t )

Sustituyendo ahora la condición x (2π ) = 2 en la solución general, se


obtiene una ecuación que permite determinar el valor de la constante de
integración C:

x (2π ) = 2
Ce 2
π cotg π
2
()
2=
sen 2π ()
Ce0
2=
1
C=2
Por tanto, la solución del problema de valor inicial es:

2etcotg(t )
x=
sen(t )

Observación. Las ecuaciones lineales completas no son, en general, ecuaciones de


variables separadas, por lo que se precisarán métodos de integración específicos para su
resolución. Sin embargo, las ecuaciones lineales completas de la forma
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 89

a1 (x) y ′ + a0 (x ) y = b(x ) , tales que a0 (x) = b(x ) , constituyen un caso particular de


ecuaciones lineales completas que, al igual que las homogéneas tienen integración
directa, pues se trata de ecuaciones de variables separadas. El siguiente ejemplo
presenta un problema clásico de caída de cuerpos cuya modelización se realiza
mediante una ecuación de este tipo.

Ejemplo: Caída libre de un cuerpo. Un cuerpo de masa 50 kg cae hacia la superficie


terrestre partiendo del reposo desde gran altura. El objeto es frenado por
la resistencia del aire con una fuerza Fr igual al doble de la velocidad de
descenso (figura 2.5). Hallar la velocidad y la distancia recorrida por el
cuerpo en función del tiempo. ¿A qué altitud comienza la caída el cuerpo
si llega al nivel del mar en 25 segundos?

Fr

Figura 2.2: Caída libre de un cuerpo.

Velocidad de caída del cuerpo: Sea v = f1(t ) la velocidad de caída del


cuerpo en un instante t. La fuerza F que produce el movimiento de caída
del cuerpo es igual a su masa por la aceleración de la gravedad, es decir,
F = (50)(9'8) = 490 N. Dado que la fuerza Fr que opone la resistencia del
aire es el doble de la velocidad de descenso, puede ponerse Fr = 2v . Por
tanto, de acuerdo con la segunda ley de Newton, se debe cumplir que:
ma = F − Fr

50a = 490 − 2v
donde a es la aceleración de la caída del cuerpo. De aquí, como a = v ′ ,
resulta la ecuación lineal completa de coeficientes constantes y término
independiente constante que determina la velocidad de descenso del
cuerpo:
50v ′ = 490 − 2v
además, dado que el cuerpo parte del reposo se trata de resolver el
problema de valor inicial:
50v ′ = 490 − 2v ; v (0) = 0
La solución general de la ecuación es:
dv
25 = 245 − v
dt
−1 −1 dt
dv = 25
245 − v
90 Ecuaciones Diferenciales

−1
∫ 245 − v dv = ∫ 25
−1 dt

(
ln 245 − v = ) −1 t
25
+C

−1 t
245 − v = eCe 25
−1 t
245 − v = Ce 25
−1 t
v = 245 − Ce 25
Sustituyendo ahora la condición v (0) = 0 en la solución general, se
obtiene el valor de la constante de integración C:
v (0) = 0
−1 ( 0 )
0 = 245 − Ce 25
0 = 245 − C
C = 245
Por tanto, la velocidad de caída del cuerpo en cada instante t es:

 −1 t 
v = 2451 − e 25 
 

Distancia recorrida por el cuerpo: Sea x = f2 (t ) la distancia recorrida por el


cuerpo en un instante t. Como v = x ′ , de la solución del problema de valor
inicial anterior resulta la ecuación de variables separadas que determina
la distancia recorrida por el cuerpo:

 −1 t 
x′ = 2451 − e 25 
 
y considerando que la distancia recorrida en el instante inicial es nula se
trata de resolver el problema de valor inicial:

 −1 t 
x′ = 2451 − e 25  ; x (0) = 0
 
La solución general de la ecuación es:

dx  −1 t 
= 2451 − e 25 
dt  

 −1 t 
dx = 2451 − e 25 dt
 

 −1 t 
∫ dx = ∫ 
2451 − e 25 dt

 

 −1 t 
x = 245 t + 25e 25  + C
 
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 91

Sustituyendo ahora la condición x (0) = 0 en la solución general, se


obtiene el valor de la constante de integración C:
x (0) = 0

 −1 ( 0 ) 
0 = 245 0 + 25e 25  + C
 
0 = 6125 + C
C = −6125
Por tanto, la distancia recorrida por el cuerpo en cada instante t es:

 −1 t 
x = 245 t + 25e 25  − 6125
 

Altitud inicial: Dado que el cuerpo tarda en llegar al nivel del mar 25
segundos, la altitud inicial coincide con la distancia total recorrida por el
cuerpo en 25 s, es decir:

( )
x = 245 25 + 25e−1 − 6125 ≅ 2253'26 m

2.6.3. Resolución de ecuaciones lineales completas

Solución general de la ecuación lineal completa

El siguiente teorema afirma que la solución general de una ecuación lineal completa
y ′ + G ( x) ⋅ y = H (x ) puede obtenerse sumando a la solución general de la correspondiente
ecuación lineal homogénea y ′ + G (x) ⋅ y = 0 una solución particular de la ecuación lineal
completa.

Teorema (Solución general de la ecuación lineal completa). Dada una ecuación lineal
completa y ′ + G (x) ⋅ y = H (x ) , si y = f (x, C) es la solución general de la correspondiente
ecuación lineal homogénea y ′ + G (x) ⋅ y = 0 , e y = g(x ) es una solución particular no
trivial de la ecuación lineal completa, entonces la solución general de la ecuación lineal
completa es:
y = f (x, C) + g(x)

Demostración:

Sea y = f (x, C) la solución general de la ecuación lineal homogénea, e y = g(x) una solución particular de la
ecuación lineal completa. Se trata de probar que la solución general de la ecuación lineal completa es
y = f (x, C) + g (x) . Derivando esta expresión, se obtiene:

y ′ = f ′(x, C) + g ′(x)

Si se llevan las expresiones de la solución general de la ecuación lineal completa y de su derivada a la ecuación
lineal completa y ′ + G (x ) ⋅ y = H (x ) , resulta:

f ′(x, C) + g ′(x) + G (x)(f (x, C) + g(x) ) = H (x )

f ′( x, C) + g ′(x ) + G (x )f (x, C) + G (x )g( x) = H (x )


92 Ecuaciones Diferenciales

(f ′(x,C) + G(x)f (x,C)) + (g ′(x) + G(x)g(x)) = H(x)


ahora bien, como y = f (x, C) es solución general de la ecuación lineal homogénea, se cumple que
f ′(x, C) + G (x)f (x,C) = 0 , y como y = g(x) es solución particular de la ecuación lineal completa, se cumple
también que g ′(x ) + G (x )g(x ) = H (x ) , obteniéndose finalmente:

0 + H (x ) = H (x)

H (x) = H (x )

es decir, la familia uniparamétrica de funciones y = f (x, C) + g (x) verifica la ecuación lineal completa, luego
constituye su solución general, como se quería probar.

Ejemplo: Integrar la siguiente ecuación lineal completa, sabiendo que una solución
particular es y = x 2 :

( )
xy ′ + x 2 − 1 y = x 4 + x 2

La solución de la ecuación lineal homogénea correspondiente es:

(
xy ′ + x 2 − 1 y = 0 )
x
dy
dx
( )
+ x2 − 1 y = 0

xdy = (1 − x )ydx 2

1 1 − x2
dy = dx
y x

1 1 − x2
∫ y dy = ∫ x
dx

2
1 
( ) ∫ 1 −xx
ln y = dx = ∫  − x dx = ln x − ( ) 1
2
x2 + C
x 

( )
∴ ln y = ln x − ( ) 1
2
x2 + C

− 1 x2
y = eC x e 2

− 1 x2
y = ±eC xe 2

− 1 x2
y = Cxe 2

Por tanto, la solución general de la ecuación lineal completa es:


− 1 x2
y = Cxe 2 + x2

Observación. El teorema anterior no es demasiado útil en la práctica para la resolución


de ecuaciones lineales completas, ya que encontrar una solución particular no trivial de
la ecuación es equivalente a encontrar la solución general de la misma. Se estudian a
continuación otros dos métodos prácticos para resolver estas ecuaciones: el método del
factor integrante (ya conocido), y el método de Lagrange o de variación de constantes.
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 93

Resolución por el método del factor integrante

Las ecuaciones lineales completas son, en general, ecuaciones no exactas, pero pueden
convertirse en exactas e integrarse mediante el método del factor integrante. Además, en
este caso el cálculo del factor integrante es sencillo, ya que siempre es un factor que
depende sólo de la variable independiente de la ecuación, lo que se prueba en la
siguiente proposición.

Proposición (Factor integrante de la ecuación lineal completa). La ecuación lineal


completa y ′ + G (x) ⋅ y = H (x ) admite un factor integrante µ dependiente sólo de la variable
x de la forma:

µ(x) = Ce ∫ G ( x )dx

Demostración:

En efecto, la ecuación lineal completa y ′ + G (x ) ⋅ y = H (x ) expresada en forma diferencial resulta:

dy
+ G (x ) ⋅ y = H (x )
dx

(G(x) ⋅ y − H(x))dx + dy = 0
luego, de acuerdo con la condición necesaria y suficiente para que una ecuación diferencial ordinaria admita
un factor integrante dependiente sólo de x, se tiene que:


(G(x) ⋅ y − H(x)) − ∂ (1)
∂y ∂x G (x ) − 0
= = G ( x ) = f (x )
1 1

es decir, toda ecuación lineal completa admite un factor integrante de la forma µ(x ) . Para hallar el factor
integrante, basta aplicar la fórmula:

µ(x ) = Ce ∫ f ( x )dx = Ce ∫ G (x )dx

Observación. Nótese que la expresión del factor integrante de una ecuación lineal
completa, µ(x) = Ce ∫ G ( x )dx , coincide con la inversa de la solución general de la ecuación
lineal homogénea correspondiente, y = Ce ∫ −G ( x )dx .

Procedimiento (Resolución de una EDO lineal completa mediante el método del factor
integrante). Para resolver una ecuación lineal completa y ′ + G (x) ⋅ y = H (x ) mediante un
factor integrante se procede de la siguiente forma:
1. Se calcula el factor integrante dependiente sólo de x asociado a la ecuación lineal
mediante la fórmula:

µ(x) = Ce ∫ G ( x )dx

2. Se multiplica la ecuación lineal original y ′ + G ( x) ⋅ y = H (x ) por el factor integrante


µ(x) , obteniéndose la EDO exacta:

µ(x )(G (x ) ⋅ y − H (x) )dx + µ(x)dy = 0


94 Ecuaciones Diferenciales

3. Se resuelve la EDO exacta, cuya solución general es equivalente a la de la


ecuación lineal original, mediante el procedimiento conocido de integración-
derivación-integración:

• F (x, y ) = ∫ µ(x)dy = µ(x) ⋅ y + ϕ(x)

∂F (x, y ) dϕ
• = µ ′(x) ⋅ y + = µ(x)(G (x) ⋅ y − H (x, y ) )
∂x dx

• ϕ(x) = ∫ (µ(x)(G (x) ⋅ y − H(x, y ) ) − µ ′(x) ⋅ y )dx

Ejemplo: Integrar la siguiente EDO lineal completa por el método del factor
integrante:
y
y′ + + x3 = 0
x
Expresando la ecuación en forma estándar, queda:
1
y′ + y = −x 3
x
1
donde G (x ) = , H (x) = −x 3 . Esta EDO lineal completa admite un factor
x
integrante dependiente de la variable x, de la forma:
1
∫ x dx ( )
ln x
µ( x ) = Ce ∫ G ( x )dx = Ce = Ce = C x = {C = 1} = x

Se multiplica la EDO lineal por el factor integrante µ(x ) = x , obteniéndose


la EDO exacta:

xy ′ + y = − x 4

dy
x + y = −x 4
dx

(x 4
)
+ y dx + xdy = 0

Se integra µ(x) = x parcialmente respecto a y:

F (x, y ) = ∫ xdy = xy + ϕ(x)

A continuación, se deriva F(x,y) parcialmente respecto a x, y se iguala a


µ(x)(G (x) ⋅ y − H (x ) ) = x 4 + y :

∂F (x, y )
=

(xy ) + dϕ = y + dϕ
∂x ∂x dx dx

y+ = x4 + y
dx

Por último, se despeja en la ecuación formada y se integra respecto a
dx
la variable x:

= x4
dx
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 95

ϕ′(x ) = x 4
4
∫ ϕ′(x)dx = ∫ x dx

ϕ(x) = 1
5
x5 + C

Por tanto, la solución general de la ecuación diferencial es:

xy + 1
5
x5 + C = 0

C
y = − 15 x 4 −
x

Ejemplo: Resolver el siguiente problema de valor inicial asociado a una EDO lineal
completa por el método del factor integrante:

(x 2
)
+ 1 y ′ + 4 xy = x ; y (0 ) = −1

Dividiendo la ecuación por x 2 + 1 , queda en forma estándar:


4xy x
y′ + =
2 2
x +1 x +1
4x x
donde G (x) = , H (x ) =
. Esta EDO lineal completa admite un
2 2
x +1 x +1
factor integrante dependiente de la variable x, de la forma:
4x
∫ x2 +1 dx ( 2 +1) = C(x 2 + 1)2 = {C = 1} = (x 2 + 1)2
µ( x ) = Ce ∫ G ( x )dx = Ce = Ce 2 ln x

Se multiplica la EDO lineal por el factor integrante µ(x) = x 2 + 1 , ( )


2

obteniéndose la EDO exacta:

(x 2
)
2
( ) (
+ 1 y ′ + 4x x 2 + 1 y = x x 2 + 1 )
(x + 2x + 1) dy
4
dx
+ (4x y + 4xy ) = x + x
2 3 3

(4x y + 4xy − x − x )dx + (x + 2x + 1)dy = 0


3 3 4 2

Se integra µ(x) = x 4 + 2x 2 + 1 parcialmente respecto a y:

∫ (x ) ( )
4
F (x, y ) = + 2x 2 + 1 dy = x 4 + 2x 2 + 1 y + ϕ(x )

A continuación, se deriva F(x,y) parcialmente respecto a x, y se iguala a


µ(x)(G (x) ⋅ y − H (x) ) = 4x 3 y + 4 xy − x 3 − x :

∂F (x, y )
∂x
=

∂x
((
x 4 + 2x 2 + 1 y +

dx
= 4x 3 + 4x y +

dx
)) ( )

4 x 3 y + 4 xy + = 4x 3 y + 4 xy − x 3 − x
dx

Por último, se despeja en la ecuación formada y se integra respecto a
dx
la variable x:
96 Ecuaciones Diferenciales


= −x 3 − x
dx

ϕ′(x ) = − x 3 − x

∫ ϕ′(x)dx = ∫ (− x )
3
− x dx

ϕ(x) = −1 x 4 − 1 x2 + C
4 2

La solución general de la ecuación diferencial es:

(x 4
)
+ 2x 2 + 1 y − 1
4
x4 − 1
2
x2 + C = 0

1 +C
1 4
y= −
4
(x 2
+1 )
2

Sustituyendo ahora la condición y (0) = −1 en la solución general, se


obtiene una ecuación que permite determinar el valor de la constante de
integración C:
y (0 ) = −1
1 +C
1 4
− = −1
4
((0) 2
+1 )
2

1 − 1 − C = −1
4 4

C =1
Por tanto, la solución del problema de valor inicial es:
5
1 4
y= −
4
(x 2
+1 )
2

Resolución por el método de Lagrange o de variación de constantes

Este método, de nombre cuando menos curioso (variación de constantes), fue ideado por
Joseph Louis Lagrange partiendo de una hipótesis en cierta forma similar a la planteada
en el teorema de la solución general de la ecuación lineal completa. La hipótesis consiste
en suponer que la ecuación lineal completa y ′ + G ( x) ⋅ y = H (x ) tiene una solución general
y = f ( x,C) de la misma forma que la solución general de la correspondiente ecuación
lineal homogénea y ′ + G (x ) ⋅ y = 0 , es decir, y = f (x,C) = Ce ∫ −G ( x )dx , pero tomando la
constante de integración C como una función que depende de la variable x.
Desarrollando esta idea, la proposición siguiente establece la expresión de la solución
general de la ecuación lineal completa.
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 97

Proposición (Fórmula de variación de constantes para la ecuación lineal completa).


Dada una ecuación lineal completa y ′ + G (x) ⋅ y = H (x ) , con ecuación lineal homogénea
asociada y ′ + G (x) ⋅ y = 0 de solución general y = Ce ∫ −G ( x )dx , la solución general de la
ecuación lineal completa es:

y = C(x )e ∫ −G ( x )dx = (∫ H(x)e )


∫ G ( x )dxdx e ∫ −G ( x )dx

Demostración:

Se trata de probar que la solución general de la ecuación lineal completa y ′ + G ( x ) ⋅ y = H (x ) es igual a la

solución general y = Ce ∫ −G ( x )dx de la correspondiente ecuación lineal homogénea, pero tomando la constante
de integración C como una función que depende de la variable x, es decir, una expresión de la forma
y = C (x )e ∫ −G ( x )dx . Derivando esta expresión, se obtiene:

y ′ = C ′(x )e ∫ −G (x )dx − C (x )G ( x )e ∫ −G ( x )dx

Si se llevan las expresiones de la solución general de la ecuación lineal completa y de su derivada a la ecuación
lineal completa y ′ + G (x ) ⋅ y = H (x ) , resulta:

C′(x )e ∫ −G ( x )dx − C (x )G ( x )e ∫ −G ( x )dx  + G (x ) C(x )e ∫ −G ( x )dx  = H (x )


   

C′(x )e ∫ −G (x )dx − C (x )G (x )e ∫ −G ( x )dx + C(x )G (x )e ∫ −G ( x )dx = H (x )

C′(x )e ∫ −G ( x )dx = H (x )

H (x)
C′(x ) =
e ∫ −G ( x )dx

C′(x ) = H (x )e ∫ G ( x )dx

e integrando esta igualdad, se obtiene:

C (x ) = ∫ H (x )e ∫ G ( x )dxdx

En conclusión, la solución general de la ecuación lineal completa es:

y = C(x )e ∫ −G ( x )dx =  ∫ H (x )e ∫ G ( x )dxdx e ∫ −G ( x )dx


 

Procedimiento (Resolución de una EDO lineal completa mediante el método de


variación de constantes). Para resolver una ecuación lineal completa y ′ + G ( x) ⋅ y = H (x )
mediante variación de constantes se procede de la siguiente forma:
1. Se calcula la solución general de la ecuación lineal homogénea asociada
y ′ + G (x ) ⋅ y = 0 :

y = Ce ∫ −G ( x )dx

2. En la solución general de la ecuación lineal homogénea se toma la constante de


integración como una función dependiente de x, y se obtiene la expresión de su
derivada:

y = C(x )e ∫ −G ( x )dx

y ′ = C ′(x )e ∫ −G ( x )dx − C(x)G (x)e ∫ −G ( x )dx


98 Ecuaciones Diferenciales

3. Se llevan ambas expresiones a la ecuación lineal completa y se despeja C'(x):

(C′(x)e ∫ −G ( x )dx ) (
− C (x )G ( x)e ∫ −G ( x )dx + G (x ) C (x )e ∫ −G ( x )dx = H (x ) )
C ′(x ) = H (x )e ∫ G ( x )dx
4. Se integra C'(x), y se lleva el factor resultante C(x) a la solución general de la
ecuación lineal:

y = C(x )e ∫ −G ( x )dx = (∫ H(x)e )


∫ G ( x )dxdx e ∫ −G ( x )dx

Ejemplo: Integrar la siguiente EDO lineal completa por el método de variación de


constantes:
y ′ + 2xy = 4x

En primer lugar, se resuelve la ecuación lineal homogénea asociada:


y ′ + 2xy = 0

dy
+ 2xy = 0
dx
dy = −2xydx

1
dy = −2xdx
y

1
∫ y dy = ∫ − 2xdx

ln( y ) = − x2 + C
2
y = eCe − x
2
y = Ce− x

Se deriva esta solución general tomando la constante de integración C


como una función dependiente de x:
2
y = C(x )e− x
2 2
y ′ = C′(x )e− x − 2xC (x )e− x

Se llevan ambas expresiones a la ecuación lineal completa y se despeja


C'(x):
2 2 2
C′( x)e− x − 2xC ( x)e− x + 2xC ( x)e− x = 4x
2
C ′( x ) = 4 xe x
Se integra C'(x), y se lleva el resultado a la solución general de la ecuación
lineal:
x2
∫ C′(x)dx = ∫ 4xe dx
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 99

2
C (x ) = 2e x + C

 2  2 2
y =  2e x + C e− x = 2 + Ce − x
 

Ejemplo: Intensidad de la corriente en un circuito eléctrico serie RL. Sea un circuito


eléctrico RL conectado a una fuente de alimentación senoidal (figura 2.6).
La fuente de alimentación proporciona una fuerza electromotriz variable
E (t ) = 12sen(500t ) , y los valores de resistencia e inducción del circuito
son R = 10 Ω y L = 20 mH. Hallar la intensidad de la corriente i del
circuito en función del tiempo.

E
−+

L R

Figura 2.3: Circuito eléctrico serie RL.

Sea i = f (t ) la intensidad de corriente del circuito en un instante t. La


corriente provoca caídas de voltaje en los elementos de resistencia e
inducción al circular a través de ellos, de magnitudes ER = 10 ⋅ i ,
EL = 0'02 ⋅ i′ respectivamente. De acuerdo con la ley de Kirchhoff, se debe
cumplir que:
ER + EL = E

de donde, sustituyendo ER, EL y E por sus respectivas expresiones, se


obtiene la ecuación lineal completa que determina la intensidad de
corriente en función del tiempo:

10 ⋅ i + 1 ⋅ i′ = 12sen(500t )
50

di
+ 500 ⋅ i = 600sen(500t )
dt
En primer lugar, se resuelve la ecuación lineal homogénea asociada:
di
+ 500 ⋅ i = 0
dt
1
di = −500dt
i
1
∫ i di = ∫ − 500dt
ln( i ) = −500t + C

i = eCe−500t
100 Ecuaciones Diferenciales

i = Ce−500t
Se deriva esta solución general tomando la constante de integración C
como una función dependiente de t:

i = C (t )e −500t

di
= C′(t )e −500t − 500C (t )e−500t
dt
Se llevan ambas expresiones a la ecuación lineal completa y se despeja
C'(x):

C′(t )e−500t − 500C(t )e−500t + 500C(t )e−500t = 600sen(500t )

C′(t ) = 600sen(500t )e500t


Se integra C'(x), y se lleva el resultado a la solución general de la ecuación
lineal:
500t
∫ C′(t )dt = ∫ 600sen(500t )e dt

 u = sen(500t ) dv = 600e500t dt 
C(t ) = ∫ 600sen(500t )e500t dt =  =
du = 500 cos ( 500t) v = 65 e500t 

= 6 sen(500t )e500t − ∫ 600 cos(500t )e500t dt =


5

 u = cos(500t ) dv = 600e500t dt 
= =
du = −500sen( 500t) v = 65 e500t 

= 6 sen(500t )e500t − 6 cos(500t )e500t − ∫ 600sen(500t )e500t dt


5 5

∴ C (t ) = ∫ 600sen(500t )e500t dt = 3
5
e500t (sen(500t ) − cos(500t )) + C

i= 3 (sen(500t ) − cos(500t )) + Ce −500t


5

2.6.4. Ecuación de Bernoulli

Definición (Ecuación de Bernoulli). Una ecuación diferencial ordinaria de primer orden


es una ecuación de Bernoulli si y sólo si se expresa en la forma:

y ′ + G (x ) ⋅ y = H (x) ⋅ y p (p ∈ o)

donde G y H son expresiones analíticas elementales.

Ejemplo: La ecuación diferencial ordinaria de incógnita y (x ) :

y ′ + y = xy 4

es una ecuación de Bernoulli con G (x) = 1 , H ( x) = x , p = 4 .


2. Ecuaciones Diferenciales Ordinarias de Primer Orden 101

Ejemplo: La ecuación diferencial ordinaria de incógnita y (x ) :

y ′ + sen(x ) y = (x + 1) y 3

es una ecuación de Bernoulli con G (x ) = sen(x ) , H (x ) = x + 1 , p = 3 .


2

Ejemplo: La ecuación diferencial ordinaria de incógnita x (t ) :

dx x
+ = tx 2
dt t
1
es una ecuación de Bernoulli con G (t ) = , H (t ) = t , p = 2 .
t

Observación. La ecuación lineal puede considerarse un caso particular de la ecuación de


Bernoulli, ya que:
• Si p = 0 , la ecuación de Bernoulli es lineal completa:

y ′ + G (x ) ⋅ y = H (x) ⋅ y 0 ⇔ y ′ + G (x ) ⋅ y = H (x)

• Si p = 1 , la ecuación de Bernoulli es lineal homogénea:

y ′ + G (x ) ⋅ y = H (x) ⋅ y 1 ⇔ y ′ + (G (x) − H (x)) ⋅ y = 0

Proposición (Cambio de variable en la ecuación de Bernoulli). Dada una ecuación de


Bernoulli y ′ + G ( x) ⋅ y = H ( x) ⋅ y p , siendo p ∈ o − {0,1} , con el cambio de variable u = y1− p
se transforma en una ecuación lineal completa de la forma:
 u = y 1− p 
y ′ + G (x) ⋅ y = H (x ) ⋅ y p ⇒  − p  ⇒ u′ + (1 − p)G (x) ⋅ u = (1 − p)H (x )
du = (1 − p) y dy 

Demostración:

En efecto, sea la ecuación de Bernoulli:

y′ + G (x) ⋅ y = H (x) ⋅ y p

dy
+ G (x) ⋅ y = H (x ) ⋅ y p
dx

dy =  H (x ) ⋅ y p − G (x) ⋅ y dx
 

(1 − p) y − pdy = (1 − p) y − p  H (x ) ⋅ y p − G (x ) ⋅ y dx
 

(1 − p) y − p dy =  (1 − p)H (x ) − (1 − p)G (x) ⋅ y 1− p dx


 

 u = y1− p 
Si a esta ecuación se le aplica el cambio   , resulta:
du = (1 − p) y − pdy 

du = ((1 − p)H (x ) − (1 − p)G (x ) ⋅ u)dx

du
+ (1 − p)G (x ) ⋅ u = (1 − p)H (x)
dx
102 Ecuaciones Diferenciales

u′ + (1 − p)G (x ) ⋅ u = (1 − p)H (x )

que es una ecuación lineal completa, como se quería demostrar.

Procedimiento (Resolución de la ecuación de Bernoulli). Para resolver una ecuación de


Bernoulli y ′ + G (x ) ⋅ y = H (x) ⋅ y p , siendo p ∈ o − {0,1} , se procede de la siguiente manera:

1. Se aplica el cambio de variable u = y1− p , obteniéndose una ecuación lineal


completa:
 u = y 1− p 
y ′ + G (x ) ⋅ y = H (x) ⋅ y p ⇒  − p  ⇒ u′ + (1 − p)G (x ) ⋅ u = (1 − p)H (x )
du = (1 − p) y dy 

2. Se resuelve la ecuación lineal completa empleando el método del factor


integrante o el método de variación de constantes.

3. Se deshace el cambio en la solución obtenida, sustituyendo u = y1− p .

Ejemplo: Hallar la solución general de la ecuación de Bernoulli:

y − y2
y′ − =
x x3
−1 −1
Se trata de una EDO de Bernoulli con G (x ) = , H (x ) = , p = 2.
x x3
 u = y −1 
Aplicando el cambio   , resulta la ecuación lineal completa:
du = − y − 2dy 

u 1
u′ + =
x x3

Esta EDO lineal completa admite un factor integrante dependiente de la


variable x, de la forma:
1
∫ dx ( )
µ(x) = Ce ∫
− G ( x )dx ln x
= Ce x = Ce = C x = {C = 1} = x

Se multiplica la EDO lineal por el factor integrante µ(x ) = x , obteniéndose


la EDO exacta:
1
xu′ + u =
x2
du 1
x +u=
dx x2
 1 
 u − dx + xdu = 0
 x2 
Se integra µ(x ) = x parcialmente respecto a u:

F (x, u) = ∫ xdu = xu + ϕ(x)


2. Ecuaciones Diferenciales Ordinarias de Primer Orden 103

A continuación, se deriva F(x,u) parcialmente respecto a x, y se iguala a


1
µ(x)(− G (x) ⋅ u + H (x ) ) = u − :
x2
∂F (x, u)
=

(xu) + dϕ = u + dϕ
∂x ∂x dx dx
dϕ 1
u+ =u−
dx x2

Se despeja en la ecuación formada y se integra respecto a la variable
dx
x:
dϕ −1
=
dx x 2

−1
ϕ′(x ) =
x2
−1
∫ ϕ′(x)dx = ∫ x 2 dx
1
ϕ( x) = +C
x
La solución general de la ecuación lineal completa es:
1
xu + +C = 0
x

Finalmente, se deshace el cambio de variable sustituyendo u = y −1 , y se


obtiene la solución general de la ecuación de Bernoulli:
x 1
+ +C =0
y x

x
y=
−1
−C
x

Ejemplo: Hallar la solución general de la ecuación de Bernoulli:

y ′ + y = xy 3

Se trata de una EDO de Bernoulli con G (x ) = 1 , H ( x ) = x , p = 3 . Aplicando


 u = y −2 
el cambio   , resulta la ecuación lineal completa:
du = −2 y − 3dy 

u′ − 2u = −2x
En este caso, se aplica el método de variación de constantes. Primero, se
resuelve la ecuación lineal homogénea asociada:
u′ − 2u = 0
du
− 2u = 0
dx
104 Ecuaciones Diferenciales

du = 2udx
1
du = 2dx
u
1
∫ u du = ∫ 2dx
( )
ln u = 2x + C

u = eC e 2x

u = Ce 2x
Tomando la constante de integración C como una función dependiente de
x, se deriva esta solución general:

u = C(x )e 2x

u′ = C ′(x )e 2x + 2C(x )e 2x

Se llevan ambas expresiones a la ecuación lineal completa y se despeja


C'(x):

C ′(x )e 2x + 2C(x )e 2x − 2C(x )e 2x = −2x

C ′(x ) = −2xe −2x

Se integra C'(x), y se lleva el resultado a la solución general de la ecuación


lineal:
−2x
∫ C ′(x)dx = ∫ − 2xe dx

 u = x dv = −2e −2x dv 
C(x ) = ∫ − 2xe − 2x dx =  =
du = dx e − 2x 

= xe −2x − ∫ e −2x dx = xe −2x − 12 e −2x + C

(
u = xe −2x − 12 e −2x + C e 2x )
u = x − 12 + Ce 2x

Finalmente, se deshace el cambio de variable sustituyendo u = y −2 , y se


obtiene la solución general de la ecuación de Bernoulli:
1
=x− 1
2
+ Ce2x
2
y
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 105

2.6.5. Ecuación de Riccati

Definición (Ecuación de Riccati). Una ecuación diferencial ordinaria de primer orden es


una ecuación de Riccati si y sólo si se expresa en la forma:

y ′ = a( x ) y 2 + b(x ) y + c(x )

donde a, b, y c son expresiones analíticas elementales.

Ejemplo: La ecuación diferencial ordinaria de incógnita y (x ) :

y′ = y2 + y + 3

es una ecuación de Riccati con a(x ) = 1 , b(x ) = 1 , c(x ) = 3 .

Ejemplo: La ecuación diferencial ordinaria de incógnita y (x ) :

y ′ = xy 2 + x 2 y + 3x 3

es una ecuación de Riccati con a(x ) = x , b(x ) = x 2 , c(x ) = 3x .

Ejemplo: La ecuación diferencial ordinaria de incógnita x (t ) :

dx
= t ⋅ x 2 + (t − 3)x + t 3
dt

es una ecuación de Riccati con a(t ) = t , b(t ) = t − 3 , c(t ) = t 3 .

Observación. La ecuación de Riccati puede considerarse una generalización de las


ecuaciones lineal y de Bernoulli, ya que:
• Si a( x ) = 0 , la ecuación de Riccati es una ecuación lineal completa:

y ′ = 0 ⋅ y 2 + b(x ) y + c(x ) ⇔ y ′ − b(x ) y = c(x )

• Si c(x ) = 0 , la ecuación de Riccati es una ecuación de Bernoulli con p = 2 :

y ′ = a(x ) y 2 + b(x ) y + 0 ⇔ y ′ − b(x ) y = a( x ) y 2

Observación. La ecuación de Riccati sólo puede resolverse si se conoce una solución


particular, en cuyo caso se aplicará un cambio de variable específico que transforma la
ecuación de Riccati en una EDO lineal completa.

Proposición (Cambio de variable en la ecuación de Riccati). Dada una ecuación de


Riccati y ′ = a( x ) y 2 + b(x ) y + c(x ) , si y = g (x ) es una solución particular no trivial de la
1
ecuación, entonces con el cambio de variable y = g (x ) + se transforma en una
u
ecuación lineal completa de la forma:

 1 
 y = g(x) +
 u 
y ′ = a(x ) y 2 + b(x ) y + c(x ) ⇒   ⇒ u′ + (2a(x )g(x ) + b(x ) )u = −a(x )
1
dy = g ′(x )dx − du
 u 2 
106 Ecuaciones Diferenciales

Demostración:

En efecto, sea la ecuación de Riccati:

y′ = a( x) y 2 + b( x) y + c(x )

 1 
 y = g(x ) + 
u
Si a esta ecuación se le aplica el cambio  1  , donde y = g(x) es una solución particular de
dy = g ′(x)dx − du
 u 2 
la misma, resulta:

2
u′  1  1
g ′(x ) − = a(x ) g(x ) +  + b(x) g(x ) +  + c(x )
2  u  u
u

u′  1 1   1
g ′(x ) − = a(x ) g 2 (x ) + 2g(x ) +  + b(x ) g(x ) +  + c(x )
2 u 2  u
u  u 

u′  1 1  1
g ′(x ) − = a(x)g 2 (x ) + b( x)g( x) + c(x ) + a(x ) 2g(x ) +  + b(x )
u2  u u2  u

y como y = g(x) es solución particular de ecuación se cumple que g ′(x ) = a(x)g 2 (x ) + b(x)g (x) + c(x ) , luego:

u′  1 1  1
− = a(x ) 2g(x ) +  + b(x ) 
u2  u u2  u

y multiplicando la ecuación por − u2 :

u′ = a(x)(−2g (x )u − 1) − b(x)u

u′ + (2a(x )g(x ) + b( x) )u = −a(x)

que es una ecuación lineal completa, como se quería demostrar.

Procedimiento (Resolución de la ecuación de Riccati). Para resolver una ecuación de


Riccati y ′ = a( x ) y 2 + b(x ) y + c(x ) , siendo y = g (x ) una solución particular no trivial de la
ecuación, se procede de la siguiente manera:
1
1. Se aplica el cambio de variable y = g (x ) + , obteniéndose una ecuación lineal
u
completa:

 1 
 y = g(x ) + 
u
 ⇒ u′ + (2a( x )g( x ) + b(x ) )u = −a(x )
2
y ′ = a(x ) y + b(x ) y + c(x ) ⇒ 
1
dy = g ′(x )dx − du
 u2 

2. Se resuelve la ecuación lineal completa empleando el método del factor


integrante o el método de variación de constantes.
1
3. Se deshace el cambio en la solución obtenida, sustituyendo u = .
y − g( x)

Ejemplo: Hallar la solución general de la ecuación de Riccati:

y ′ = y 2 − 2xy + 1 + x 2

sabiendo que posee una solución particular de la forma y = mx .


2. Ecuaciones Diferenciales Ordinarias de Primer Orden 107

Se trata de una EDO de Riccati con a(x ) = 1 , b(x ) = −2x , c( x ) = 1 + x 2 . En


primer lugar, se debe obtener la solución particular, para lo cual se deriva
y = mx , y se llevan y e y' a la ecuación para así despejar m:

y′ = m

m = m2 x 2 − 2mx 2 + 1 + x 2

( )
m2 x 2 + m − 2x 2 − 1 + 1 + x 2 = 0

m=
2x 2 + 1 ± (− 2x − 1)
2 2
(
− 4x 2 1 + x 2 )=
2
2x
 1
2x 2 + 1 ± 4 x 4 + 4x 2 + 1 − 4x 2 − 4x 4 2x 2 + 1 ± 1 1 +
= = = x2
2x 2 2x 2 1

Por tanto, m = 1 y la solución particular de la ecuación de Riccati es


 1 
y=x+
 u 
y = x . A continuación, se aplica el cambio de variable  ,
1
dy = dx − du
 u 2 
obteniéndose la EDO lineal:
u′ + (2x − 2x)u = −1

u ′ = −1
Se trata de una EDO de variables separadas, de integración directa:
du
= −1
dx
du = −dx

∫ du = ∫ − dx
u = −x + C
Finalmente, se deshace el cambio en la solución obtenida, sustituyendo
1
u= :
y−x

1
= −x + C
y−x

1
y=x+
− x +C

Ejemplo: Hallar la solución general de la ecuación de Riccati:

(1 + x )y ′ + 2xy
3 2
+ x2y + 1 = 0

sabiendo que posee una solución particular de la forma y = mx .

Despejando y' en la ecuación, se obtiene:


108 Ecuaciones Diferenciales

− 2x x2 1
y′ = y2 − y−
3 3
1+ x 1+ x 1 + x3

−2x − x2 −1
que es una EDO de Riccati con a(x ) = , b(x ) =
. , c(x) =
3 3
1+ x 1+ x 1 + x3
En primer lugar, se debe obtener la solución particular, para lo cual se
deriva y = mx , y se llevan y e y' a la ecuación para así despejar m:

y′ = m

− 2x 3 x3 1
m= m2 − m−
3 3
1+ x 1+ x 1 + x3

2x 3 2x 3 + 1 1
m2 + m+ =0
3 3
1+ x 1+ x 1+ x3

(
2x 3 m2 + 2x 3 + 1 m + 1 = 0 )
m=
− 2x 3 − 1 ± (2x 3
+1 )2
− 8x 3
=
− 2x 3 − 1 ± 4x 6 − 4x 3 + 1
=
4x 3 4x 3

=
− 2x 3 − 1 ± (2x 3
−1 )2
=
− 2x 3 − 1 ± 2x 3 − 1 ( ) = 2−x1 3
4x 3 4x 3 − 1

Por tanto, m = −1 y la solución particular de la ecuación de Riccati es


y = −x . A continuación, se aplica el cambio de variable
 1 
 y = −x + u 
 1  , obteniéndose la EDO lineal:
dy = −dx − du
 u 2 

 4x 2 x2 
u′ +  − u = 2x
1 + x3 1 + x3  1 + x3
 

3x 2 2x
u′ + u=
3
1+ x 1 + x3
En este caso, se aplica el método de variación de constantes. Primero, se
resuelve la ecuación lineal homogénea asociada:

3x 2
u′ + u=0
1 + x3

du 3x 2
+ u=0
dx 1 + x 3

1 − 3x 2
du = dx
u 1 + x3

1 − 3x 2
∫u du = ∫ 1 + x 3 dx
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 109

ln u = − ln 1 + x 3  + C
( )  

eC
u =
1 + x3

C
u=
1 + x3
Tomando la constante de integración C como una función dependiente de
x, se deriva esta solución general:
C( x )
u=
1 + x3

C ′(x) 3x 2C(x )
u′ = −
1 + x3 (1 + x ) 3 2

Se llevan ambas expresiones a la ecuación lineal completa y se despeja


C'(x):

C ′(x ) 3x 2C(x) 3x 2C(x) 2x


− + =
1+ x 3
(1 + x ) (1 + x )
3 2 3 2 1 + x3

C ′(x) = 2x

Se integra C'(x), y se lleva el resultado a la solución general de la ecuación


lineal:

∫ C ′(x)dx = ∫ 2xdx
C( x ) = x 2 + C

x2 + C
u=
1 + x3
1
Finalmente, se deshace el cambio de variable sustituyendo u = , y se
y+x
obtiene la solución general de la ecuación de Riccati:

1 x2 + C
=
y + x 1 + x3

1 + x3
y= −x
x2 + C
110 Ecuaciones Diferenciales

2.7. ECUACIÓN DE LAGRANGE

2.7.1. Definición de la ecuación de Lagrange

Definición (Ecuación de Lagrange). Una ecuación diferencial ordinaria de primer orden


es una ecuación de Lagrange si y sólo si se expresa en la forma:
y = xf ( y ′) + g ( y ′)

donde f y g son expresiones analíticas elementales.

Ejemplo: La ecuación diferencial ordinaria de incógnita y (x ) :

y = xtg(2 y ′) + cos( y ′ + π)

es una ecuación de Lagrange con f ( y ′) = tg(2 y ′) , g ( y ′) = cos( y ′ + π) .

Ejemplo: La ecuación diferencial ordinaria de incógnita x (t ) :

x′
x = tx ′ 2 +
x′ + 1
x′
es una ecuación de Lagrange con f (x ′) = x ′2 , g (x′) = .
x′ + 1

Ejemplo: La ecuación diferencial ordinaria de incógnita y (x ) :

y = xy ′2 + xy ′ + 3 y ′

es una ecuación de Lagrange con f ( y ′) = y ′2 + y ′ , g ( y ′) = 3 y ′ .

2.7.2. Resolución de la ecuación de Lagrange

La resolución de la ecuación de Lagrange y = xf ( y ′) + g ( y ′) se realiza mediante el cambio


de variable y ′ = t , que permite transformarla en una ecuación lineal completa de variable
dependiente x y variable independiente t. Integrando esta ecuación lineal es posible
expresar la solución general de la ecuación de Lagrange en forma de ecuaciones
paramétricas x = h(t,C ) , y = h(t,C)f (t ) + g (t ) .

Proposición (Cambio de variable en la ecuación de Lagrange). Dada una ecuación de


Lagrange y = xf ( y ′) + g ( y ′) , con el cambio de variable y ′ = t se transforma en una
ecuación lineal completa de variable dependiente x y variable independiente t de la
forma:
dx f ′(t ) g ′(t )
y = xf ( y ′) + g ( y ′) ⇒ {y ′ = t } ⇒ − x=
dt t − f (t ) t − f (t )

Demostración:

En efecto, sea la ecuación de Lagrange:


y = xf ( y ′) + g ( y ′)
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 111

Si a esta ecuación se le aplica el cambio {y ′ = t } , resulta:

y = xf (t ) + g (t )

y derivando esta igualdad respecto a x se elimina la variable y de la ecuación:

dt dt
y ′ = t = f (t ) + xf ′(t ) + g ′(t )
dx dx

dt
t = f (t ) + (xf ′(t ) + g ′(t ) )
dx

dt
de donde despejando , se obtiene:
dx

dt t − f (t )
=
dx xf ′(t ) + g ′(t )

e invirtiendo esta ecuación, lo que equivale a tomar x como variable dependiente y t como variable
independiente, se obtiene:

dx xf ′(t ) + g ′(t )
=
dt t − f (t )

dx f ′(t ) g ′(t )
− x=
dt t − f (t ) t − f (t )

que es una ecuación lineal completa, como se quería demostrar.

Procedimiento (Resolución de la ecuación de Lagrange). Para resolver una ecuación de


Lagrange y = xf ( y ′) + g ( y ′) , se procede de la siguiente manera:

1. Se aplica a la ecuación de Lagrange el cambio de variable y ′ = t , y se deriva


respecto a x la ecuación resultante, obteniéndose una ecuación lineal completa
de variable dependiente x y variable independiente t:
dx f ′(t ) g ′(t )
y = xf ( y ′) + g ( y ′) ⇒ {y ′ = t } ⇒ − x=
dt t − f (t ) t − f (t )
2. Se resuelve la ecuación lineal completa empleando el método del factor
integrante o el método de variación de constantes, obteniéndose una solución
general de la forma x = h(t,C ) .

3. Se lleva la expresión x = h(t,C ) a la ecuación inicial y = xf (t ) + g (t ) , con lo que se


obtienen las ecuaciones paramétricas solución general de la ecuación de
Lagrange:
x = h(t,C)

 y = h(t,C)f (t ) + g (t )
4. Finalmente, en caso de que h sea una función inyectiva y en la ecuación
x = h(t,C ) pueda despejarse t en función de x, será posible eliminar el parámetro
t de las ecuaciones y obtener la expresión explícita de la solución general de la
ecuación de Lagrange:

( ) (
y = xf h−1( x,C) + g h−1(x,C) )
Ejemplo: Hallar la solución general de la ecuación de Lagrange:

y = xy ′2 + y ′ + 5
112 Ecuaciones Diferenciales

Se trata de una EDO de Lagrange con f ( y ′) = y ′2 , g ( y ′) = y ′ + 5 . En primer


lugar, se aplica el cambio de variable y ′ = t , y se deriva respecto a x la
ecuación resultante, con lo que se obtiene una EDO lineal completa:

y = xt 2 + t + 5

dt dt
y ′ = t = t 2 + 2xt +
dx dx
dt
t = t 2 + (2xt + 1)
dx

t − t2 dt
=
2xt + 1 dx
dx 2xt + 1
=
dt t − t2
dx 2 1
− x=
dt 1 − t t − t2
En este caso, se aplica el método de variación de constantes. Primero, se
resuelve la ecuación lineal homogénea asociada:
dx 2
− x=0
dt 1 − t
1 2
dx = dt
x 1− t
1 2
∫ x dx = ∫ 1 − t dt
ln( x ) = −2ln(1 − t ) + C

eC
x =
(1 − t )2

C
x=
(1 − t )2

Tomando la constante de integración C como una función dependiente de


t, se deriva esta solución general:
C (t )
x=
(1 − t )2

dx C′(t ) 2C(t )
= 2
+
dt (1 − t ) (1 − t )3

Se llevan ambas expresiones a la ecuación lineal completa y se despeja


C'(t):
C′(t ) 2C (t ) 2C(t ) 1
2
+ 3
− 3
=
(1 − t ) (1 − t ) (1 − t ) t − t2
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 113

1− t
C′(t ) =
t
Se integra C'(t), y se lleva el resultado a la solución general de la ecuación
lineal:
1− t
∫ C′(t )dt = ∫ t
dt

C(t ) = ln( t ) − t + C

x=
( )
ln t − t + C
(1 − t )2

Se lleva la expresión de x a la ecuación inicial y = xt 2 + t + 5 , con lo que


se obtienen las ecuaciones paramétricas solución general de la ecuación
de Lagrange:
 ln( t ) − t + C
x =
 (1 − t )2

 y = t ln( t ) − t + Ct + t + 5
2 3 2

 (1 − t )2

Ejemplo: Hallar la solución general de la ecuación de Lagrange:

y = x (1 + y ′) + y ′2

Se trata de una EDO de Lagrange con f ( y ′) = 1 + y ′ , g( y ′) = y ′2 . En primer


lugar, se aplica el cambio de variable y ′ = t , y se deriva respecto a x la
ecuación resultante, con lo que se obtiene una EDO lineal completa:

y = x (1 + t ) + t 2

dt dt
y′ = t = 1 + t + x + 2t
dx dx
dt
0 = 1 + (x + 2t )
dx
−1 dt
=
x + 2t dx
dx
= −x − 2t
dt
dx
+ x = −2t
dt
En este caso, se aplica el método de variación de constantes. Primero, se
resuelve la ecuación lineal homogénea asociada:
dx
+x=0
dt
1
dx = −dt
x
114 Ecuaciones Diferenciales

1
∫ x dx = ∫ − dt
( )
ln x = −t + C

x = eC e −t

x = Ce −t
Tomando la constante de integración C como una función dependiente de
t, se deriva esta solución general:

x = C(t )e −t

dx
= C ′(t )e − t − C(t )e − t
dt
Se llevan ambas expresiones a la ecuación lineal completa y se despeja
C'(t):

C ′(t )e −t − C(t )e −t + C(t )e −t = −2t

C ′(t ) = −2te t

Se integra C'(t), y se lleva el resultado a la solución general de la ecuación


lineal:
t
∫ C ′(t )dt = ∫ − 2te dt
 u = −2t dv = e t dt 
C(t ) = ∫ − 2te t dt =  t t
 = −2te + ∫ 2e dt =
du = −2dt v = e t 

= −2te t + 2e t + C = (−2t + 2)e t + C

(
x = (−2t + 2)e t + C e −t )
x = −2t + 2 + Ce −t

Se lleva la expresión de x a la ecuación inicial y = x(1 + t ) + t 2 , con lo que


se obtienen las ecuaciones paramétricas solución general de la ecuación
de Lagrange:
x = −2t + 2 + Ce −t

( )
 y = − 2t + 2 + Ce − t (1 + t ) + t 2
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 115

2.8. ECUACIÓN DE CLAIRAUT

2.8.1. Definición de la ecuación de Clairaut

Definición (Ecuación de Clairaut). Una ecuación diferencial ordinaria de primer orden es


una ecuación de Clairaut si y sólo si se expresa en la forma:
y = xy ′ + g (y ′)

donde g es una expresión analítica elemental.

Ejemplo: La ecuación diferencial ordinaria de incógnita y (x ) :

y = xy ′ + arctg( y ′)

es una ecuación de Clairaut con g ( y ′) = arctg( y ′) .

Ejemplo: La ecuación diferencial ordinaria de incógnita y (x ) :

y = xy ′ − 2 y ′

es una ecuación de Clairaut con g ( y ′) = −2 y ′ .

Ejemplo: La ecuación diferencial ordinaria de incógnita x(t ) :

1
x = tx ′ +
x′
1
es una ecuación de Clairaut con g (x′) = .
x′

Observación. La ecuación de Clairaut no es más que un caso particular de la ecuación


de Lagrange. En efecto, si en la ecuación de Lagrange y = xf ( y ′) + g ( y ′) se hace f ( y ′) = y ′
se obtiene la ecuación de Clairaut y = xy ′ + g (y ′) .

2.8.2. Resolución de la ecuación de Clairaut

La resolución de la ecuación de Clairaut y = xy ′ + g (y ′) se realiza mediante el mismo


cambio de variable y ′ = t de la ecuación de Lagrange, con lo que se transforma en una
ecuación diferencial de variables separadas. A partir de esta ecuación es posible
expresar la solución general de la ecuación de Clairaut como una familia de rectas
y = Ct + g (C) , más una solución singular en forma paramétrica x = − g ′(t ) ,
y = −g ′(t ) ⋅ t + g (t ) .

Proposición (Cambio de variable en la ecuación de Clairaut). Dada una ecuación de


Clairaut y = xy ′ + g (y ′) , con el cambio de variable y ′ = t se transforma en una ecuación
diferencial de variables separadas de la forma:
dt
y = xy ′ + g ( y ′) ⇒ {y ′ = t } ⇒ (x + g ′(t ) ) =0
dx
116 Ecuaciones Diferenciales

Demostración:

En efecto, sea la ecuación de Clairaut:


y = xy ′ + g ( y ′)

Si a esta ecuación se le aplica el cambio {y ′ = t } , resulta:

y = xt + g (t )

y derivando esta igualdad respecto a x se elimina la variable y de la ecuación:

dt dt
y′ = t = t + x + g ′(t )
dx dx

dt dt
0=x + g ′(t )
dx dx

y de aquí, resulta:

(x + g′(t )) dt =0
dx

Procedimiento (Resolución de la ecuación de Clairaut). Para resolver una ecuación de


Clairaut y = xy ′ + g (y ′) , se procede de la siguiente manera:

1. Se aplica a la ecuación de Clairaut el cambio de variable y ′ = t , y se deriva


respecto a x la ecuación resultante, obteniéndose una ecuación diferencial de
variables separadas:
dt
y = xy ′ + g ( y ′) ⇒ {y ′ = t } ⇒ (x + g ′(t ) ) =0
dx
2. Se resuelve esta ecuación diferencial, teniendo en cuenta dos casos:
dt
• Si = 0 , entonces t = C , de donde, llevando esta expresión a la
dx
ecuación inicial y = xt + g (t ) , se obtiene como solución general de la
ecuación de Clairaut la familia uniparamétrica de rectas y = Cx + g(C) .

• Si x + g ′(t ) = 0 , entonces x = −g ′(t ) , de donde, llevando esta expresión a la


ecuación inicial y = xt + g (t ) , se obtienen las ecuaciones paramétricas de
la solución singular de la ecuación de Clairaut:
x = −g ′(t )

 y = −g ′(t ) ⋅ t + g (t )

Ejemplo: Hallar la solución general de la ecuación de Clairaut:

y = xy ′ + y ′ − y ′2

Se trata de una EDO de Clairaut con g( y ′) = y ′ − y ′2 . En primer lugar, se


aplica el cambio de variable y ′ = t , y se deriva respecto a x la ecuación
resultante, con lo que se obtiene una EDO de variables separadas:

y = xt + t − t 2

dt dt dt
y′ = t = t + x + − 2t
dx dx dx
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 117

dt
(x + 1 − 2t ) =0
dx
y de aquí, pueden considerarse dos casos:
dt
• = 0 ⇒ t = C ⇒ y = Cx + C − C 2 (Solución general).
dx
x = 2t − 1 x = 2t − 1
• x + 1 − 2t = 0 ⇒ x = 2t − 1 ⇒  2 ⇒  2
 y = (2t − 1)t + t − t y = t
(Solución singular).

2.9. APLICACIONES GEOMÉTRICAS DE LAS ECUACIONES


DIFERENCIALES ORDINARIAS DE PRIMER ORDEN

2.9.1. Trayectorias ortogonales a una familia de curvas

En el estudio diferencial de las funciones reales de una variable se define el ángulo entre
dos curvas funcionales derivables y = f (x ) , y = g(x) en un punto de corte (x, y ) entre
ambas, como el ángulo θ que forman entre sí las rectas tangentes a f y g en dicho punto,
el cual viene determinado por la expresión:
θ = α − β = arctg(f ′(x ) ) − arctg(g ′(x ) )

y = g(x)

θ
y y = f (x)

α
β
X
O x

Figura 2.4: Ángulo θ entre las curvas y = f (x ) , y = g(x) .

Supóngase que y = f (x,C) es una familia uniparamétrica de curvas planas, y que


y ′ = F (x, y ) es la EDO de primer orden asociada a dicha familia. De acuerdo con el
significado geométrico de la derivada, es evidente que F(x,y) determina la pendiente de
la recta tangente a la curva de la familia que pasa por el punto (x,y). Esta recta tangente
forma un ángulo α respecto al eje OX dado por:
tg (α) = y ′ = F (x, y )

Por tanto, la recta tangente a una trayectoria y = g (x ) que corte perpendicularmente a la


curva y = f (x,C) en el punto (x,y) (trayectoria ortogonal) tendrá una pendiente igual a la
tangente del ángulo α + π , es decir:
2
118 Ecuaciones Diferenciales

(
tg α + π
2
) = F (−x1, y )

y = g (x )

π
2

y y = f ( x,C)

α+ π
α 2

O x X

Figura 2.5: Trayectoria y = g (x ) ortogonal a cada curva y = f (x,C ) .

Por definición de derivada se cumple que g ′(x ) = tg α + ( π


2
) = F (−x1, y ) , de donde se obtiene
la EDO de primer orden en forma normal que define a la familia uniparamétrica de
trayectorias y = g(x, K ) que verifican la condición de ortogonalidad planteada:

−1
y′ =
F (x, y )

Ejemplo: Trayectorias ortogonales a una familia de circunferencias. Dada la familia


uniparamétrica de circunferencias de centro el origen x2 + y 2 = C2 ,
determinar la ecuación de las trayectorias ortogonales a cada
circunferencia de la familia.
La EDO de primer orden asociada a la familia de circunferencias
x2 + y 2 = C2 puede obtenerse derivando la ecuación de la familia y
eliminando el parámetro C del sistema formado:

x2 + y 2 = C2

2x + 2 yy ′ = 0

−x
y′ =
y
Según esto, la recta tangente a la circunferencia de la familia que pasa
−x
por un punto arbitrario (x,y) tiene como pendiente y ′ = en ese punto.
y
Por tanto, la EDO de primer orden que define a la familia uniparamétrica
de trayectorias y = g(x, K ) que verifican la condición de ortogonalidad
planteada es:
−1 y
y′ = =
F (x, y ) x
Se trata de una EDO de primer orden de variables separadas:
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 119

y
y′ =
x
dy y
=
dx x
dy dx
=
y x

1 1
∫ y dy = ∫ x dx
( )
ln y = ln x + K( )
y =e
( )
ln x + K

y = ±e K x

y = Kx

Como se observa, la familia uniparamétrica de trayectorias ortogonales a


la familia de circunferencias x2 + y 2 = C2 es una familia de rectas que
pasan por el origen del plano: y = Kx .

y = Kx
2 2 2
x + y =C π
2
y

α+ π
2
α

O x X

Figura 2.6: Recta y = Kx ortogonal a cada circunferencia x2 + y 2 = C2 .

2.9.2. Trayectorias oblicuas a una familia de curvas

Al igual que en el caso anterior, supóngase que y = f (x,C ) es una familia uniparamétrica
de curvas planas, y que y ′ = F (x, y ) es la EDO de primer orden asociada a dicha familia.
De acuerdo con el significado geométrico de la derivada, es evidente que F(x,y)
determina la pendiente de la recta tangente a la curva de la familia que pasa por el
punto (x,y). Esta recta tangente forma un ángulo α respecto al eje OX dado por:
tg (α) = y ′ = F (x, y )

Por tanto, la recta tangente a una trayectoria y = g (x ) que corte a la curva y = f (x,C ) en
el punto (x,y) con un ángulo θ ≠ π (trayectoria oblicua) tendrá una pendiente igual a la
2
tangente del ángulo suma α + θ , es decir:
120 Ecuaciones Diferenciales

tg(α) + tg (θ) F (x, y ) + tg (θ)


tg (α + θ) = =
1 − tg (α)tg(θ) 1 − F (x, y )tg(θ)

y = g (x )

θ
y y = f ( x,C)

α+θ
α
O x X

Figura 2.7: Trayectoria y = g (x ) oblicua a cada curva y = f (x,C ) con ángulo constante
θ≠ π .
2

F (x, y ) + tg(θ)
Por definición de derivada se cumple que g ′(x ) = tg(α + θ) = , de donde se
1 − F (x, y )tg(θ)
obtiene la EDO de primer orden en forma normal que define a la familia uniparamétrica
de trayectorias y = g(x, K ) que cortan a cada curva y = f (x,C ) con ángulo constante
θ≠ π :
2

F (x, y ) + tg(θ)
y′ =
1 − F (x, y )tg(θ)

Ejemplo: Trayectorias oblicuas a una familia de rectas. Dada la familia uniparamétrica


de rectas y = Cx , determinar la ecuación de la trayectoria que corta a
cada recta de la familia formando un ángulo constante θ ≠ π .
2

La EDO de primer orden asociada a la familia de rectas y = Cx puede


obtenerse derivando la ecuación de la familia y eliminando el parámetro C
del sistema formado:
y = Cx

y′ = C

 y = Cx  y
C = y
 ⇒ x ⇒ y′ =

y = C C = y ′ x

Según esto, la recta tangente a la recta de la familia y = Cx que pasa por


un punto arbitrario (x,y) (que en este caso coincide con la propia recta
y
y = Cx ) tiene como pendiente y ′ = en ese punto. Por tanto, la EDO de
x
primer orden que define a la familia uniparamétrica de trayectorias
y = g(x, K ) que cortan a cada recta de la familia formando un ángulo
constante θ ≠ π es:
2
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 121

y
+ tg(θ)
F (x, y ) + tg(θ)
y′ = = x
1 − F (x, y )tg (θ) y
1 − tg(θ)
x
Se trata de una EDO de primer orden homogénea:
y
+ tg(θ)
y′ = x
y
1 − tg(θ)
x
 y = tx 
Para resolverla, se aplica el cambio   , obteniéndose:
dy = xdt + tdx 
y
+ tg(θ)
dy
= x
dx 1 − y tg(θ)
x
y
+ tg(θ)
dy = x dx
y
1 − tg(θ)
x
tx
+ tg(θ)
xdt + tdx = x dx
tx
1 − tg(θ)
x
t + tg(θ)
xdt + tdx = dx
1 − t ⋅ tg(θ)

( )
(1 − t ⋅ tg(θ))xdt + t − t 2 ⋅ tg(θ) dx = (t + tg(θ))dx

(1 + t )⋅ tg(θ)dx = (1 − t ⋅ tg(θ))xdt
2

A continuación, se integra la ecuación de variables separadas obtenida:


1 1 − t ⋅ tg(θ)
dx = dt
x (
1 + t 2 ⋅ tg(θ) )
1 1 − t ⋅ tg(θ)
∫ x dx = ∫ (1 + t 2 )⋅ tg(θ) dt
1 1 2t
∫ x dx = 1
tg ( θ) ∫ 1 + t 2 dt − 12 ∫ 1 + t 2 dt
( )
ln x = 1 arctg (t )
tg (θ)
( )
− 12 ln 1 + t 2 + K

( )
1 arctg (t ) − 1 ln 1+t 2 +K
x = e tg ( θ) 2

1 arctg (t )
± eK e tg (θ)
x=
1 + t2
122 Ecuaciones Diferenciales

1 arctg (t )
Ke tg (θ)
x=
1+ t2
Finalmente, se deshace el cambio de variable:
1 arctg  y 
tg ( θ) x
Ke
x=
y2
1+
x2

Y
1 arctg  y 
tg ( θ) x
Ke
x=
y2 y = Cx
1+
x2 θ
y

α+θ
α
O x X

1 arctg  y 
tg ( θ) x
Ke
Figura 2.8: Trayectoria x = oblicua a cada recta y = Cx con
y2
1+
x2
ángulo constante θ ≠ 2π .
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 123

EJERCICIOS

1. Aplicando el teorema básico de existencia y unicidad, determinar el dominio D


del plano en cada uno de cuyos puntos (a, b) ∈ D existe y es única la solución de
las siguientes EDO:

a. y′ = x4 + y 4 . Solución: D = o2

b. y′ = 3 x + y .
Solución: D = {(x, y ) ∈ o 2 | y ≠ x }

c. y′ = x y .
Solución: D = {(x, y ) ∈ o 2 | y > 0}

dx 3
d. = tx .
dt
Solución: D = {(x, y ) ∈ o 2 | x ≠ 0}

2. Verificar si se cumplen las condiciones del teorema básico de existencia y


unicidad para los siguientes problemas de valor inicial:

a. ( )
y ′ = ln 1 + x2 y 2 ; y (0) = 0 .

b. y′ = 3
xy 2 ; y (0) = 0 .

c. y′ = y 3 ; y (0) = 0 .

d. y ′ = 3 x 3 + y 3 ; y (0) = 0 .

3. Hallar la solución general de las siguientes ecuaciones de variables separadas:

a. ( )
(x − 1) y 2 − 1 dx + x2dy = 0 . Solución:


1
x

y = th ln( x ) + + C 

b. y ′ = x ln( x ) cos2 ( y ) .
Solución: (
y = arctg 12 x 2 ln( x ) − 1
4
x2 + C )
1 − cos (x)
c. y ′sen(x) = y . Solución: y=C
1 + cos(x )

1
d. x 2 ydy = 1 − y 2 dx . Solución: 1 − y2 = +C
x

4. Resolver los siguientes problemas de valor inicial asociados a ecuaciones de


variables separadas:

a. dx = e2t cos(t )dt ; x(0) = 1 .


Solución: x= 1 e 2t (sen(t ) + 2 cos(t )) + 35
5
124 Ecuaciones Diferenciales

1 x − 1 sen(2x )
b. y ′ = sen2 (x )(1 + y ) ; y (0) = 0 . Solución: y = e2 4 −1

c. − e− xdt + (t 2 + 2t + 2)dx = 0 ; x(0) = ln 2π . ()


Solución: (
x = ln arctg(t + 1) + 4π )
dy ch( y )   1 
d. = ; y (1) = 0 . Solución: y = argsh tg1 −  
dx x2   x 

5. Verificar si las siguientes ecuaciones diferenciales son homogéneas:

x2
a. (y 2
)
− xy dx +
x+y
dy = 0 .

1 5
b. dx + dy = 0 .
4 3
x +x y x + y2
2

 x − 3y 
c. y ′ = cos  .
 2x + 2y 

3 2
dy 2x − xy
d. = .
dx x + 2y

6. Hallar la solución general de las siguientes ecuaciones diferenciales homogéneas:


2
a. (y − x )dx + x dy = 0 . Solución: y=
x
y ln( x ) − C

y
4 xy − y 2 x
b. y′ = . Solución: x =C
2x 2 y
2−
x

y 2 − 4x 2 C
c. y′ = . Solución: x=
2xy 2
y
+4
x2
C
d. 2xyy ′ = x 2 − y 2 . Solución: x=
3
3y 2
−1
x2

7. Resolver los siguientes problemas de valor inicial asociados a ecuaciones


diferenciales homogéneas:

 y 
a.  xtg  + y dx − xdy = 0 ; y (1) = 6π . Solución: y = x ⋅ arcsen 12 x ( )
 x 

b. ( x+y + )
x − y dx + ( x−y − )
x + y dy = 0 ; y (1) = 1 .
1
Solución: x=
y2
1+ 1−
x2
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 125

x
e
−x t
2 t x
dx (t + x) e + x2 − e 1+
c. = ; x (1) = 1 . Solución: t= e2e t
dt tx

x2 + y 2
− 5
y
xy 2e 2 xe
d. y′ = ; y (1) = 2 . Solución: x=
x + y x2 + y 2
2 y

8. Integrar las siguientes ecuaciones diferenciales reducibles a homogéneas:


x− y −3
a. y′ = .
x + y −1
C
Solución: x −2 =
2( y + 1) ( y + 1) 2
1− −
x −2 (x − 2) 2

b. (x − 2y + 4)dy = (x − y + 1)dx .
C
Solución: x −2 =
2( y − 3) 2( y − 3)2
1− +
x −2 (x − 2) 2

2x − 2y + 4
c. y′ = .
2x − 2 y − 1
5
 2( y − x) − 4 
−5  − 8
Solución: x= ln 1 − +
16  2( y − x) + 1  2
   2( y − x ) − 4 
1 − 
 2( y − x) + 1 

  5
5 ln 1 + 2( y − x ) − 4
+ 16 − 8 +C
 2( y − x ) + 1  2( y − x) − 4
  1+
2( y − x ) + 1

dx t +x+2
d. = ; x(1) = 1 .
dt 2t + 2x − 1
Solución: t = 23 (t + x) − 5 ln
9
(3t + 3x + 1 ) + 59 ln(7) − 13
9. Verificar si las siguientes ecuaciones diferenciales son exactas:

 1  ln(x ) 
a.  2xe y + dx +  x 2e y − dy = 0 .
 xy   y 2 

b. (ln(x) − 3 x − y )dx + (sen( y) + 3 x − y )dy = 0 .


3x
c. y ′ = xy − 3tg( y ) .
cos 2 ( y )

d. (xch(t ) + x )dt + (sh(t ) + 2tx )dx = 0 .


2
126 Ecuaciones Diferenciales

10. Hallar la solución general de las siguientes ecuaciones diferenciales exactas:

 2 y   1 y
a.  y − 2 dx +  2xy + dy = 0 . Solución: xy 2 + +C = 0
 x   x x

1 1 y   x 
b.  + − dx +  y cos( y ) + dy = 0 .
 x x2 x2 + y 2   x2 + y 2 
  
1 x
Solución: ln( x ) − − arctg  + ysen( y ) + cos( y ) + C = 0
x y

c. (x 3
) (
− yexy dx + 2y − xexy dy = 0 . ) Solución: 1 x4
4
− e xy + y 2 + C = 0

 
 1  dρ θ
d.  − θ  dθ = ρ − e .
2
 − 3 + 4ρ − ρ 
Solución: eθ − θρ + arcsen(ρ − 2) + C = 0

11. Resolver los siguientes problemas de valor inicial asociados a ecuaciones


diferenciales exactas:

− ex − y
a. y′ = ; y (0) = 1 .
2 + x + ye y
Solución: ex + xy + 2 y + ( y − 1)e y − 3 = 0

 
b.
 1− y
 2
+ cos( x ) − 2xy
 2
(
dy = y + ysen(x ) dx ; y (0) = 1 . )
 1− y 
Solución: − xy 2 + y cos(x ) + arcsen( y ) + 1 − y 2 − 1 − π
2
=0

 
c. (x + y )2dx +  x2 + 2xy + 2 y − 21 dy = 0 ; y (1) = 1 .
 1+ y 
Solución: 1 (x
3
+ y )3 − 1
3
( )
y 3 + ln 1 + y 2 − arctg( y ) + π
4
− ln(2) − 7
3
=0

d.



( )
 1 − x2 + y 3 dx + 3xy 2 + 1 dy = 0 ; y (0) = 0 .

Solución: x3 y + y + 1 arcsen( x ) + 1 x
2 2
1 − x2 = 0

12. Hallar la solución general de las siguientes ecuaciones diferenciales mediante


factores integrantes:

a. (x 2
) (
− y 2 − x − y dx + y 2 − x2 dy = 0 , ) sabiendo que admite un factor
integrante dependiente de x+y.
Solución: 1 x2 − xy − x + 1 y2 + C = 0
2 2

b. (1 − x y )dx + x (y − x )dy = 0 ,
2 2
sabiendo que admite un factor integrante
dependiente sólo de x.
1
Solución: 1
2
y 2 − xy − +C = 0
x
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 127

y2
c. y′ = , sabiendo que admite un factor integrante dependiente de
x2 + xy
x2y.
y
Solución:
x
( )
+ ln y + C = 0

d. (3y 2
) ( )
− x dx + 2 y 3 − 6xy dy = 0 , sabiendo que admite un factor integrante
dependiente de x+y2.
− 2y 2 1
Solución: + +C = 0
(x + y ) 2 2 x + y2

13. Hallar la solución general de las siguientes ecuaciones diferenciales mediante


factores integrantes:

a. (
6xydx + 4 y + 9x2 dy = 0 . ) Solución: 3x 2 y 3 + y 4 + C = 0

b. (x 4
) (
ln( x ) − 2xy 3 dx + 3x2 y 2 dy = 0 . )
Solución: y= 3
− x 3 ln(x) + x 3 − Cx 2

c. (− xysen(x ) + 2 y cos(x ))dx + 2x cos(x )dy = 0 .


−C
Solución: y=
2x cos( x)

1
−3
2
y2 + −1
d. xy ′ = x .
y
Solución: 1
2
x3y2 + 1
3
x3 − 1
2
x2 + C = 0

14. Resolver los siguientes problemas de valor inicial mediante factores integrantes:

a. (2x 3
) (
y + 3x 2 dx + 4x 2 y 2 + x 4 dy = 0 ; y (0) = 1 .)
2
Solución: x y + 3x + 4
3
y3 − 4
3
=0

3xy 2
b. (x 2
− xy + y 2 dx + ) x+y
dy = 0 ; y (1) = 1 .

5
Solución: y= 3− 1
4
x3 +
4x

 y 
c.  xych(x ) + sh( y ) dx − ych( y )dy = 0 ; y (1) = 0 .
 x 
Solución: y = argsh(xsh(x ) − xsh(1) )

x  π 
d. 2arctg( y )dx + 2
dy = 0 ; y (1) = 1 . Solución: y = tg  
1+ y  4x 2 
128 Ecuaciones Diferenciales

−4 y 2x + 2y
15. Dada la EDO dx + dy = 0 , deducir la solución general a partir de
2
(x − y ) (x − y ) 2
un factor integrante dependiente sólo de y, y un factor integrante dependiente de
x−y.

x−y
Solución: =C
y

16. Integrar las siguientes ecuaciones lineales homogéneas:


4 x 3 − 1 x + 1 sen(2x )
a. (
y ′ + sen2 (x) − 4x2 y = 0 . ) Solución: y = Ce 3 2 4

dy
b. = y 1 + x2 .
dx
1 x 1+ x2 + 1 argsh( x )
Solución: y = Ce 2 2

1 + 2x − x2
c. y′ + y = 0 .
x
 
1+ 2x − x2 − arcsen 1 x − 1 
Solución: y = Ce  2 2

1
1−
−e x
d. x 2 y ′ = (1 − x ) y ; y (1) = −1 . Solución: y=
x

17. Integrar las siguientes ecuaciones lineales completas:


y
a. y′ − = sen(x ) + 1 , sabiendo que una solución particular es
cos( x )
y = − cos( x ) .
1 + sen(x )
Solución: y =C − cos (x)
1 − sen(x )

b. y ′ − y = e x , sabiendo que una solución particular es y = xe x .


Solución: y = Ce x + xe x

1 + 3x 3 −1
c. y ′ − 3x 2 y = , sabiendo que una solución particular es y = .
2 x
x
3 1
Solución: y = Ce x −
x
d. y ′ − 2th(x ) y = −sh(x ) , sabiendo que una solución particular es y = ch(x ) .
Solución: y = C ⋅ ch2 (x ) + ch(x )

18. Hallar la solución general de las siguientes ecuaciones lineales completas


mediante el método del factor integrante:
C
a. y ′ + 2xy = 2x 3 . Solución: y = x2 − 1 − 2
ex
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 129

2xy
b. y′ − = x + 1.
1 + x2
Solución: (
y = 1 + x2 1
2
)( ( )
ln 1 + x 2 + arctg(x ) − C )
dy 3x2 − 1
c. − 3 y = x+3.
dx x −x
Solución: ( )( ( ) (
y = x 3 − x − 3 ln x + ln x + 1 + 2 ln x − 1 − C ) ( ) )
dx tg2 (t ) tg(t ) − t − C
d. +x = . Solución: x=
dt et et

19. Resolver los siguientes problemas de valor inicial asociados a ecuaciones lineales
completas mediante el método del factor integrante:
y
a. y′ − = x2 ; y (1) = 0 . Solución: y = 12 x 3 − 12 x
x
1
b. y ′ + ytg(x ) = ; y (0) = 1 . Solución: y = sen(x ) + cos( x )
cos(x )

c. xy ′ − y = ln( x ) ; y (1) = 1 . Solución: y = − ln(x ) + 2x − 1

dx
d. − x = sen(2t ) ; x (0) = 0 .
dt
Solución: x = −51 (sen(2t ) + 2 cos(2t ) ) + 2
5
et

20. Hallar la solución general de las siguientes ecuaciones lineales completas


mediante el método de variación de constantes:
a. y ′ + y = sen(x ) .
Solución: y = 12 (sen(x ) − cos( x )) + Ce − x

y
b. 1 − x2 ⋅ y ′ − =1.
arcsen(x )
Solución: ( ( )
y = ln arcsen(x ) + C arcsen(x ) )
dy 3
c. + 3x 2 y = x 2 . Solución: y = 13 + Ce − x
dx

d.
dy
+
y
dx 2 x
−4 = 0. Solución: ( )
y = 8 x − 1 + Ce − x

21. Resolver los siguientes problemas de valor inicial asociados a ecuaciones lineales
completas mediante el método de variación de constantes:
y
a. y′ − = xarctg(x ) ; y (1) = 1 .
x
Solución: ( ) (
y = x 2arctg (x ) − 12 x ln 1 + x 2 + 1 + 12 ln(2) − 4π x )
(2x + 1) y
b. y′ + = e −2x ; y (1) = 1
2
.
x
−2x
Solución: (
y = 12 xe − 2x + 12 e 2 − 1 )e x
130 Ecuaciones Diferenciales

dy 2
c. + 2xy = x ; y (0) = 2 . Solución: y = 12 + 3
2
e−x
dx
dy y 1
d. + = ; y (1) = 1 . Solución: y =1
dx x2 x2

22. Integrar las siguientes ecuaciones diferenciales reducibles a lineales:


y x C
a. y′ + = 3. Solución: y 4 = x2 −
2x y x2

dx 4
b. t + x = t 3x 3 ; x (1) = 4 . Solución: x=
dt t (t − 2) 2

c. y ′ = (1 − x ) y 2 + (2x − 1) y − x , sabiendo que una solución particular es de la


forma y = m .
1
Solución: y =1+
x − 2 + Ce − x
6
d. y ′ = −2 y 2 + , sabiendo que una solución particular es de la forma
x2
m
y= .
x
2 1
Solución: y= +
x −2
7
x + Cx 8

23. Integrar las siguientes ecuaciones de Lagrange y Clairaut:

 −t −C
x =
a. y = 2xy ′ + ln( y ′) . Solución:  t2
 y = − 2t − 2C + ln(t )
 t

 −1 C
3 x = 2 +
b. y = − xy ′ − . Solución:  t t
y′ y = − 2 − C t
 t

c. y = xy ′ + 1 + y ′ .
 −1
x =
Solución: y = Cx + 1 + C (general);  2 1+ t (singular)
−t
y = + 1+ t
 2 1+ t

d. y = xy ′ − 2 − y ′ .
x = 1
Solución: y = Cx − 2 − C (general);  (singular)
 y = −2

24. Obtener las trayectorias ortogonales a las siguientes familias de curvas planas:
1 x2 + 1 y 2
a. Familia de elipses C 2 x 2 + y 2 = 1 . Solución: y = Ke 2 2
2. Ecuaciones Diferenciales Ordinarias de Primer Orden 131

K
b. Familia de hipérbolas x 2 − y 2 = C . Solución: y=
x

25. Obtener las trayectorias oblicuas que cortan a las siguientes familias de curvas
planas con un ángulo constante θ = 4π :

a. Familia de parábolas y = Cx 2 .
3 arctg  4 y + 1 
 
Ke 17  17 x 17 
Solución: x=
2y 2 y
+ +1
2 x
x

b. Familia de circunferencias x 2 + y 2 = C .
K
Solución: x=
y
arctg  
y2 x
+1 ⋅e
x2

También podría gustarte