Taller EDO 2
Taller EDO 2
Trabajo Evaluativo
Estadistica Inferencial
Autores:
Erick Andres Báez Cruz
Juan Sebastian López Urbano
Jhonnatan David Mendez Celeita
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
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
M (x, y) = 2y 2 + 3x ⇒ My (x, y) = 4y
N (x, y) = 2xy ⇒ Nx (x, y) = 2y
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
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
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
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
(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
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
1
u=
y3
1
y3 =
u
1
y= √ 3
u
1
y = u− 3
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.
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
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
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
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 · 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:
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)
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
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
Factorizando
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
8±
m =
√4
8 ± −64
m =
4
8 8i
m = ±
4
4
m = 2 ± 2i
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
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
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
2m3 + 2m2 − 5m − 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
4.4
6
m2 − 8m + 20 = 0
√ √
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 = 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
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
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