TALLER CORTE 3. MÉTODOS NUMÉRICOS.
DOCENTE: LEANDRO MANUEL PARDO CALDERÓN
Nota: Todos los ejercicios deben incluir la tabla de datos y la solución analítica en la tabla
de comandos, y las gráficas superpuestas de la solución analítica y la numérica.
1. Usando el método de Euler con incremento ℎ = 0.1 para determinar un valor
aproximado de la solución en el rango hasta [0,1] para el problema de valor
inicial. Compruebe la solución dada.
𝑑𝑦
= 1 + 𝑡 − 𝑦; 𝑦(0) = 0; 𝑦(𝑡) = 𝑡
𝑑𝑡
2. Use el método de Euler con incremento ℎ = 0.025 para determinar un valor
aproximado de la solución en el rango hasta [0,1] para el problema de valor
inicial. Compruebe la solución dada.
𝑑𝑦 𝑡
= 𝑡𝑒 −𝑦 + ; 𝑦(0) = 0; 𝑦(𝑡) = 𝑙𝑛(1 + 𝑡 2 )
𝑑𝑡 1 + 𝑡2
3. Aplique el método de Heun ℎ = 0.1 y ℎ = 0.025 para determinar un valor
aproximado de la solución en 0 ≤ 𝑡 ≤ 2 para el problema de valor inicial.
Compare con el valor dado de la solución.
𝑑𝑦
= 1 + 𝑦2 − 𝑡 2; 𝑦(0) = 0; 𝑦(𝑡) = 𝑡
𝑑𝑡
4. Use el método de Runge-Kutta orden 2, variación método del punto medio, con
𝜋
ℎ = 40 para determinar un valor aproximado de la solución del problema de
valor inicial
𝑑𝑦
= 2𝑠𝑒𝑐 2 𝑡 − (1 + 𝑦 2 ); 𝑦(0) = 0
𝑑𝑡
en 0 ≤ 𝑡 ≤ 1. Compare los resultados con la solución 𝑦(𝑡) = tan 𝑡.
5. Repita los cálculos del ejercicio anterior con el método de Runge-Kutta orden 2,
𝜋
variación método de Ralston y ℎ = 160.
6. Aplique el método Runge-Kutta de cuarto orden con incrementos ℎ = 0.1 para
determinar un valor aproximado de la solución en 0 ≤ 𝑡 ≤ 4 para el problema
de valor inicial. Compare con el valor dado de la solución.
𝑑𝑦 𝑦2
= −1 + 2𝑡 + ; 𝑦(0) = 1; 𝑦(𝑡) = 1 + 𝑡 2
𝑑𝑡 (1 + 𝑡 2 )2
7. Resuelva el sistema de ecuaciones mediante el método de Runge-Kutta, en el
rango 0 ≤ 𝑡 ≤ 2, con ℎ = 0.1.
𝑑𝑥
= −𝑥 + 𝑦
𝑑𝑡
𝑑𝑦
= 2𝑥
𝑑𝑡
𝑥(0) = 0, 𝑦(0) = 1
Compare con la solución analítica del sistema
1 1
𝑥 = − 𝑒 −2𝑡 + 𝑒 𝑡
3 3
1 2
𝑦 = 𝑒 −2𝑡 + 𝑒 𝑡
3 3
8. Resuelva la EDO de segundo orden en 0 ≤ 𝑡 ≤ 2
𝑦 ′′ − 3𝑦 ′ − 4𝑦 = 0; 𝑦(0) = 1, 𝑦′(0) = 0. Use ℎ = 0.02
Compare con la solución:
1
𝑦 = (𝑒 4𝑥 + 4𝑒 −𝑥 )
5
9. Use el método de diferencias finitas para resolver la EDP:
𝜕2 𝑢 𝜕𝑢
0.75 𝜕𝑥 2 = ; 𝑢(0, 𝑡) = 120, 𝑢(10, 𝑡) = 70, 𝑢(𝑥, 0) = 50.
𝜕𝑡
Tome ∆𝑥 = 1 y ∆𝑡 tal que el sistema sea estable. Tome 𝑡𝑚𝑎𝑥 tal que se visualicen
20 incrementos temporales.
10. Resuelva la ecuación anterior mediante el esquema implícito.