Programación de Métodos Numéricos - PMN
INGENIERÍA Semestre 2023 - V
PRACTICA 5 SOLUCIÓN
Hora de inicio: 2:00 pm
Duración: 1hrs45min
Viernes, 17 de Enero de 2023
EL ORDEN Y LA LIMPIEZA INFLUIRÁN EN LA CALIFICACIÓN
SOLO CALCULADORA SIMPLE
Pregunta 1: (10 ptos)
Resolver la ecuación diferencial usando el método de RUNGE KUTTA de 4to orden para la ecuación
diferencial mostrada, considerar los valores iniciales de y’’(2)=-4, y’(2)=3, y(2)=7 un paso de 0.1 y
un tiempo final de 2.5 s. Utilizar una tabla para el desarrollo de sus resultados.
𝑠𝑒𝑛(𝑡)𝑦 ′′′ + cos(𝑡𝑦) + 𝑠𝑒𝑛(𝑡 2 + 𝑦 ′′ ) + (𝑦′)3 = 2𝐿𝑜𝑔10 (𝑡)
t h*X2 X1 h*X3 X2 K1 K2 K3 K4 X3
2 7 3 -4
2.1000 0.3000 7.3000 -0.4000 2.6000 -2.9182 -2.2404 -2.2581 -1.7424 -6.2762
2.2000 0.2600 7.5600 -0.6276 1.9724 -1.7430 -1.1640 -1.1390 -0.7235 -7.4550
2.3000 0.1972 7.7572 -0.7455 1.2269 -0.7276 -0.3812 -0.3614 -0.1448 -7.8479
2.4000 0.1227 7.8799 -0.7848 0.4421 -0.1483 -0.0073 -0.0003 0.0807 -7.8617
2.5000 0.0442 7.9241 -0.7862 -0.3441 0.0797 0.1372 0.1382 0.2104 -7.7216
Pregunta 2: (10 ptos)
Utilice el método implícito para resolver la siguiente EDP a continuación. Utilizar una tabla para
el desarrollo de sus resultados. x=y=1.
𝜕 2 𝑢(𝑥, 𝑦) 𝜕 2 𝑢(𝑥, 𝑦)
+ =0
𝜕𝑥 2 𝜕𝑦 2
Condiciones de frontera:
Página 1 de 4 Tema del documento
Universidad de Piura Programación de Métodos Numéricos - PMN
Programación de Métodos Numéricos - PMN
Página 2 de 4 Tema del documento
Universidad de Piura Programación de Métodos Numéricos - PMN
Programación de Métodos Numéricos - PMN
Página 3 de 4 Tema del documento
Universidad de Piura Programación de Métodos Numéricos - PMN
Programación de Métodos Numéricos - PMN
FORMULARIO
Método Runge Kutta 4to Orden
𝑘1 = 𝑓(𝑥𝑖 , 𝑦𝑖 )
1 1
𝑘2 = 𝑓(𝑥𝑖 + ℎ, 𝑦𝑖 + ℎ𝑘1 ) ℎ
2 2 𝑦𝑖+1 = 𝑦𝑖 + (𝑘1 + 2𝑘2 + 2𝑘3 + 𝑘4 )
1 1 6
𝑘3 = 𝑓(𝑥𝑖 + ℎ, 𝑦𝑖 + ℎ𝑘2 )
2 2
𝑘4 = 𝑓(𝑥𝑖 + ℎ , 𝑦𝑖 + ℎ𝑘3 )
Método de Adams 4to Orden
ℎ
𝑦̂𝑖+1 = 𝑦𝑖 + [55𝑦 ′ 𝑖 − 59𝑦 ′ 𝑖−1 + 37𝑦 ′ 𝑖−2 − 9𝑦 ′ 𝑖−3 ]
24
ℎ
𝑦𝑖+1 = 𝑦𝑖 + [9𝑦̂ ′ 𝑖+1 + 19𝑦 ′ 𝑖 − 5𝑦 ′ 𝑖−1 + 𝑦 ′ 𝑖−2 ]
24
Aproximación diferencia finita
𝑑𝑢 𝑢𝑖+1,𝑗 − 𝑢𝑖,𝑗
=
𝑑𝑥 𝛥𝑥
𝑑 2 𝑢 𝑢𝑖−1,𝑗 − 2𝑢𝑖,𝑗 + 𝑢𝑖+1,𝑗
=
𝑑𝑥 2 (𝛥𝑥)2
Página 4 de 4 Tema del documento