0% encontró este documento útil (0 votos)
63 vistas22 páginas

Taller EDO 2

Este documento presenta el resumen de un trabajo evaluativo sobre estadística inferencial realizado por tres autores. En menos de 3 oraciones, resume los principales puntos tratados en el documento, como la solución de una ecuación diferencial exacta mediante el método de factor integrante y el cambio de una ecuación diferencial a la forma estándar para su resolución.

Cargado por

no se
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
63 vistas22 páginas

Taller EDO 2

Este documento presenta el resumen de un trabajo evaluativo sobre estadística inferencial realizado por tres autores. En menos de 3 oraciones, resume los principales puntos tratados en el documento, como la solución de una ecuación diferencial exacta mediante el método de factor integrante y el cambio de una ecuación diferencial a la forma estándar para su resolución.

Cargado por

no se
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

Universidad Coopetiva de Colombia

Facultad de Ingeniería Sistemas

Trabajo Evaluativo

Estadistica Inferencial

Autores:
Erick Andres Báez Cruz
Juan Sebastian López Urbano
Jhonnatan David Mendez Celeita

9 Septiembre del 2022


2.4
25

(y 2 cos x − 3x2 y − 2x)dx + (2y sin x − x3 + ln y)dy = 0, y(0) = e

M (x, y) = y 2 cos x − 3x2 y − 2x ⇒ My (x, y) = 2y cos x − 3x2


3
N (x, y) = 2y sin x − x + ln y ⇒ Nx (x, y) = 2y cos x − 3x2

Como vemos si es exacta, comenzamos el procedimiento para una EDO exacta.

fx (x, y) = M (x, y)
fx (x, y) = y 2 cos x − 3x2 y − 2x
Z Z
fx (x, y)dx = (y 2 cos x − 3x2 y − 2x)dx

f (x, y) = y 2 sin x − x3 y − x2 + ϕ(y)


fy (x, y) = 2y sin x − x3 + ϕ′ (y)
fy (x, y) = N (x, y)
3
+ ϕ′ (y) = ( − x3 + ln y
(( ((
2y
( (sin
( x(− x 2y(sin
(x(
ϕ′ (y) = ln y
Z Z
ϕ′ (y)dy = ln ydy

Sustitucion por partes.

u = ln y, dv = 1
1
du = , v = ydy
x
Z
1
y ln y − y dy
y
y ln y − y

ϕ(y) = y ln y − y
f (x, y) = y 2 sin x − x3 y − x2 + y ln y − y

y 2 sin x − x3 y − x2 + y ln y − y = C

1
31

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

M (x, y) = 2y 2 + 3x ⇒ My (x, y) = 4y
N (x, y) = 2xy ⇒ Nx (x, y) = 2y

Se observa que My (x, y) ̸= Nx (x, y), se procede convertir la ecuación a exacta.

My − Nx
N

4y − 2y
2xy

yy(4−2)
2xy

2
2x

1
x

1
R
e x dx = eln x = x

 
x (2y 2 + 3x)dx + (2xy)dy
(2xy 2 + 3x2 )dx + (2x2 y)dy

Comprobamos si es exacta

M (x, y) = 2xy 2 + 3x2 ⇒ My = 4xy


N (x, y) = 2xy ⇒ Nx = 4xy

Como vemos si es exacta, comenzamos el procedimiento para una EDO exacta.

fx (x, y) = M (x, y)
fx (x, y) = 2xy 2 + 3x2
Z Z
fx (x, y)dx = (2xy 2 + 3x2 )dx

f (x, y) = x2 y 2 + x3 + ϕ(y)
fy (x, y) = 2x2 y + ϕ′ (y)
fy (x, y) = N (x, y)
2x y + ϕ′ (y)

2
 =  2
2x y

ϕ (y) = 0
Z Z
ϕ′ (y)dy = 0dy

ϕ(y) = 0
f (x, y) = x2 y 2 + x3 + 0

x2 y 2 + x3 = C

2
36
(y 2 + xy 3 )dx + (5y 2 − xy + y 3 sin y)dy = 0

M (x, y) = y 2 + xy 3 ⇒ My (x, y) = 2y + 3xy 2


2 3
N (x, y) = 5y − xy + y sin y ⇒ Nx (x, y) = −y

Se observa que My (x, y) ̸= Nx (x, y), se procede convertir la ecuación a exacta.

Nx − My
M

−y − 2y − 3xy 2
y 2 + xy 3

−3y−3xy 2
y 2 +xy 3
−3y (1+xy)

y 2
(1+xy)

−3 yy2
− y3

−3
− y3 1
R R
e = e−3 y = e−3 ln y = eln y = y− 3

y −3 (y 2 + xy 3 )dx + (5y 2 − xy + y 3 sin y)dy


 

(y −1 + x)dx + (5y −1 − xy −2 + sin y)dy


Comprobamos si es exacta

M (x, y) = y −1 + x ⇒ My = −y −2
N (x, y) = 5y −1 − xy −2 + sin y ⇒ Nx = −y −2
Como vemos si es exacta, comenzamos el procedimiento para una EDO exacta.

fx (x, y) = M (x, y)
fx (x, y) = y −1 + x
Z Z
fx (x, y)dx = (y −1 + x)dx

x2
f (x, y) = y −1 x + + ϕ(y)
2
fy (x, y) = −y −2 x + ϕ′ (y)
fy (x, y) = Nx (x, y)
−2 ′
−xy + ϕ (y) = 5y −1 − xy −2 + sin y
ϕ′ (y) = 5y −1 + sin y
Z Z
ϕ′ (y)dy = (5y −1 + sin y)dy

ϕ(y) = 5 ln y − cos y
x2
f (x, y) = y −1 x + + 5 ln y − cos y
2

1
xy −1 + x2 + 5 ln y − cos y = C
2

3
2.5
18

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

Para cambiarla a la forma estandar:

(1 + x)
y′ − y = y2
x
Para la solucion observe que n=2

u = y −1
dy dy du
=
dx du dx
dy 1 du
=− 2
dx u dx
Efectuando la sustitucion:

1 du (1 + x) 1 1
− 2
− = 2
u dx x u u
multiplicamos a ambos lados de la igualdad por −u2

du (1 + x)
+ u = −1
dx x
Resolucion mediante factor integrante:

1+x
R R
P (x)dx x dx
e =e = eln x+x = eln x ex = xex

du (1 + x)
xex [ + u] = −(xex )
dx x
du
xex + (1 + x)ex u = −(xex )
dx
d
[xex u] = −xex
dx Z
xex u = − xex dx

Metodo de integracion por partes:

u = x
du = dx
dv = ex
v = ex

Z
x x
xe u = −(xe − ex dx)
xex u = −xex + ex + c
x x
1 xe e c
=−  x − 
 
+ x
y xe
 xex xe

4
1
y=− 1 c
1+ x + xex

20

dy
3(1 + t2 ) = 2ty(y 3 − 1)
dt
Cambiarla a la forma estandar

3(1 + t2 )y ′ + 2ty = 2ty 4


2t 2t
y′ + y= y4
3(1 + t2 ) 3(1 + t2 )

1
u=
y3
1
y3 =
u
1
y= √ 3
u
1
y = u− 3

Multiplicar ambos lados de la igualdad por 1


3(1+t2

1 4 du 2t −1 2t 14
− u− 3 + u 3 = u− 4
3 dx 3(1 + t2 ) 3(1 + t2 )
1 4 du 2t −1 2t 4
− u− 3 + u 3 = u− 3
3 dx 3(1 + t2 ) 3(1 + t2 )
du 6t 6
+ u =
dx 3(1 + t2 ) 3(1 + t2 )
2t 2t
u′ − u = −
1 + t2 1 + t2
Resolucion mediante factor integrante


R 2t
dt
R dw 1 1 1
e 1+t2 = e− w = eln w = =
w 1 + t2
Reemplazando

1 2t 2t
2
(u′ − 2
u=− )
1+t 1+t 1 + t2
1 1 1 2t
2
u′ − 2
( 2
)u = −
1+t 1 + t 2t + t (1 + t2 )2
d 1 2t
[ 2
u] = −
dx 1 + t (1 + t2 )2
Z
1 1
u=− dt
1 + t2 (1 + t2 )2

v = 1 + t2 dv = 2t

5
Z
1 du
u = −
1 + t2 u2
1
u = v −1 + c
1 + t2
1 1 1
2 3
= +c
1+t y 1 + t2
1
= 1 + (1 + t2 )c
y3

s
1
y= c
1 + (1 + t2 )

38

Crecimiento de la población En el estudio de la población dinámica uno de los más famosos modelos
para un crecimiento poblacional limitado es la ecuación logística

dP
= P (a − bP ),
dt
Crecimiento de la población En el estudio de la población dinámica uno de los más famosos modelos
para un crecimiento poblacional limitado es la ecuación logística.

3.1 y 3.2 Resumen

1. Crecimiento de bacterias
dP
− kP = 0
dt
2. Vida media
dA
= kA, A(0) = A0
dt
3.Ley de enfriamiento de newton
dT
= K (T − T m)
dt
4. Mezclas
dA
= (Razondeentradadesal) − (Razondesalidadesal) = Rentra − Rsale
dt
5. Circuito en serie LR
di
L + Ri = E(t)
dt
6. Circuito en serie RC
1
Ri + q = E(t)
C
7. Dinamica Poblacional
dP/dt
P
8. Hipotesis de dependencia de densidad
dP/dt
= f (P )
P

6
o
dP
= P f (P )
dt
Ecuacion logistica
9
dP  r 
=P r− P
dt K

dP
= P (a − bP )
dt
11.Solucion de ecuacion logistica
aP
Pt =
bP0 + (a − bPo ) ∗ e−at

aP
Pt =
bP0 + (a − bPo ) ∗ e−at
12. Modificacion de la ecuacion logistica
dP
= P (a − bP ) − h
dt
y

dP
= P (a − bP ) + h
dt

dP
= P (a − blnP )
dt
15. Reaccion quimicas
M
a− X
M +N
y
N
− X
M +N

  
dX M N
α a− X b− X
dt M +N M +N

dX
= K (α − X) (β − X)
dt
3.1
2

Suponga que se sabe que la población de la comunidad del problema 1 es de 10 000 después de tres
años. ¿Cuál era la población inicial P0 ? ¿Cuál será la población en 10 años? ¿Qué tan rápido está
creciendo la población en t = 10?
Solucion Punto 2
Para la solucion del ejercicio nos vamos al ejemplo 1 del libro que se encuentra como ’Ejemplo 1
crecimiento de bacterias’ donde nos dan las siguientes ecuaciones:
dP
i) = kP
dt

7
ii) P (t) = P0 · ekt
donde se nos presenta que k es una constante de proporcionalidad, t es el tiempo, en este caso, en años,
P es la poblacion y P0 es la poblacion inicial. En el ejercicio nos dicen que este esta relacionado con el
problema anterior el cual dice que ’la Poblacion inicial P0 se duplico en 5 años’, entonces tenemos que

P (5) = 2 · P0

Al tener esto en cuenta miramos la segunda ecuacion dada por el libro y notamos que las poblaciones
iniciales se cancelan por lo que la unica variable desconocida seria la constante de proporcionalidad k
por lo que la despejariamos de la siguiente manera:

P (5) = P0 · e5k

2 · P0 = P0 · e5k
P0
2 = e5k
P0
2 · 1 = e5k
cancelamos euler con su funcion inversa

ln 2 = ln e5k

ln 2 = 5k
ln 2
=k
5
Al conocer el valor de la constante de proporcionalidad k, ya podriamos hallar la Poblacion inicial P0 ,
teniendo en cuenta que,

P (3) = 10000

Hallamos la poblacion inicial P0 de la siguiente manera:


ln 2
P (3) = P0 · e 5 ·(3)

3
10000 = P0 · e 5 ln 2
3
10000 = P0 · eln 2 5
3
10000 = P0 · 2 5
√5
10000 = P0 · 23
10000

5
= P0
23
racionalizamos 10000

5 3
2
√5
10000 23

5
· √5
23 23
√5
10000 · 23

5
26
√5
√5
10000 · 22 · 21

5
√5
26 · 21
√5
√5
10000 · 22 · 2 1

5
√5 
6
2 · 2 1

10000 √ 5
· 22
2
√5
5000 · 4
Una vez racionalizado tenemos que la poblacion inicial es:

5
P0 = 5000 · 4 = 6597,55

Una vez hallada la poblacion inicial podemos hallar la poblacion en 10 años de la 2 maneras:

8
1. Primeramente podriamos multiplicar por 4 la poblacion inical P0 por que si en si años se duplico,
en 10 habra quadruplicado su poblacion inicial P0

4 · P0

4 · 6597,55 = 26390,2

2. Segundamente con la ecuacion P (t) = P0 · ekt y quedaria:


ln 2
P (10) = 6597,55 · e 5 (10)

10
P (10) = 6597,55 · eln 2 5
P (10) = 6597,55 · 4 = 26390,2

Para saber que tan rapido crece la poblacion en t = 10 usamos la primera ecuacion dada por el libro
dP
dt = kP y despajamos la constante de proporcionalidad k y la poblacion por la que habra en 10 años
es decir 26390,2 por lo que quedaria:
dP
= kP
dt
dP ln 2
= · 26390,2
dt 5
dP 26390,2
= · ln 2
dt 5
dP 26390,2
= ln 2 5
dt
dP
= 3658,45
dt
R: La Poblacion inicial de la comunidad era de 6597,55 habitantes, en 10 años fueron
26390,2 habitantes y durante t = 10 crecieron 3658,45 habitantes por año

14

Un termómetro se lleva de una habitación hasta el ambiente exterior, donde la temperatura del aire es
5° F. Después de 1 minuto, el termómetro indica 55° F y después de 5 minutos indica 30° F. ¿Cuál era
la temperatura inicial de la habitación?
Solucion punto 14
La ley de enfriamiento de Newton establece que ’la tasa de perdida de calor de un cuerpo es
proporcional a la diferencia de temperatura entre el cuerpo y sus alrededores’ y tiene como ecuacion:
dT
= k(T − Tm )
dt
la cual tiene como solucion:

T (t) = Tm + C · ekt

y tenemos que k es la constante de proporcionalidad, C es una constante, Tm es la temperatura de


ambiente, T es la temperatura y t es el tiempo, en este caso, en minutos. en el ejercicio nos dan las
siguientes temperaturas :

Tm = 5◦ F

T (1) = 55◦ F
T (5) = 30◦ F
se define a Tm = 5F debido a que es la unica temperatura presentada sin tiempo por lo que inferimos
que es la temperatura del ambiente, Para hallar las constantes debemos remplazar los valores y hacer
un sistema de ecuaciones 2x2.

T (1) = 5 + Cek·1

55 = 5 + Cek

9
55 − 5 = Cek
50 = Cek
T (5) = 5 + Cek·5
30 = 5 + Ce5k
30 − 5 = Ce5k
25 = Ce5k
Despejamos la constante C en cualquiera de las 2 y la remplamzamos en la ecuacion opuesta
25 = Ce5k
25
C = 5k
e
remplazamos en
50 = Cek
 
25
50 = ek
e5k
 
25
50 = ek
5k
e

50 = 25 · e−4k
50
= e−4k
25
ln 2 = −4k
ln 2
− =k
4
una vez hallada la constante de proporcionalidad k la remplazamos y despejamos C en cualquiera de la
2 ecuaciones
25
C = 5k
e
25
C=
5(− ln42 )
e

25
C= 5
e− 4 ln 2
25
C= −5
eln 2 4
25
C = −5
2 4
5
C = 25 · 2 4
√4
C = 25 · 25

4

4
C = 25 · 24 · 2

4
C = 25 · 2 · 2

4
C = 50 2
C = 59,46
Hallada la constante C ahora en la funcion de Temperatura colocamos t = 0 para hallar la temperatura
inicial y tendriamos:
T (t) = Tm + Cekt

T (0) = 5 + 59,46 · e(− )·(0)


ln 2
4

T (0) = 5 + 59,46 · e0
T (0) = 5 + 59,46
T (0) = 64,46◦ F
R: La temperatura incial de la habitacion es de 64,46 ◦ F

10
31

Se aplica una fuerza electromotriz de 100 volts a un circuito en serie RC, en el que la resistencia es de
200 ohms y la capacitancia es de 10−4 farads. Determine la carga q(t) del capacitor, si q(0) = 0.
Encuentre la corriente i(t).
Solucion punto 31
la ecuacion diferencial para los circuitos RC se define como:
dq 1
R + q = E(t)
dt Ca
teniendo en cuenta que, R es la resistencia, Ca es la capacitancia y E(t) es el voltaje, en el ejercicio
estas variables estan definidas de la siguiente manera:

R = 200

Ca = 10−4
E(t) = 100
remplazamos en la ecuacion diferencial
dq 1
200 + −4 = 100
dt 10
dq
200+ 104 q = 100
dt
Pasamos la ecuacion a forma estandar
200 dq 10000 100
+ q = 
200 dt 200 200
dq 1
+ 50q =
dt 2
hallamos el factor integrante
R R
P (x)dx P (t)dt
e =e
R R
50dt
e = e50 dt
= e50t
La forma de la solucion de la ecuacion es:

q(t) = qtr(t) + qs(t)


R
qtr(t) = Ce P (t)dt
Z R
1
qs(t) = R · e P (t)dt · f (t)dt
e P (t)dt
Donde qtr es la carga transitoria del capacitor en el circuito en serie RC y qs es la carga estacionaria del
capacitor
remplazamos qtr, qs y resolvemos
R
P (t)dt
qtr(t) = Ce = Ce50t
Z
1 R
P (t)dt
qs(t) = R
P (t)dt
· e · f (t)dt
e
Z
1
qs(t) = · e50t · E(t)dt
e50t
Z
1
1
qs(t) = e50t · dt ·
e50t2
Z
1 1
qs(t) = 50t · e50t dt
e 2
1 1 1
qs(t) = · · e50t
e50t 2 50

11
1 1 50t
qs(t) = e·
e50t
100
1 1 50t
qs(t) = 50t · e

e
  100
1
qs(t) =
100
Vemos como va quedando la solucion de q(t)

q(t) = qtr(t) + qs(t)


1
q(t) = Ce50t +
100
Vemos que tenemos la constante C por lo que para hallar su valor tenemos que usar el Valor inicial del
ejercicio que es q(0) = 0
1
q(0) = Ce50(0) +
100
1
0 = Ce0 +
100
1
0=C+
100
1
− =C
100
Y conociendo la constante C ya tenemos completa la funcion q(t)

q(t) = qtr(t) + qs(t)


1
q(t) = Ce50t +
100
1 50t 1
q(t) = − e +
100 100
Como en el ejercicio nos piden hallar la corriente i(t) y esta es igual a i(t) = dq
dt asi que resolvemos la
siguiente ecuacion y despues despejamos dq dt

dq 1
+ 50q =
dt 2
dq 1
+ 50q =
dt 2
 
dq 1 50t 1 1
+ 50 − e + =
dt 100 100 2
dq 1 50t 1 1
− e + =
dt 2 2 2
dq 1 50t 1 1
− e = −
dt 2 2 2
dq 1 50t
− e =0
dt 2
dq 1
= e50t
dt 2
dq 1
i(t) = = e50t
dt 2
1 50t 1
R: se determina que la carga q(t) cuando q(0) = 0 es q(t) = − 100 e + 100 y tambien se
dq 1 50t
determina que i(t) = dt = 2 e

3.2
5

12
dp
= P (a − bP )
dt
dp
= aP − bP 2
dt
dp
− aP = −bP 2
dt

1 dP dP du dP du
u = p−1 P = = = −u−1
u dt dt dt dt dt

du
+ (−au) = b
dt R
eP (x) = ea dt = eat
du
eat + aeat = beat
dt Z
d  at 
e u = b eat
dt
beat
eat u = +c
a
at
be + ac
u =
aeat
1 b c
= + at
p a e
1
p = b c
a + eat
eat a
p =
beat + ac
4.3
7

12y ′′ − 5y ′ − 2y = 0
La convertimos a le ecuación característica

12m2 − 5m − 2 = 0
Hallar m1 y m2


−(−5) ± 25 − 96
m =
2(12)

−(−5) + 25 − 96
m1 =
2(12)

−(−5) − 25 − 96
m2 =
2(12)

−(−5) + −71
m1 =
2(12)

−(−5) − −71
m2 =
2(12)

5 71
m1 = + i
24 √24
5 71
m1 = − i
24 24

13
como vemos m1 y m2 son complejos conjugados así que aplicamos la siguiente formula

√ √
5 71 71
y = e [C1 cos(
24 x) + C2 sin( x]
24 24

14

2y ′′ − 3y ′ + 4y = 0

La convertimos a le ecuación característica

2m2 + 3m + 4 = 0

Hallar m1 y m2

p
−(−3) ± 32 − 4(2)(4)
m =
p2a
−(−3) ± 32 − 4(2)(4)
m =
2(2)

3 ± 9 − 32
m =
√4
3 ± 23i
m =
4√
3 23i
m1 = +
4 √4
3 23i
m2 = −
4 4
como vemos m1 y m2 son complejos conjugados así que aplicamos la siguiente formula

√ √ !
3 23 23
y=e 4x C1 cos x + C2 sin x
4 4

28

d5 x d4 x d3 x d2 x
2 5
− 7 4 + 12 3 + 8 2 = 0
ds ds ds ds
Pasandolo a ecuacion caracteristica

2m5 − 7m4 + 12m3 + 8m2 = 0

Factorizando

m2 (2m3 − 7m2 + 12m + 8) = 0

Por analisis:

m1 = 0

y mediante ruffini:

14
2 −7 12 8
1
− 2 −1 4 −8
2 −8 16 0

se obtiene la ecuacion:

1
(m + )(2m2 − 8m + 16) = 0
2
De la ecuacion cuadratica

p
−(−8) ± (−8)2 − 4(2)(−8)
m =
2(2)

64 − 128

m =
√4
8 ± −64
m =
4
8 8i
m = ±
4
 4
m = 2 ± 2i

Dejando como soluciones:

m1 = 0
m2 = 0
m3 = 2 + 2i
m4 = 2 − 2i
1
m5 = −
2
Solucion implicita m1 = m2

C1 + C2 x

m3
1
c3 e − 2 x

m4 y m5 complejos

e2x [C4 cos 2x + C5 sin 2x]

Total

1
C1 + C2 x + c3 e− 2 x + e2x [C4 cos 2x + C5 sin 2x]

34

y ′′ − 2y ′ + y = 0, y(0) = 5, y ′ (0) = 10

m2 − 2m + 1 = 0

15

−b ± b2 − 4ac
m =
2a p
−(−2) ± (−2)2 − 4(1)(1)
m =
2(1)

2± 4−4
m =
2
2
m =
2

m1 = m2 = 1

y = C1 ex + C2 xex

Problema de valores inicales

C1 e0 + C2 0e0 = 5
C1 = 5

y′ = C1 ex + C2 ex + C2 xex
y′ = C1 ex + C2 (ex + xex )

C1 e0 + C2 (e0 + 0e0 ) = 10
C1 + C2 = 10
5 + C2 = 10
C2 = 10 − 5
C2 = 5

C1 = 5, C2 = 5

Resultado final

y = 5ex + 5xex

36

y ′′′ + 2y ′′ − 5y ′ − 6y = 0

La convertimos a le ecuación característica

2m3 + 2m2 − 5m − 6 = 0

Se realiza Ruffini para factorizar el polinomio


1 2 −5 −6
−1 −1 −1 6
1 1 −6 0

16
(m + 1)(m2 + m − 6)

Por analisis:

m3 = −1

p
−(1) ± 12 − 4(1)(−6)
m =
2(1)
p
−(1) ± 1 − (−24)
m =
2(1)

−(−1) + 1 + 24
m1 =
2(1)

−(1) − 1 + 24
m2 =
2(1)

−1 + 25
m1 =
2√
−1 − 25
m2 =
2
−1 + 5
m1 =
2
−1 − 5
m2 =
2
4
m1 =
2
−6
m2 =
2
m1 = 2
m2 = −3

y = C1 e−x + C2 e2x + C3 e−3x

4.4
6

y ′′ − 8y ′ + 20y = 100x2 − 26xex

Pasamos la ecuacion a su alternativa

m2 − 8m + 20 = 0

Usamos la formula cuadratica y resolvemos


p
8 ± (−8)2 − 4(1)(20)
m=
2(1)

√ √
16i 16i
m1 = 4 + , m2 = 4 −
2 2
m1 = 4 + 2i, m2 = 4 − 2i
yc = e4x (C1 cos(2x) + C2 sin(2x))
yp = Ax2 + Bx + C + (Dx + F )ex => Ax2 + Bx + C + Dxex + F ex

17
Yp′ = 2Ax + B + Dex + Dxex + F ex
yp′′ = 2A + Bex + Dex + Dxex + F ex
Se considero una solucion para la derecha del igual

20Ax2 + (−16A + 20B)x + (2D + F − 8D − 8F + 20F )ex + (D − 8D + 20D)xex ...


+(2A − 8B + 20C) = 100x2 − 26xex

∗20Ax2 = 100x2
100
A=
20
A=5
∗(−16A + 20B)x = 0
−16(5) + 20B = 0
20B = 80
80
B=
20
B=4

x x
∗(D − 8D + 20D)
xe = −26
xe
13D = −26
26
D=−
13
D = −2
∗(2D + F − 8D − 8F + 20F )ex = 0
−6D + 13F = 0
12
13F = −12F = −
13
∗2A − 8B + 20C = 0
2(5) − 8(4) + 20C = 0
20C = 22
22 11
C= =
20 10
Ya con A,B,C,D,F hallados remplazamos en yp

11 12
yp = 5x2 + 4x + − 2xe − ex
10 13
y = yc + yp
11 12
R: y = e4x (C1 cos(2x) + C2 sin(2x)) + 5x2 + 4x + − 2xex − ex
10 13
22

y ′′′ − 2y ′′ − 4y ′ + 8y = 6xe2x

18
1 −2 −4 8
2 2 0 −8
1 0 −4 0

m3 − 2m2 − 4m + 8 = 0
(m − 2)(m2 − 4) = (m + 2)
m1 = 2
m2 = −2
m3 = 2
Y c = C1 e−2x C2e2x − C3xe2x
Y p = (Ax + B)e2x = Axe2x + Be2x
Y p′ = Ae2X + ZAxe2x + 2Be2x
Y p′′ = 2Ae2X + ZAxe2x + 4Axe2x + 4Be2x
Y p′′′ = 4Ae2X + 4Ae2X + 4Ae2x + 8Axe2x + 8Be2x
3(4Ae2x ) + 8Axe2x + 8Be2x − 2(2(2Ae2x ) + 4Axe2x + 4Be2x ) + 8(Axe2x + Be2x ) = 6xe2x
12Ae2x ) + 8Axe2x + 8Be2x − 8Ae2x − 8Axe2x − 8Be2x ) + 8Axe2x + 8Be2x = 6xe2x
(12A + 8B − 8A − 8B + 8B)e2X + (8A + 8A + 8A)xe2x = 6xe2x
4A + 8B = 0
 
3
4 + 8B = 0
4
8B = −3
 
3
B= −
8
8A = 6
 
6
A=
8
 
3
A=
4
 
3 2x 3 2x
Yp= xe − e
4 8
Y =Yc+Yp
3 3
Y = C1e−2x ∗ C2e−2x C3 e−2x + xe2x − e2x
4 8

30

19
y ′′ + 4y ′ + 4y = (3 + x)e−2x
m2 + 4m + 4 = (3 + x)e−2x

(m + 2)(m + 2) = 0
m1 = −2
m2 = −2

yc = C1 e−2x + C2 xe−2x

yp = (Ax3 + Bx2 )e−2x


yp′ = (3Ax2 + 2Bx)e−2x − 2(Ax3 + Bx2 )e−2x
= (3Ax2 + 2Bx − 2Ax3 − 2Bx2 )e−2x

yp′′ = (6Ax + 2B − 6Ax2 − 4Bx)e−2x − 2(3Ax2 + 2Bx − 2Ax3 − 2Bx2 )e−2x


= (6Ax + 2B − 6Ax2 − 4Bx − 6Ax2 − 4Bx + 4Ax3 + 4Bx2 )e−2x
= (6Ax + 2B − 12Ax2 − 8Bx + 4Ax3 + 4Bx2 )e−2x

y ′′ + 4y ′ + 4y = (3 + x)e−2x
(6Ax + 2B − 12Ax2 − 8Bx + 4Ax3 + 4Bx2 )e−2x ...
... + (12Ax2 + 8Bx − 8Ax3 − 8Bx2 )e−2x ...
... + (4Ax3 + 4Bx2 )e−2x = (3 + x)e−2x

(6A − 8B + 8B)x + (−12A + 4B + 12A − 8B + 4B)x2 ...


... + (4A − 8A + 4A)x3 + 2B = 3 + x

2B = 3
3
B=
2

(6A − 8B + 8B)x = x
6A = 1
1
A=
6

1 3
yp = ( x3 + x2 )e−2x
6 2

y = yc + yp
1 3
y = C1 e−2x + C2 xe−2x + ( x3 + x2 )e−2x
6 2

20
y(0) = 2
1 3
y = C1 e−2x + C2 xe−2x + ( x3 + x2 )e−2x
6 2
1 3
2 = C1 e−2(0) + C2 (0)e−2(0) + ( (0)3 + (0)2 )e−2(0)
6 2
C1 = 2

y ′ (0) = 5
1 3
y = C1 e−2x + C2 xe−2x + ( x3 + x2 )e−2x
6 2
′ −2x −2x −2x (3x2 e−2x − 2x3 e−2x )(6) (6xe−2x − 6x2 e−2x )(2)
y = −2C1 e + C2 e − 2C2 xe + +
62 22
((0) e − 2(0) e )(6) (6(0)e − 6(0)2 e0 )(2)
2 0 3 0 0
y ′ = −2(−2)e0 + C2 e0 − 2C2 (0)e0 + +
62 22
5 = −2(2) + C2
C2 = 9

21

También podría gustarte