FOLLETOEDO1 V 1
FOLLETOEDO1 V 1
2018
Índice general
Sección 1:
Elementos básicos y teorema de existencia
unicidad en Ecuaciones Diferenciales
Ordinarias.
1 — CONTENIDO 1: Elementos básicos y Teorema de existencia y unicidad en EDO1
1.0.2 Contenidos
Definición de ecuación diferencial.
Orden y grado de una ecuación diferencial.
Clasificación de ecuaciones diferenciales ordinarias.
Solución y tipos de solución en una ecuación diferencial ordinaria.
Problema de valor inicial (PVI) y problema de valor frontera (PVF).
Teorema de existencia y unicidad en ecuaciones diferenciales ordinarias de primer orden.
1.1 Conceptos Básicos en Ecuaciones Diferenciales
A continuación se muestran los aspectos teóricos relacionados con los elementos básicos en las ecuaciones
diferenciales.
Definición 1.1
Una ecuación diferencial es una ecuación que relaciona una función desconocida (la variable
dependiente), las variables de las que depende (variables independientes) y sus derivadas respecto de
estas variables independientes.
En particular si dicha ecuación contiene derivadas respecto a una sola variable independiente se le
denomina Ecuación Diferencial Ordinaria (EDO) y si contienen derivadas respecto a dos o más
variables independientes se le llama Ecuaciones diferenciales en derivadas parciales (EDP).
1.1
dy
a) + 5y = ex : EDO de variable dependiente y, variable independiente x.
dx
00
b) y + y − 2x = 0 : EDO de variable dependiente y, variable independiente x.
dx dy
c) + = 2x + y : EDO de variable dependiente x, y, variable independiente t.
dt dt
Una ecuación diferencial ordinaria (EDO) puede contener más de una variable dependiente como en el
ejemplo anterior c).
1.2
∂u ∂u
a) =− : EDP de variable dependiente u, variable independiente y, x
∂y ∂x
∂ 2u ∂ 2u
b) + = 0 : EDP de variable dependiente u, variable independiente x, y
∂ x 2 ∂ y2
∂ 2u ∂ 2u ∂u
c) 2
= 2 −2 : EDP de variable dependiente u, variable independiente x,t
∂x ∂t ∂t
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
Definición 1.2
El orden de una EDO se define por el orden de la más alta derivada en la ecuación dada. En general,
el n−ésimo orden de una EDO puede simbolizarse por
1.3
Considere las siguientes ecuaciones diferenciales ordinarias y determinar que orden tienen cada una
de ellas.
00
a) y − 5y0 + 3y3 = 0.
d3y d2y
b) + 2 + y = x.
dx3 dx2
dy
c) M (x, y) + N (x, y) = 0.
dx
De acuerdo a cada ecuación diferencial dada se tiene lo siguiente:
00
(a) la ecuación diferencial y − 5y0 + 3y3 = 0 es de orden 2.
d3y d2y
(b) la ecuación diferencial + 2 + y = x es de orden 3.
dx3 dx2
dy
(c) toda ecuación de la forma M (x, y) + N (x, y) = 0 es de orden 1.
dx
Definición 1.3
El grado de una ecuación diferencial ordinaria se define por la potencia sobre y(i) , i = 1, . . . , n,
que tiene la más alta derivada en la ecuación dada.
1.4
donde a0 (x) , a1 (x) , ..., an (x) , g (x) son funciones reales continuas que dependen sólo de la
variable independiente x y definidas en un intervalo abierto común I, con an 6= 0 para toda x ∈ I.
Una ecuación diferencial que no puede escribirse de la forma 1.2 es No lineal.
ecuaciones diferenciales lineales con coeficientes constantes, cuando todos las funciones
ai (x) ∀i = 1, ..., n son funciones constantes.
1.5
Considere las siguientes ecuaciones diferenciales ordinarias y clasificarlas en EDO lineal (con
coeficientes contantes o variables) y no lineal.
00
(a) y − 2y0 + y = 0
d2y dy
(b) 2 + 3 = ex .
dx dx
(c) x3 y000 − x2 y00 + 3xy0 + 5y = x.
d2y dy
(d) 2 + (x + 1) + 2 = 0.
dx dx
(e) ex y00 + 2xy = 0.
(f) yy00 − 2y0 = x.
d2y
(g) 2 + sen y = 0.
dx
(h) y000 + y2 = cos x.
De acuerdo a las definiciones dadas anteriormente, se tiene que las ecuaciones diferenciales (a) y
(b) son lineales con coeficientes constantes. Las ecuaciones diferenciales (c), (d) y (e) son lineales
con coeficientes variables. Las ecuaciones diferenciales ( f ), (g) y (h) son no lineales.
Definición 1.5
Si una ecuación diferencial de orden n dada en la forma F x, y, y0 , ..., y(n) = 0 puede expresarse de
manera que la derivada de orden n aparezca despejada, es decir:
y(n) = G x, y, y0 , ..., y(n−1) (1.3)
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
1.6
Para cada una de las ecuaciones diferenciales ordinarias dadas, escribir en forma normal las ecuacio-
nes diferenciales:
a) yy00 − 2y0 = x.
d2y
b) + sen x = 0.
dx2
c) y000 + y2 = cos x.
Se despeja la derivada más alta en cada una de las ecuaciones diferenciales del ejemplo y obtenemos
su forma normal:
y0 x
(a) yy00 − 2y0 = x, por tanto y00 = 2 + con y 6= 0.
y y
d2y d2y
(b) + sen x = 0, por tanto = − sen x.
dx2 dx2
(c) y000 + y2 = cos x, por tanto y000 = −y2 + cos x.
00
Resolver una ecuación diferencial ordinaria, F x, y0 , y , y000 , ... = 0, es hallar una expresión para la función
y(x) que satisfaga la relación de igualdad que determina dicha ecuación, ∀x ∈ I un intervalo común.
Definición 1.6
00
Se llama solución de una ecuación diferencial ordinaria F x, y0 , y , y000 , ... = 0 en un intervalo
J a una función ϕ(x) definida en J que, sustituida en la ecuación junto con sus derivadas, la
verifica en dicho intervalo, es decir:
Es importante rescatar que una solución puede ser definida en alguno de los siguientes intervalos J: ]α, β [,
]α, β ], [α, β [, [α, β ], [α, ∞[, ]α, ∞[, ] − ∞, β ], ] − ∞, β [ y ] − ∞, ∞[= R; donde α, β ∈ R y α < β .
Se puede comprobar que la función ϕ(x) = e5x es solución de la ecuación diferencial ordinaria de
primer orden y0 = 5y en el intervalo ] − ∞, ∞[.
En efecto, derivando ϕ(x) = e5x respecto de x, se obtiene ϕ(x)0 = 5e5x , luego sustituyendo en la
ecuación diferencial, se verifica la ecuación diferencial: ϕ(x)0 = 5e5x = 5ϕ (x) , ∀x ∈ R.
1.8
1
Considere la función ϕ(x) = x2 − y verificar es solución de la ecuación diferencial ordinaria de
x
00 2
segundo orden y − 2 y = 0, ∀x 6= 0.
x
1
Derivando dos veces la función ϕ(x) = x2 − respecto de x, se tiene:
x
1 2
ϕ(x)0 = 2x + 2
, ϕ(x)00 = 2 − 3 .
x x
2 2 2 1 2 2
2− 3 − 2 x − = 2 − 3 − 2 + 3 = 0, ∀x 6= 0.
x x x x x
1
Por tanto, ϕ(x) = x2 − es solución en el intervalo ]−∞, 0 [ ∪ ] 0, +∞[.
x
1.9
Se puede comprobar que toda función de la forma ϕ(x) = C1 e−x +C2 e2x es solución de la ecuación
diferencial y00 − y0 − 2y = 0 para cualquier valor de las constantes C1 y C2 .
En efecto, derivando dos veces ϕ(x) = C1 e−x +C2 e2x , se tiene que:
C1 e−x + 4C2 e2x +C1 e−x − 2C2 e2x −C1 e2x − 2C2 e2x = 0. ∀x ∈ R.
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
1.10
c dy
Verifique que la función φ (x) = + 1, x > 0 es una solución de la ecuación x + y = 1.
x dx
dy
En efecto, si sustituimos φ en el lugar de y en la ecuación x + y = 1:
dx
d c c
x +1 + +1
dx x x
, luego simplificando obtenemos:
c c
x − 2 + +1 = 1
x x
c
que es igual al lado derecho de la ecuación. Por lo que φ (x) = + 1 es solución de la ecuación
x
dy
x +y = 1
dx
.
Observe que el intervalo en el cual la función φ está definida I = (0, +∞) excluye el punto de
c
discontinuidad, igual podemos ver que la función definida φ (x) = + 1 pero para x < 0 es solución
x
también.
1.11
La ecuación diferencial (y0 )2 + 1 = 0 no tiene solución real, pues no existe una función real y
derivable que la satisfaga, es decir, no es posible se cumpla que exista y(x) tal que (y0 )2 sea -1.
1.12
2x + 2yy0 = 0,
es decir,
dy
2x + 2y = 0.
dx
dy x
Por tanto = − (y 6= 0) y de esta forma se verifica la ecuación diferencial para −2 < x < 2,
dx
√ y
pues y = ± 4 − x2 .
sen x + y ln y = C · y
la y es una función implícita de x (y > 0), adémas C una constante arbitraria. Compruebe que ésta es
solución implícita de la ecuación diferencial
y cos x
y0 =
sen x − y
.
Si en la familia de curvas dada sen x + y ln y = C · y despejamos la constante C se tiene
sen x + y ln y
C= [1]
y
Ahora derivando [1] implicítamente con respecto a la variable x, tenemos
Como se observa en los ejemplos, las soluciones de una ecuación diferencial pueden venir dadas como
una función o una expresión que verifica la ecuación. Por ello, las llamamos soluciones explícitas si la
solución es una expresión de la forma y = y(x) y llamamos soluciones implícitas si la solución es una
expresión de la forma g(x, y) = 0.
1.14
2xy
La relación y − ln y = x2 + 1 es una solución implícita de la EDO y0 = para y > 0 y x ∈ R,
y−1
y0
pues, derivando la relación dada con respecto a x, se tiene y0 − = 2x, es decir, y0 y − y0 = 2xy.
y
2xy
Luego como y0 = y sustituyendo en y0 y − y0 = 2xy se verifica dicha igualdad para y > 0 y
y−1
x ∈ R.
1.15
3x2
Se puede verificar que y2 − x3 + 8 = 0 es solución implícita de la ecuación diferencial y0 =
2y
en ]2, +∞[ .
3x2
y0 =
2y
y sustituyendo en la ecuación diferencial, se tiene
3x2 3x2
=
2y 2y
y esto√se verifica ∀x ∈ ]2, +∞[ , pues si despejamos y en la expresión de la solución implícita,
y = ± x3 − 8, vemos que esta definida para x ≥ 2.
1.16
Z x
sen(t)
Considere la función definida por y = x dt y comprobar es solución de la ecuación
0 t
diferencial xy0 = y + x sen(x).
En ejemplos anteriores se ha resuelto ejercicios sobre comprobación o verificación de si cierta función, una
curva o familia de estas es o no solución ya sea explícita o implícita de una ecuación diferencial dada, en
próximos capitulos nos encargaremos del problema de hallar la solución general de una ecuación diferencial.
Ahora se aborda el problema inverso; a saber, el de hallar la ecuación diferencial a partir de su solución
general. Para ello se ilustra tal proceso a través de los siguientes ejemplos.
1.18
Determine una ecuación diferencial cuya solución general viene dada en forma implícita por xy2 =
C
x + , con C constante arbitraria. Note que:
x
C
xy2 = x + ⇒ x2 y2 = x2 +C
x
Aplicando derivación implícita con respecto a x a ambos lados le lo obtenido en paso anterior, se
tiene que la ecuación diferencial correspondiente es:
1.19
Determine una ecuación diferencial para la siguiente familia de curvas dada por y = A(x − B)3 .
y
Si despejamos la constante A y derivamos implícitamente con respecto a x la curva dada = A,
(x − B)3
de donde
y0 (x − B)3 − 3y(x − B)2 y0 3y
= 0 ⇒ − =0
(x − B)6 (x − B)3 (x − B)4
3y
Multiplicando por (x − B)4 , obtenemos y0 (x − B) = 3y ⇒ x − B − = 0. Derivando de nuevo
y0
implícitamente con respecto a x
3(y0 )2 − 3yy00
1=
(y0 )2
⇒ (y0 )2 = 3(y0 )2 − 3yy00
Halle la ecuación diferencial cuya solución viene dada por la curva H : ecx + y = x.
Observación 1.1 Note que en estos dos últimos ejemplos, la idea consiste básicamente en derivar las
veces según cantidad de constantes que posea la solución general,así como efectuar ciertas maniobras
algebraicas hasta llegar a tener una ecuación diferencial libre de dichas constantes.
1.21
Por ejemplo y(x) = C1 e−x + C2 e2x es solución 2-paramétrica (solución general) de la ecuación
diferencial ordinaria y00 − y0 − 2y = 0 para C1 y C2 constantes arbitrarias.
x2
La familia de curvas solución y = + C es una solución 1-paramétrica (solución general) de la
0
2
ecuación diferencial y = x, las cuales se muestran en la figura 1.1 para ciertos valores de C.
y
3
2
C=-2
1 C=-1
C=0
-2 -1 1 2 x
C=1
-1 C=2
-2
-3
x2
Figura 1.1: gráficas de la familia 1-paramétrica y = +C, para y0 (x) = x.
2
Fuente: elaboración propia con M ATHEMATICA 9.0.
Solución particular: es una solución de la ecuación diferencial que no contiene constantes arbitrarias y
que se obtiene dando valores numéricos a las constantes de la familia n-paramétrica de soluciones. Por
x2
ejemplo para y = +C solución 1-paramétrica de la ecuación diferencial y0 = x, si C = 0, C = −1 se
2
x2 x2
obtienen las soluciones particulares y = , y = − 1 respectivamente y mostradas en la figura 1.2.
2 2
y
3 2
y= x
2
-2 -1 1 2 x
2
-1 y= x -1
2
-2
-3
x x
Figura 1.2: gráficas de soluciones particulares y = , y y= − 1.
2 2
Fuente: elaboración propia con M ATHEMATICA 9.0.
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
Solución singular: es una solución de la ecuación diferencial que no contiene constantes arbitrarias y no
esta contenida en solución general. No siempre existen; si existe, se trata de la curva llamada envolvente de
la familia de curvas integrales definida por la familia n-paramétrica de soluciones.
1.23
x2
La función y(x) = es una solución singular de la ecuación (y0 )2 − xy0 + y = 0 en J = R, ya
4
que satisface dicha EDO y no puede obtenerse de sustituir C por algún valor numérico de la solución
general y(x) = Cx −C2 .
y
2
Solución 1-Paramétrica
1
y = C x-C 2
-2 -1 1 2 x
-1 Solución Singular
y HxL=
-2 x2
4
-3
x2
Figura 1.3: gráfica de la solución singular y(x) = y familia 1-paramétrica y(x) = Cx −C2 .
4
Fuente: elaboración propia con M ATHEMATICA 9.0.
satisface dicha EDO y√no puede obtenerse de sustituir C por algún valor numérico de la solución
general y(x) = xC + 1 +C2 .
.
−x
En efecto, note que al derivar implícitamente y con respecto a x se obtiene 2x + 2yy0 = 0 ⇔ y0 = .
y
Luego al sustituir este resultado en la ecuación diferencial se verifica que
s 2 s
−x2 x2
−x −x
y=x + 1+ ⇔ y= + 1+ 2
y y y y
s
−x2 x 2 + y2
⇔ y= +
y y2
s
−x2 1
⇔ y= +
y y2
−x2 1
⇔ y= +p ; y>0
y y2
−x2 1
⇔ y= +
y y
2
−x + 1
⇔ y=
y
y 2
⇔ y=
y
⇔ y = y.
√
Por otro lado, al sustituir y(x) = xC + 1 +C2 en x2 + y2 = 1 para encontrar un C en común, se
determina que
p p
x2 + (xC + 1 +C2 )2 = 1 ⇔ (xC + 1 +C2 )2 = 1 − x2
p
⇔ x2C2 + 2xC 1 +C2 + 1 +C2 = 1 − x2 .
Sin embargo, al cambiar C por cualquier valor numérico se determina que x debe ser fija,
√ lo cual
contradice el hecho de que x es variable. En la figura
√ 1.4 se muestra como y = 1 − x2 es
2
envolvente de la familia 1-paramétrica y(x) = C.x + 1 +C .
Observación 1.2 En este caso, la solución singular está restringida al dominio [−1, 1]. De esta manera,
J ⊆ [−1, 1] donde J es cualquier subconjunto de restricción del dominio.
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
y
2
Solución 1-Paramétrica
1
y = Cx+ 1 + C2
-2 -1 1 2 x
-1 Solución Singular
-2
1 - x2
-3
√
Figura 1.4: gráfica de la solución singular x2 + y2 = 1 y la solución general y(x) = C.x + 1 +C2 .
Fuente: elaboración propia con M ATHEMATICA 9.0.
Solución Total de una ecuación diferencial ordinaria de orden n : es la que contiene todas las solucio-
nes de la ecuación. Esta formada por la solución general y las soluciones singulares en caso de que hallan o
se puedan determinar.
Observación 1.3 Por lo general para efectos prácticos cuando resolvemos una ecuación diferencial,
suponemos están las funciones y derivadas de ella bien definidas para evitar analizar casos de donde se
puede desprender algunas soluciones singulares, con ello, se pide solo solución general y no total.
Definición 1.9
La gráfica de una solución de una ecuación diferencial se denomina curva integral de la ecuación
diferencial.
1.25
√ x +C
S= y= ,C ∈ R ∪ {y = 0} .
2
Hx + 2L
5
1
C = -2 C = -1 y =
2
C = -1,5 C =-0.5
4
C=0 C = 0.5
C=1 C = 1.5
3
C=2
Hx - 2L
Solució n particular
1
y =
2
0
Solució n singular y = 0
-1
-2
-4 -2 0 2 4
Video1: https://youtu.be/kF_O6DpwBvA
Video2: https://youtu.be/sB6mXqvfibc
Video3: https://youtu.be/UYIP3z5HlHo
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
(a) x2 + exy + xy dx + y2 + x dy = 0.
(a) Orden 1, no lineal. (b) Orden 2, lineal. (c) Orden 1, no lineal. (d) Orden 4, lineal. (e) Orden 6,
no lineal.
Respuestas
dy x2 + exy + xy
(a) =− .
dx x2 + y
5 0 8 ex cos x
(b) y00 = − y + 2y+ .
3x 3x x2
dy x 2 + y2
(c) =− x .
dx 3e − 2y
d4y d2y
dy
(d) = − tan x − ex + xy + e3x .
dx4 dx 2 dx
5 3 x ln x
(e) yvi = y000 − 2y00 − yy0 + .
4 4 4 (x + 1)
1. Comprobar si las siguientes funciones o relaciones son solución de la ecuación diferencial dada:
x4 1
b) La función y(x) = solución de y0 (x) − x · y 2 = 0 en ] − ∞, +∞[.
16
x 1 ln x ln2
e) La función ϕ(x) = + − + , para la ecuación diferencial x3 y000 + 6x2 y00 + 7xy0 + y = x.
8 x x x
f ) La relación x2 − 2y2 = 3, para la ecuación x − 2yy0 = 0.
Respuestas : (a) Sí. (b) Sí. (c) ϕ1 (x). sí, ϕ2 (x) sí. (d) No. (e) Sí, (f) Sí . (g) No.
2. Sustituir y = erx en la ecuación diferencial 5y00 − 3y = 0 y determina todos los valores de r para
los que ϕ(x) = erx es solución de dicha ecuación.
q
Respuesta: r = ± 35
Respuestas:
4. Sabiendo que las relaciones (o ecuaciones) dadas definen a y como función implícita de x, com-
pruebe que cada una de éstas es solución implícita de la correspondiente ecuación diferencial.
y
a) ln(xy) − yey = C · y, ED: y0 = .
xyey (1 + y) − x
b) C · y2 = xy +C, ED: (xy2 + x) · y0 = y3 − y
−1
5. La expresión y(x) = define una familia 1-paramétrica de soluciones de la ecuación diferen-
C + 3x
0 2
cial y = 3y ¿Hay algún valor de C para el que y(0) = 0?
Respuesta
No
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
Respuestas
No; Si y(x) = 0.
7. Compruebe que (x2 − 2x + 2)(y2 + 1)e2 arctan(y) = C, donde C es una constante arbitraria, es
solución general en forma implícita de la ecuación (xy2 − y2 + x − 1)dx + (x2 y − 2xy + x2 + 2y −
2x + 2)dy = 0.
y
8. Si curva con ecuación y2 − xy = A es una solución implícita de la ecuación diferencial y0 = ,
2y − x
halle una solución particular que satisfaga la condición y(1) = 2.
Respuesta : y2 − xy = 2
9. Sabiendo que las siguientes funciones (relaciones o ecuaciones) corresponden a la solución general
de una determinada ecuación diferencial,(A, B, C constantes arbitrarias) encuentre la ecuación
diferencial respectiva.
a) y = Ax + Bx2
b) y = Ax3 + Bx2 +Cx
c) y = Ax2 + Bex + x − 1
d) y = x + A sen x
e) y = Ax2 + A2
f) (x − A)2 + y2 = 4
g) y = Aex + Be−2x
Respuestas
(a) x2 y00 − 2xy0 + 2y = 0; (b) x3 y000 − 3x2 y00 + 6xy0 − 6y = 0 (c) falta (d) y0 = 1 + (y − x) cot(x)
(e) (y0 )2 + 2x3 y0 − 4x2 y = 0 (f) y2 (y0 )2 + y2 = 4, (g) y00 + y0 − 2y = 0
En algunos casos es importante no solo determinar la solución general de una ecuación diferencial sino
también se busca una solución y(x) de una ecuación diferencial de modo que y(x) satisfaga condiciones
adicionales prescritas, es decir, condiciones impuestas en la y(x) desconocida o sus derivadas. Para ello, se
expone la siguiente definición.
dny
= f (x, y, y0 , · · · , y(n−1) ) con , y(x0 ) = y0 , y0 (x0 ) = y1 , · · · , y(n−1) (x0 ) = yn−1 (1.5)
dxn
donde y0 , y1 , · · · , yn−1 son constantes reales especificadas de manera arbitraria, se llama problema
de valores iniciales (PVI). Los valores de y(x) y sus primeras n − 1 derivadas en un solo punto x0 :
El problema que se expone en 1.5 se llama también problema de valores iniciales de n-ésimo orden. Por
ejemplo, resolver
dy
= f (x, y) sujeta a y(x0 ) = y0 (1.6)
dx
y
d2y
= f (x, y, y0 ) sujeta a y(x0 ) = y0 , y0 (x0 ) = y1 (1.7)
dx2
son problemas de valores iniciales de primer y segundo orden, respectivamente. En el caso de la ecuación
1.6, se busca una solución y(x) de la ecuación diferencial y0 = f (x, y) en el intervalo I que contiene a x0 de
modo que su gráfica pase por el punto (x0 , y0 ), como se muestra en la figura 1.6.
En el caso de 1.7 se desea encontrar una solución y(x) de la ecuación diferencial y00 = f (x, y, y0 ) en un
intervalo I que contiene a x0 , de manera que su gráfica pase no solo por (x0 , y0 ), sino que la pendiente de la
curva en este punto sea el número y1 , como se muestra en la figura 1.7.
Las palabras condiciones iniciales derivan de sistemas físicos donde la variable independiente es el tiempo
t y donde y(t0 ) = y0 y y0 (t0 ) = y1 representan la posición y la velocidad, respectivamente, de un objeto en
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
La estrategia para resolver un PVI para una EDO se resume, (como se muestra en la ecuación 1.5) por
lo general, en hallar primero una familia de soluciones n paramétricas (solución general) de la ecuación
diferencial dada y luego usar las n condiciones iniciales en x0 , para determinar valores numéricos de las n
constantes en la familia. La solución particular resultante se define en algún intervalo I que contiene al
punto inicial x0 .
Para efectos de esta sección, se hará énfasis en los PVI para EDO de primer orden. Por ende, la estrategia
del párrafo anterior se resume en resolver la ecuación diferencial en términos de un parámetro, despejando
luego este con la condición inicial dada.
1.26
Al resolver el problema de valor inicial: y0 (x) = y(x), y (0) = 3, es fácil comprobar que y(x) = Cex ,
con C una constante arbitraria, representa la solución general de la ecuación diferencial y0 (x) = y(x)
en el intervalo ] − ∞, +∞[.
La figura 1.8, muestra la gráfica de la familia de curvas solución generadas a partir de dicha solución
general. Luego para hallar la solución particular cuya curva integral pasa por el punto (0, 3), se
sustituye la condición inicial y (0) = 3 en la solución general, de donde:
3 = Ce0 ⇔ C = 3,
y así de esta forma la solución particular buscada es y = 3ex , la cual se muestra en la figura 1.8.
y y
y = 3 ãx y = 3 ãx
4 4
(0,3)
2 2
-4 -2 2 4
x -4 -2 2 4
x
-2 -2
-4 -4
,
Figura 1.8: gráfica de familia de curvas solución de y0 (x) = y(x) ysolución del P.V.I y0 = y, y (0) = 3
Fuente: elaboración propia con M ATHEMATICA 9.0.
Se puede comprobar que la función ϕ(x) = senx − cosx es solución del problema de valor inicial:
y00 + y = 0 con y (0) = −1, y0 (0) = 1.
En efecto, note que para ϕ(x) = sen x − cos x, se tiene que ϕ 0 (x) = cos x + sen x y que ϕ 00 (x) =
− sen x + cos x. Luego, sustituyendo en la ecuación diferencial se comprueba que − sen x + cos x +
sen x − cos x = 0; por tanto, ϕ(x) es solución de la ecuación diferencial. Pero además: ϕ(0) =
sen 0 − cos 0 = −1 y ϕ 0 (0) = cos 0 + sen 0 = 1, así por tanto, se verifican las condiciones iniciales
dadas al P.V.I, su representación gráfica se muestran en figuras 1.9.
y y
2 2
1 1
-4 -2 2 4
x -4 -2 2 4
x
-1 -1
-2 -2
Figura 1.9: gráfica de familia curvas solución y = A cos(x) + B sen(x) y gráfica de solución del P.V.I :y00 + y = 0 con
y (0) = −1, y0 (0) = 1.
Fuente: elaboración propia con M ATHEMATICA 9.0.
1.28
x4
Un PVI puede tener varias soluciones. Por ejemplo, las funciones y = 0; y = satisfacen la
16
dy 1
ecuación diferencial = xy 2 , y la condición inicial y(0) = 0, por lo cual, el problema de valor
dx
dy 1
inicial = xy , y(0) = 0, tiene al menos dos soluciones, como se ilustra en la figura 1.10, donde
2
dx
las gráficas de ambas funciones pasan por el mismo punto (0, 0).
Este ejemplo nos manifiesta que al tratar de resolver un problema de valores iniciales es posible que
existan más de una solución que satisfaga la condición inicial dada. En consecuencia, es deseable saber por
adelantado si existe una solución y , cuando es así, si la solución es única. Para ello, es necesario establecer
un teorema que dé las condiciones suficientes de existencia y unicidad de una solución de un problema de
valor inicial de primer orden, como se muestra en apartado más adelante.
d2y
= f (x, y, y0 ) sujeta a y(a) = y0 , y(b) = y1 (1.8)
dx2
se llama un problema de valores en la frontera (PVF). Los valores necesarios, y(a) = y0 ; y(b) = y1 , se
denominan condiciones en la frontera.
1.29
dT (t)
= β (T (t) − Ta )
dt
Al resolver dicha ecuación se necesita conocer la temperatura del objeto en dos instantes diferentes,
ya que hay dos constantes por determinar: la constante de proporcionalidad β y la constante de
integración. De esta manera un PVF que lo permite resolver es
dT
= β (T (t) − Ta ) sujeta a T (0) = T0 , T (t1 ) = T1 (1.9)
dt
La solución del problema de valor de frontera permite obtener la Ley de Variación de la temperatura
en función del tiempo ( esto es, una ecuación para T (t)).
a) y( π4 ) = 1 ; y( π2 ) = 2.
√
b) y(0) = 1 ; y( π8 ) = 2
Ante un problema de valor inicial surgen dos cuestiones fundamentales, ¿ existe solución al problema?, ¿es
única?. Antes de resolver un problema es interesante saber si tiene solución y si es única. Geometricamente,
equivale a preguntarse si de toda la familia de curvas integrales existe alguna que pase por el punto definido
por la condición inicial y si por dicho punto pasa una única curva.
Poder asegurar la existencia de una solución no implica que se sea capaz de hallarla, aunque sí es importante
conocer su existencia y unicidad; para ello se utilizará el siguiente teorema que nos garantiza la existencia
y unicidad de soluciones en un P.V.I para una ecuación diferencial ordinaria de primer orden.
dy
= f (x, y), y(x0 ) = y0 . (1.10)
dx
1.30
2y
Considerando la ecuación diferencial xy0 = 2y se observa que f (x, y) = de donde por TEU sólo
x
garantiza la existencia de soluciones en la parte del plano donde x 6= 0 y además por otra parte que
∂f 2y
como = − 2 la unicidad también se garantiza en la misma parte.
∂y x
1.31
1
En la ecuación diferencial y0 = p se tiene que f (x, y) = √ 21 2 es continua en
2
x −y 2 x −y
2 2 2 2
(x, y) : x > y . (x, y) : x = y está formado por las dos rectas que pasan por el origen: y = x,
y = −x la región del plano buscada es R = {(x, y) : |y| < x}
1.32
Considere el problema de valor inicial y0 (x) = y(x) con y (0) = 3. Se quiere determinar mediante el
teorema de existencia y unicidad para EDO1 si ¿existe solución única?.
∂f
Para ello, observe que tanto f (x, y) = y como = 1 son funciones continuas en R2 , por tanto
∂y
son continuas en un entorno de (0, 3), luego por el Teorema de existencia y unicidad (TEU) se puede
asegurar que existe un entorno de x0 = 0 donde existe solución y además cumple ser única. De hecho,
en ejemplo 3.19 se determinó que dicha solución es y = 3ex y también mostrada en la figura 1.8.
1.33
Al considerar el problema de valor inicial, y0 (x) = y(x) con y(0) = 0, se puede verificar que no
p
En efecto,
p en el ejemplo 3.20, se determinó que la solución total de la ecuación diferencial
0
p
y (x) = y(x) esta formada por la solución general y(x) = 12 (x + C) y la solución singular
y(x) = 0. Asimismo mediante la figura 1.5 se observa que por el punto (0, 0) pasan dos curvas
solución, y por tanto no hay unicidad.
√ 1 ∂f 1
y0 =
p
y, es decir, f (x, y) = y(x) = y(x) 2 , y derivando se tiene = p .
∂y 2 y(x)
∂f
Las funciones f (x, y) y son continuas en el semiplano y > 0, pero no son continuas en un
∂y
rectángulo que contenga a (0, 0). En cambio, sí podríamos asegurar que ∀ (x0 , y0 ) con y0 > 0
existe un intervalo centrado en x0 en el que el problema dado tiene solución única. Lo obtenido
√
Figura 1.12: región Rectangular R, para y0 = y con y(0) = 0.
Fuente: elaboración propia con Mathematica 9.0.
1.34
√
Considere el problema de valor inicial y0 = −3 + 3x + y con y(a) = b; a, b ∈ R.
a) Determinar la región R del plano real xy para la cual la ecuación diferencial dada admite solución
única que pase por el punto de la región R.
x +C 2
b) Sabiendo que la solución general de la ecuación diferencial propuesta es y = − 3x ¿ se
2
puede hallar una solución para y(−1) = 3? En caso afirmativo, ¿ es esta solución única? Justifique.
Solución
Por lo tanto, sí hay solución para el punto (−1, 3), la cual está dada por
x+1 2
y= − 3x.
2
√
Dicha solución NO es única, pues en y0 = −3 + 3x + y se puede verificar que tiene otra solución
que pasa por el punto (−1, 3) ; a saber, la función y(x) = −3x, la cual es una solución singular. De
√
donde finalmente, se determina que la ecuación diferencial y0 = −3 + 3x + y tiene al menos dos
soluciones que satisfacen la condición inicial y(−1) = 3, a saber
2
x+1
y= − 3x ∨ y = −3x.
2
Por otra parte, la gráfica dada en la figura 1.13 muestra los resultados obtenidos en dicho ejemplo.
√
Figura 1.13: región Rectangular R, para y0 = −3 + 3x + y con y(−1) = 3.
Fuente: elaboración propia con Mathematica 9.0.
√
2. Considere la ecuación diferencial y0 = x y con y(a) = b y mediante página interactiva
√
Teorema existencia y unicidad en y0 = x y con y(a) = b, ingrese diferentes valores para
a y b asimismo la región R del plano xy en donde se cumplen las hipótesis del TEU para visualizar
y comprobar lo que sucede en puntos del plano xy que están dentro y fuera dicha región R. Además
verifique lo reflejado en el ejemplo 3.21.
Respuestas: (a) No. (b) Sí. (c) No. (d) Sí. (e) Sí.
a) Determine y gráfique la mayor región del plano xy en donde el teorema de existencia y unicidad
(T.E.U) garantiza solución única.
√
b) Justifique si T.E.U garantiza que existe una solución que pase por el punto (1, 3).
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
p
7. Considere la ecuación diferencial y0 = 2 y − x2 + 2x.
a) Determine y gráfique la mayor región del plano xy en donde el teorema de existencia y unicidad
(T.E.U) garantiza solución única.
y(x0 ) = y0
a) Determine y gráfique la mayor región del plano xy para la cual la ecuación diferencial admite
solución única por el punto (x0 , y0 ) de la región R.
x +C 2
b) Sabiendo que la solución general en su forma explícita esta dada por y = 2x − .
2
¿ Se puede hallar una solución para y(−1) = −2?. En caso afirmativo ¿ Es esta solución única?.
Justifique.
y0 = 2 +
p
y − 2x + 3
tiene solución única para cada punto (x0 , y0 ) de la región. La ecuación anterior posee una solución
singular. Determínela. ? Afecta la existencia de tal solución a la unicidad establecida en el punto
anterior?
entonces
a) Determine la región R del plano xy sobre la cual podemos garantizar con certeza la existencia
y unicidad de soluciones. Haga una representación gráfica de la región R.
Sección 2:
Ecuaciones Diferenciales Ordinarias en
variables separables y homogéneas.
2 — Ecuaciones diferenciales variables separables y homogéneas
2.0.2 Contenidos
Las formas 2.1 y la dada por 2.2 son equivalentes tomando N(x, y) = 1 y M(x, y) = − f (x, y), pues
dy
M(x, y) + N(x, y) = 0 puede escribirse como
dx
dy −M(x, y)
= = f (x, y),
dx N(x, y)
con N(x, y) 6= 0.
dy
Por otra parte se tiene que M(x, y) + N(x, y) = 0 puede escribirse de forma equivalente como
dx
dy
pues multiplicando M(x, y) + N(x, y) = 0 por dx, se obtiene
dx
dy
M(x, y) + N(x, y) dx = 0 ⇔ M(x, y)dx + N(x, y)dy = 0.
dx
Esta última forma 2.3, se utiliza ampliamente y tiene ventaja de ser simétrica en x e y en el sentido de que
cualquiera de ellas puede ser tratada como variable independiente, según las circunstancias.
Definición 2.1
Se dice que una ecuación diferencial de primer orden y0 = f (x, y), es de variables separables si
puede ser escrita de la forma y0 = g(x)h(y), donde g(x) y h(y) son funciones continuas de x e y
respectivamente en cierto intervalo común I.
Al dividir entre la función h(y), una ecuación separable se puede escribir en la forma
dy
p(y) = g(x) (2.4)
dx
1
con donde, por comodidad, p(y) representa a (h(y) 6= 0).
h(y)
Si y = φ (x) representa una solución de 2.4, se debe cumplir
Z Z
p(φ (x))φ 0 (x)dx = g(x)dx (2.5)
Z Z
p(y)dy = g(x)dx (2.6)
De donde se tiene que solución general para este tipo de ecuaciones diferenciales viene dada por
Z
F(y) = g(x)dx +C
.
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
2.1
dy
Considere la ecuación diferencial y2 = x, determine si es de variable separable y encuentre la
dx
solución general.
Note que dicha ecuación diferencial es de variable separable pues esta dada de acuerdo a la forma
y3
2.4. Así luego, dado que una antiderivada de q(y) = y2 es se cumple que
3
3
d[ y3 ]
= x,
dx
de donde al integrar a ambos miembros con respecto a variable independiente x se obtiene que la
solución general viene dada de forma implícita por
y3 x 2
= +C.
3 2
Definición 2.2
Una ecuación diferencial de la forma diferencial: M(x, y)dx + N(x, y)dy = 0 es una ecuación sepa-
rable o de variables separables si se puede escribir de la forma:
1
Ahora si la ecuación diferencial 2.7, se multiplica por , con (g(y)h(x) 6= 0):
g(y)h(x)
f (x) k(y)
dx + dy = 0, (2.8)
h(x) g(y)
donde separando la variables: m(y)dy = −n(x)dx,, luego integrando a ambos lados de la igualdad:
Z Z
m(y)dy = −n(x)dx,
de esta forma se obtiene la solución de 2.7 en su forma implícita viene dada por :
Observación 2.1 En el proceso de resolución se pueden perder soluciones con las manipulaciones alge-
braicas. Por ejemplo, el proceso al pasar de 2.7 a 2.8, es válido si h(x)g(y) 6= 0. Sin embargo, para los casos
cuando h(x) = 0 y g(y) = 0 se deben comprobar en la ecuación original si alguna de ellas es o no solución,
y en caso de serlo, añadirla a la solución general.
Comprobar mediante las definiciones dadas anteriormente que las siguientes ecuaciones diferenciales
son de variables separables:
dy dy 2x + xy
(a) = x2 + x2 y3 . (b) = 2 (c) y0 = 2 − xy.
dx dx y + 1.
(d) cos xey dx + (x + 1) dy = 0. (e) (xy2 − y2 + x − 1)dx + (x2 y − 2xy + x2 + 2y − 2x + 2)dy = 0.
En efecto, fácilmente utilizando las definiciones aportadas para ecuaciones diferenciales en variables
separable se comprueba que:
dy
= x2 1 + y3 , luego sí es una ecuación separable.
(a) Esta ecuación se puede reescribir como:
dx
dy x(3 − y)
(b) Esta ecuación se puede reescribir como = 2 , por tanto, el término de la derecha se
dx y +1
dy 3−y
puede escribir como el producto de una función de x por una función de y: = x· 2 luego sí
dx y +1
es una ecuación separable.
2.3
dy x2 + 1
Encontrar la solución general de = y determinar la solución particular para la cual y = 4
dx 2−y
cuando x = −3.
Separando las variables, se puede escribir la ecuación dada en la forma (2 − y)dy = (x2 + 1)dx,
luego integrando la ecuación anterior se obtiene fácilmente la solución general de la siguiente forma
y2 x 3
Z Z
(2 − y)dy = (x2 + 1)dx ⇔ 2y − = + x +C. (2.10)
2 3
Posteriormente, sustituyendo x = −3 y y = 4 en 2.10 se obtiene C = 12, por lo que la solución
particular requerida es
−x3 y2
− x + 2y − = −12. (2.11)
3 2
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
Observación 2.2 Se debe notar que para llegar a la solución general de la EDO del ejemplo anterior, se
−x3 y2
debe considerar que y 6= 2. Gráficamente, las curvas solución en − x + 2y − = C se observan en la
3 2
figura 2.1 y asimismo que tienen una pendiente vertical en el punto donde y = 2.
y
8
Valores para C
C=- 5 C=3
4
C=- 4 C=4
C=- 3 C=5
y=2 C=- 2
2
C=- 1
C=0
C=1
-2 2
x C=2
-2
-4
dy x2 + 1
Figura 2.1: algunas soluciones particulares de la EDO = .
dx 2−y
Fuente: elaboración propia con M ATHEMATICA 9.0.
y
6
H-3, 4L 4
Funciones
1
6- 6 - 30 - 3 x - x3
3
2 1
6+ 6 - 30 - 3 x - x3
3
2
-4 -3 -2 -1
x
-2
dy x2 + 1
Figura 2.2: solución particular de la EDO que pasa = por (−3, 4).
dx 2−y
Fuente: elaboración propia con M ATHEMATICA 9.0.
2.4
dy
Resuelva el PVI xy = yx + x + y + 1, con y(1) = 0, x 6= 0.
dx
y x+1
Z Z
dy = dx
y+1 x
y − ln|y + 1|+ C1 = x + ln|x|+ C2
y − ln|y + 1| = x + ln|x|+ C
y = x + ln|(y + 1)x|+ C (∗)
2.5
dy x2
= , y(0) = 1
dx 1 − y2
1
Z
(1 − y2 )dy = y − y3 + C1
3
1
Z
x2 dx = x3 +C2
3
Por lo que la solución está dada por:
1 1
y − y3 − x3 = C(∗)
3 3
donde C = C1 −C2 es una constante arbitraria. Ahora si utilizamos en en (*) la condición y(0) = 1
2
para obtener la solución del problema de valor inicial. obtenemos C = .
3
Por lo que la solución del problema de valor inicial viene dada por :
1 1 2
y − y3 − x 3 =
3 3 3
Observe que aunque el TEU no garantiza una solución al problema de valores iniciales en el punto
(0, 1) (pues el denominador se anula), esta sí existe.
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
2.6
dy
Resuelva la ecuación diferencial x = 2y hasta obtener su solución general.
dx
Si en la ED dada separamos variables (y 6= 0), e integramos a ambos lados se tiene
ln |y| = 2 ln |x| + K,
y
⇒ ln 2 = K
y x
⇒ 2 = eK ,
x
⇒ |y| = x2 eK
Por lo que (*) representa la solución general de la ED dada. Por otro lado si analizamos el caso
cuando y(x) = 0 se tiene que esta se desprende de (*) para la constante C = 0, esto indica no es
solución singular. Además se aprecia que y(0) = 0 para todas las soluciones, lo cual se puede apreciar
en la figura 2.3.
dy
Figura 2.3: soluciones de la EDO x = 2y por (0, 0).
dx
Fuente: elaboración propia con M ATHEMATICA 11.0.
|c1 |− sen x − 1
ey =
sen x + 1
es decir,
K − sen x
y = ln (∗)
sen x + 1
con K = |c1 |−1. Se quiere ahora una solución que cumple con y(π) = 0, entonces
K − sen 0
0 = ln = ln K
sen 0 + 1
de donde
K=1
Sustituyendo en (*) obtenemos la solución del problema de valor inicial viene dada por
1 − sen x
y(x) = ln
sen x + 1
2.8
x−1 y+1
dx + dy = 0
x2 − 2x + 2 y2 + 1
x−1 y+1
2
dx = − 2 dy
x − 2x + 2 y +1
1 1
⇒ ln |x2 − 2x + 2| = − ln |y2 + 1| + arctan y + k
2 2
1 1
⇒ ln(x2 − 2x + 2) + ln(y2 + 1) + arctan y = k
2 2
⇒ ln[(x2 − 2x + 2)(y2 + 1)] + 2 arctan y = 2k
2 −2x+2)(y2 +1)]+2 arctan y
⇒ eln[(x = e2k
⇒ (x2 − 2x + 2)(y2 + 1) · e2 arctan y = e2k
2.9
t +1 2t 2 + 2t
Z Z
· 2tdt = dt
t2 + 1 t2 + 1
2t 2
Z
2t
= + dt
t2 + 1 t2 + 1
Z
2 2t
= 2− 2 + 2 dt
t +1 t +1
= 2t − 2 arctant + ln(t 2 + 1) +C1
Esta última expresión representa la solución general de la ecuación diferencial en forma implícita.
2.10
dy
Resolver x − y = 2x2 y y determinar el dominio máximo de las soluciones de la EDO dada.
dx
dy
Al multiplicar por dx la ecuación x − y = 2x2 y se tiene
dx
xdy − ydx = 2x ydx ⇔ (2x2 y + y)dx − xdy = 0,
2
2x2 + 1
1 1 1
y(2x2 + 1)dx − xdy = (0) ⇔
dx − dy = 0.
xy xy x y
En tal caso, x 6= 0 y y 6= 0. Luego, se tiene que la solución está dada por
2x2 + 1
Z Z Z
1 1 1
Z
dx − dy = C ⇔ 2x + dx − dy = C,
x y x y
⇔ x2 + ln |x| − ln |y| = C. (2.12)
Esto también puede escribirse en una forma libre de logaritmos escribiendo 2.12 como
2 2
x
x + ln |x| − ln |y| = C ⇔ x + ln = C
y
x
⇔ ln = C − x2
y
x
⇔ = eC−x2
y
x
⇔ = eC e−x2
y
x 2
⇔ = ±eC e−x
y
2
⇔ x = ±yeC e−x
2
⇔ xex = ±yeC
2
⇔ ± e−C xex = y
2
⇔ y = Axex . (2.13)
Note que el caso y = 0 esta contenido en la solución general para el caso A = 0.
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
2.11
Resolver sen(x + y)dx + sen y(csc xdy − cos xdx) = 0, en un intervalo I que satisfaga el teorema de
existencia y unicidad.
Al considerar la ecuación diferencial sen(x + y)dx + sen y(csc xdy − cos xdx) = 0 y efectuar
simplificaciones previamente se tiene
⇔ (sen x cos y + cos x sen y)dx + (sen y csc xdy − sen y cos xdx) = 0
⇔ (sen x cos y + cos x sen y − sen y cos x)dx + (sen y csc x)dy = 0
En algunos casos se pueden hacer cambios de variable apropiados que permiten transformar una ecuación
dy
de la forma = f (x, y), en otra de variables separables, como se ilustra en el siguiente ejemplo.
dx
2.12
dy p
= 3 3 y − 3x. (2.14)
dx
La raíz cúbica en la ecuación diferencial sugiere cambiar la variable dependiente, por ejemplo y a z,
de modo tal que la raíz pueda evitarse.
dy dz
= 3 + 3z2 ,
dx dx
de donde sustituyendo en la ecuación 2.14 se tiene
dz
3 + 3z2 = 3z,
dx
es decir,
dz
z2 = z − 1,
dx
al separar variables
z2
dz = dx,
z−1
posteriormente integrando a ambos lados
z2
+ z + ln | z − 1 |= x + K.
2
√
Finalmente como z = 3
y − 3x, se concluye que
q p p
3
(y − 3x)2 + 2 3 y − 3x + 2 ln | 3 y − 3x − 1 |= 2x +C. (C = 2K)
2.13
dy y − xy2
Verifique que la ecuación diferencial = se puede transformar en una ecuación diferencial
dx x + x2 y
de variables separables al aplicar la transformación zx−1 = y (no resuelva esta última ecuación).
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
dz −1 dy
Observe que: zx−1 = y ⇒ x − x−2 z = , luego de esta forma
dx dx
dy y − xy2 dz −1 zx−1 − x · (zx−1 )2
= ⇒ x − x−2 z =
dx x + x2 y dx x + x2 · zx−1
dz −1 zx−1 + x−1 z
⇒ x =
dx x + xz
dz −1 2zx−1
⇒ x =
dx x(1 + z)
1+z 2dx
⇒ dz = −→ Ecuación diferencial en variables separables z(x)
z x
2.14
√
Resolver la siguiente ecuación diferencial 2 xy − y dx − xdy = 0, mediante un cambio de
variable apropiado.
La raíz nos sugiere hacer el cambio de variable u2 = xy, así, u y y son funciones que dependen de x.
du dy
luego derivando con respecto a x se tiene que 2u = y+x ,
dx dx
dy dy 2u du −y
despejando se obtiene = dx , (x 6= 0).
dx dx x
√
Ahora sustituyendo en la ecuación 2 xy − y dx − xdy = 0, pero expresada en su forma normal
√
2u du
dx − y 2 u2 − y
= ,
x x
du
es decir, 2u = 2u, de donde separando variable e integrando a ambos miembros de la ecuación
dx
u = x + c,
finalmente devolviendo el cambio de variable u2 = xy se obtiene por solución general en su forma
implícita
√
± xy = x + c.
dy
a) y = e−4x + e−2x
dx
6x2 + 1
b) y0 =
2y + 3
c) − csc(y)dx + (x2 + 1)dy = 0
d) y0 = x2 (y2 + 1)
y cos x
(c) y0 =
1 − 2y
(d) e−y x2 + x2 + 2 yy0 = 0
dy
(g) (1 + x2 + y2 + x2 y2 ) =1
dx
(h) y0 + y2 − y = 0
(j) e−x + y0 = √ 1
x2 +1
+ 6x con y(0) = e
Respuestas
(f) y = arctan x +C
1
(g) y =
1 −Cx
C · ex
(h) y =
1 +C · ex
−e3x + 4
(i) y =
3
(j) y = sen h−1 + 3x2 + e−x + e + 1
dy y 2x2 + 1
2
b) = , y (1) = 1 SP:: y = xe(x −1)
dx x
4. Muestre que en la ecuación diferencial ydx + (1 + y2 e2x )dy = 0 el cambio de variables u = yex
separa variables, luego obtenga la solución general de dicha ecuación.
u3 −1
Respuesta: ED en variables separables u0 = 2
. SG: + ln |yex | = x +C. SS: y = 0.
1+u 2y2 e2x
dy
5. Considere la ecuación diferencial = f (ax + by + c); b 6= 0.
dx
a) Demuestre que esta ecuación se puede transformar a una de variables separables mediante la
sustitución z = ax + by + c.
1) y0 = sen2 (x − y + 1)
2) y0 = 2xy − x2 − y2 .
1
3) y0 = +1
ln(2x + y + 3)
Respuesta
2k+1
1. arctan(x + y) = x +C; 2. SG: tan(x − y + 1) = x +C SS: x − y + 1 = 2 π con K ∈ Z
a puede ingresar de manera directa mediante el link https://tecdigital.tec.ac.cr/revistamatematica/Libros/
Existen algunas ecuaciones diferenciales de variables NO separables las cuales mediante un cambio ade-
cuado de variables, se pueden transformar en ecuaciones de variables separables, como las que se muestran
en las sesiones siguientes.
Definición 2.3
Una ecuación diferencial en la forma m(x, y)dx + n(x, y)dy = 0 es homogénea si y solo si las
funciones m(x, y) y n(x, y) son homogéneasa y del mismo grado.
a Una función f (x, y) es homogénea de grado n, con n ∈ R si cumple la propiedad f (tx,ty) = t n f (x, y).
2.15
La ecuación (x2 + y2 )dx + (x2 − xy)dy = 0 es homogénea de grado dos, pues tanto la función
m(x, y) = x2 + y2 como la función n(x, y) = (x2 − xy) son funciones homogéneas y de grado dos
ambas, ya que m(tx,ty) = t 2 (x2 + y2 ) = t 2 m(x, y) y n(tx,ty) = t 2 n(x, y).
Definición 2.4
Una ecuación diferencial y0 = f (x, y), es homogénea si la función f (x, y) es homogénea de grado
cero, es decir, si cumple que f (tx,ty) = f (x, y).
2.16
2x + 3y
La ecuación diferencial y0 = es homogénea de grado cero pues,
x−y
2tx + 3ty t(2x + 3y) 0 2x + 3y 2x + 3y
f (tx,ty) = = =t = = f (x, y).
tx − ty t(x − y) x−y x−y
Definición 2.5
dy y dx x
Una EDO de primer orden es homogénea si tiene la forma =F o la forma =F . En
dx x dy y
x y
estas se cumple que tanto F como F son funciones homogéneas y de grado cero.
y x
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
Observación 2.3 Cuando la ecuación diferencial esta escrita de la forma normal m(x, y)dx + n(x, y)dy =
0, suponiendo que n(x, y) 6= 0 conviene escribirla en su forma diferencial
dy m(x, y)
=− .
dx n(x, y)
Algunas operaciones algebraicas tienen validez solo bajo ciertas condiciones, por ejemplo, pasar de
dy xy + y2
(x2 − y2 )dy + (xy + y2 )dy = 0 a =− 2 es válido si y solo si x2 − y2 6= 0. Los casos x2 − y2 =
dx x − y2
0 ⇔ y = ±x deben estudiarse por aparte, pues pueden generar soluciones singulares.
Teorema 2.1
Si la ecuación diferencial de primer orden y0 = f (x, y) es homogénea, entonces el cambio de variable
y
u = o equivalentemente y = u · x la reduce a una ecuación diferencial en variables separables.
x
Demostración
2.17
p
La siguiente ecuación diferencial xy0 = x2 − y2 + y, se transforma en una ecuación de variables
y
separables al utilizar cambio de variable u = .
x
Primero dividiendo por x, suponiendo x 6= 0, se tiene
r y 2 y
0
y = 1− + ,
x x
y
luego, haciendo la sustitución u = , se obtiene
x
p du dx
u0 x + u = 1 − u2 + u ⇔ √ = .
1−u 2 x
Integrando a ambos miembros de la ecuación
arc sen(u) = ln |x| + ln |C|,
En términos generales, si y0 = f (x, y) representa una EDO y (x, y) es un punto en el dominio de f (x, y),
entonces se puede mostrar que las siguientes propiedades son equivalentes:
2. f (tx,ty) = t 0 f (x, y) = f (x, y) para todo t tal que (tx,ty) pertenezca al dominio de f . En este caso, se
dice que la EDO de primer orden es homogénea de grado 0.
y
3. f (x, y) = F o f (x, y) = G xy para alguna función F o G de variable real.
x
2.18
dy y2 + yx y 2 y
=− = +
dx x2 x x
y du dy
de esta forma tomamos u = ⇒ x + u = , luego sustituyendo en la ED dada:
x dx dx
du
x + u = u2 + u
dx
du
x = u2
dx
−u−1 = ln|x|+C
y
como u =
x
x
− = ln|x|+C (∗)
y
Por lo tanto la solución general de ED viene dad por (∗) y además note que las soluciones no están
definidas para x = 0 (ni para y = 0).
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
2.19
y
Resuelva la ecuación diferencial x · y0 − y ln − y = 0 hasta obtener su solución general.
x
y
Tomando = u ⇒ y0 = u + u0 x. Luego sustituyendo en la ED dada se tiene que
x
u + u0 x − u ln u − u = 0
u0 1
⇒ =
u ln u x
du
Z
⇒ = ln|x|+ ln|C|
u lnyu
⇒ ln ln = ln|C · x|
x
2.20
y
Escribiendo la ecuación diferencial de la forma y0 = F , es decir,
x
y y 2
dy e− x + x
= y
dx x
y dy dv
de donde se toma el cambio de variable v = ⇔ y = xv ⇒ = x +v
x dx dx
y sustituyendo en ED dada
dv e−v + v2
x +v =
dx v
Separando variables
dx
vev dv =
x
e integrando
vev − ev = ln |x| + c
finalmente retomando las variables originales, se tiene que la solución en forma implícita viene dada
por :
y y
ye x − xe x − x ln |x| = cx
Mostrar que (x2 − y2 )dx + (xy + y2 )dy = 0 es una EDO homogénea y hallar su solución.
En primer lugar, como las funciones m(x, y) = (x2 − y2 ) y n(x, y) = (xy + y2 ) son ambas funciones
homogéneas del mismo grado, se cumple que dicha ecuación diferencial es homogénea. Además
por otro lado si se expresa la ecuación diferencial en su forma normal, es decir, se encuentra f (x, y).
Esto es:
y
1. Reescribir f (x, y) como F , siendo F una función real de variable real. En efecto
x
y2 −x2 y2 y 2
y2 − x 2 x2 x2
−1 x −1 y
f (x, y) = = xy+y2
= 2 = 2
= F (x 6= 0).
xy + y2 y
+ y2 y
+ y x
2x x x x x
x
2. Reescribir f (x, y) como G siendo G una función real de variable real. En efecto
y
2
y2 −x2 2
1 − xy
y2 − x2 y2
1 − xy2
x
f (x, y) = = xy+y2
= x = =G (y 6= 0).
xy + y2 y +1
x
+ 1 y
y2 y
3. Mostrando que f (tx,ty) = t 0 f (x, y) = f (x, y) para todo t tal que (tx,ty) pertenezca al dominio
de f .
Con cualquiera de los tres alternativas, se muestra que la EDO dada es homogénea.
y2 − x 2
Para resolver la ecuación diferencial dada, considérela de la forma y0 = , suponiendo
xy + y2
xy + y2 6= 0 . Ahora si dividimos por x2 6= 0 se obtiene la ecuación en la forma
( xy )2 − 1
y0 = y y 2,
x +(x)
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
y
luego tomando la sustitución u = ⇒ y0 = u0 x + u y sustituyendo en la ecuación anterior se tiene
x
u2 − 1 du (1 + u3 )
u0 x + u = ⇔ x = − ,
u + u2 dx u(u + 1)
de donde resolviendo y separando variables se tiene
−u 1
du = dx,
u2 − u + 1 x
integrando a ambos miembros de la ecuación
tan −1 2u−1
√
1 3
− log u2 − u + 1 −
√ = C + log(x),
2 3
finalmente retomando el cambio de variable y = ux, se tendría la solución general en su forma
implícita
2y
2 tan−1 x√−1
1 y y 3
− log 2 − + 1 − √ = C + log(x).a
2 x x 3
En la figura 2.4 se muestran gráficas de curvas solución de dicha ecuación diferencial.
a Utilice página interactiva Homogéneas2 o ingrese de manera directa mediante el link https://tecdigital.tec.
ac.cr/revistamatematica/Libros/LibrosCDF/EDO-NOviedo/CDFEDO1-Internet-Noviedo2018.pdf o en su
defecto escaneando el código QR adjunto para verificar los resultados.
y
y = -x
4
2 Valorespara C
C=1
y=0
C=2
-4 -2 2 4 x C=3
C=4
C=5
-2
-4
y2 − x2
Figura 2.4: algunas soluciones particulares de la EDO y0 = .
xy + y2
Fuente: elaboración propia con M ATHEMATICA 11.0.
Determine el valor de a de tal forma que la ecuación 2xz3a dx + a(x2 z3a−1 − za−1 )dz = 0 se
convierta en homogenéa. Luego resuelva hasta obtener una solución general.
Observe que:
N(x, z) = a(x2 z3a−1 − za−1 ) ⇒ N(tx,tz) = a(t 1+3a x2 z3a−1 − t a−1 za−1 ) [2]
De [1] y [2] se tiene que para la ED dada sea homogénea debe darse
t 1+3a = t a−1 ⇔ 1 + 3a = a − 1 ⇔ a = −1
Ahora si sustituimos a = −1 en la ED se tendría 2xz−3 dx − (x2 z−4 − z−2 )dz = 0 una ecuacion
diferencial homogénea que puede ser escrita como
dz 2xz−3 2xz 2 · xz
= 2 −4 −2 = 2 2 =
dx x z − z x −z 1 − ( xz )2
z
Tome de esta forma u = ⇒ z = ux ⇒ z0 = u0 x+u, luego sustituyendo en la ecuación anterior se tiene
x
2u 0 2u 0 u3 + u
u0 x + u = ⇒ u x = − u ⇒ u x = , luego separando variables se tendría
1 − u2 1 − u2 1 − u2
1 − u2 1
3
du = dx
u +u x
⇒ ln(u) − ln u2 + 1 = ln x + lnC (suponga u > 0, x > 0, C > 0)
u
⇒ ln 2 = lnC.x
u +1
u
⇒ 2 = C.x
u +1
z
De donde devolviéndose a las variables originales z, x con u = se obtiene por solución general
x
zx
=Cx
z2 + x 2
2.23
Existen algunos casos donde la ecuación diferencial puede transformarse a una homogénea con una
sustitución apropiada. En tales circunstancias, los ejemplos son específicos u obedecen a una estructura.
Para ilustrar este hecho, considere la EDO
dy ax + by + k
=G (2.17)
dx Ax + By + K
Primer Caso
u = x−α aα + bβ + k = 0
con α, β , solución del sistema (2.18)
v = y−β Aα + Bβ + K = 0
transformará la EDO en homogénea. En efecto, a partir de 2.18 se despejan las variables x y y obteniendo
x = u+α
(2.19)
y = v+β
dx = du
(2.20)
dy = dv
Por otro lado, note que en el paso 2.21 se tiene la condición de que ambas expresiones se anulen, lo
cual es posible resolviendo el sistema
aα + bβ + k = 0
(2.23)
Aα + Bβ + K = 0
o equivalente el sistema
aα + bβ = −k a b α −k
⇔ · = (2.24)
Aα + Bβ = −K A B β −K
Luego, mediante la Regla de Cramer, el sistema de ecuaciones tiene una única solución para α y β dada
por
−k b a −k
−K B −kB + Kb A −K −aK + Ak
α= = y β= = . (2.25)
a b aB − Ab a b aB − Ab
A B A B
Segundo Caso
Si aB − Ab = 0 y G es una función real continua en una región abierta, es posible el despeje de a o b en los
casos en que A 6= 0 o B 6= 0. De esta manera, se puede considerar la sustitución
v = Ax + By, (2.26)
en donde
dv dy 1 dv A dy
= A+B ⇔ · − = . (2.27)
dx dx B dx B dx
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
En primer lugar, se encuentra f (x, y). Esto es se transforma dicha ecuación diferencial en su forma
normal:
Se sigue que 2.39 representa una EDO homogénea, la cual se transforma a una de variables separables
mediante la sustitución
v
V= , (2.40)
u
esto permite cambiar v en términos de V manteniendo a u como variable independiente. En conse-
cuencia
dv
dV du u − v dV dv v dv dV
= ⇒u = − ⇔ =u +V. (2.41)
du u2 du du u du du
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
dv 2 − uv dV 2 −V
= v ⇔u +V = . (2.42)
du 1 − 2 u du 1 − 2V
dV 2 −V dV 2 −V
u +V = ⇔ u = −V
du 1 − 2V du 1 − 2V
dV 2V 2 − 2V + 2
⇔ u =
du 1 − 2V
u dV 2
V −V + 1
⇔ =
2 du 1 − 2V
1 − 2V 2
⇔ 2 dV = du
V −V + 1 u
1 − 2V 2
Z Z
⇒ 2
dV − du = ln A
V −V + 1 u
⇔ − ln(V 2 −V + 1) − 2 ln |u| = ln A
⇔ − ln(V 2 −V + 1) = ln u2 + ln A
⇔ ln(V 2 −V + 1)−1 = ln(Au2 )
1
⇔ ln 2 = ln(Au2 )
V −V + 1
1
⇒ 2 = Au2 . (2.43)
V −V + 1
Recuperando la sustitución realizada en 2.40, se tiene ahora que
1 1
= Au2 ⇔ = Au2 ⇔ v2 − uv + u2 = C. (2.44)
V 2 −V +1 v 2
− uv + 1
u
1
con C = . Se realiza luego el cambio realizado en 2.38, para dar la solución obtenida en 2.44 en
A
términos de las variables iniciales x, y, lo cual solución general es
Observación 2.4 Puede verificar resultados obtenidos del ejemplo anterior, adémas complementar
teoría mediante página interactiva ecuación diferencial especial o ingresando de forma directa a la
dirección electrónica https://tecdigital.tec.ac.cr/revistamatematica/Libros/LibrosCDF/
EDO-NOviedo/CDFEDO1-Internet-Noviedo2018.pdf o en su defecto escaneando el código QR ad-
junto que lleva a página respectiva.
xy2 + 2x
a) y0 =
xy
b) xy2 dy − (x3 + y3 )dx = 0
y y
c) y0 = + cos2
x x
dy
d) −(x + y ) + (3x + y2 ) = 0
3 3
dx
2y
dy
e) y = e x (x + y)
x
2. Resuelva de manera algebraica las siguientes ecuaciones diferenciales, en caso de no estar la ED
dada en su forma normal conviene tenerla en dicha forma y0 = f (x, y). Posteriormente compruebe
lo obtenido mediante página interactiva Homogéneas2 o de forma directa escaneando el código
QR adjunto para comprobar dichos resultados en dicha página.
3xy + y2
a) y0 =
2x2
dy y3 + x3 + xy2
b) = .
dx xy2
y
0
c) xy = y + 2xe x .
d) (x − y)dx + (x + y)dy = 0.
y2 − x 2
e) y0 = .
xy + y2
y dy
f ) − y + x tan +x = 0 .
x dx
3. Resuelva las siguientes ecuaciones diferenciales.
−x
a) (y2 + xy)dx − x2 dy = 0. SG: = ln |x| +C
y
b) y(ln x − ln y)dx = (x ln x − x ln y − y)dy SG: x(ln x − 1) + (y − x) ln y = Cy
c) (x2 + 2xy)y0 = −2y2 − 3xy SG: x2 y2 + x3 y = C
d) (2xy + x2 + 3y2 )y0 + (y2 + 2xy + 3x2 ) = 0 SG: (y + x)(y2 + 2xy + 3x2 ) = C
x x x 1
e) x0 = + sec2 ( ). SG: + sen(2x) = ln |y| +C
y y 2 4
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
2x
dx −2x ey
f ) y = x + 4ye y SG: = 4 ln |y| +C
dy 2 s r
√ y
g) (2 xy − y)dx − xdy = 0. . SG: C · x = 1 +
x
p p dy 3
h) (x3 + y2 x2 + y2 ) − xy x2 + y2 · =0 SG: x3 ln(Cx3 ) = (x2 + y2 ) 2
dx
i) (2x + 2y − 1) + y0 (x + y − 2) = 0 (Sug. tome u = x + y) SG: 3 ln |x + y +
1| −2x −
y = C
y 1
j) x2 y0 = xy + x2 e x ; x > 0, SG: y = −x ln ln +C
x
4. Considere la ED 2x4 yy0 + y4 = 4x6 , si al aplicar el cambio de variable y = zn para cierto n la
convierte en ED homogénea.sDetermine el valor de n y solución general de la misma.
3 3 4C + x5
,R/ n = , SG: y = −x 2
2 x5 −C
0 ax + by + h
5. Considere la ecuación diferencial y = f .
dx + cy + k
a) Si ac − bd 6= 0, las rectas ax + by + h = 0 y dx + cy + k = 0 se intersecan en un único punto
(x0 , y0 ). Muestre que, en este caso , la sustitución u = x − x0 y v = y − y0 transforma la
ecuación en homogénea en las variables u y v.
6. Aplique ejercicio anterior para resolver las siguientes ecuaciones diferenciales y posteriormente
verifique resultados obtenidos por medio de la página interactiva ecuación diferencial especial
o de forma directa escaneando el código QR adjunto para comprobar dichos resultados en dicha
página.
6x + 3y − 3
a) y0 =
4x + 8y − 2
b) (2x − y + 2)dx + (4x − 2y − 1)dy = 0
x−y−1
c) y0 =
x + 3y − 5
d) (2x − y − 4)dx − (x + y − 1)dy = 0
e) (x − y + 1)dy − (x + y − 1)dx = 0
3x − y + 2
f ) y0 =
6x − 2y
0 12x + 5y − 9
g) y =
−5x − 2y + 3
h) (x + y − 4)dx + (x − y + 2)dy = 0.
dy x+y+7
i) =− .
dx −2(x + y) + 1
Sección 3:
Ecuaciones Diferenciales Ordinarias exactas y
con factor integrante.
3 — EDO Exactas y con factor integrante.
3.0.2 Contenidos
Definición 3.1
Una ecuación diferencial de primer orden
dy
M(x, y) + N(x, y) o equivalentemente M(x, y)dx + N(x, y)dy = 0 (3.1)
dx
es exacta en un rectángulo R, si existe una función F(x, y) tal que el diferencial total de F(x, y)
denotado dF(x, y) a sea igual a M(x, y)dx + N(x, y)dy, es decir, si existe una función F(x, y) tal que
∂ F(x, y) ∂ F(x, y)
dF(x, y) = dx + dy = M(x, y)dx + N(x, y)dy,
∂x ∂y
∂ F(x, y) ∂ F(x, y)
de donde se extrae, en tal caso que = M(x, y) y = N(x, y), ∀(x, y) perteneciente
∂x ∂x
a la región rectangular R.
a En cálculo diferencial, se tiene que si z = f (x, y) es una función con primeras derivadas parciales continuas en una
∂ f (x, y) ∂ f (x, y)
región R del plano xy, su diferencial (que también llamamos diferencial total) es dz = dx + dy
∂x ∂y
De acuerdo a lo anterior, como dF(x, y) = 0, la solución de la ecuación diferencial que se obtiene por
integración, viene dada en forma implícita por F(x, y) = C, con C una constante arbitraria.
El siguiente teorema nos da una condición necesaria y suficiente para conocer cuándo una ecuación es
exacta y su demostración nos proporciona un método para obtener la solución general:
F(x, y) = C.
Teorema 3.1
Sean M(x, y) y N(x, y) funciones continuas con derivadas parciales de primer orden continuas en un
rectángulo R. Entonces, la ecuación M(x, y)dx + N(x, y)dy = 0 es exacta si y sólo si se verifica:
∂ M (x, y) ∂ N (x, y)
= ; ∀(x, y) perteneciente a la región R. (3.2)
∂y ∂x
Demostración.
(⇒) Supongamos que la ecuación M(x, y)dx + N(x, y)dy = 0 es exacta. Entonces, existe una función F(x, y)
tal que:
∂ F (x, y) ∂ F (x, y)
= M(x, y) y = N(x, y);
∂x ∂y
por tanto:
Puesto que las primeras derivadas parciales de M y N son continuas en R, también lo son las derivadas
parciales segundas cruzadas de F ; por tanto, éstas son iguales y se tiene que:
∂ M(x, y) ∂ N(x, y)
= , ∀(x, y) ∈ R.
∂y ∂x
(⇐) Dada la ecuación M(x, y)dx + N(x, y)dy = 0, suponga que My (x, y) = Nx (x, y). Hay que demostrar que
∂ F(x, y) ∂ F(x, y)
existe una función F tal que dx + dy = M(x, y)dx + N(x, y)dy, osea, F debe satisfacer
∂x ∂y
simultáneamente dos condiciones a saber:
∂ F(x, y) ∂ F(x, y)
= M(x, y) y = N(x, y).
∂x ∂y
∂ F(x, y)
Si = M(x, y) entonces:
∂x
Z
F(x, y) = M(x, y)dx = G(x, y) + ϕ(y), (3.3)
donde G(x, y) es una primitiva de M(x, y) respecto de x. Ahora, derivamos parcialmente respecto de y la
expresión obtenida para F y la igualamos a N(x, y) :
∂ F(x, y) ∂ G(x, y)
= + ϕ 0 (y) = N(x, y),
∂y ∂y
Si esta expresión sólo depende de y, podremos integrar y obtener ϕ(y) que, sustituida en 3.3 nos dará la
expresión de F(x, y).
Queda por demostrar que 3.4 sólo depende de y. Para ello, comprobamos que su derivada parcial respecto
de x es cero:
∂ ∂ G(x, y) ∂N ∂ ∂G ∂N ∂ ∂G ∂N ∂M
N(x, y) − = − = − = − = 0.
∂x ∂y ∂x ∂x ∂y ∂x ∂y ∂x ∂x ∂y
Método de resolución
donde F(x, y) se obtiene siguiendo los pasos vistos en la demostración del teorema 4.2, con C es una
constante arbitraria.
Solución
∂ 2 ∂
(x + y2 + 2x) = 2y = (2xy + 3y2 ).
∂y ∂x
Por tanto, la solución general viene dada por F(x, y) = C, siendo F una función tal que:
∂F
= x2 + y2 + 2x = M(x, y), (3.5)
∂x
y
∂F
== 2xy + 3y2 = N(x, y). (3.6)
∂y
Se calcula F(x, y). Integrando la ecuación 3.5 respecto de x se tiene de esta manera:
x3
Z
F(x, y) = (x2 + y2 + 2x)dx = + y2 x + x2 + ϕ(y).
3
Derivando esta expresión respecto de y, se tiene que:
∂F
= 2xy + ϕ 0 (y).
∂y
Igualando esta ecuación a la ecuación 3.6:
Se puede tomar k = 0, ya que en la solución final de la ecuación diferencial esta constante quedará
suprimida en la constante C. Por tanto, la función F buscada es:
x3
F(x, y) =+ y2 x + x2 + y3 ,
3
y la solución general de la ecuación diferencial es:
x3
+ y2 x + x2 + y3 = C.
3
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
∂ ∂
(cos x sen x − xy2 ) = −2xy = (y(1 − x2 )).
∂y ∂x
Por tanto, la solución general viene dada por F(x, y) = C, siendo F una función tal que:
∂F
= cos xsenx − xy2 , (3.7)
∂x
y
∂F
= y(1 − x2 ). (3.8)
∂y
En este caso, es más sencillo comenzar integrando la ecuación 3.8 respecto de y:
(1 − x2 )y2
Z
F(x, y) = y(1 − x2 )dy = + Ψ(x).
2
Derivando ahora esta expresión respecto de x, se tiene
∂F
= −xy2 + Ψ0 (x).
∂x
Igualando esta ecuación a la ecuación 3.7 se tiene:
(1 − x2 )y2 cos2 x
F(x, y) = − ,
2 2
y la solución general de la ecuación diferencial es:
(1 − x2 )y2 cos2 x
− = C, o bien, (1 − x2 )y2 − cos2 x = C1 .
2 2
Imponiendo la condición inicial, y(0) = 2, se tiene que 4 = C1 y sustituyendo este valor de C1 , se
obtiene la solución del problema de valor inicial dado:
4 − cos2 x
(1 − x2 )y2 − cos2 x = 4 o bien, y2 = con 1 − x2 6= 0.
(1 − x2 )
Para comprobar y visualizar paso a paso dichos resultados de los ejemplos ante-
riores, puede ingresar a página interactiva Exactas o de forma directa a la direc-
ción electrónica https://tecdigital.tec.ac.cr/revistamatematica/Libros/LibrosCDF/
EDO-NOviedo/CDFEDO1-Internet-Noviedo2018.pdf o en su defecto escaneando el código QR
adjunto que lleva a página respectiva.
Observe que:
x−y −1
M(x, y) = 2
⇒ My = 2
x x
1
N(x, y) = 1 + f (x) · ln x ⇒ Nx (x, y) = ln x f 0 (x) + f (x)
x
Luego para que la ED dada sea exacta debe cumplirse
1 −1
ln x f 0 (x) + f (x) = 2
x x
1 1
es decir, f 0 (x) + f (x) = − 2 : ED Lineal en f (x)
x ln x x ln x
1
Z
1 1
ln x · f 0 (x) + f (x) = − 2
x x
1
Z Z
⇒ (ln x · f (x))0 dx = − dx
x2
1
ln x · f (x) = +C
x
1 C 1
f (x) = + haciendo C = 0 puede tomar f (x) =
x ln x ln x x ln x
2. Sustituya la f(x) obtenida en el punto anterior y resuelva la ED exacta F
1
Cambiando f (x) = en F se tiene que resolver ahora la ED exacta:
x ln x
1 y 1
− dx + 1 + dy = 0 (A)
x x2 x
∂F ∂F
Note existe una función potencial F(x,y) tal que la ED se escribe dx + dy = 0
∂x ∂y
(B)
∂F 1 y
= − 2 (1)
∂x x x
∂F 1
= 1+ (2)
∂y x
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
y
Integrando (1) conrespecto a x a ambos lados, se tiene F(x, y) = ln x + + P(y) para alguna
x
función P, luego si la derivamos parcialmente con respecto a y se tiene
∂F 1
= + P0 (y) e igualando a (2)
∂y x
1 1
1+ = + P0 (y) ⇒ 1 = P0 (y) ⇒ P(y) = y + K
x x
y
de donde la solución generala de F viene dada por ln x + + y + K = C, es decir,
x
y
ln x + + y = B, con B = C − K
x
3.3
(yexy + 2x − 1) dx + (xexy − 2y + 1) dy = 0
⇒ exy − y2 + y + x2 − x = C
3.4
por lo que la ecuación es exacta y de esta forma existe una función f (x, y) tal que
∂f ∂f 1
= −y(y + sen x) y = + cos x − 2xy
∂x ∂y 1 + y2
∂f
considerando = −y(y + sen x), luego integrando con respecto a x, se tiene
∂x
f (x, y) = −y2 x + y cos x + h(y)
∂ f (x, y)
de donde derivando parcialmente f (x, y) con respecto a y, es decir, = −2yx + cos x + h0 (y),
∂y
∂f 1
posteriormente usando que = + cos x − 2xy, se obtiene
∂y 1 + y2
1
h0 (y) = ⇒ h(y) = arctan y +C1
1 + y2
y de esta forma la familia de soluciones implícitas de la ecuación viene dada por :
1
−y2 x + y cos x + arctan y = 1 + π
4
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
3.5
En efecto , patra para que la ecuación diferencial dada 3x2 − y3 dx + 3xy2 g(x)dy = 0 se exacta
2. Determine una función N(x, y) de tal forma que la ED (y sen x + x2 y − x)dx + N(x, y)dy = 0 sea
exacta y luego resuelvala.
1. Verifique que las siguientes ecuaciones diferenciales son exactas, luego resolverlas. En caso de no
estar la ED dada en la forma m(x, y)dx + n(x, y)dy = 0 escribala de dicha forma antes de resolverla
hasta obtener solución general. Posteriormente ingrese a página interactiva Exactas o en su defecto
escanear el código QR adjunto que lleva a página respectiva para comprobar paso a paso dichos
resultados.
Respuestas:
(a) x2 y3 − 2x3 y2 = C.
(b) e2y x − sen xy + y2 = C.
(c) xyex − 2y + x = C.
(d) x2 y − y2 − tan x = C
(e) xy2 + cos x − tan y = C
(f) xy3 − 5y + 2x√
3 = −2.
y 2y 1
d) cos x + 2 sen x + sen xy0 = 0
x x px √
2
e) 2x ln ydx + xy + y y2 + 1 dy = 0, y(0) = 2
f) (x + y2 )dx − 2xydy = 0
g) (2xy + y4 )dx + (3x2 + 6xy3 )dy = 0
h) (2xy4 ey + 2xy3 + y)dx + (x2 y4 ey − x2 y2 − 3x)dy = 0
i) (y + y cos(xy))dx + (x + x cos(xy))dy = 0
Respuestas
a) e2y x − sen xy + y2 = C
b) xyex + 2y + x = C
c) x2 y + y2 − tan x = C
d) x2 y sen x = C
1 3 √
e) x2 ln y + (y2 + 1) 2 = 3
3
−y2
f) + ln |x| = C
x
g) x y + xy6 = 0
2 3
x2 x
h) x2 ey + + 3 = C
y y
i) xy + sen(xy) = C
3. Determine el valor de k de tal modo que las siguientes ecuaciones diferenciales sean exactas.
1
M(x, y)dx + (xexy + 2xy + )dy = 0,
x
y
sea exacta. Respuesta: M(x, y) = exy y − − y2 + φ (x)
x2
5. Determine una función g(y) de tal forma que la siguiente ecuación diferencial
6. Considere la ecuación diferencial (3ya + 10xy2 )dx + (6xa−1 y − 2 + 10x2 y)dy = 0 con a ∈ R. De-
termine el valor de a de tal forma dicha ecuación sea exacta y luego resolverla.
Respuesta: A = 3 y B = 2
8. Obtener una función M(x, y) de modo tal que sea exacta la ecuación diferencial:
1 1 y3
Respuesta: f (x) = 1 con Solución General: −xy2 + 3x + ex sen(x) − ex cos(x) − = C
2 2 3
xy
11. Halle la función P(y) de tal forma que (yP(y) + ln x) dx + dy = 0 sea una ecuación diferencial
y−1
exacta.
3.7
pero si se multiplica toda la ecuación por x se obtiene una ecuación diferencial exacta.
∂ ∂
(2ey + 3x sen y) = 2ey + 3x cos y 6= xey + x2 cos y = ey + 2x cos y.
∂y ∂x
∂ ∂ 2 y
2xey + 3x2 sen y = 2xey + 3x2 cos y = x e + x3 cos y .
∂y ∂x
Al factor µ(x, y) tal que, al multiplicar una ecuación diferencial por dicho factor, ésta se convierte
en exacta, se le llama factor integrante. Puede abrir página interactiva Factor integrante o en
su defecto escaneando el código QR adjunto que lleva a página respectiva para comprobar dichos
resultados.
Observación 3.1 . Ambas ecuaciones tienen esencialmente las mismas soluciones, pero al multiplicar
por el factor integrante µ(x, y) es posible ganar o perder soluciones. En el siguiente ejemplo vemos cómo
determinar si se ha ganado o perdido alguna solución en el proceso.
3.8
Al resolver la ecuación 3.10, se tendría que su solución será de la forma F(x, y) = C, siendo F una
función tal que
∂F
= 2xy3 − 6x2 y2 , (3.11)
∂x
y
∂F
= 3x2 y2 − 4x3 y. (3.12)
∂y
Integrando 3.11 respecto de x se tiene que
Z
F(x, y) = (2xy3 − 6x2 y2 )dx = x2 y3 − 2x3 y2 + ϕ(y).
∂F
= 3x2 y2 − 4x3 y + ϕ 0 (y).
∂y
Igualando esta ecuación con la ecuación 3.12:
Por tanto, F(x, y) = x2 y3 − 2x3 y2 y la solución general de la ecuación diferencial 3.10 esta dada por
x2 y3 − 2x3 y2 = C.
Como se ha comentado, al multiplicar 3.9 por el factor integrante µ(x, y) es posible ganar o perder
soluciones. En este caso, al multiplicar la ecuación 3.9 por xy2 , se ha obtenido y = 0 como solución
de 3.10, pero no lo es de 3.9.
Observación 3.2 Puede ingresar a página interactiva Exactas o en su defecto escanear el código QR
adjunto que lleva a página respectiva para comprobar paso a paso dichos resultados.
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
Se ha visto que µ(x, y) es un factor integrante de M(x, y)dx + N(x, y)dy = 0 si se cumple que la ecuación
∂µ ∂M ∂µ ∂N
M +µ =N +µ ,
∂y ∂y ∂x ∂x
es decir
∂µ ∂µ ∂N ∂M
M −N = − µ. (3.13)
∂y ∂x ∂x ∂y
Obtener µ de esta ecuación no es sencillo puesto que a veces se llega a una ecuación en derivadas parciales
cuya resolución es más complicada que la ecuación inicial. Sin embargo, esta ecuación se simplifica si se
busca un factor de la forma µ(x), es decir µ sólo depende de la variable x, o de la forma µ(y), es decir que
sólo depende de y.
Suponga que se quiere hallar un factor integrante que sólo dependa de x. Entonces la ecuación 3.13 se
transforma en:
0 ∂N ∂M
−Nµ (x) = − µ(x),
∂x ∂y
de donde,
∂M
µ 0 (x) ∂y − ∂∂Nx
= .
µ(x) N
Si la expresión de la derecha sólo depende de x y es continua, entonces existe el factor integrante µ(x) y se
obtiene integrando esta ecuación, es decir:
R ∂∂M
y
− ∂N
∂x
dx
µ(x) = e N . (3.14)
Si ahora se tiene un factor integrante que solo dependa de y, entonces la ecuación 3.13 queda:
0 ∂N ∂M
Mµ (y) = − µ(y).
∂x ∂y
De donde ahora:
Método de resolución
∂M ∂N
Dada la ecuación diferencial en la forma M(x, y)dx + N(x, y)dy = 0, al calcular y :
∂y ∂x
∂M ∂N
Si 6= , la ecuación no es exacta.
∂y ∂x
∂M
∂y − ∂∂Nx
Si depende sólo de x, entonces un factor integrante es
N
R ∂∂M
y
− ∂N
∂x
dx
µ(x) = e N
∂N
∂x − ∂∂My
Si depende sólo de y, entonces un factor integrante es:
M
R ∂∂Nx − ∂∂M
y
dy
µ(y) = e M
Al encontrar el factor integrante, se multiplica toda la ecuación diferencial por dicho factor y, puesto
que la ecuación obtenida es ahora exacta, se resuelve hallando la solución F(x, y) = C según los
procedimientos ya anteriormente descritos para este tipo de ecuaciones diferenciales.
3.9
En efecto, considerando las ecuaciones dadas y utilizando la teoría antes descrita se verifica lo
siguiente para cada una de ellas:
∂ ∂ 2
(2x2 + y) = 1 6=
x y − x = (2xy − 1).
∂y ∂x
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
3.9
∂M
∂y − ∂∂Nx 1 − (2xy − 1) −2xy + 2 2(1 − xy) 2
En este caso se tiene que = 2
= = =− ,
N (x y − x) −x(1 − xy) −x(1 − xy) x
.
R −2
µ(x) = e x dx = x−2 .
Puede ingresar a la página interactiva Factor integrante o en su defecto escanear el código QR
adjunto que lleva a página respectiva para verificar dichos resultados obtenidos.
∂F y
= 2+ 2 (3.16)
∂x x
∂F 1
= y− (3.17)
∂y x
Si se integra a ambos lados la ecuación 3.17, con respecto a y entonces:
y2 y
Z
1
F(x, y) = y− dy = − + Ψ(x).
x 2 x
Luego derivando esta expresión respecto de x :
∂F y
= 2 + Ψ0 (x).
∂x x
Igualando esta ecuación con la ecuación 3.16 se obtiene:
y y
2
+ Ψ0 (x) = 2 + 2 ,
x x
Ψ(x) = 2x,
y de esta forma se obtiene la solución general
y2 y
− + 2x = C.
2 x
∂ ∂
(y(x + y + 1)) = x + 2y + 1 6= (x(x + 3y + 2)) = 2x + 3y + 2.
∂y ∂x
En este caso:
∂M
∂y − ∂∂Nx x + 2y + 1 − (2x + 3y + 2) −x − y − 1
= = ,
N x(x + 3y + 2) x(x + 3y + 2)
no depende sólo de x, por tanto, no existe un factor de la forma µ(x).
∂F
= y2 (x + y + 1), (3.18)
∂x
y
∂F
= xy(x + 3y + 2). (3.19)
∂y
Para hallar F, se integra la ecuación 3.18 con respecto a x:
x2 y2
Z
y2 x + y3 + y2 dx = + y3 x + ϕ(y).
F(x, y) = (3.20)
2
Para hallar ϕ(y), se deriva esta expresión respecto de y
∂F
= x2 y + 3y2 x + 2yx + ϕ 0 (y),
∂y
e igualando a la expresión dada con 3.19 se tendría:
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
Puesto que se puede elegir cualquier valor constante para k, se elige de esta forma ϕ(y) = 0 y se
sustituye en 3.20 para obtener la expresión de F, de donde finalmente la solución general es:
x 2 y2
+ y3 x + y2 x = C.
2
3.10
∂M
∂y − ∂∂Nx 2x + 4y3 − 6x − 6y3
= (2)
N 3x2 + 6xy3
−(2y3 + 4x)
=
3x2 + 6xy3
La expresión en (2) no es una función exclusivamente de x, por lo que no admite factor integrante
µ(x), sin embargo, si uno µ(y), como se muestra a continuación :
∂N
∂x − ∂∂My 6x + 6y3 − 2x − 4y3
= (3)
M 2xy + y4
2(y3 + 2x)
=
y(y3 + 2x)
2
=
y
R 2
2
e y dy = e2 ln y = eln y = y2
Así: µ(y) = y2 . Puede ingresar a la página interactiva Factor integrante o en su defecto escanear
el código QR adjunto que lleva a página respectiva para verificar dichos resultados obtenidos.
Ya que se conoce un factor integrante de (1), multiplicamos dicha ecuación por y2 , se tiene :
x2 y3 + xy6 = C
3.11
R dx
µ(x) = e x = eln x = x
luego multiplicamos (1) por µ(x) y obtenemos la ecuación:
y
+ 2x dx + (1 + ln(xy))dy = 0
x
la cual es una ecuación diferencial exacta y al resolverla se obtiene que por solución general la curva
definida implícitamente en la ecuación
y ln(xy) + x2 = C
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
3.12
Esta ecuación no es exacta pero admite factor integrante de la forma µ(x, y) = x · f (y). Halle función
f (y) y luego determine la solución general de dicha ecuación diferencial.
Para que la ecuación obtenida sea exacta, se debe cumplir que My = Nx , es decir:
f 0 (y) 2
⇒ =
f (y) y
Z 0
f (y) 2
Z
⇒ dy = dy
f (y) y
⇒ ln( f (y)) = lny2
⇒ f (y) = y2
∴ µ(x, y) = xy2
My = 8xy3 + 12x3 y2 = Nx
∂F ∂F
Por lo que existe una función F(x, y) tal que, = M(x, y) [1] y = N(x, y) [2]. De esta
∂x ∂y
forma usando [1] entonces:
Z
F(x, y) = (2xy4 + 4x3 y3 )dx + h(y) ⇒ F(x, y) = x2 y4 + x4 y3 + h(y)
F(x, y) = x2 y4 + x4 y3 + K
= x2 y4 + x4 y3 (K = 0)
Finalmente de esta forma la solución general viene dada de forma implícita por
x2 y4 + x4 y3 = C
∂M
∂y − ∂∂Nx R
f (x)dx
P1: = f (x) µ(x) = e
N
∂N
∂x − ∂∂My R
f (y)dy
P2: = f (y) µ(y) = e
M
∂M
∂y − ∂∂Nx R
f (xy)d(xy)
P3: = f (xy) µ(xy) = e tome u = xy
y·N −x·M
∂N
− ∂∂My y y y y
= e f ( x )d( x ) tome u =
∂x
R
y
P4: · x2 = f µ x
y·N +x·M x x
∂M
− ∂∂Nx
x x x x
= e f ( y )d( y ) tome u =
R
∂y 2 x
P5: ·y = f µ y
y·N +x·M y y
3.13
Considere la ecuacion (x2 y3 + 2y)dx + (2x − 2x3 y2 )dy = 0, verifique que no es exacta y que no
admite factor integrante dependiente de solo de x o solo de y, sino admite factor integrante una
función f (xy). Luego resuelvala y determine su solución general.
My − Nx 9xy2 My − Nx 9x2 y
Como = no aplica P1 y ademas como = − no aplica P2.
N 2 − 2x2 y2 −M x2 y2 + 2
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
My − Nx −3 R −3 R −3
Por otro lado como = si aplica P3. Así µ = e xy d(xy) =e u du = u−3 = x−3 y−3
yN − xM xy
es un factor integrante de la ED.
Al multiplicar (x2 y3 + 2y)dx + (2x − 2x3 y2 )dy = 0 por µ = x−3 y−3 se obtiene la ecuación diferencial
exacta
(x−1 + 2x−3 y−2 )dx + (2x−2 y−3 − 2y−1 )dy = 0 [1]
Dicha ecuación diferencial [1] al resolverla se obtiene por solución general
3.14
Encuentre los valores de n y m para que xn ym sea un factor integrante de y2 dx + (xy − x3 )dy = 0 y
verifique para los valores obtenidos, que efectivamente es exacta.
Para que la ecuación diferencial resultante [1] sea exacta debe ahora cumplirse:
∂M ∂N
De [2] M(x, y) = x−3 y−2 ⇒ = −2x−3 y−3 ; N(x, y) = x−2 y−3 − y−4 ⇒ = −2x−3 y−3
∂y ∂x
Por lo tanto, de esta forma se comprueba que la ecuación diferencial dada por [2] es exacta.
Puede visualizar explicación y solución general bajo los valores de n y m encontrados de la ecuación
diferencial asociada. Eso mediante el video a través del link: https://youtu.be/UjYmutO_KfM o
en su defecto escaneando el código QR adjunto que lleva a página respectiva para comprobar paso a
paso dichos resultados.
−ydx + xdy = 0
Es fácil comprobar que My 6= Nx y por tanto cumple no ser una ED exacta. Luego en teoría si se multi-
plica −ydx + xdy = 0 por un factor integrante de la forma η = η(x2 + y2 ) esta seria exacta, es decir,
∂ ∂
Por lo tanto, deberá cumplir (−y · η(x2 + y2 )) = (x · η(x2 + y2 )), es decir,
∂y ∂y
−η(x2 + y2 ) − yη(x2 + y2 ) · 2y = η(x2 + y2 ) + xη(x2 + y2 ) · 2x
de donde agrupando en términos de η y η 0 se tiene
η 0 (x2 + y2 ) 1
⇔ =− 2
η(x2 + y2 ) x + y2
Importante notar que todos los términos tienen argumento x2 + y2 , de lo contrario nos se puede
esperar un factor integrante de la forma supuesta.
k k
|η| = ⇔ η(v) = ±
v v
1 1
Como es suficiente un factor integrante se puede escoger η(v) = , es decir, η(x2 + y2 ) = 2 .
v x + y2
1
finalmente para resolver −ydx + xdy = 0 se multiplica por el factor integrante hallado 2 y se
x + y2
obtiene la ED exacta
y x
− 2 dx + 2 dy = 0
x + y2 x + y2
x
y cuya solución general al resolver se obtiene − arctan = C.
y
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
1. Considere la ecuación diferencial y · dx + (x2 − xy)dy = 0. Halle una función f (y) de tal forma
f (y)
que ED tenga un factor integrante de la forma µ = 2 .
x
1
a) (y ln y + yex )dx + (x + y cos y)dy = 0 con µ(y) = .
y
b) (2ey + 3x sen y)dx + (xey + x2 cos y)dy = 0 con µ(x) = x.
1
c) (x − y)dx + xdy = 0 con µ(x) = .
x2
2
d) −xdx + x2 y + y3 dy = 0 con µ(y) = e−y .
e) (2x2 + 2y2 − x)dx + (x2 + y2 − y)dy = 0 no admite factor integrante solo de x ni solo de y, sino
1
µ(x, y) = x2 +y 2.
2. Encuentre el factor integrante µ apropiado en cada una de las siguientes ecuaciones diferenciales,
posteriormente resuelvalas hasta obtener su solucin general.
y4
a) ydx + (3 + 3x − y)dy = 0. FI: µ(y) = y2 SG :y3 + xy3 − =C
4
2 x2 cos(2x) y2 1
b) (y sen(2x) + xy2 )dx + (y3 − sen2 x)dy = 0. FI: µ(y) = SG : − + + =C
y 2 2y 2 2y
c) (3x5 tan y − 2y3 )dx + (x6 sec2 y + 4x3 y3 + 3xy2 )dy = 0
1 3 y3
. FI: µ(x) = SG :x tan(y) + + y4 = C
x3 x2
2
d) (2xy + y4 )dx + (3x2 + 6xy3 )dy = 0 ,FI: µ(y) = SG: x2 y3 + xy6 = C
y
1
Respuesta: a = −1, b = −2. FI: µ(x, y) = x−1 y−2 SG: ln(x) − =C
2x2 y2
4. Halle un factor integrante de la forma µ(x, y) = xa yb para la ecuación diferencial
x4 y2 − y dx + x2 y4 − x dy = 0 y resuelvala
Respuesta: Se toman a = b = −2. El factor integrante es ρ(x, y) = x−2 y−2 . Solución General es
x3 1 y3
+ + =C
3 xy 3
1
Respuesta: m = n = −1, FI: = , SG: y + x + ln y ln x = C
xy
6. Halle el valor de α para que la ecuación diferencial (6xy3 + cos y)dx + (αx2 y2 − x sen y)dy = 0 sea
exacta y resuélvala en dicho caso.
7. Halle una función N(y) tal que h(x, y) = 3xy2 sea un factor integrante de la ecuación diferencial
x 1
− dx + N(y)dy = 0.
y2 xy
−1
Respuesta: N(y) = SG: x3 − 3xy = C.
y2
∂M ∂N ∂M ∂N 1
8. Sí = , y =− , verifique que es un factor integrante
∂x ∂y ∂y ∂x [M(x, y)] + [N(x, y)]2
2
de la ecuación diferencial M(x, y)dx + N(x, y)dy = 0.
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
9. Encuentre la solución general de la ecuación diferencial, usando un factor integrante µ(x, y), de la
forma indicada.
Respuesta
a 2x2 y2 + x3 y + xy3 = C
d) Para la ecuación (x +x4 +2x2 y2 +y4 )dx +ydy = 0, se necesita que a = b y el factor integrante
es
1
µ(x, y) = 2 2
b (x + y2 )2
y la solución general es
−1 x
+ =C
2b2 (x2 + y2 ) b2
e) Para la ecuación (3x + 2y + y2 )dx + (x + 4xy + 5y2 )dy = 0, el factor integrante es µ(x, y) =
x + y2 y la solución general es
1
µ(x, y) =
(1 − x2 + y2 )2
y la solución general es
x
=C
1 − x 2 + y2
ey
11. Encuentre todas las funciones f (x, y) que cumplan f (2018, y) = y tales que (x + y)dx +
2018
y 2x + f (x, y) dy = 0 tiene por factor integrante a µ(x, y) = x. Luego resuelva la ecuación según
f (x, y)
My − Nx
12. Demuestre que si es una función g(z) donde z = x + y entonces se cumple que µ =
Z N −M
g(z)dz
e es un factor integrante para la ecuación
Respuesta
Respuesta
a) Para la ecuación y(x2 y2 + 2) + x(2 − 2x2 y2 )y0 = 0 tenemos es factor integrante µ(x, y) =
1
siempre que xy 6= 0, y la solución general es
3(xy)3
1 1 2
ln |x| − 2 2 − ln |y| = c
3 3x y 3
b) Para la ecuación y(2xy + 1)dx + x(1 + 2xy − x3 y3 )dy = 0 tenemos el factor integrante
1
µ(x, y) =
(xy)4
Respuesta
x cos(xy) 1 x cos(2011x) 1
g(x, y) = sen(xy) + − 2; g(x, 2011) = sen(2011x) + − .
y y2 y 2011 (2011)2 (2011)2
16. Considere la siguiente ecuación diferencial f (t)dy + (t 2 + y)dt = 0. Se sabe que tal ecuación posee
un factor integrante de la forma µ(t, y) = t. Encuentre todas las posibles funciones f (t). Luego
resuelva la ecuación para las f que encontró.
Respuesta
t t4 t2
Las funciones están dadas por f (t) = + ct −1 . La solución General es + y +Cy = D.
2 4 2
17. Supongamos que µ(x, y) = ey sin x es un factor integrante de la ecuacion diferencial
yF(x)dx + x2 G(y)dy, = 0.
Respuesta
F(x) = 2x + x2 cot x
a)
G(y) = 1 + y
b) La solución General es x2 y sen(x)ey = C
dy
18. La ecuación diferencial et sec y − tan y + = 0 posee un factor integrate de la forma µ(t, y) =
dx
−at
e cos y para alguna constante a. Encuentre a y luego de acuerdo a ello resuelva dicha ED hasta
obtener la solución general.
Sección 4:
Ecuaciones Diferenciales Ordinarias de
primer orden lineales y de Bernouilli
4 — EDO Lineal y Bernoulli
4.0.2 Contenidos
Una clase de ecuaciones diferenciales que aparecen con frecuencia en las aplicaciones es la constituida por
las ecuaciones lineales. Aunque en el tema siguiente se habla de las ecuaciones diferenciales lineales en
general, se verá ahora cómo resolver las de primer orden mediante factores integrantes.
Definición 4.1
Una ecuación lineal de primer orden es aquélla que tiene la forma:
Suponga que a1 (x) , a0 (x) y b (x) son funciones continuas en un intervalo común I y con a1 (x) 6= 0.
Dividiendo por a1 (x) se puede reescribir esta ecuación en forma canónica:
donde P(x) y Q(x) son funciones continuas en dicho intervalo I. Si expresamos la ecuación 4.2 en
forma diferencial, se tiene:
Método de resolución
Como al multiplicar la ecuación 4.3 por µ(x) se obtiene una ecuación exacta, se tiene que µ(x)P(x) = µ 0 (x),
por tanto la ecuación anterior queda:
es decir,
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
d
(µ (x) y) = µ(x)Q(x),
dx
luego integrando respecto de x se tiene
Z
µ(x)y = µ(x)Q(x)dx +C.
Por tanto la solución general esta dada por:
Z
−1
y = µ(x) µ(x)Q(x)dx +C. (4.5)
Teorema 4.1
Si P(x) y Q(x) son funciones continuas en un intervalo ]a, b[ que contiene al punto x0 , para cualquier
valor inicial y0 existe una única solución y(x) en ]a, b[ del problema de valor inicial:
4.1
2
Considere la ecuación diferencial lineal de primer orden y0 + 2xy − 2xe−x = 0 y halle la solución
general.
2
En esta ecuación lineal P(x) = 2x y Q(x) = 2xe−x . Luego el factor integrante es:
2
R
2xdx
µ(x) = e = ex .
2 2
Entonces al multiplicar la ecuación y0 + 2xy = 2xe−x por dicho factor integrante µ(x) = ex , se
tiene
2 2
ex y0 + ex 2xy = 2x,
d x2
es decir, e y = 2x, de donde
dx
Z
2
ex y = 2xdx = x2 +C,
Las curvas solución se aprecian en la figura 4.1. De donde claramente se puede apreciar que para
valores muy grandes positivos o negativos, sus soluciones tienden a y = 0.
y = ã-x Ix2 + CM
Solución General
2
-4 -2 2 4
-2
-4
2
Figura 4.1: curvas solución de la ecuacion y0 + 2xy = 2xe−x
Fuente: elaboración propia con M ATHEMATICA 9
4.2
es decir,
d 2x
e y = 3e3x ,
dx
luego,
Z
e2x y = 3e3x dx = e3x +C,
y = ex +C · e−2x
Puesto que P(x) = 2 y Q(x) = 3ex son funciones continuas en toda la recta real, la solución es válida
en todo R.
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
En la figura 4.2 se muestra la familia de curvas solución de la ecuación diferencial dada. En ellas
se aprecia que para valores de x grandes positivos estas curvas solución tiende a la curva solución
y = ex , como claramente se obtuvo por medio de la solución algebraica.
4
Solución General
y = ãx + ã-2 x C
2
Solución Particular
y=ãx
-4 -2 2 4
-2
-4
4.3
1 2 π
Resolver el problema de valor inicial y0 − 2 y = x cos x con la condición inicial y = 3.
x x 2
En primer lugar, se expresa la ecuación en forma canónica:
2
y0 − y = x2 cos x (4.6)
x
−2
entonces de esta forma se tiene que P(x) = y Q(x) = x2 cos x y el factor integrante
R −2 x
µ(x) = e x dx = x−2 .
d −2
x−2 y0 − 2x−3 y =
(x y) = cos x,
dx
de donde se obtiene que al integrar a ambos lados con respecto a x :
Z
x−2 y = cos xdx = sen x +C,
y = x2 (sen x +C).
π2 2 12
Sustituyendo ahora la condición inicial y π2 = 3, se tiene 3 = 4 +C π4 ⇔ C = π2
− 1,
6
4
-4 -2 p 2 4
-4 -2 2 4 2
-2
-2
-4 -4
1 2 π
Figura 4.3: curvas solución de y0 − 2 y = x cos x y la solución particular obtenida de la condición inicial y =3
x x 2
Fuente: elaboración propia con M ATHEMATICA 9
Note que la figura 4.3 muestra las gráficas de las curvas solución general y la gráfica de la solución
particular obtenida a partir de la condición inicial dada, además es importante rescatar que en el punto (0, 0)
se visualiza claramente no existe solución única, sino infinitas y asimismo no se contradice el teorema de
existencia y unicidad en dicho punto.
4.4
x+2
Resolver la ecuación diferencial (x + 2)y0 + 2y = .
x2 + 4
Divida la ecuación dada por x + 2 6= 0, esto para escribirla en su forma
2 1
y0 + y= 2 . (4.7)
x+2 x +4
R 2
De esta forma, el factor integrante viene dado por µ = e x+2 dx = e2 ln |x+2| = (x + 2)2 (constante de
integración se escoge como 0).
(x + 2)2
(x + 2)2 y0 + 2(x + 2)y = ,
x2 + 4
es decir,
d (x + 2)2
((x + 2)2 y) = 2 ,
dx x +4
de donde al integrar con respecto a x en ambos lados y resolviendo se obtiene por solución general
de la forma implícita
(x + 2)2 y = x + 2 ln(x2 + 4) +C.
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
Puede verificar resultados obtenidos del ejemplo anterior, adémas complementar teoría mediante
página interactiva LINEAL o en su defecto escaneando el código QR adjunto que lleva a página
respectiva.
En algunos casos puede que en una ecuación diferencial dada se deba invertir las variables dependientes
e independientes, es decir, si la ED es en y(x) pueda no resolverse en variable dependiente y, variable
independiente x pero al invertirlas y tomarlas como ED en x(y) esto si sea posible, como lo mostramos a
continuación por medio del siguiente ejemplo:
4.5
1
Resuelva la siguiente ecuación diferencial y0 = [1] hasta obtener la solución general.
ey − x
Claramente la ED dada por [1] no representa una ecuación diferencial lineal en variables depen-
diente y e independiente x, es decir, en y(x), sin embargo, de acuerdo con [1] se tiene que
dy 1 dx
= y ⇒ = ey − x
dx e − x dy
es decir, al invirtir las variables en [1] se llega a una ED de la forma
Ahora resolviendo [2] se tiene en ella µ(y) = ey que al multiplicarla por dicho factor integrante se
tiene
e2y
ey · x = +C
2
De esta forma se tiene por solución general
ey
x(y) = +C · e−y
2
con C una constante arbitraria.
1. Resuelva cada una de las siguientes ecuaciones diferenciales, luego verifique sus procedimientos
y resultados obtenidos paso a paso mediante página interactiva LINEAL o de forma directa
escaneando el código QR adjunto para comprobar dichos resultados en dicha página.
a) xy0 − 4y = x6 ex .
b) cos x · y0 + y · sen x = 1
1 − e−2x
c) y0 + y = −x
e + ex
d) y0 + y · cotx = 2 cos x
e) y0 − y · tan x = e− sen x .
1
f) xy0 + 3y = ex , y (1) = 1.
x √
g) (x2 + 1)y0 + xy = (1 − 2x) x2 + 1.
h) y2 x0 + xy = 2y2 + 1 .
2. Muestre que la ecuación y0 + P(x) = Q(x)e−y puede transformarse en una ecuación lineal con la
sustitución u = ey . De acuerdo a ello resuelva xy0 + 3 = 4xe−y , con y(2) = 0.
24
Respuesta SP : ey = 4 − 3
x
3. Verifique que la ecuación lineal y0 + P(x)y = Q(x) expresada en la forma M(x, y)dx + N(x, y)dy =
0, no es exacta y admite factor integrante solo en x.
Definición 4.2
Se llama ecuación de Bernoulli a una ecuación diferencial de primer orden que se puede expresar de
la forma:
Método de resolución
transformándose de esta forma la ecuación 4.8 en una lineal para las variables (x, v). Para ello, se divide
primero la ecuación 4.8 por yn , luego,
y0 · y−n + P(x)y1−n = Q(x). (4.9)
dv dy
De este modo, tome v = y1−n y = (1 − n)y−n , luego sustituyendo en 4.9 se tiene
dx dx
1 dv
+ P(x)v = Q(x) ⇒ v0 + (1 − n)P(x)v = (1 − n)Q(x),
1 − n dx
que es una ecuación diferencial lineal de variable independiente x y variable dependiente v.
Observación 4.1 Se debe tener en cuenta que al dividir por yn se pierde la solución y = 0, luego habrá
que ver si esta solución está contenida en la solución general y si no lo está especificar que también lo es.
Una vez resuelta la ecuación para v(x) se deshace el cambio.
4.6
Al considerar las ecuaciones diferenciales dadas, se tiene en cada una de ellas lo siguiente:
1 ln x 2
(a) Expresando la ecuación en forma normal, se tiene y0 + y = y
x x
que es una ecuación de Bernoulli, ya que es de la forma 4.8, con n = 2. Al multiplicar esta
ecuación por y−2 se tiene
1 ln x
y−2 y0 + y−2 y = (4.10)
x x
dv dy
luego al hacer el cambio de variable v = y1−2 = y−1 , se tiene que = −y−2 , donde realizando
dx dx
las sustituciones correspondiente en 4.10
dv 1 ln x
− + v= ,
dx x x
es decir,
dv 1 ln x
− v=− ,
dx x x
1 ln x
que es una ecuación lineal para v con P(x) = − y Q (x) = , entonces, la solución viene dada
x x
por
Z
−1
v = µ (x) µ (x) Q (x) dx +C ,
1
y= ,
ln x + 1 +Cx
además de la solución y = 0.
Puesto que P(x) y Q(x) son continuas en el intervalo ]0, +∞[, la solución obtenida es válida en dicho
intervalo.
(b) La ecuación diferencial y0 + y = ex y−2 es una ecuación de Bernoulli con n = −2. Por tanto, se
multiplica la ecuación por y2 .
y2 y0 + y3 = ex , (4.11)
dv dy
y al hacer el cambio v = y3 ⇒ = 3y2 .
dx dx
Luego al realizar las sustituciones correspondientes en 4.11 se obtiene
1 dv
+ v = ex ,
3 dx
es decir,
dv
+ 3v = 3ex ,
dx
ecuación lineal cuya solución es
e4x
3
Z
−3x −3x
v(x) = e 4x
3e dx = 3e +C = ex +Ce−3x .
4 4
Deshaciendo el cambio inicial v = y3 , se obtiene la solución general
3
y3 = ex +Ce−3x .
4
En este caso, y = 0 no era solución de la ecuación dada. Como P(x) y Q(x) son continuas en todo el
intervalo real, la solución obtenida es válida ∀x ∈ R.
Puede verificar resultados obtenidos de los ejemplos anteriores, adémas complementar teoría me-
diante página interactiva BERNOULLI o en su defecto escaneando el código QR adjunto que lleva
a página respectiva.
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
4.7
Resolver la ecuación de Bernoulli en la variable x dada por 2y3 x0 − y2 x = −3x−1 con y > 0.
Sea z = x2 , de esta forma z0 = 2xx0 , luego sustituyendo en 4.12 se tiene y3 z0 − y2 z = −3, es decir,
z
z0 = − = −3y−3 . (4.13)
y
R dy
1
Multiplicando 4.13 por el factor integrante µ(y) = e− y = e− ln y = se obtiene
y
d −1
Z
−1 0 −2 −4
y z − y z = −3y ⇒ [zy ] = −3y−4 ⇒ zy−1 = − 3y−4 dy +C,
dy
,es decir,
zy−1 = y−3 +C ⇒ z = y−2 +C.
Como z = x2 se tiene que la solución general viene dada por la solución implícita x2 = y2 +Cy.
4.8
dy 1 dx
= ⇒ = x sen(y) + x2 sen(2y)
dx x sen(y) + x2 sen(2y) dy
de donde se tiene que
x0 (y) − x(y) sen(y) = sen2 (2y) · x2 (4.15)
esta última ecuación diferencial 4.15 una de Bernouilli en x(y) con n = 2. Ahora si en 4.15 la
dividimos entre x2 se tiene
Esta última ecuación diferencial 4.17 una lineal en v(y), la cual al resolverla se tiene por solución
Finalmente retomando las variables originales mediante el cambio de variable inicial v = x−1 se
tiene por solución general
1. Resuelva cada una de las siguientes ecuaciones diferenciales, luego verifique sus procedimientos y
resultados obtenidos paso a paso mediante página interactiva BERNOULLI o de forma directa
escaneando el código QR adjunto para comprobar dichos resultados en dicha página.
1 + 2x2 2 1 2
a) y0 + 2xy = y . SG: = C · ex + 1
x2 y
2 3 1 9
b) y0 − y = 2 y4 . SG: 3 = − +C · x−6
x x y 5x
2 1
c) x2 y0 + 2x3 y = y2 (1 + 2x2 ). SG: y−1 = C · ex +
x
C
d) (1 + x2 )y0 − xy = xy2 . SG: y−1 = √ −1
x2 + 1
2 4
e) 3xy0 − 2y = x3 y−1 . SG: y2 = x3 +C · x 3
5
1 31
f) xy0 + 2y + x10 y5 = 0 con y(−1) = 2 SP: 4 = x8 (2x2 − )
y 16
g) 2y3 x0 − y2 x = −3x−1 , (y > 0) (Ecuación Bernoulli en variable x).
Tecnlógico de Costa Rica Profesor: M.Sc. Norberto Oviedo Ugalde
y x(x + ln(x))
a) xy0 + = SG: y2 ln2 (x) = x2 + 2x ln(x) − 2x +C
ln(x) y ln2 (x)
dy y 4y
b) = SG: x =
dx x + y3 x2 C − y4
dy 1
c) 3(1 + x2 ) = 2xy(y3 − 1) SG: y = p
dx 3
1 +C · (1 + x)
x
d) xy0 = xy2 + y + xy SG: = −x + 1 +C · e−x
y
3. Encuentre la solución del siguiente problema de valores iniciales:
(y2 + 2y + x2 )y0 + 2x = 0
.
y(1) = 0
Respuesta
Se debe considerar x = x(y) para que quede una ecuación de Bernoulli cuya solución general es
x2 = −y2 +Ce−y y la solución particular que satisface y(1) = 0 es x2 = −y2 + e−y .
dy
f 0 (y) + p(x) f (y) = q(x)
dx
a) Haga un cambio de variable que convierta la ecuación anterior en una ecuación lineal.
Respuesta
Respuesta
a) Demostrativa.
Respuesta
2x 4 arctan x √
a) La ecuación después del cambio es y0 − 2
y− √ y = 0.
1+x x2 + 1
√ 2
p p
2 + c 1 + x2
b) La solución general de la ecuación anterior es y = (arctan x) 1 + x
√ p p
c) La solución de la ecuación original es xu = (arctan2 x) 1 + x2 + c 1 + x2
Respuesta
x2 ln y − x y0 = y
Respuesta
4.3 Bibliografía
Abell, Martha L. y Braselton James P. (2004), Differential Equations with Mathematica, Elsevier
Science & Technology Books.
Mora Walter (2013). Plantilla del formato y diseño de Revista Matemática ITCR.
Murray R. Spiegel (1983). Ecuaciones diferenciales aplicadas. Primera edición. México. Prentice-
Hall Hispanoamericana, S.A.
http://demonstrations.wolfram.com/