AREA DE CONOCIMIENTO: ANÁLISIS FUNCIONAL NRC: 3207

Unidad III (Derivación, Integración y Ecuaciones Diferenciales

Versión 1.1
Ordinarias)

GRUPO 05: Montalvo Kevin, Solano Erick, Villafuerte Jéssica

EMAIL: [email protected]

Tema:
Método de Euler Predictor - Corrector

Análisis Matemático
2.1. Demostraciones Preliminares
Preliminar 1: Polinomio de Lagrange

AFIRMACIÓN JUSTIFICACIÓN

𝑛𝑛 Es simplemente una reformulación del

𝑓𝑓𝑛𝑛 (𝑥𝑥 ) = � 𝐿𝐿𝑖𝑖 (𝑥𝑥)𝑓𝑓(𝑥𝑥) polinomio de Newton que evita el cálculo
𝑖𝑖=0 por diferencias divididas.

𝑛𝑛
𝑥𝑥 − 𝑥𝑥𝑗𝑗
𝐿𝐿𝑖𝑖 (𝑥𝑥 ) = � Donde pi designa el “producto de”. Por
𝑥𝑥𝑖𝑖 − 𝑥𝑥𝑗𝑗 ejemplo, la versión lineal (n=1)
𝑗𝑗=0
𝑗𝑗≠1
 El resultado de la formulación de Lagrange
se puede captar de manera directa al darse
cuenta que cada termino 𝐿𝐿𝑖𝑖 (𝑥𝑥) será 1 en 𝑥𝑥 =
𝑥𝑥𝑖𝑖 y 0 en todos los otros puntos de la
muestra. De esta manera, cada producto
𝑥𝑥 − 𝑥𝑥𝑖𝑖 𝑥𝑥 − 𝑥𝑥0 𝐿𝐿𝑖𝑖 (𝑥𝑥 )𝑓𝑓(𝑥𝑥𝑖𝑖 ) toma el valor de 𝑓𝑓(𝑥𝑥𝑖𝑖 ) en el
𝑓𝑓1 (𝑥𝑥 ) = 𝑓𝑓 (𝑥𝑥0 ) + 𝑓𝑓(𝑥𝑥1 )
𝑥𝑥0 − 𝑥𝑥1 𝑥𝑥1 − 𝑥𝑥0 punto de muestra 𝑥𝑥𝑖𝑖 . En consecuencia, la
sumatoria de todos los productos designados
para la ecuación de Lagrange es el único
polinomio de n-ésimo orden que pasa de
manera exacta a través de todos los n+1
puntos.
Puesto que el polinomio de interpolación
𝑛𝑛
correspondiente a un polinomio p de Pn es
el propio polinomio, las componentes de p
𝑝𝑝(𝑥𝑥 ) = � 𝑙𝑙𝑖𝑖 (𝑥𝑥 )𝑝𝑝(𝑥𝑥𝑖𝑖 ) en esta base, son los valores del polinomio
𝑖𝑖=0 en los nodos. Es decir, su expresión en
términos de las funciones de L, es :
lo que prueba que L es un sistema de
generadores de Pn. Puesto que hay n + 1
funciones de forma y la dimensión de Pn es
𝑛𝑛 n + 1, se concluye que L es una base de Pn.
𝑥𝑥 = � 𝑙𝑙𝑙𝑙 (𝑥𝑥 )𝑝𝑝(𝑥𝑥𝑖𝑖𝑗𝑗 )
𝑗𝑗 De la unicidad de las componentes de un
𝑖𝑖=0 vector en una base se desprende la unicidad
de L. Si se expresan los elementos de la base
{1, x, ..., xn} de Pn en términos de la base
de Lagrange, se obtiene:

1 1…1 1
𝑙𝑙0
𝑥𝑥0 𝑥𝑥1 … 𝑥𝑥𝑛𝑛 𝑥𝑥 En forma matricial, las relaciones anteriores
�⋮ ⋮ ⋱ ⋮� �𝑙𝑙 1 � = � �
⋮ ⋮ se escriben como
𝑥𝑥0𝑛𝑛 𝑥𝑥1𝑛𝑛 … 𝑥𝑥𝑛𝑛𝑛𝑛 𝑙𝑙𝑛𝑛 𝑥𝑥 𝑛𝑛

La matriz de coeficientes A de este sistema

es la matriz de Vandermonde del conjunto
|𝐴𝐴𝑖𝑖 |
𝑙𝑙𝑖𝑖 (𝑥𝑥 ) = {,x_1,...,x_n} y por consiguiente es
|𝐴𝐴|
invertible. Los elementos de la base de
Lagrange están dados por;
 donde Ai representa la matriz de coeficientes en la
𝑙𝑙𝑖𝑖 �𝑥𝑥𝑗𝑗 � = 𝛿𝛿𝑖𝑖𝑖𝑖 que se ha sustituido la columna i-ésima por el
término independiente. Así pues, li verifica

Donde δij representa la delta de Kronecker

(1 si i = j y 0 si i ≠ j). Estas relaciones pueden
(𝑥𝑥 − 𝑥𝑥0 )(𝑥𝑥 − 𝑥𝑥1 ) … (𝑥𝑥 − 𝑥𝑥𝑖𝑖−1 )(𝑥𝑥 − 𝑥𝑥𝑖𝑖−1 ) … (𝑥𝑥 − 𝑥𝑥𝑛𝑛 ) ser consideradas en si mismas como un problema
𝑙𝑙𝑖𝑖 (𝑥𝑥 ) =
(𝑥𝑥𝑖𝑖 − 𝑥𝑥0 )(𝑥𝑥𝑖𝑖 − 𝑥𝑥1 ) … (𝑥𝑥𝑖𝑖 − 𝑥𝑥𝑖𝑖−1 )(𝑥𝑥𝑖𝑖 − 𝑥𝑥𝑖𝑖−1 ) … (𝑥𝑥𝑖𝑖 − 𝑥𝑥𝑛𝑛 ) de interpolación y consecuentemente determinan
unívocamente las funciones básicas. De hecho, la
siguiente expresión de las funciones básicas

𝑛𝑛 El conocimiento explıcito de la base de Lagrange

𝑝𝑝(𝑥𝑥 ) = � 𝑙𝑙𝑖𝑖 (𝑥𝑥 )𝑝𝑝(𝑥𝑥𝑖𝑖 ) permite, dada una función cualquiera f, construir un
𝑖𝑖=0 polinomio de grado menor o igual que n

Que interpole a f en n + 1 puntos x_0,x_1,...,x_n.

Conviene recordar que la solución del problema de
𝑛𝑛
interpolación es la solución del problema de mejor
〈𝑓𝑓, 𝑔𝑔〉 = � 𝑓𝑓(𝑥𝑥𝑖𝑖 )𝑔𝑔(𝑥𝑥𝑖𝑖 )
aproximación de una función en el espacio de los
𝑖𝑖=0
polinomios de grado n respecto a la seminormal
definida por el producto discreto

𝑛𝑛

〈𝑙𝑙𝑖𝑖 , 𝑙𝑙𝑗𝑗 〉 = � 𝛿𝛿𝑖𝑖𝑖𝑖 𝛿𝛿𝑘𝑘𝑘𝑘 = 𝛿𝛿𝑖𝑖𝑖𝑖 La base de Lagrange es ortonormal respecto a 〈 , 〉 .

𝑖𝑖=0 En efecto, se cumple que

Fig.1: Polinomio de Lagrange

 2.2. Demostraciones Específicas

Especifica 1: Método de Euler

AFIRMACIÓN JUSTIFICACIÓN

𝑦𝑦′ = 𝑓𝑓(𝑥𝑥, 𝑦𝑦) (1)

Disponemos de las condiciones (1) y (2)
𝑦𝑦(𝑥𝑥𝑜𝑜 ) = 𝑦𝑦𝑜𝑜 (2)

Curva aproximada, donde 𝑦𝑦𝑜𝑜 y 𝑥𝑥𝑜𝑜 son las

condiciones iniciales y h la distancia entre cada
punto en el eje x

Fig 1. Curva solución exacta y curva solución del método numérico

ℎ = 𝑥𝑥𝑖𝑖+1 − 𝑥𝑥𝑖𝑖 Definimos el tamaño del paso

𝑑𝑑𝑑𝑑 La derivada por definición es la pendiente en un

𝑦𝑦 ′ = (3)
𝑑𝑑𝑑𝑑 punto

𝑑𝑑𝑑𝑑
= 𝑓𝑓(𝑥𝑥, 𝑦𝑦) Igualando (3) y (1)
𝑑𝑑𝑑𝑑

∆𝑦𝑦 𝑦𝑦𝑖𝑖+1 − 𝑦𝑦𝑖𝑖 Por otro lado tenemos la variación entre y y x

= (4)
∆𝑥𝑥 𝑥𝑥𝑖𝑖+1 − 𝑥𝑥𝑖𝑖
 𝑦𝑦𝑖𝑖+1 − 𝑦𝑦𝑖𝑖 Igualando las ecuaciones (4) y (1) y sustituyendo la
𝑓𝑓(𝑥𝑥, 𝑦𝑦) = (5)
ℎ variación en x por h

Despejando 𝑦𝑦𝑖𝑖+1 , finalmente obtenemos la formula

𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + ℎ𝑓𝑓(𝑥𝑥, 𝑦𝑦)
general del Método de Euler

Método Euler
Calculo Valores en Calculo valores en
x y
x1=xo+h y1=yO+hf(xo,yo)
x2=x1+h y2=y1+hf(x1,y1)
x3=x2+h y3=y2+hf(x2,y2)
x4=x3+h y4=y3+hf(x3,y3) Aplicación del método, como se observa primero
x5=x4+h y5=y4+hf(x4,y4) se debe calcular el valor de x y luego el valor de y
el número de iteraciones dependerá de ha y en
x6=x5+h y6=y5+hf(x5,y5) intervalo dado
x7=x6+h y7=y6+hf(x6,y6)
x8=x7+h y8=y7+hf(x7,y7)
. .
. .
. .
Tabla 1. Ejemplo de aplicación del método con valores x e y

Especifica 2: Método de Euler Predictor

AFIRMACIÓN JUSTIFICACIÓN

𝑓𝑓(𝑡𝑡, 𝑦𝑦(𝑡𝑡))
Polinomio que pasa por “n”
(𝑡𝑡𝑖𝑖 , 𝑓𝑓𝑖𝑖 ), (𝑡𝑡𝑖𝑖−1 , 𝑓𝑓𝑖𝑖−1 ),…,(𝑡𝑡𝑖𝑖−𝑛𝑛+1 , 𝑓𝑓𝑖𝑖−𝑛𝑛+1 ) puntos

𝑦𝑦 ′ = 𝑓𝑓(𝑡𝑡, 𝑦𝑦) Resuelve EDO

𝑡𝑡𝑖𝑖+1 𝑡𝑡𝑖𝑖+1

� 𝑦𝑦′(𝑡𝑡)𝑑𝑑𝑑𝑑 = � 𝑓𝑓�𝑡𝑡, 𝑦𝑦(𝑡𝑡)�𝑑𝑑𝑑𝑑

𝑡𝑡𝑖𝑖 𝑡𝑡𝑖𝑖
Integrando desde t i hasta
𝑡𝑡𝑖𝑖+1 t i+1
𝑦𝑦𝑖𝑖+1 − 𝑦𝑦𝑖𝑖 = � 𝑓𝑓�𝑡𝑡, 𝑦𝑦(𝑡𝑡)�𝑑𝑑𝑑𝑑
𝑡𝑡𝑖𝑖
 𝑡𝑡𝑖𝑖+1

𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � 𝑓𝑓�𝑡𝑡, 𝑦𝑦(𝑡𝑡)�𝑑𝑑𝑑𝑑

Despejando 𝑦𝑦𝑖𝑖+1
𝑡𝑡𝑖𝑖

Aproximando 𝑓𝑓�𝑡𝑡, 𝑦𝑦(𝑡𝑡)�

mediante el polinomio de
interpolación que pasa por
(𝑡𝑡𝑖𝑖 , 𝑓𝑓𝑖𝑖 ),
(𝑡𝑡𝑖𝑖 − 𝑡𝑡 )𝑓𝑓𝑖𝑖−1 + (𝑡𝑡 − 𝑡𝑡𝑖𝑖−1 )𝑓𝑓𝑖𝑖 los puntos
𝑃𝑃(𝑡𝑡 ) = (𝑡𝑡𝑖𝑖−1 , 𝑓𝑓𝑖𝑖−1 ), donde 𝑓𝑓𝑖𝑖−1 =
ℎ 𝑓𝑓(𝑡𝑡𝑖𝑖−1 , 𝑦𝑦(𝑡𝑡𝑖𝑖−1 ))
Y 𝑓𝑓𝑖𝑖 = 𝑓𝑓(𝑡𝑡𝑖𝑖 , 𝑦𝑦(𝑡𝑡𝑖𝑖 )) el
polinomio interpolante esta
dado por:
𝑡𝑡𝑖𝑖+1

𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � 𝑃𝑃(𝑡𝑡)𝑑𝑑𝑑𝑑

𝑡𝑡𝑖𝑖

Reemplazando:
𝑡𝑡𝑖𝑖+1
(𝑡𝑡𝑖𝑖 − 𝑡𝑡)𝑓𝑓𝑖𝑖−1 + (𝑡𝑡 − 𝑡𝑡𝑖𝑖−1 )𝑓𝑓𝑖𝑖
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � � 𝑑𝑑𝑑𝑑
ℎ
𝑡𝑡𝑖𝑖

𝑡𝑡𝑖𝑖+1
(𝑡𝑡𝑖𝑖 )𝑓𝑓𝑖𝑖−1
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � � 𝑑𝑑𝑑𝑑
ℎ
𝑡𝑡𝑖𝑖
𝑡𝑡𝑖𝑖+1 𝑡𝑡𝑖𝑖+1 𝑡𝑡𝑖𝑖+1
(𝑡𝑡)𝑓𝑓𝑖𝑖−1 (𝑡𝑡)𝑓𝑓𝑖𝑖 (𝑡𝑡𝑖𝑖−1 )𝑓𝑓𝑖𝑖
−� � � 𝑑𝑑𝑑𝑑 + � � � 𝑑𝑑𝑑𝑑 − � � � 𝑑𝑑𝑑𝑑
ℎ ℎ ℎ
𝑡𝑡𝑖𝑖 𝑡𝑡𝑖𝑖 𝑡𝑡𝑖𝑖

(𝑡𝑡𝑖𝑖 )𝑓𝑓𝑖𝑖−1 𝑓𝑓 Integrando:

𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � (𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 ) − � 𝑖𝑖−1 � (𝑡𝑡𝑖𝑖+1 2 − 𝑡𝑡12 )
ℎ 2ℎ
𝑓𝑓 (𝑡𝑡𝑖𝑖−1 )𝑓𝑓𝑖𝑖
+ � 𝑖𝑖 � (𝑡𝑡𝑖𝑖+1 2 − 𝑡𝑡12 ) − � � (𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 )
2ℎ ℎ
 (𝑡𝑡𝑖𝑖 ) 1
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + �𝑓𝑓𝑖𝑖−1 � �� � (𝑡𝑡 2 − 𝑡𝑡12 )�
� (𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 ) − �
ℎ 2ℎ 𝑖𝑖+1
1 (𝑡𝑡𝑖𝑖−1 )
+ (𝑓𝑓𝑖𝑖 ) �� � (𝑡𝑡𝑖𝑖+1 2 − 𝑡𝑡12 ) − � � (𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 )�
2ℎ ℎ

2𝑡𝑡𝑖𝑖 (𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 ) − 𝑡𝑡𝑖𝑖+1 2 + 𝑡𝑡12

𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + �𝑓𝑓𝑖𝑖−1 � � �
2ℎ
𝑡𝑡𝑖𝑖+1 2 − 𝑡𝑡12 − 2𝑡𝑡𝑖𝑖−1 𝑡𝑡𝑖𝑖+1 + 2𝑡𝑡𝑖𝑖−1 𝑡𝑡𝑖𝑖
+ (𝑓𝑓𝑖𝑖 ) � �
2ℎ

𝑦𝑦𝑖𝑖+1
Agrupando términos
= 𝑦𝑦𝑖𝑖
��𝑓𝑓𝑖𝑖−1 �(2𝑡𝑡𝑖𝑖 𝑡𝑡𝑖𝑖+1 −𝑡𝑡𝑖𝑖+1 2 − 𝑡𝑡12 ) + (𝑓𝑓𝑖𝑖 )(𝑡𝑡𝑖𝑖+1 2 − 𝑡𝑡12 − 2𝑡𝑡𝑖𝑖+1 𝑡𝑡𝑖𝑖−1 + 2𝑡𝑡𝑖𝑖−1 𝑡𝑡𝑖𝑖 )�
+
2ℎ

1
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � ��𝑓𝑓𝑖𝑖−1 �(2𝑡𝑡𝑖𝑖 𝑡𝑡𝑖𝑖+1 −𝑡𝑡𝑖𝑖+1 2 − 𝑡𝑡12 )
2ℎ
+ (𝑓𝑓𝑖𝑖 )(𝑡𝑡𝑖𝑖+1 2 − 𝑡𝑡12 − 2𝑡𝑡𝑖𝑖+1 𝑡𝑡𝑖𝑖−1 + 2𝑡𝑡𝑖𝑖−1 𝑡𝑡𝑖𝑖 )�

1 Sumamos y restamos 𝑡𝑡𝑖𝑖−1

𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � �−�𝑓𝑓𝑖𝑖−1 �(𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 )2 + (𝑓𝑓𝑖𝑖 )((𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖−1 )2 − (𝑡𝑡𝑖𝑖−1 − 𝑡𝑡𝑖𝑖 )2 )�
2ℎ para completar el trinomio
cuadrado perfecto donde sea
necesario
𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 = ℎ

𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖−1 = 2ℎ
𝑡𝑡𝑖𝑖−1 − 𝑡𝑡𝑖𝑖 = −ℎ

1 Reemplazando:
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � �−�𝑓𝑓𝑖𝑖−1 �(ℎ)2 + (𝑓𝑓𝑖𝑖 )((2ℎ)2 − (−ℎ)2 )�
2ℎ

ℎ
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � �−�𝑓𝑓𝑖𝑖−1 � + (3𝑓𝑓𝑖𝑖 )� Sacando factor común ℎ2
2
 ℎ
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � �(3𝑓𝑓𝑖𝑖 ) − �𝑓𝑓𝑖𝑖−1 �� Organizando nos queda la
2 solución:

Curva aproximada, donde

𝑦𝑦𝑖𝑖−1 y 𝑡𝑡𝑖𝑖−1 son las
condiciones iniciales y h la
distancia entre cada punto
en el eje t

Especifica 3: Método de Euler Corrector

AFIRMACIÓN JUSTIFICACIÓN

𝑓𝑓(𝑡𝑡, 𝑦𝑦(𝑡𝑡))
Polinomio que pasa por “n”
(𝑡𝑡𝑖𝑖+1 , 𝑓𝑓𝑖𝑖+1 ), (𝑡𝑡𝑖𝑖 , 𝑓𝑓𝑖𝑖 ), ,…,(𝑡𝑡𝑖𝑖−𝑛𝑛+1 , 𝑓𝑓𝑖𝑖−𝑛𝑛+1 ) puntos

Se dice que la derivada de la

𝑦𝑦 ′ = 𝑓𝑓(𝑡𝑡, 𝑦𝑦) función y es igual a otra
función f(t,y)
𝑡𝑡𝑖𝑖+1 𝑡𝑡𝑖𝑖+1

� 𝑦𝑦′(𝑡𝑡)𝑑𝑑𝑑𝑑 = � 𝑓𝑓�𝑡𝑡, 𝑦𝑦(𝑡𝑡)�𝑑𝑑𝑑𝑑

𝑡𝑡𝑖𝑖 𝑡𝑡𝑖𝑖

𝑡𝑡𝑖𝑖+1 Integrando desde t i hasta

t i+1
𝑦𝑦𝑖𝑖+1 − 𝑦𝑦𝑖𝑖 = � 𝑓𝑓�𝑡𝑡, 𝑦𝑦(𝑡𝑡)�𝑑𝑑𝑑𝑑
𝑡𝑡𝑖𝑖

𝑡𝑡𝑖𝑖+1

𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � 𝑓𝑓�𝑡𝑡, 𝑦𝑦(𝑡𝑡)�𝑑𝑑𝑑𝑑

Despejando 𝑦𝑦𝑖𝑖+1
𝑡𝑡𝑖𝑖
 Aproximando 𝑓𝑓�𝑡𝑡, 𝑦𝑦(𝑡𝑡)�
mediante el polinomio de
interpolación que pasa por
(𝑡𝑡𝑖𝑖+1 − 𝑡𝑡 )𝑓𝑓𝑖𝑖 + (𝑡𝑡 − 𝑡𝑡𝑖𝑖 )𝑓𝑓𝑖𝑖+1 los puntos (𝑡𝑡𝑖𝑖+1 , 𝑓𝑓𝑖𝑖+1 ),
𝑃𝑃(𝑡𝑡 ) = (𝑡𝑡𝑖𝑖 , 𝑓𝑓𝑖𝑖 ), donde 𝑓𝑓𝑖𝑖+1 =
ℎ 𝑓𝑓(𝑡𝑡𝑖𝑖+1 , 𝑦𝑦(𝑡𝑡𝑖𝑖+1 ))
Y 𝑓𝑓𝑖𝑖 = 𝑓𝑓(𝑡𝑡𝑖𝑖 , 𝑦𝑦(𝑡𝑡𝑖𝑖 )) el
polinomio interpolante esta
dado por:
𝑡𝑡𝑖𝑖+1

𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � 𝑃𝑃(𝑡𝑡)𝑑𝑑𝑑𝑑

𝑡𝑡𝑖𝑖

Reemplazando:
𝑡𝑡𝑖𝑖+1
(𝑡𝑡𝑖𝑖+1 − 𝑡𝑡)𝑓𝑓𝑖𝑖 + (𝑡𝑡 − 𝑡𝑡𝑖𝑖 )𝑓𝑓𝑖𝑖+1
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � � 𝑑𝑑𝑑𝑑
ℎ
𝑡𝑡𝑖𝑖

𝑡𝑡𝑖𝑖+1
(𝑡𝑡𝑖𝑖+1 )𝑓𝑓𝑖𝑖
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � � 𝑑𝑑𝑑𝑑
ℎ
𝑡𝑡𝑖𝑖
𝑡𝑡𝑖𝑖+1 𝑡𝑡𝑖𝑖+1 𝑡𝑡𝑖𝑖+1
(𝑡𝑡)𝑓𝑓𝑖𝑖 (𝑡𝑡)𝑓𝑓𝑖𝑖+1 (𝑡𝑡𝑖𝑖 )𝑓𝑓𝑖𝑖+1
−� � � 𝑑𝑑𝑑𝑑 + � � � 𝑑𝑑𝑑𝑑 − � � � 𝑑𝑑𝑑𝑑
ℎ ℎ ℎ
𝑡𝑡𝑖𝑖 𝑡𝑡𝑖𝑖 𝑡𝑡𝑖𝑖

(𝑡𝑡𝑖𝑖+1 )𝑓𝑓𝑖𝑖 𝑓𝑓 Integrando:

𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � (𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 ) − � 𝑖𝑖 � (𝑡𝑡𝑖𝑖+1 2 − 𝑡𝑡12 )
ℎ 2ℎ
𝑓𝑓 (𝑡𝑡𝑖𝑖 )𝑓𝑓𝑖𝑖+1
+ � 𝑖𝑖+1 � (𝑡𝑡𝑖𝑖+1 2 − 𝑡𝑡12 ) − � � (𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 )
2ℎ ℎ

(𝑡𝑡𝑖𝑖+1 ) 1
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + �𝑓𝑓𝑖𝑖 � �� � �𝑡𝑡 2 − 𝑡𝑡𝑖𝑖2 ��
� (𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 ) − �
ℎ 2ℎ 𝑖𝑖+1
1 (𝑡𝑡𝑖𝑖 )
+ (𝑓𝑓𝑖𝑖+1 ) �� � �𝑡𝑡𝑖𝑖+1 2 − 𝑡𝑡𝑖𝑖2 � − � � (𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 )�
2ℎ ℎ
Agrupando términos

2𝑡𝑡𝑖𝑖+1 (𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 ) − 𝑡𝑡𝑖𝑖+1 2 + 𝑡𝑡𝑖𝑖2

𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + �𝑓𝑓𝑖𝑖 � � �
2ℎ
𝑡𝑡𝑖𝑖+1 2 − 𝑡𝑡𝑖𝑖2 − 2𝑡𝑡𝑖𝑖 (𝑡𝑡𝑖𝑖+1 ) + 2𝑡𝑡𝑖𝑖2
+ (𝑓𝑓𝑖𝑖+1 ) � �
2ℎ
 ��𝑓𝑓𝑖𝑖 ��𝑡𝑡𝑖𝑖+1 2 − 2(𝑡𝑡𝑖𝑖+1 )(𝑡𝑡𝑖𝑖 ) + 𝑡𝑡2𝑖𝑖 � + (𝑓𝑓𝑖𝑖+1 )�𝑡𝑡𝑖𝑖+1 2 + 𝑡𝑡𝑖𝑖2 − 2𝑡𝑡𝑖𝑖+1 𝑡𝑡𝑖𝑖 ��
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 +
2ℎ

1
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � �(𝑓𝑓𝑖𝑖 )�𝑡𝑡𝑖𝑖+1 2 − 2(𝑡𝑡𝑖𝑖+1 )(𝑡𝑡𝑖𝑖 ) + 𝑡𝑡𝑖𝑖2 �
2ℎ
+ (𝑓𝑓𝑖𝑖+1 )�𝑡𝑡𝑖𝑖+1 2 − 2𝑡𝑡𝑖𝑖+1 𝑡𝑡𝑖𝑖 + 𝑡𝑡𝑖𝑖2 ��

1
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � ��𝑓𝑓𝑖𝑖 �(𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 )2 + (𝑓𝑓𝑖𝑖+1 )((𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 )2 )� Factoramos el trinomio
2ℎ cuadrado perfecto donde sea
necesario
𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖 = ℎ

𝑡𝑡𝑖𝑖+1 − 𝑡𝑡𝑖𝑖−1 = 2ℎ
𝑡𝑡𝑖𝑖−1 − 𝑡𝑡𝑖𝑖 = −ℎ

1 Reemplazando:
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � �𝑓𝑓𝑖𝑖 (ℎ)2 + (𝑓𝑓𝑖𝑖+1 )(ℎ)2 �
2ℎ

ℎ2
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � �𝑓𝑓𝑖𝑖+1 + 𝑓𝑓𝑖𝑖 �
2ℎ Sacando factor común ℎ2

ℎ
𝑦𝑦𝑖𝑖+1 = 𝑦𝑦𝑖𝑖 + � � �𝑓𝑓𝑖𝑖+1 + 𝑓𝑓𝑖𝑖 � Organizando nos queda la
2 solución:
 Curva aproximada, donde
𝑦𝑦𝑖𝑖+1 y 𝑡𝑡𝑖𝑖+1 son las
condiciones iniciales y h la
distancia entre cada punto
en el eje t

2.3. Casos Especiales:

Caso Especial 1: Método de Euler Mejorado

AFIRMACIÓN JUSTIFICACIÓN
𝒅𝒅𝒅𝒅
= 𝑭𝑭(𝒙𝒙, 𝒚𝒚) 𝒚𝒚(𝒙𝒙𝟎𝟎 ) = 𝒚𝒚𝟎𝟎
𝒅𝒅𝒅𝒅 Dada una EDO y su condición inicial

𝑢𝑢𝑛𝑛+1 = 𝑦𝑦𝑛𝑛+1 Planteamos una igualdad

𝑢𝑢𝑛𝑛+1 = 𝑦𝑦𝑛𝑛 + ℎ𝐹𝐹(𝑥𝑥𝑛𝑛 , 𝑦𝑦𝑛𝑛 )

Después de escoger el paso h podemos
𝑢𝑢𝑛𝑛+1 = 𝑦𝑦𝑛𝑛 + ℎ𝑘𝑘1
plantear una aproximación
𝒙𝒙𝒏𝒏+𝟏𝟏 = 𝒙𝒙𝒏𝒏 + 𝒉𝒉
𝑘𝑘2 = 𝐹𝐹(𝑥𝑥𝑛𝑛+1 , 𝑢𝑢𝑛𝑛+1 )
Tomamos una segunda estimación de la
𝑘𝑘1 + 𝑘𝑘2
𝑦𝑦𝑛𝑛+1 = 𝑦𝑦𝑛𝑛 + ℎ𝑘𝑘 𝑘𝑘 = pendiente K2 para promediar con la
2
primera.
𝑢𝑢𝑛𝑛+1 = 𝑦𝑦𝑛𝑛 ∗ ℎ ∗ 𝐹𝐹(𝑥𝑥𝑛𝑛 , 𝑦𝑦𝑛𝑛 )

𝟏𝟏 Se obtienen la ecuación que va a converger

𝒚𝒚𝒏𝒏+𝟏𝟏 = 𝒚𝒚𝒏𝒏 + 𝒉𝒉 ∗ (𝑭𝑭(𝒙𝒙𝒏𝒏 , 𝒚𝒚𝒏𝒏 ) + 𝑭𝑭(𝒙𝒙𝒏𝒏+𝟏𝟏 , 𝒖𝒖𝒏𝒏+𝟏𝟏 )
𝟐𝟐 más rápido que el método normal.
 2.4. Demostraciones de la Teoría de Errores:

Error 1: Error

AFIRMACIÓN JUSTIFICACION

𝑒𝑒𝐿𝐿 = |𝑦𝑦𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 𝑜𝑜 𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟 − 𝑦𝑦𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 | Expresado en fórmula matemática

Como se muestra en la figura el error local se

define como la variación existente entre cada
valor real con el valor aproximado

𝑛𝑛 𝑛𝑛

𝑒𝑒𝑔𝑔 = � 𝑒𝑒𝑖𝑖 = � |𝑦𝑦𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣 𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 𝑜𝑜 𝑟𝑟𝑟𝑟𝑟𝑟𝑟𝑟 − 𝑦𝑦𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣𝑣 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 | Error global expresado en fórmula matemática
𝑖𝑖=1 𝑖𝑖=1
 Se observa en la gráfica que el error global es la
suma del error local en todos los puntos.

3. Caracterización:

Caracterización 1: Ventajas / Desventajas

VENTAJAS DESVENTAJAS

Caracterización 2: Conclusiones / Recomendaciones

CONCLUSIONES RECOMENDACIONES
 5. Referencias Bibliográficas

[1] Nakamura Shoichiro; “Métodos Numéricos Aplicados con Software”; Pearson Ed.;
México; 1992.

[2] Ecuaciones Diferenciales Ordinarias, Universidad de las Fuerzas Armadas-ESPE,

(Agosto - Diciembre), 2013, Cuaderno de Materia Ing. Ximena Hidalgo.

[3] Chapra, S. C., Canale, R. P., & Del Valle Sotelo, J. C. (2008). Métodos numéricos para
ingenieros (sexta edición). México, etc: McGraw-Hill.