Integrador de Ecuaciones Carlos
Integrador de Ecuaciones Carlos
Integrador de Ecuaciones Carlos
David Sarumeño
30 De Noviembre Del 2016
Ecuaciones Diferenciales
Literal a
C1 = 1 0 = −C1eπ − C2eπ
C2 = −1 C1 + C2 = 0
Literal b
y(0) = 1 y(π) = −1
1
C1 = 1 Reemplazamos el C1 = −1 en −1 = C1eπ (−1) + C2eπ (0)
−1 = −C1eπ
Literal c
y(0) = 1
C1 = 1
• y( π2 ) = 1
π π
• 1 = C1e 2 (0) + C2e 2 (1)
π
1 = C2 e 2
π
C2 = e− 2
π
y = ex cos x + e− 2 ex sin x
Literal c
y(0) = 0 y(π) = 0
C1 = 0 0 = C1 eπ
La respuesta del C1 es = 0
y = C2ex sin x
2
Verifique si las funciones dadas determinan un conjunto funda-
mental de soluciones de la Ecuación Diferencial en el intervalo
indicado.Forme una solución general.
00 0
y − 2y + 5y = 0 ex cos 2x , ex sin 2x (−∞, ∞)
w = ex cos (2x)(ex sin (2x) +2ex cos (2x)) −ex sin (2x)(ex cos (2x)
−2ex sin (2x))
w = e2x(cos (2x) sin (2x) + 2(cos (2x))2) −e2x(cos (2x) sin (2x) −
2(sin (2x))2)
w=0
Ejericio.33
00 0
y − 4y + 4y = 2e2x + 4x − 12
Yc = C1e2x + C2xe2x
m2 − 4m + 4 = 0
m=2
Yc = C1e2x + C2xe2x
y = Ce2x + Dx + E
0
y = 2Ce2x + D
00
y = 4Ce2x
3
4Ce2x − 4(2Ce2x + D) + 4Ce2x + Dx + E = 2e2x + 4x − 12
4Ce2x − 8Ce2x − 4D + 4Ce2x + 4Dx + 4E = 2e2x + 4x − 12
4Dx − 4D + 4E == 2e2x + 4x − 12
C = 2x
4D = 4 D = 44
D = 1 −4D + 4E = −12
− 4 + 4E = −12
4E = 8
E=2
Yp = x − 2
Y = Yc + Yp = C1e2x + C2xe2x + x − 2
00 0 2x 2x
9y −12y +4y = 0 Y1 = e 3 y = u(x)e 3 Sacos la primera
2x
derivada de Y1 = e 3
0 2x 2x
y = 23 e 3 u + e 3 u
00 2x 00 2x 0 2x
Sacamos la segunda derivada y = e 3 u + 43 e 3 u + 49 e 3 u
2x 00 2x 0 2x 2x 2x 2x 2x 00
9e 3 u + 43 e 3 u + 49 e 3 u − 12 23 e 3 u + e 3 u + 4u(x)e 3 = 9e 3 u = 0
2x 00 2x 0 2x 2x 2x 0 2x 2x 00
9e 3 u +12e 3 u +4e 3 u−8e 3 u+12e 3 u +4u(x)e 3 = 9e 3 u = 0
2x 00 2x 0 2x 2x 2x 00
9e 3 u + 24e 3 u − 4ue 3 + 4u(x)e 3 = 9e 3 u = 0
2x 00
9e 3 u = 0
00 2x 2x
u =0 −4ue 3 + 4u(x)e 3
2x
ue 3 (4 + 4x)
u(C1x − C2)
u = C1 x + C2
2x
Y2 = xe 3
Ejericio12
00
4x2y + y = 0 y = xm
4
0
y = m(xm−1)
00
y = m(m − 1)xm−2
4x2m(m − 1)xm−2 + xm = 0
xm(4m(m − 1) + 1) = 0
4m2 − 4m + 1 = 0
b = −4, c = 1
a = 4, √
4+ (−4)2 −4(4)(1)
m1 = √ 82
4− (−4) −4(4)(1)
m2 = 8
1
m1 = 2
1
m2 = 2
1 1
Y2 = C1x 2 + C2x 2 ln (x)
00 0
•2y + √ 2y + y = 0
−2 4−8
m= 4
+
m1 = (− 12 α 21 β)
−
x
y = e− 2 (C1 cos( x2 ) + C2 sin( x2 ))
Ejercicio30
d2 y
dθ2
+y
=0
0
y( π3 ) =
0, y ( πe ) = 2
00
y +y =0
m2 + 1 = 0
m2 = −1
+
m= i
−
α=0
+
β= =i
−
5
y = e0(C1 cos(θ) + C2 sin(θ))$ y = C1 cos(θ) + C2 sin(θ)
0 = C1 cos( π3 ) + C2 sin( π3 )
0.5C1√+ 0.86C2 = 0
C1 3
2 + 2 C2
0
y = −C1 sin(θ) + C2 cos(θ)
2 = −C √1
sin( π3 ) + C2 cos( π3 )
2 = − 23 C1 + 12 C2
En la grafica representa la solucion de la Ecuacion Diferencial Or-
dinaria de Segundo orden. Explique porque coinciden.
Ejercicio 45
00 0
y + 2y + 2y = 0
m2 + 2m +
√2 = 0
−2+ 4−4(2)(1)
m1 = √2
−2− 4−4(2)(1)
m2 = 2
−2
m1 = 2 + 2i2 = −1 + 1i
m2 = −2 −2i
2 − 2 = −1 − 1i
y = ex(C1 cos(x) + C2 sin(x))
•Existe Oscilacion entre sin y cos
•Es una funcion decreciente
•α < 0
6
•Existen Raices imaginarias con α, β
Ejercicio 54
Considere el problema de los valores de frontera.
00
• y + xy = 0
m1 = xi m2 = xi
m2 + λ = 0
y = C1 cos (λ x) +C2 sin (λ x)
y (0) = 0
0 = C1 cos (λ 0) +C2 sin (λ 0)
0 = C1 cos (0) +C2 sin (0)
C1 = 0
y( π2 ) = 0
0 = C1 cos (λ π2 ) +C2 sin (λ π2 )
0 = 0 + C2 sin (λ π2 )
C2 = 0
La Solucion de Y es y = 0
C1 + C2 = 0
Condiciones :
0 0
t = 0 x(0) = 1 x = 0 x (0) = 0
7
x =1 t = 1
1 1
C1 = 2 C2 = 2 t = 2.07
Ejercicio 30
00 0
y + 4y + 4 y = (3 + x)e−2x
0
y(0) = 2 y (0) = 5
m2 + 4m + 4 = 0
(m + 2) = 0
(m + 2) = 0
m1 = −2
m2 = −2
2B = 3
B = 32
A = 1
yn = C1 e−2x + C2 x e−2x
yp = (Ax + B)
0
yp = (2(B + Ax)x − 2(B + Ax)x2)e−2x
2(B + Ax) e−2x = (3 + x)e−2x
Reemplazamos en esta f ormula
yp = (Ax + B)
8
yp = (x + 32 )x2 e−2x
y = C1 e−2x + C2 e−2x + (x + 3 2
2 )x e−2x
(−A cos (x) − B sin (x)) + 4(A cos (x) + B sin (x)) = sin (x)
− A cos (x) − B sin (x) + 4 A cos (x) + 4B sin (x) = sin (x)
3 A cos (x) + 3 B sin (x) = sin (x)
3 (A cos (x) + B sin (x)) = sin (x)
(A cos (x) + B sin (x)) = sin3(x)
sin (x)
yp = 3
+√
0 0 − 4(1)(4)
−
2
00
y + 4 y = 0
m2 + 4 = 0
9
1 = C1 cos (2(0)) + C2 sin (2(0))
C1 = 1
0
y = C1 sin (x) + C 2 cos (2x)
0
y (0) = 2
2 = 2 C1sin (0) − 2 C2 cos (0)
2 C2 = 2
C2 = 1
yn = cos (2x) + sin (2x)
yn = cos (2x) + sin (2x) + sin3(x)
3.5
Ejercicio 12
00
Encuentre la solucion particular de la ecuacion diferencial. y −
0 x
2y + y = 1 +e x2
m2 − 2m + 1 = 0
(m − 1) = 0
(m − 1) = 0
m1 = 1
m2 = 1
0
y2 = x ex + ex
0
y2 = ex (x + 1)
x
x ex
e
w = x x = e2x (x + 1) −
e e (x + 1)
x e2x = e2x (x + 1 + x) = e2x
x
0 x e −x e2x
w1 = e
x x
=
1 + x2
2 e (x + 1)
1 x+ x
e 0 e2x
w2 = x
e x = 1 + x2
e 1 + x2
10
−x e2x −x
U1 = 1 + x2
= 1 + x2
U = x2 du =
e2x
2x.dx
du
dx = 2x
−x −x du
R R
U1 = 1 + x2
dx = 1 + 4 . 2x
− ln (1 + x2 )
U1 = 2
e2x 1
U2 = 1 + x2
= 1 + x2
e 2x
1
R
U2 = 1 + x2
dx = arctan (x)
x x x ln (1 + x2 )
y = C1e + C2xe + e 2 + xex arctan (x)
x
e2 0
w2 = 1 x2 x = x ex
2e x e
2
1 x.dx = x2
R
U1 =R
x ex
U2 = −1 = − ex(x − 1)
11
x
x2 e 2 x
yp = 2 − ex(x − 1)e− 2
3.6
Ejercicio 23
Resuelva la ecuacion diferencial por variacion de parametros.
00 0
x2y − xy − y = 0 m(m − 1) + (m − 1) = 0
m2 − m + m − 1 = 0
m2 −1 = 0
yc = P
√
m 0 − 4(1)(−1)
m = 2 = 1
y1 = x y2 = x ln (x)
yc = C1 x + C2 x ln (x)
−x x ln (x)
w = = x[ln (x) + 1] −
1 ln (x) + 1
x ln (x) = x
0 x ln (x)
w1 = = −[x ln (x). ln (x)] =
ln (x) ln (x) + 1
−x [ln (x)]2
12
x 0
w2 = = x ln (x)
1 ln (x)
−x [ln (x)]2
w1 =
Rx
[ln (x)]2.dx
R
du w1 =
x ln (x)
w2 = x
R R
du w2 = ln (x) ; U2 = x ln (x) − x
y (1) = 5
0
y (1) = 3
00 0
x2 y − 3xy + 4y = 0
2
m + (−3 − 1)m + 4 = 0
m2 − 4m + 4 = 0
(m − 2) (m − 2) = 0
13
y2 = 2 C1 x + C2 (2x ln (x) + x)
3 = 2 C1 + C2
C2 = −7
2
yp = 5x2 − 72 x2 ln (x)
3.8
Ejercicio 9
Una masa que pesa 8lb se sujeta un resorte cuando se pone en
movimiento el sistema de resorte masa exhibe un movimiento ar-
monico simple. Determine la ecuacion de movimiento si la constante
de resorte es de libra sobre pie y la masa se libera inicialmente desde
un punto que esta 6in por debajo de la posicion de equilibrio con
velocidad descendente de 23 pies
in
d2 x
dt2
+ w2x = 0
w 2
= mk
8
m = 32 = 14
d2 x
dt2
+ 4x = 0
x(t) = C1 cos (2t) + C2 sin (2t)
1
2 = C1 cos (0) + C2 sin (0)
x(t) = C1 cos (2t) + C2 sin (2t)
3
2 = −C1 sin (0) 2 + C2 cos (0)2
2 C2 = 23 C2 = 34
x(t) = 12 cos (2t) + 34 sin (2t)
x(t) = qA sin (wt + 0)
2 2
A = ( 12 ) + ( 43 )
A = 0.90
tan = θ = CC12
14
1
2
θ = arctan ( )3
4
θ = 33, 69
md2 x
dt + B dx
dt + kx = 0
d2 x B dx K
dt +m dt + m = 0
2
√
m + 4 2m + 8 = 0
m1=m=−2√2 √
√
−2 2 −2 kt
x = C1 e + C2 e t
x (0) = 0
C1=0
ecuacion de la velocidad
√ −2√2t
x = 5 − 10 2t e
5dt
x= s
15
t = 0.35
Una masa de un slug se sujeta a un resorte cuya constante es de
5lb/pie. En inicio la masa libera desde un pie por debajo de
la posicion de equilibrio a velocidad descendente de 5pies/s y el
movimiento subsiguiente se desarrolla en un medio que ofrece una
fuerza amortiguarte numericamente igual al doble de la velocidad
instantanea
a) Encuentre la ecuación del movimiento y la masa es conducido
por una fuerza externa igual a f (t) = 12 cos (2t) + 3sen (2t)
b) grafique las soluciones transitorias y permanente en los mismo
ejes de coordenadas
c) grafique la solución de movimiento x00 + 2x0 + 5x = 12cos (2t) +
3sen (2t) x(0) = 1 x2(0) = 5
x = 3sen(2t)
A = 15013
100
B = 13
qp (t) = 150 100
13 cos (t) + 13 sen (t)
ip (t) = − 150 100
13 sen (t) + 13 cos (t)
16
¿Cual es la carga en la que capacitor despues un largo tiempo ?
1 00
2q + 10q + 150
−10t 3
q (t) = e C1cos (10t) + 2sen (10t) + 2
q (0) = 1 y q 0 (0) = 0 C1 y C2 = − 21
3
q(t) = 2
17