CÁLCULO NUMÉRICO

EDOs-Resolución de PVI

1
¿QUE ES UNA ECUACIÓN DIFERENCIAL?
Es una ecuación en la que intervienen derivadas de una o más funciones desconocidas.

 Clasificación por tipos:

➢ Ecuaciones diferenciales ordinarias

𝑑𝑦
+ 5𝑦 = 𝑒 𝑥
𝑑𝑥
➢ Ecuaciones en derivadas parciales

𝜕2𝑦 𝛿2𝑢
+ 2 ∗ 𝑥 ∗ 𝑦 ∗ +𝑢 =!
𝜕𝑥 2 𝛿𝑦 2

2
SOLUCIÓN EDOs
 Solución general: una solución de tipo genérico, expresada con una o más constantes. Tiene
un orden de infinitud de acuerdo a su cantidad de constantes.

 Solución particular: Si fijando cualquier punto P(𝑥0 , 𝑦0 ) por donde debe pasar necesariamente
la solución de la ecuación diferencial, existe un único valor de C, y por lo tanto de la curva
integral que satisface la ecuación, éste recibirá el nombre de solución particular de la
ecuación en P(𝑥0 , 𝑦0 ), que recibe el nombre de condición inicial.

 Solución singular: una función que verifica la ecuación, pero que no se obtiene
particularizando la solución genera

3
PROBLEMA DE VALOR INICIAL
Un problema de valor inicial es una ecuación diferencial ordinaria que tiene un
prerrequisito inicial sujeto con la misma. Este prerrequisito es la salida de la función indefinida
para algún valor que se encuentra dentro del dominio de la ecuación diferencial dada.
𝑑𝑛 𝑦
= 𝑓 𝑥, 𝑦, 𝑦 ′ , … , 𝑦 𝑛−1
𝑑𝑛 𝑥

Con condiciones:
𝑦 𝑥0 = 𝑦0 , 𝑦 ′ 𝑥0 = 𝑦1 , … , 𝑦 𝑛−1 𝑥0 = 𝑦𝑛−1

4
MÉTODOS NUMÉRICOS
 Estos procedimientos utilizados para encontrar aproximaciones numéricas a las soluciones
de ecuaciones diferenciales ordinarias (EDO) son lo que llamamos Métodos Numéricos

 Métodos de un paso:
Los métodos de un paso tienen por objetivo obtener una aproximación de la solución
de un problema bien planteado de valor inicial en cada punto, basándose en el
resultado obtenido para el punto anterior.

 Métodos de pasos múltiples:

Los métodos de varios pasos intentan obtener eficiencia manteniendo y utilizando la
información de los pasos anteriores, en lugar de descartarla. Por consiguiente, se refieren
a distintos puntos anteriores y a los valores de sus derivadas.

5
MÉTODOS DE UN PASO
Euler
4to Orden

METODOS DE 3er Orden Taylor

RUNGE-KUTTA
Ralston Heun

Euler
Mejorado /
Punto medio

6
 MÉTODO EULER

Suponemos que el PVI tiene solución

y ′ = 𝑓 𝑥, 𝑦 , 𝑦 𝑥0 = 𝑦0

Usamos la linealización de una solución incógnita y(x) en 𝑥 = 𝑥0

𝐿 𝑥 = 𝑦0 +𝑓(𝑥0 , 𝑦0 ) 𝑥 − 𝑥0

Sea h un incremento positivo del eje x tal que 𝑥1 = 𝑥0 + ℎ

𝑦1 = 𝑦0 +𝑓(𝑥1 , 𝑦1 )ℎ
𝑦1 ≈ 𝑦 (𝑥1 )

7
 MÉTODO EULER

Repitiendo el procedimiento usando una segunda recta tangente en (𝑥1 , 𝑦1 ), obtenemos una
aproximación 𝑦2 ≈ 𝑦 (𝑥2 )

Entonces, Se puede definir recursivamente

𝑦𝑛+1 = 𝑦𝑛 + ℎ𝑓 𝑥𝑛 , 𝑦𝑛

Donde
𝑥𝑛 = 𝑥0 + 𝑛ℎ

Este procedimiento de uso sucesivo de rectas tangentes es el que denominamos “método de Euler

8
MÉTODO EULER: INTERPRETACIÓN GEOMÉTRICA

9
VINCULO: https://www.geogebra.org/m/dvrjbvwy
 𝑦 𝑡

ANIMACION PARA REPRESENTACION GRAFICA (NO PROPIA)

 𝑦 𝑡

𝑦 ′ = 𝑓 𝑡, 𝑦
𝑦 𝑎 =𝛼

ANIMACION PARA REPRESENTACION GRAFICA (NO PROPIA)

 𝑦 𝑡

𝑦 ′ = 𝑓 𝑡, 𝑦
𝑦 𝑎 =𝛼

ANIMACION PARA REPRESENTACION GRAFICA (NO PROPIA)

 𝑦 𝑡

𝑦 ′ = 𝑓 𝑡, 𝑦
𝑦 𝑎 =𝛼

Recta 𝑦 = 𝑥𝑓 𝑎, 𝛼 − 𝑎𝑓 𝑎, 𝛼 + 𝛼

ANIMACION PARA REPRESENTACION GRAFICA (NO PROPIA)

 𝑦 𝑡

𝑦 ′ = 𝑓 𝑡, 𝑦
𝑦 𝑎 =𝛼

Recta 𝑦 = 𝑥𝑓 𝑎, 𝛼 − 𝑎𝑓 𝑎, 𝛼 + 𝛼

𝑎 𝑡1

ANIMACION PARA REPRESENTACION GRAFICA (NO PROPIA)

 𝑦 𝑡

𝑦 ′ = 𝑓 𝑡, 𝑦
𝑦 𝑎 =𝛼

Recta 𝑦 = 𝑥𝑓 𝑎, 𝛼 − 𝑎𝑓 𝑎, 𝛼 + 𝛼
𝑤1
𝛼

𝑎 𝑡1

ANIMACION PARA REPRESENTACION GRAFICA (NO PROPIA)

 𝑦 𝑡
𝑦 ′ = 𝑓 𝑡, 𝑦
𝑦 𝑎 =𝛼
𝑤N
…

𝑤2
𝑤1
𝛼

𝑎 𝑡1 𝑡2 … 𝑡N = 𝑏

ANIMACION PARA REPRESENTACION GRAFICA (NO PROPIA)

 𝑦 𝑡
𝑦 ′ = 𝑓 𝑡, 𝑦
𝑦 𝑎 =𝛼
𝑤N
…

𝑤2
𝑤1
𝛼

𝑎 𝑡1 𝑡2 … 𝑡N = 𝑏

ANIMACION PARA REPRESENTACION GRAFICA (NO PROPIA)

ANALISIS PARA EL ERROR EN EULER
La solución numérica de las EDO implica dos tipos de error

1. Errores de redondeo causados por el número limitado de cifras significativas que una
computadora puede retener.

2. Errores de truncamiento originados por la naturaleza de las técnicas empleadas para

aproximar los valores de y. se componen en 2 partes:
I. un error de truncamiento local, resulta de una aplicación del método considerado, en un solo paso.

II. un error de truncamiento propagado, resulta de las aproximaciones producidas durante los pasos
previos.

La suma de ambos es el error de truncamiento global o total

18
 ERRORES DE TRUNCAMIENTO PARA EL MÉTODO
DE EULER

 Para deducir una fórmula para el error de truncamiento local del método de Euler, se usa la
fórmula de Taylor con residuo.

 Si una función y(x) tiene k+1 derivadas que son continuas en un intervalo abierto que contiene 𝑎 a
y a 𝑥, entonces:

𝑘
𝑥−𝑎
′
𝑥−𝑎 𝑥 − 𝑎 𝑘+1
𝑦 𝑥 =𝑦 𝑎 +𝑦 𝑎 + ⋯ + 𝑦𝑘 𝑎 +𝑦 𝑘+1
𝑐
1! 𝑘! (𝑘 + 1)!

Donde 𝑎 ≤ 𝑐 ≤ 𝑥
 ERRORES DE TRUNCAMIENTO PARA EL MÉTODO
DE EULER

Si establecemos 𝑘 = 1 , 𝑎 = 𝑥𝑛 , 𝑥 = 𝑥𝑛+1 = 𝑥𝑛 + ℎ
Se obtiene
ℎ ℎ2
𝑦 𝑥𝑛+1 = 𝑦 𝑥𝑛 + 𝑦 ′ 𝑥𝑛 +𝑦 ′′ 𝑐
1! 2!
ℎ2
𝑦 𝑥𝑛+1 = 𝑦𝑛 + ℎ𝑓 𝑥𝑛 , 𝑦𝑛 + 𝑦 ′′ 𝑐 ቋ 𝑀É𝑇𝑂𝐷𝑂 𝐷𝐸 𝐸𝑈𝐿𝐸𝑅
2!
El error de truncamiento local en 𝑦𝑛+1 es
ℎ2
𝑦 ′′ 𝑐 𝑑𝑜𝑛𝑑𝑒 𝑥𝑛 < 𝑐 < 𝑥𝑛+1
2!
 𝑒(ℎ) es de orden ℎ𝑛 , denotado con 𝑂 ℎ𝑛 , si existe una constante 𝐶 y un entero positivo 𝑛 tal que 𝑒 ℎ ≤ 𝐶ℎ𝑛
para ℎ suficientemente pequeña.
 Por lo que el error de truncamiento local para el método de euler es 𝑂(ℎ2 )
 MÉTODOS DE UN PASO

4to Orden Euler

METODOS DE 3er Orden Taylor

RUNGE-KUTTA
Ralston Heun

Euler
Mejorado /
Punto medio

21
 MÉTODO DE TAYLOR

Si 𝑓 𝑒𝑠 𝑢𝑛𝑎 𝑓𝑢𝑛𝑐𝑖𝑜𝑛 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑎𝑚𝑒𝑛𝑡𝑒 𝑑𝑖𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎𝑏𝑙𝑒 𝑒𝑛 𝑎, 𝑏 entonces se tiene que:

𝑦 ′′ 𝑡𝑖 𝑡𝑖+1 −𝑡𝑖 2 𝑦 ′′′ 𝑡𝑖 𝑡𝑖+1 −𝑡𝑖 3 𝑦 𝑛 𝑡𝑖 𝑡𝑖+1 −𝑡𝑖 𝑛

𝑦 𝑡𝑖+1 = 𝑦 𝑡𝑖 + 𝑦′(𝑡𝑖 )(𝑡𝑖+1 − 𝑡𝑖 ) + + + …+
2 6 𝑛!

Se aproxima 𝑦𝑖 como:

𝑤0 = 𝑦0
2 ′
ቐ ℎ 𝑓 𝑡𝑖 , 𝑤𝑖 ℎ𝑛 𝑓 (𝑛−1) 𝑡𝑖 , 𝑤𝑖
𝑤𝑖+1 = 𝑤𝑖 + ℎ 𝑓 𝑡𝑖 , 𝑤𝑖 + +⋯+ ; 𝑖 = 0, 1, 2, … , 𝑛
2 𝑛!

𝜕 𝑓 (𝑡,𝑤) 𝜕 𝑓 (𝑡,𝑤) 𝑑𝑤
Donde 𝑦 ′ = 𝑓 𝑡, 𝑤 𝑦 𝑦 ′′ = 𝑓 ′ 𝑡, 𝑤 = + ∙ 22
𝜕𝑡 𝜕𝑤 𝑑𝑡
 MÉTODOS DE RUNGE-KUTTA

 Los métodos de Runge-Kutta son generalizaciones de la fórmula básica de Euler en que la función pendiente f se
reemplaza por un promedio ponderado de pendientes en el intervalo x x xn

𝑦𝑖+1 = 𝑦𝑖 + 𝜃 𝑥𝑖 , 𝑦𝑖 , ℎ ℎ

 Donde 𝜃 𝑥𝑖 , 𝑦𝑖 , ℎ se conoce como función incremento, la cual puede interpretarse como una pendiente
representativa en el intervalo. La función incremento se escribe en forma general como.
𝜃 = 𝑎1 𝑘1 + 𝑎2 𝑘2 + ⋯ + 𝑎𝑛 𝑘𝑛
Donde las a son constantes y las k
𝑘1 = 𝑓 𝑥𝑖 , 𝑦𝑖
𝑘2 = 𝑓 𝑥𝑖 + 𝑝1 ℎ, 𝑦𝑖 + 𝑞11 𝑘1 ℎ
𝑘3 = 𝑓 𝑥𝑖 + 𝑝2 ℎ, 𝑦𝑖 + 𝑞21 𝑘1 ℎ+𝑞22 𝑘2 ℎ
𝑘𝑛 = 𝑓 𝑥𝑖 + 𝑝𝑛−1 ℎ, 𝑦𝑖 + 𝑞𝑛−1 𝑘1 ℎ+𝑞𝑛−1 ,2 𝑘2 ℎ + ⋯ +𝑞𝑛−1 ,𝑛−1 𝑘𝑛−1 ℎ
23
Donde las p y las q son constantes
TÉCNICAS DE RUNGE-KUTTA: PRIMER ORDEN

▪ Al tomar 𝑚 = 1, 𝑤1 = 1 y 𝑘1 = 𝑓(𝑥𝑛 , 𝑦𝑛 )

▪ Obtenemos la conocida como fórmula de

Euler

𝑦𝑛+1 = 𝑦𝑛 + ℎ𝑓 𝑥𝑛 , 𝑦𝑛
 MÉTODOS DE RUNGE-KUTTA: 2do ORDEN
El método de RUNGE-KUTTA en versión de segundo orden es

𝑦𝑖+1 = 𝑦𝑖 + 𝑎1 𝑘1 , 𝑎2 𝑘2 ℎ
Donde:

𝑘1 = 𝑓 𝑥𝑖 , 𝑦𝑖
𝑘2 = 𝑓 𝑥𝑖 + 𝑝1 ℎ, 𝑦𝑖 + 𝑞11 𝑘1 ℎ

Los valores a1,a2,p1 y q11 se evalúan al igualar la expansión de la serie de Taylor hasta el termino de segundo orden. Al hacerlo
desarrollamos 3 ecuaciones con 4 incógnitas:

𝑎1 + 𝑎2 = 1
1
𝑎2 𝑝1 =
2
1
𝑎2 𝑞11 =
2

Decidimos darle un valor a 𝑎2 para así sacar los valores de las otras incógnitas. Dependiendo del valor que le damos a 𝑎2 es el método
que vamos a utilizar para resolver
25
MÉTODOS DE UN PASO

4to Orden Euler

METODOS DE 3er Orden Taylor

RUNGE-KUTTA
Ralston
Heun
Euler
Mejorado /
Punto medio

26
MÉTODOS DE RUNGE-KUTTA: 2do ORDEN
 Método de Heun cuando 𝑎2 = 1/2:
Suponiendo que 𝑎1 = 1/2 y que 𝑝1 = 𝑞11 = 1 sustituyendo en la ecuación anterior:

1 1
𝑦𝑖+1 = 𝑦𝑖 + 𝑘 , 𝑘 ℎ
2 1 2 2
Donde

𝑘1 = 𝑓 𝑥𝑖 , 𝑦𝑖
𝑘2 = 𝑓 𝑥𝑖 + 𝑝1 ℎ, 𝑦𝑖 + 𝑞11 𝑘1 ℎ

Apreciando el grafico a continuación

27
 Gráfico runge-kutta de orden 2

VINCULO
GEOGEBRA

28
Métodos numéricos para ingenieros - 5ª edición-PAGINA N° 25
29
MÉTODOS DE UN PASO

4to Orden Euler

METODOS DE 3er Orden Taylor

RUNGE-KUTTA
Ralston Heun

Euler
Mejorado /
Punto medio
30
MÉTODOS DE RUNGE-KUTTA: 2do ORDEN
 Método Euler Mejorado o Método del punto medio cuando 𝑎2 = 1:

Suponiendo que 𝑎1 = 0 y que 𝑝1 = 𝑞11 = 1/2 sustituyendo en la ecuación anterior:

1
𝑦𝑖+1 = 𝑦𝑖 + 0𝑘1 , 𝑘2 ℎ
2
Donde
𝑘1 = 𝑓 𝑥𝑖 , 𝑦𝑖
1 1
𝑘2 = 𝑓 𝑥𝑖 + ℎ, 𝑦𝑖 𝑘1 ℎ
2 2
Se lo define como método de punto medio que podemos apreciarlo gráficamente a
continuación

31
GRAFICO DE MéTODO DE PUNTO MEDIO

VINCULO
GEOGEBRA

32
MÉTODOS DE UN PASO

4to Orden Euler

METODOS DE 3er Orden Taylor

RUNGE-KUTTA
Heun

Ralston
Euler
Mejorado /
Punto medio

33
MÉTODOS DE RUNGE-KUTTA: 2do ORDEN
 Método de Ralston cuando 𝑎2 = 2/3:

Suponiendo que 𝑎1 = 1/3 y que 𝑝1 = 𝑞11 = 3/4 sustituyendo en la ecuación anterior:

1 2
𝑦𝑖+1 = 𝑦𝑖 + 𝑘 , 𝑘2 ℎ
3 1 3
Donde
𝑘1 = 𝑓 𝑥𝑖 , 𝑦𝑖
𝑘2 = 𝑓 𝑥𝑖 + 𝑝1 ℎ, 𝑦𝑖 + 𝑞11 𝑘1 ℎ
Se lo define como método de Ralston que podemos apreciarlo gráficamente a continuación

34
MÉTODOS DE UN PASO
Euler

4to Orden

Taylor
METODOS DE 3er Orden

RUNGE-KUTTA

Heun
Ralston

Euler
Mejorado /
Punto medio 35
Comparación de la solución
verdadera con soluciones
numéricas usando tres
métodos RK de segundo
orden y el método de Euler

36
Métodos numéricos para ingenieros - 5ª edición-CAPITULO N°25
MÉTODOS DE UN PASO
4to Orden

Euler

METODOS DE 3er Orden Taylor

RUNGE-KUTTA
Ralston Heun

Euler
Mejorado /
Punto medio

37
MÉTODOS DE RUNGE-KUTTA: 4to ORDEN
 Un procedimiento de Runge-Kutta de cuarto orden consiste en determinar
parámetros de modo que la fórmula

𝑦𝑖+1 = 𝑦𝑖 + 𝑤1 𝑘1 + 𝑤2 𝑘2 + 𝑤3 𝑘3 + 𝑤4 𝑘4 ℎ
Donde
𝑘1 = 𝑓 𝑥𝑖 , 𝑦𝑖
𝑘2 = 𝑓 𝑥𝑖 + 𝛼1 ℎ, 𝛽1 𝑘1 ℎ
𝑘3 = 𝑓 𝑥𝑖 + 𝛼2 ℎ, 𝛽2 𝑘1 ℎ + 𝛽3 𝑘2 ℎ
𝑘4 = 𝑓 𝑥𝑖 + 𝛼3 ℎ, 𝛽4 𝑘1 ℎ + 𝛽5 𝑘2 ℎ + 𝛽6 𝑘3 ℎ

Concuerda con el polinomio de Taylor de grado 4, donde nos da un

sistema de ecuaciones de 11 ecuaciones con 13 incógnitas

38
ECUACIONES DE RK 4orden
ℎ
𝑦𝑖+1 = 𝑦𝑖 + 𝑘1 + 𝑘2 + 𝑘3 + 𝑘4
6

𝑘1 = 𝑓 𝑥𝑖 , 𝑦𝑖
1
𝑘2 = 𝑓 𝑥𝑖 + ℎ, 𝑦𝑖 + 1/2𝑘1 ℎ
2
1
𝑘3 = 𝑓 𝑥𝑖 + ℎ, 𝑦𝑖 + 1/2𝑘2 ℎ
2
1
𝑘4 = 𝑓 𝑥𝑖 + ℎ, 𝑦𝑖 + 1/2𝑘3 ℎ
2

Mientras que las otras formulas de cuarto orden se deducen con facilidad, este es
el algoritmo resumido
39
Gráfico runge-kutta de orden 4

40
VINCULO GEOGEBRA
MÉTODOS DE PASOS MÚLTIPLES
Guardan información de pasos anteriores para obtener de manera más efectiva la trayectoria de
la solución

 Método de Heun sin autoinicio: En el caso de Heun sin autoinicio, utiliza la información de 2
puntos anteriores, para la obtención del predictor. El corrector es el mismo de Heun de 1
paso.

 Método de Adams-Bashforth: La idea fundamental del método de Adams-Bashforth de n

pasos es usar un polinomio de interpolación de f(t, y(1)) que pasa por los n puntos: (t i ,f),(t i 1
,f 1 ) ,...,(t n+1 ,f i-n+1 ) .

 Método de Adams-Moulton: La idea fundamental del método de Adams-Moulton de n pasos

es usar un polinomio de interpolación de f(t, y(t)) que pasa por los n + 1 puntos. (t 1,1 ,1 1,1 ) (t
i f j ),...,(t n+1 ,f i-n+1 ) .

41
BIBLIOGRAFIA
 “Análisis Numérico” BURDEN, Richard L.& FAIRES, J. Douglas Grupo Editorial Iberoamericana
1998 6ta ed.
 “Métodos Numéricos con Mathematica” GARCIA RAFFI, L. et. all Alfaomega Grupo Editor.
2005 1ra ed.
 "Ecuaciones diferenciales con Aplicaciones de Modelado ZILL, Dennis G. & CULLEN, Michael R.
McGraw-Hill 2.009
 “Métodos Numéricos y Computación” CHENEY, W. & KINCAID, D. Cengage Learning 2.011
9na ed.
 "Métodos Numéricos Aplicados con Software" NAKAMURA, SHOICHIRO Pearson Education
1.992 1ra ed.
 CalcLabs With Mathematica HOLLYS S. Cengage Learning 2.012 5ta ed.

42