Tarea de Valores de Frontera: Métodos
númericos
Barrales Ordoñez Denise
Universidad Nacional Autónoma de México
ENES Unidad Juriquilla
Licenciatura de Ingeniería en Energías Renovables
26 de marzo de 2022
Ejercicio 1
En un laboratorio experimentan como disminuir una nueva sepa viral con una vacuna la
disminución del virus se representa por la siguiente ecuación:
x2 y 00 + xy 00 + y = x2
Donde y representa el tiempo después de aplicar la vacuna en meses y x la cantidad de
antiviral en mililitros aplicada se sabe que a partir del primer mes de la cantidad de virus
es de 0 y al tercer sigue siendo 0. Aplicando el método de disparo con método de Runge
Kutta-Euler (con paso de 0.5) y cambiando la variable encuentra el valor para z donde se
cumplen las condiciones de tiempo anterior.
Solución:
y(1)=0 y(3)=0
−z y
zx = − 2 +1
x x
1
Primera suposición:
z(1) = 1
k x k yk f y
1 1 0 1 0 1
2 1,5 0,5 1 0,111 1
≈ y(3) = 0
3 2 1 1,056 0,22 1,056
4 2,5 1,52 1,666 0,28 1,16
5 3 2,11 1,311
Segunda suposición:
z(1) = -2
k xk yk f y
1 1 0 −2 3 −2
2 1,5 −1 −0,5 1,77 −0,5
3 2 −1,25 0,33 1,11 0,38
4 2,5 1,055 0,94 0,78 0,94
5 3 −0,58 1,34
Se interpola para encontrar el valor de z:
−2−1
zsolución = 1 + −0,58−2,11
(0 − 2,11) = −1,35
Ejercicio 2
Considera la ecuación lineal:
y 00 + x2 y 0 − 2xy = 3 + x
Defina en el intervalo [0,2].En cada apartado se darán unas condiciones para los extremos
del intervalo; se pide encontrar los valores aproximados de la solución de la ecuación bajo
esas condiciones en los nodos x1 = 0,5, x2 = 1 y x3 = 1,5 utilizando el método de diferencias
finitas.
Solución:
1
Se define un número de puntos intermedios y tamaño de paso: n = 3, y h = 2
. Sustituimos en la ecuación:
1 p(xk ) 2 1 p(xk )
( 2
− )wk−1 + (− 2 + q(xk ))wk + ( 2 + )wk−1 = g(xk ), k = 1, 2, 3
h 2h h h 2h
(-4-p(xk ))wk−1 + (−8 + q(xk ))wk + (4 + p(xk ))wk+1 = g(xk ), k = 1, 2, 3
2
Tenemos que:
p(x) = x2 q(x) = −2x g(x) = 3 + x
Entonces:
1 3
x1 = 2
x2 = 1 x3 = 2
1 9
p(xk ) 4
1 4
q(xk ) −1 −2 −3
7 9
g(xk ) 2
4 2
Ejercicio 3
d2 y
Sea la EDO dx2
+ 4y = 8 cos(2x) con valores de frontera y(0) = 0, y(10) = 18,26
Solución:
Sustituimos la 2da derivada:
yi−1 + (4∆x2 − 2)yi + yi+1 = 8∆x2 cos (2x)
Comenzamos las iteraciones:
Primera iteración en i= 1, x=2
y0 + (4 ∗ 22 − 2)y1 + y2 = 8 ∗ 22 cos (2 ∗ 2)14y1 + y2 = −20,9166
Segunda iteración
y1 + (4 ∗ 22 − 2)y2 + y3 = 8 ∗ 22 cos (2 ∗ 4)y1 + 14y2 + y3 = −4,6560
Tercera iteración en i=3, x=6
y2 + (4 ∗ 22 − 2)y3 + y4 = 8 ∗ 22 cos (2 ∗ 6)y2 + 14y3 + y4 = −27,0003
Tercera iteración en i=4, x=8
y3 + (4 ∗ 22 − 2)y4 + y5 = 8 ∗ 22 cos (2 ∗ 8)y3 + 14y4 = −12,3862
3
Si x = 10
y5 = −18,2589
−20,9166
14 1 0 0 y1
1 14 1 0 y −4,656
2
=
0 1 14 1 y3 −27,0003
0 0 1 14 y
4 −12,3962
Solucionando el sistema de ecuaciones:
y1 = −1,4874 y2 = −0,0229 y3 = −1,8683 y4 = −0,7513