This is page i

Printer: Opaque this

Ecuación diferencial de Bernoulli

Pedro Portillo Diaz Alejandro Ortiz Rivera

David Naverrete

29 de Mayo del 2018

ii

ABSTRACT Estudio de la ecuación diferencial de Bernoulli

 This is page iii
Printer: Opaque this

Contents
0.1 CRITERIO DE IDENTIFICACION . . . . . . . . . . . . . v
0.2 DEDUCCION DE LA SOLUCIÓN . . . . . . . . . . . . . . v
0.3 EJEMPLOS . . . . . . . . . . . . . . . . . . . . . . . . . . . vii
0.4 EJEMPLO DE BERNOULLI CON CONDICIÓN INICIAL x
0.5 EJERCICIO EN CLASE POR EQUIPO . . . . . . . . . . . xii
0.6 SOLUCIÓN DEL EJERCICIO EN CLASE POR EQUIPO xii
0.7 EJERCICIOS DE BERNOULLI . . . . . . . . . . . . . . . xiii
0.8 SOLUCIONES DE LOS EJERCICIOS DE BERNOULLI . xiii
0.9 EJERCICIOS DE BERNOULLI CON CONDICIÓN INICIAL xiv
0.10 SOLUCIÓN DE LOS EJERCICIOS DE BERNOULLI CON
CONDICIÓN INICIAL . . . . . . . . . . . . . . . . . . . . xiv
0.11 PROBLEMAS DE ED DE BERNOULLI QUE VIENEN EN
EXAMENES GLOBALES Y DE RECUPERACIÓN . . . . xiv
0.12 SOLUCIÓN DE LOS EJERCICIOS USANDO SCIENTIFIC
WORKS PLACE . . . . . . . . . . . . . . . . . . . . . . . . xv
0.12.1 INTRODUCCIÓN . . . . . . . . . . . . . . . . . . . xv
0.12.2 SOLUCIÓN DE LOS EJEMPLOS USANDO SWP . xvi
0.12.3 USANDO SWP CON CONDICIONES INICIALES xviii
0.12.4 EJEMPLO DE CONDICIÓN INICIAL CON SWP . xviii
0.13 EJERCICIOS DE BERNOULLI USANDO SWP . . . . . . xx
0.14 EJERCICIOS DE BERNOULLI CON CONDICIÓN INI-
CIAL USANDO SWP . . . . . . . . . . . . . . . . . . . . . xxi
0.15 SOLUCIÓN DE LOS EJERCICIOS USANDO MATHEMAT-
ICA . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xxi
0.15.1 INTRODUCCIÓN . . . . . . . . . . . . . . . . . . . xxi
0.15.2 SOLUCIÓN DE LOS EJEMPLOS USANDO MATH-
EMATICA . . . . . . . . . . . . . . . . . . . . . . . xxii
0.15.3 USANDO MATHEMATICA CON CONDICIONES
INICIALES . . . . . . . . . . . . . . . . . . . . . . . xxv
0.15.4 EJEMPLOS DE CONDICIÓN INICIAL CON SWP xxv
0.16 EJERCICIOS DE LINEALES USANDO MATHEMATICA xxvii
0.17 EJERCIOS DE LINEALES CON CONDICIÓN INICIAL
USANDO MATHEMATICA . . . . . . . . . . . . . . . . . xxviii
iv
 This is page v
Printer: Opaque this

ECUACION DIFERENCIAL
DE BERNOULLI
0.1 CRITERIO DE IDENTIFICACION
La ecuación diferencial de Bernoulli es una ED de la siguiente forma:

dy
+ P (x)y = Q(x)y n
dx

donde n 6= 0; 1; si n = 1 la ecuación es Lineal de primer orden, si n =

0 la ecuación es lineal homogenea de primer orden. En nuestros cursos,
generalmente se toma a n como ENTERO, aunque el método es aplicable
cuando n es un REAL.

0.2 DEDUCCION DE LA SOLUCIÓN

i). En la mayoría de los textos de ecuaciones diferenciales, la ecuación
diferencial de Bernoulli se resuelve transformándola en una ED lineal.
Para hacer esto:
dy
l.- Se multiplica dx + P (x)y = Q(x)y n por y n

dy
+ P (x)y = Q(x)y n y n
dx

n dy
y + P (x)y 1 n
= Q(x)
dx
2. - El factor que multiplica a P (x) nos sugiere el siguiente cambio de
variable:
u = y1 n

que nos permite transformar la ecuación de Bernoulli en una lineal en u.

dy
3.- Para poder hacer esto, es necesario que el término y n dx se escriba
en términos de u, lo cuál se logra usando la regla de la cadena

du du dy
=
dx dy dx
vi

al usar el cambio de variable, obtenemos

du d 1 ndy
= (y )
dx dy dx
dy
= (1 n) y n
dx
n dy
de aquí,despejando y dx encontramos

n dy 1 du
y =
dx 1 n dx
y sustituyendo en

n dy
y + P (x)y 1 n
= Q(x)
dx
encontramos la siguiente ecuación diferencial en u

1 du
+ P (x)u = Q(x)
1 n dx
la cuál al normalizar nos queda

du
+ (1 n)P (x)u = (1 n)Q(x)
dx
haciendo

(1 n)P (x) = P1 (x); (1 n)Q(x) = Q1 (x)

encontramos
du
+ P1 (x)u = Q1 (x)
dx
que es lineal en u
4. Se resuelve la ED lineal en u con el procedimiento que ya estudiamos
para este tipo de ED.
5. Para encontrar la solución de la ED de Bernoulli, en la solución del
paso anterior se deberá hacer el cambio de variable u = y 1 n
RESUMEN
i) La ecuación diferencial de Bernoulli

dy
+ P (x)y = Q(x)y n
dx
se transforma a la ED lineal en u:
du
+ P1 (x)u = Q1 (x)
dx
 vii

donde
P1 (x) = (1 n)P (x); Q1 (x) = (1 n)Q(x))
al hacer el cambio de variable

u = y1 n

ii) En una forma más rápida, podemos pasar de la ED de Bernoulli a la

lineal usando la siguiente tabla

u = y1 n
se elimina y n
se cambia y por u
dy du
dx + P (x)y = Q(x)y n ! dx + (1 n)P (x)u = (1 n)Q(x)
P (x) ! (1 n)P (x)
Q(x) ! (1 n)Q(x)

es decir, tanto P (x) como Q(x) en la lineal, quedan multiplicadas por el

exponente del cambio de variable y 1 n , además de que cambiamos y por u
y se elimina y n .

0.3 EJEMPLOS
Resolver las siguientes ecuaciones diferenciales
1.
dy
+ y = xy 3
dx
SOLUCIÓN
Tenemos una ecuación diferencial de Bernoulli, la cuál se convierte en
lineal al seguir el siguiente procedimiento
i. Se multiplica la ecuación diferencial por y 3 , y nos queda

3 dy 2
y +y =x
dx
ii. De esta expresión encontramos que el cambio de variable es
2
u=y

usando la regla de la cadena

du d dy
2
= (y )
dx dy dx
dy
= 2y 3
dx
viii

3 dy
despejando y dx
1 du 3 dy
= y
dx 2 dx
y obtenemos la siguiente ecuación lineal
1 du
+u=x
2 dx
iii) normalizando, encontramos la siguiente ED lineal en u
du
2u = 2x
dx
iv) para a resolver esta ED lineal en u, procedemos a
a) encontrar su factor de integración
R
2 dx 2x
FI = e =e

b) multiplicamos la ED lineal por el FI

d 2x 2x
(e u) = 2xe
dx
c) integrando Z
2x 2x
e u= 2 xe dx

d) calculando las integrales:

calculamos la integral por partes la integral
Z
xe 2x dx
w = x; dw = dx
2x 1 2x
dv = e dx; v = e
2
encontramos que
Z Z
2x 1 2x 1 2x
xe dx = xe + e dx
2 2
integrando
Z
2x 1 1 2x
xe dx = xe 2x e
2 4
1 2x
= e (2x + 1)
4
sustituyendo en Z
2x 2x
e u= 2 xe dx
 ix

encontramos

2x 1 2x
e u = 2 e (2x + 1) + c
4
1 2x
= e (2x + 1) + c
2
e) despejando u
1 2x
2e (2x + 1) + c
u =
e 2x
1
= (2x + 1) + ce2x
2
1
= x + + ce2x
2
que es la solución de la ED lineal en u.
v) Finalmente, la solución de la ED de Bernoulli la encontramos cuando
usamos el cambio de variable

2
u=y
en la expresión anterior

2 1
y =x+ + ce2x
2
2.
dy
+ y = y 2 sin x
dx
SOLUCIÓN
La ecuación diferencial es de Bernoulli con n = 2, usando la tabla, ten-
emos que el exponente del cambio de variable se encuentra al cambiar el
signo de 2 y sumarle 1, por lo que 1 n = 1, entonces, podemos pasar de
la ED de Bernoulli a la lineal, usando la siguiente tabla:

v=y 1
elimina y 2
dy dv
dx + y = y 2 sin x ! dx v= sin x
se cambia y por v
P (x) = 1 ! ( 1)(1) = 1
Q(x) = sin x ! sin x

dv
v= sin x
dx
que es una ED lineal en v, que se resuelve de la siguiente forma
x

i) encontramos su factor de integración FI

R
dx
FI = e
x
= e
ii) multiplicamos la ED por el FI
d x
(e v) = ex sin x
dx
iii) integrando, encontramos
Z
x x
e v= e sin xdx

iv) calculemos la integral

Z
x
e sin xdx

integrando por partes

Z
1
e x sin xdx = e x
(cos x + sin x)
2
por lo que
x 1 x
e v=
e (cos x + sin x)
2
v) despejando v, encontramos la solución de la ED lineal en v
1
v= (cos x + sin x + c)
2
1
vi) usando el cambio de variable v = y encontramos la solución de la
ED de Bernoulli
1 1 1
y = = (cos x + sin x + c)
y 2
o en términos de x
2
y=
cos x + sin x + c

0.4 EJEMPLO DE BERNOULLI CON

CONDICIÓN INICIAL
Resuelva el siguiente problema de valor inicial
dy
y 1=2 + y 3=2 = 1
dx
y(0) = 4
 xi

SOLUCIÓN
I. Resolvemos la ED
dy
y 1=2 + y 3=2 = 1
dx
normalizando
dy
+ y = y 1=2
dx
la cúal es una ED Bernoulli, la cuál se puede convertir en lineal por el
siguiente procedimiento
cambio de variable v = y 3=2
se elimina y 1=2
dy dv
dx +y =y 1=2
! dx + 32 v = 3
2
se cambia y por v
P (x) = 1 ! ( 32 )(1) = 23
Q(x) = 1 ! ( 32 )(1) = 32
entonces tenemos la siguiente ED lineal en v
dv 3 3
+ v=
dx 2 2
que se resuelve por el siguiente procedimiento
a. Construimos su factor de integración FI
R 3 3
FI = e 2 xdx = e2x

b. multiplicamos la ED por el FI
d 3x 3 3
[e 2 v] = e 2 x
dx 2
c. integramos ambos lados de esta expresión
Z
3 3 3
e 2 x v = e 2 x ( dx)
2
d. calculamos la integral
3 3
e 2 xv = e 2 x + c

e. despejando v
3
v = 1 + ce 2x

f. la solución de la ED de Bernoulli se obtiene al usar el cambio de variable

v = y 3=2
3 3
y 2 = 1 + ce 2 x
II. Usamos la condición inicial

y(0) = 4
xii

para calcular la constante de integración c

43=2 = 1 + c
8 = 1+c

por lo que
c=7
III. Finalmente, sustituyendo esta constante en la solución de la ED de
Bernoulli, encontramos la solución del problema de valor inicial
3 3
y 2 = 1 + 7e 2x

0.5 EJERCICIO EN CLASE POR EQUIPO

Resolver la ED
dy 1
x +y = 2
dx y
p
3
x3 +c
SOLUCIÓN: y = x

0.6 SOLUCIÓN DEL EJERCICIO EN CLASE

POR EQUIPO
Normalizando
dy 1 1 2
+ y= y
dx x x
la cuál es una ED de Bernoulli, la cuál haciendo un cambio de variable
se convierte en una lineal, de acuerdo al siguiente esquema

cambio de variable v = y 3
se elimina y 2
dy
dx + x1 y = x1 y 2
! dv
dx + x3 v = 3
x
se cambia y por v
P (x) = x1 ! (3)( x1 ) = x3
Q(x) = x1 ! (3)( x1 ) = x3

resolvemos la ED lineal
dv 3 3
+ v=
dx x x
usando el siguiente procedimiento
 xiii

i. factor de integración
R 3
FI = e x dx = e3 ln x
3
= eln x = x3
ii) multiplicamos la ED lineal por el FI
d 3
(x v) = 3x2
dx
iii) integrando ambos lados de esta ecuación
Z
x v = 3 x2 dx + c
3

iv) calculando la integral

x3 v = x3 + c
v) despejando v
c
v = 1+
x3
3
x +c
=
x3
Finalmente, la solución de la ED de Bernoulli la encontramos al sustituir
en la expresión anterior el valor de v = y 3
x3 + c
y3 =
x3
sacando raíz cúbica p
3
x3 + c
y=
x
que es la solución de este problema

0.7 EJERCICIOS DE BERNOULLI

Resuelva las siguientes ecuaciones diferenciales de Bernoulli.
dy dy
1.- x dx + y = y12 2.- dx = y xy 3 1
dy dy
3.- x2 dx + y 2 = xy 4.- 3 1 + x2 dx = 2xy y 3 1

0.8 SOLUCIONES DE LOS EJERCICIOS DE

BERNOULLI
p
3
q q
x3 +c 3 x 1
1. y = x 2. y = 3
3x+1+ce3x 3. y = c+ln x 4. y = 3
c(1+x2 )+1
xiv

0.9 EJERCICIOS DE BERNOULLI CON

CONDICIÓN INICIAL
Resuelva la ecuación diferencial dada, sujeta a la condición inicial que se
indica

dy 1
1.- x2 2xy = 3y 4 y (1) =
dx 2
dy
2.- y 1=2 + y 3=2 = 1 y (0) = 4
dx
dy
3.- xy 1 + xy 2 = 1 y (1) = 0
dx

0.10 SOLUCIÓN DE LOS EJERCICIOS DE

BERNOULLI CON CONDICIÓN INICIAL
p q 6 3 3 e2y
1 2
1. y = 4 3 5 3 576xx5 13 2. y 2 = 1 + 7e 2x 3. x = 1 2 1 2
y2 e 2 y 2e 2 y +1

0.11 PROBLEMAS DE ED DE BERNOULLI

QUE VIENEN EN EXAMENES GLOBALES
Y DE RECUPERACIÓN
1. GLOBAL 18I VESPERTINO
3.- Resolver
dy
3 1 + x2 = 2xy y 3 1
dx
2. RECUPERACIÓN 18I VESPERTINO
3.- Resolver a ecuación diferencial ordinaria
dy
= y xy 3 1
dx
3. GLOBAL 16I VESPERTINO
3.- Resolver la ecuación diferencial

xy 0 xy 2 = (1 + x)y

4. GLOBAL 14O VESPERTINO

2.- Resolver la ecuación diferencial
p p
3
1 3 y 0 + cot(x)y = sin2 (x)y
 xv

5. GLOBAL 11O VESPERTINO

1.- Resolver la ecuación diferencial

y 2 dx x2 + xy dy = 0
6. GLOBAL 11I VESPERTINO
3.- Resolver la ecuación diferencial

1 y2
y0 y=
x x2

0.12 SOLUCIÓN DE LOS EJERCICIOS

USANDO SCIENTIFIC WORKS PLACE
0.12.1 INTRODUCCIÓN
Para resolver ED de Bernoulli usando SWP procedemos de la siguiente
forma
dy
i) Pasamos la ED dx + P (x)y = Q(x)y n a la forma lineal en u, usando
el cambio de variable u = y 1 n , del cuál obtenemos la ED lineal haciendo
lo siguiente:
a) eliminamos y n
b) se cambia y ! u
c) …nalmente cambiamos P (x) ! (1 n)P (x) y Q(x) ! (1 n)Q(x), es
decir, las funciones de la variable independiente que aparecen en la función
lineal, simplemente están multiplicadas por el exponente 1 n que aparece
en el cambio de variable.
ii) Usando la información anterior, la ED asociada lineal asociada a la ED
de Bernoulli es dudx + (1 n)P (x)u = (1 n)Q(x) o haciendo (1 n)P (x) =
PL (x) y (1 n)Q(x) = QL (x) tenemos du dx + PL (x)u = QL (x) R
iii) Encontramos el Factor de Integración usando la expresión e PL (x)dx
donde para calcular la integral hacemos
controli para escribir el símbolo de la integral, y después controle para
calcular la integral, la cuál sustituimos como exponente de la función ex-
ponencial y simpli…camos usando el comando de SWP altcy
iv) Multiplicamos la ED lineal por el FI, recordando que siempre el lado
izquierdo es igual a la derivada del producto del FI y la variable dependi-
ente, donde la integral del lado derecho la calculamos insertando el simbolo
de la integral usando el comando controli y después la calculamos con el
comando controle
v) Despejando u usando el comando altcse y poniendo en el cuadro de
dialogo y y apretando enter encontramos la solución de la ED lineal
vi) Finalmente, usando el cambio de variable u = y 1 n en la solución de
la lineal, encontramos la respectiva solución de la ED de Bernoulli
xvi

v) La forma directa de encontrar la solución de la ED lineal usando SWP

es colocar el cursor en la ED y usar el comando altcoe la cuál nos da la
solución que se encontró con el procedimiento de anterior.

0.12.2 SOLUCIÓN DE LOS EJEMPLOS USANDO SWP

1.
dy
+ y = xy 3
dx
SOLUCIÓN
Tenemos una ED de Bernoulli que convertimos en lineal por medio del
cambio de variable
u=y 2
entonces la ED lineal en términos de u es
du
2u = 2x
dx

Rusando controli y controle, calculamos la integral del FI

2dx = 2x
sustituyendo como exponente de la función exponencial y simpli…cando
usando altcy
e 2x = e 2x
multiplicamos la ED lineal por el FI, recordando que siempre el lado
izquierdo es igual a la derivada del producto del FI y la variable dependiente
d 2x 2x
(e u) = 2xe
dx
integrando, tanto el lado izquierdo como el lado derecho
Z
e 2x u = 2xe 2x dx + c

donde la integral del lado derecho la calculamos insertando el simbolo

de la integral usando el comando controli y después la calculamos con el
comando
R controle
2xe 2x dx = 21 e 2x (2x + 1)
sustituyendo, encontramos

2x 1 2x
e u= e (2x + 1) + c
2
despejando u usando el comando altcse y poniendo en el cuadro de dial-
ogo u y apretando enter encontramos
e 2x u = 21 e 2x (2x + 1) + c, Solution is: e 12x c + 12 e 2x (2x + 1)
o quitando paréntesis, usando el comando altcx, encontramos
 xvii

1
e 2x c + 12 e 2x
(2x + 1) = ce2x + x + 1
2 o sea
1
u = ce2x + x +
2
2
Finalmente, usando el cambio de variable u = y , encontramos la solu-
ción de la ED de Bernoulli

2 1
y = ce2x + x +
2
La forma directa de encontrar la solución de la ED de Bernoulli usando
SWP es colocar el cursor en la ED lineal y usar el comando altcoe la cuál
nos da la solución de la lineal en u.
du
dx 2u = 2x, Exact solution is: x + 12 ; x + C34 e2x + 12
3
…nalmente, usando el cambio de variable u = y 2 encontramos la solución
de la ED de Bernoulli
1
y 2 = ce2x + x +
2
2.
dy
+ y = y 2 sin x
dx
SOLUCIÓN
Tenemos una ED de Bernoulli que convertimos en lineal por medio del
cambio de variable
u=y 1
entonces la ED lineal en términos de u es
du
u= sin x
dx
usando
R controli y controle, calculamos la integral del FI
dx = x
sustituyendo como exponente de la función exponencial y simpli…cando
usando altcy
e x=e x
multiplicamos la ED lineal por el FI, recordando que siempre el lado
izquierdo es igual a la derivada del producto del FI y la variable dependiente
d x x
(e u) = e sin x
dx
integrando, tanto el lado izquierdo como el lado derecho
Z
e xu = e x sin xdx + c

donde la integral del lado derecho la calculamos insertando el simbolo

de la integral usando el comando controli y después la calculamos con el
comando controle
xviii
R
e x sin xdx = 12 e x (cos x + sin x)
sustituyendo, encontramos

x 1 x
e u= e (cos x + sin x) + c
2
despejando u usando el comando altcse y poniendo en el cuadro de dial-
ogo u y apretando enter encontramos
e x u = 12 e x (cos x + sin x)+c, Solution is: e 1 x c + 12 e x (cos x + sin x)
o quitando paréntesis, usando el comando altcx, encontramos
1 1
e x c + 2e
x
(cos x + sin x) = 12 cos x + 21 sin x + cex o sea

1 1
u= cos x + sin x + cex
2 2
1
Finalmente, usando el cambio de variable u = y , encontramos la solu-
ción de la ED de Bernoulli

1 1 1
y = cos x + sin x + cex
2 2
La forma directa de encontrar la solución de la ED de Bernoulli usando
SWP es colocar el cursor en la ED lineal y usar el comando altcoe la cuál
nos da la solución de la lineal en u.
du
dx u = sin x, Exact solution is: 12 cos x + 21 sin x + C38 ex
…nalmente, usando el cambio de variable u = y 1 encontramos la solución
de la ED de Bernoulli

1 1 1
y = cos x + sin x + C38 ex
2 2

0.12.3 USANDO SWP CON CONDICIONES INICIALES

0.12.4 EJEMPLO DE CONDICIÓN INICIAL CON SWP
1. Resuelva el siguiente problema de valor inicial

dy
y 1=2 + y 3=2 = 1
dx
y(0) = 4

SOLUCIÓN
Normalizando la ED
dy 1=2
+y =y
dx
Es una ED de Bernoulli que convertimos en lineal por medio del cambio
de variable
u = y 3=2
 xix

entonces la ED lineal en términos de u es

du 3 3
+ u=
dx 2 2
usando
R 3 controli y controle, calculamos la integral del FI
3
2 dx = 2x
sustituyendo como exponente de la función exponencial y simpli…cando
usando altcy
3 3
e2x = e2x
multiplicamos la ED lineal por el FI, recordando que siempre el lado
izquierdo es igual a la derivada del producto del FI y la variable dependiente
d 3x 3 3
(e 2 u) = e 2 x
dx 2
integrando, tanto el lado izquierdo como el lado derecho
Z
3
x 3 3x
e u=
2 e 2 dx + c
2
donde la integral del lado derecho la calculamos insertando el simbolo
de la integral usando el comando controli y después la calculamos con el
comando
R 3 3 x controle
3
2 dx = e 2 x
2e
sustituyendo, encontramos
3 3
e 2 xu = e 2 x + c

despejando u usando el comando altcse y poniendo en el cuadro de dial-

ogo y y apretando enter encontramos
1 3
u= 3 c + e2x
e 2x

o quitando paréntesis, usando el comando altcx, encontramos

c
u = 3 +1
e2x
3
= 1 + ce 2 x
3
Finalmente, usando el cambio de variable u = y 2 , encontramos la solu-
ción de la ED de Bernoulli
3 3
y2 = 1 + c 2x

La forma directa de encontrar la solución de la ED de Bernoulli usando

SWP es colocar el cursor en la ED lineal y usar el comando altcoe la cuál
nos da la solución de la lineal en u.
xx
n 3
o
du
dx + 32 u = 23 , Exact solution is: 1
3 C13 + e 2 x
e2x
3
…nalmente, usando el cambio de variable u = y 2 encontramos la solución
de la ED de Bernoulli
3 1 3
y2 = 3 C13 + e 2 x
e 2x

o al usar el comando altcx encontramos la solución en la forma

3 3
y 2 = 1 + ce 2x

Para encontrar la constante c, usamos la condición inicial

y(0) = 4
x = 0; y = 4

sustituyendo en la solución
3 3
4 2 = 1 + ce 2 (0)

que se resuelve al usar el comando altcse en esta expresión

3 3
4 2 = 1 + ce 2 (0) , Solution is: 7
Finalmente, la solución del problema de valor inicial es
3 3
y 2 = 1 + 7e 2x

Para encontrar directamente la solución con SWP al problema de valor

inicial, insertamos una matriz de tamaño 2 1 usando el comando altix y
después eligiendo el tamaño rows 2 columns 1 y haciendo enter

después colocamos la ED y la condición inicial en esta matriz

dy 1
dx + y = y
2

y(0) = 4
usando el comando altcoe, encontramos la solución del problema de valor
inicial 1
dy
dx + y = y
2

y(0) = 4
y que en este caso el comando no funciona.

0.13 EJERCICIOS DE BERNOULLI USANDO

SWP
Resuelva las siguientes ecuaciones diferenciales de Bernoulli usando SWP
dy dy
1.- x dx + y = y12 2.- dx = y xy 3 1
dy dy
3.- x2 dx + y 2 = xy 4.- 3 1 + x2 dx = 2xy y 3 1
 xxi

0.14 EJERCICIOS DE BERNOULLI CON

CONDICIÓN INICIAL USANDO SWP
Resuelva la ecuación diferencial dada, sujeta a la condición inicial que se
indica usando SWP

dy 1
1.- x2 2xy = 3y 4 y (1) =
dx 2
dy
2.- y 1=2 + y 3=2 = 1 y (0) = 4
dx
dy
3.- xy 1 + xy 2 = 1 y (1) = 0
dx

0.15 SOLUCIÓN DE LOS EJERCICIOS

USANDO MATHEMATICA
0.15.1 INTRODUCCIÓN
Para resolver ED de Bernoulli usando MATHEMATICA procedemos de la
siguiente forma
dy
i) Pasamos la ED dx + P (x)y = Q(x)y n a la forma lineal en u, usando
el cambio de variable u = y 1 n , del cuál obtenemos la ED lineal haciendo
lo siguiente:
a) eliminamos y n
b) se cambia y ! u
c) …nalmente cambiamos P (x) ! (1 n)P (x) y Q(x) ! (1 n)Q(x), es
decir, las funciones de la variable independiente que aparecen en la función
lineal, simplemente están multiplicadas por el exponente 1 n que aparece
en el cambio de variable.
ii) Usando la información anterior, la ED asociada lineal asociada a la ED
de Bernoulli es dudx + (1 n)P (x)u = (1 n)Q(x) o haciendo (1 n)P (x) =
PL (x) y (1 n)Q(x) = QL (x) tenemos du dx + PL (x)u = QL (x) R
iii) Encontramos el Factor de Integración usando la expresión e PL (x)dx
iv) Calculamos la integral usando MATHEMATICA de la siguiente forma.
escinttesc para escribir el símbolo de la integral, y después mayusenter
para calcular la integral la cuál sustituimos como exponente de la función
exponencial y simpli…camos usando el comando de SWP Simplify[expr]
iv) Multiplicamos la ED lineal por el FI, recordando que siempre el lado
izquierdo es igual a la derivada del producto del FI y la variable dependi-
ente, donde la integral del lado derecho la calculamos insertando el simbolo
de la integral usando el comando escinttesc para escribir el símbolo de la
integral, y después mayusenter para calcular la integral
xxii

v) Despejando u usando el comando Solve[expr,vars] encontramos la

solución de la ED lineal
vi) Finalmente, usando el cambio de variable u = y 1 n en la solución de
la lineal, encontramos la respectiva solución de la ED de Bernoulli
vii) La forma directa de encontrar la solución de la ED lineal usando
MATHEMATICA es usar el comando DSolve[edo==0,y,x] que nos da la
solución general de la EDO en términos de la función y[x] y usamos DSolve[{edo==0,condinic},y,x]
para resolver la EDO en términos de la función y[x] con condiciones ini-
ciales.

0.15.2 SOLUCIÓN DE LOS EJEMPLOS USANDO

MATHEMATICA
1.
dy
+ y = xy 3
dx
SOLUCIÓN
Tenemos una ED de Bernoulli que convertimos en lineal por medio del
cambio de variable
u=y 2
entonces la ED lineal en términos de u es
du
2u = 2x
dx
usando escinttesc y mayusenter, calculamos la integral del factor de in-
tegración
In[4]:= n[Integral]-2 n[Di¤erentialD]x
Out[4]= -2 x
calculando la exponencial de este resultado usando la función Exp[expr]
In[5]:= Exp[-2 x]
Out[5]= E^(-2 x)
multiplicamos la ED lineal por el FI, recordando que siempre el lado
izquierdo es igual a la derivada del producto del FI y la variable dependiente

d 2x 2x
(e u) = 2xe
dx
integrando, tanto el lado izquierdo como el lado derecho
Z
2x
e u= 2xe 2x dx + c

usando escinttesc y mayusenter, calculamos la integral del lado derecho

In[7]:= n[Integral]-2 x E^(-2 x) n[Di¤erentialD]x
Out[7]= -2 E^(-2 x) (-(1/4) - x/2)
 xxiii

simpli…cando la expresión usando Simplify[expr]

In[8]:= Simplify[-2 E^(-2 x) (-(1/4) - x/2)]
Out[8]= 1/2 E^(-2 x) (1 + 2 x)
Por lo tanto

2x 1 2x
e u= e (2x + 1) + c
2
y ahora despejamos u como una función de x usando el comando Solve[expr,vars]
In[10]:= Solve[E^(-2 x) u == 1/2 E^(-2 x) (1 + 2 x) + c, u]
Out[10]= {{u -> 1/2 (1 + 2 c E^(2 x) + 2 x)}}
simpli…camos la expresión usando Simplify[expr]
In[11]:= Simplify[1/2 (1 + 2 c E^(2 x) + 2 x)]
Out[11]= 1/2 + c E^(2 x) + x
o sea

1
u=x+ + ce2x
2
2
Finalmente, usando el cambio de variable u = y , encontramos la solu-
ción de la ED de Bernoulli

2 1
y = ce2x + x +
2
Ahora, chequemos esta solución usando el comando DSolve[edo==0,y,x]
aplicado directamente a la ED lineal, donde tenemos que especi…car que u
es una función de x y la derivada se maneja usando primas..
In[18]:= DSolve[Derivative[1][u][x] - 2 u[x] == -2 x, u, x]
Out[18]= {{u -> Function[{x}, -2 (-(1/4) - x/2) + E^(2 x) C[1]]}}
…nalmente, usando el cambio de variable u = y 2 encontramos la solución
de la ED de Bernoulli
1
y 2 = ce2x + x + :
2
Despejando y en la expresión anterior usando Solve[expr]
In[30]:= Solve[y^-2 == -2 (-(1/4) - x/2) + E^(2 x) C[1], y]
Out[30]= {{y -> -(Sqrt[2]/Sqrt[1 + 2 x + 2 E^(2 x) C[1]])}, {y ->
Sqrt[2]/Sqrt[1 + 2 x + 2 E^(2 x) C[1]]}}
es decir r
2
y= 2
2ce2x + 2x + 1
Aplicando directamente el comando DSolve[edo==0,y,x] a la ED de Bernoulli
dy
dx + y = xy 3
In[31]:= DSolve[Derivative[1][y][x] + y[x] == x y[x]^3 , y, x]
Out[31]= {{y ->
Function[{x}, -(Sqrt[2]/Sqrt[1 + 2 x + 2 E^(2 x) C[1]])]}, {y ->
Function[{x}, Sqrt[2]/Sqrt[1 + 2 x + 2 E^(2 x) C[1]]]}}
xxiv

es decir r
2 2
y=
2ce2x + 2x + 1
que coincide con la expresión anterior.
2.
dy
+ y = y 2 sin x
dx
SOLUCIÓN
Tenemos una ED de Bernoulli que convertimos en lineal por medio del
cambio de variable
u=y 1
entonces la ED lineal en términos de u es
du
u= sin x
dx
usando escinttesc y mayusenter, calculamos la integral del factor de in-
tegración
In[19]:= n[Integral]-1 n[Di¤erentialD]x
Out[19]= -x
calculando la exponencial de este resultado usando la función Exp[expr]
In[20]:= Exp[-x]
Out[20]= E^-x
multiplicamos la ED lineal por el FI, recordando que siempre el lado
izquierdo es igual a la derivada del producto del FI y la variable dependiente
d x x
(e u) = e sin x
dx
integrando, tanto el lado izquierdo como el lado derecho
Z
x
e u= e x sin xdx + c

usando escinttesc y mayusenter, calculamos la integral del lado derecho

In[21]:= n[Integral]- E^- x Sin[x] n[Di¤erentialD]x
Out[21]= 1/2 E^-x (Cos[x] + Sin[x])
Por lo tanto

x 1 x
e u=
e (cos x + sin x) + c
2
y ahora despejamos u como una función de x usando el comando Solve[expr,vars]
In[22]:= Solve[E^- x u == 1/2 E^-x (Cos[x] + Sin[x]) + c, u]
Out[22]= {{u -> 1/2 (2 c E^x + Cos[x] + Sin[x])}}
o sea
1
u= (cos x + sin x) + cex
2
 xxv

1
Finalmente, usando el cambio de variable u = y , encontramos la solu-
ción de la ED de Bernoulli
1
1
y (cos x + sin x) + cex
=
2
Ahora, chequemos esta solución usando el comando DSolve[edo==0,y,x]
aplicado directamente a la ED lineal, donde tenemos que especi…car que u
es una función de x y la derivada se maneja usando primas..
IIn[25]:= DSolve[Derivative[1][u][x] - u[x] == -Sin[x], u, x]
Out[25]= {{u -> Function[{x}, E^x C[1] + 1/2 (Cos[x] + Sin[x])]}}
…nalmente, usando el cambio de variable u = y 1 encontramos la solución
de la ED de Bernoulli y 1 = cex + 21 (cos x + sin x).
Despejando y en la expresión anterior usando Solve[expr]
In[28]:= Solve[y^-1 == E^x C[1] + 1/2 (Cos[x] + Sin[x]), y]
Out[28]= {{y -> 2/(2 E^x C[1] + Cos[x] + Sin[x])}}
es decir
2
y=
2cex cos x + sin x
Aplicando directamente el comando DSolve[edo==0,y,x] a la ED de Bernoulli
dy
dx + y = y 2 sin x
In[27]:= DSolve[Derivative[1][y][x] + y[x] == y[x]^2 Sin[x], y, x]
Out[27]= {{y -> Function[{x}, 2/(2 E^x C[1] + Cos[x] + Sin[x])]}}
es decir
2
y= x
2ce cos x + sin x
que coincide con la expresión anterior.

0.15.3 USANDO MATHEMATICA CON CONDICIONES

INICIALES
0.15.4 EJEMPLOS DE CONDICIÓN INICIAL CON SWP
1. Resuelva el siguiente problema de valor inicial

dy
y 1=2 + y 3=2 = 1
dx
y(0) = 4

SOLUCIÓN
Normalizando
dy
+ y = y 1=2
dx
Tenemos una ED de Bernoulli que convertimos en lineal por medio del
cambio de variable
u = y 3=2
xxvi

entonces la ED lineal en términos de u es

du 3 3
+ u=
dx 2 2
usando escinttesc y mayusenter, calculamos la integral del factor de in-
tegración
In[32]:= n[Integral]3/2 n[Di¤erentialD]x
Out[32]= (3 x)/2
calculando la exponencial de este resultado usando la función Exp[expr]
In[33]:= Exp[(3 x)/2]
Out[33]= E^(3 x/2)
encontramos el FI
3
FI = e2x
multiplicando la ED por el FI, recordando que el lado izquierdo es el
producto del FI con la variable dependiente, encontramos
d 3x 3 3
(e 2 u) = e 2 x
dx 2
integrando el lado izquierdo y el lado derecho, tenemos
Z
3
x 3 3x
e u=
2 e 2 dx
2
usando escinttesc y mayusenter, calculamos la integral del lado derecho
In[34]:= n[Integral]3/2 E^(3 x/2) n[Di¤erentialD]x
Out[34]= E^(3 x/2)
entonces, sustituyendo el valor de la integral tenemos
3 3
e 2 xu = e 2 x + c

y ahora despejamos u como una función de x usando el comando Solve[expr,vars]

In[35]:= Solve[E^(3 x/2) u == E^(3 x/2) + c, u]
Out[35]= {{u -> E^(-3 x/2) (c + E^(3 x/2))}}
simpli…cando usando el comando Simplify[expr]
In[36]:= Simplify[E^(-3 x/2) (c + E^(3 x/2))]
Out[36]= 1 + c E^(-3 x/2)
3
es decir, la solución es u = 1 + ce 2 x
Ahora, chequemos esta solución usando el comando DSolve[edo==0,y,x]
aplicado directamente a la ED lineal, donde tenemos que especi…car que u
es una función de x y la derivada se maneja usando primas.
In[38]:= DSolve[Derivative[1][u][x] + 3/2 u[x] == 3/2, u, x]
Out[38]= {{u -> Function[{x}, 1 + E^(-3 x/2) C[1]]}}
3
que nos da la solución que encontramos anteriormente, u = 1 + ce 2 x .
…nalmente, usando el cambio de variable u = y 3=2 encontramos la solu-
3 3
ción de la ED de Bernoulli y 2 = 1 + ce 2 x .
 xxvii

Despejando y en la expresión anterior usando Solve[expr]

In[39]:= Solve[y^(3/2) == 1 + E^(-3 x/2) C[1], y]
Out[39]= {{y -> (E^(-3 x/2) (E^(3 x/2) + C[1]))^(2/3)}}
simpli…cando usando el comando Simplify[expr]
In[40]:= Simplify[(E^(-3 x/2) (E^(3 x/2) + C[1]))^(2/3)]
Out[40]= (1 + E^(-3 x/2) C[1])^(2/3)
es decir
3 2
y = (1 + ce 2x )3

Para encontrar c sustituimos x = 0 y y = 4 en la solución

3 2
4 = (1 + ce 2 (0) ) 3
Ahora despejamos c usando el comando Solve[expr,vars]
In[43]:= Solve[4 == (1 + c E^(-3 (0)/2) )^(2/3), c]
Out[43]= {{c -> 7}}
por lo que c = 7
Entonces la solución de este problema es
3 2
y = (1 + 7e 2x )3

Para encontrar directamente la solución con Mathematica al problema

de valor inicial, usamos el comando DSolve[{edo==0,condini},y,x] a la ED
de Bernoulli con condición inicial

dy
y 1=2 + y 3=2 = 1
dx
y(0) = 4

IIn[46]:= DSolve[{Derivative[1][y][x] + y[x] == y[x]^(-1/2),

y[0] == 4}, y, x]
Out[46]= {{y -> Function[{x}, (E^(-3 x/2) (7 + E^(3 x/2)))^(2/3)]}}
es decir
3 2
y = (1 + 7e 2x )3

que coincide con la expresión anterior.

0.16 EJERCICIOS DE LINEALES USANDO

MATHEMATICA
Resuelva las siguientes ecuaciones diferenciales de Bernoulli usando Math-
ematica
dy dy
1.- x dx + y = y12 2.- dx = y xy 3 1
dy dy
3.- x2 dx + y 2 = xy 4.- 3 1 + x2 dx = 2xy y 3 1
xxviii

0.17 EJERCIOS DE LINEALES CON

CONDICIÓN INICIAL USANDO
MATHEMATICA
Resuelva la ecuación diferencial dada, sujeta a la condición inicial que se
indica usando Mathematica

dy 1
1.- x2 2xy = 3y 4 y (1) =
dx 2
dy
2.- y 1=2 + y 3=2 = 1 y (0) = 4
dx
dy
3.- xy 1 + xy 2 = 1 y (1) = 0
dx