Facultad de Ciencias Físicas

E.A.P. Ingeniería Mecánica de fluidos

Métodos Numéricos 2
Tarea Domiciliara N°01

DOCENTE: William Wilfredo Chauca Nolasco

ALUMNO: Bazan Valles, Edson Ricardo
CODIGO: 18130110

Lima – Perú
2022
 TAREA 01 MÉTODOS NUMÉRICOS II

1. Resumen:
Se desarrolló un conjunto de fórmulas de diferencias finitas y se creó una tabla que asiste
en la construcción de la matriz de diferenciación del método de diferencias finitas.

2. Introducción:
Las ecuaciones diferenciales han sido utilizadas extensamente para modelar una gran
variedad de problemas en mecánica de sólidos y fluidos, biología, ciencia de materiales,
economía, ecología ciencias computacionales, etc. Por ejemplo, las ecuaciones de Navier-
Stokes en dinámica de fluidos, las ecuaciones de Maxwell para electromagnetismo.
Desafortunadamente, aunque las ecuaciones diferenciales pueden describir una gran cantidad
de estos problemas, solo una pequeña parte de ellos pueden ser solucionados de manera
exacta en términos de funciones elementales (polinomios, seno, coseno, logaritmos etc.) y
sus combinaciones (funciones compuestas) frecuentemente, aun cuando una ecuación
diferencial puede ser resuelta analíticamente, se requiere para esto grandes esfuerzos
utilización de teorías matemáticas avanzadas. La forma cerrada de la solución puede ser muy
complicada para ser realmente útil.
Si la solución analítica no está disponible, se desea encontrar una solución aproximada a la
ecuación diferencial. Esto se logra utilizando métodos numéricos.
Este documento muestra cómo resolver numéricamente ecuaciones diferenciales ordinarias
(EDO) y parciales (EDP) usando el método diferencias finitas (MDF).

3. Método de diferencias finitas:

El método diferencias finitas está basada en aproximaciones que permiten reemplazar
ecuaciones diferenciales por ecuaciones de diferencia.
Estas aproximaciones de diferencia finita son de forma algebraica; relacionan el valor de la
variable dependiente, en un punto dentro de la región de solución, con sus valores en algunos
puntos vecinos. En base a esto la solución por diferencia finita básicamente involucra tres
pasos:
(1) Discretización.
(2) Aproximación de la Ecuación Diferencial.
(3) Solución de las Ecuaciones de Diferencia.
3.1. Discretización (Paso 1):
Debido que los cálculos para la solución de las ecuaciones diferenciales son realizados con
la ayuda de una computadora es necesario discretizar la región de solución en un número
finito de n + 1 puntos.

Figure 1: Discretización de tipo malla uniforme.

Tal Como se muestra en la Fig.(1) para el método de diferencias finitas los n + 1 puntos están
distribuidos con una distancia uniforme h. Esta discretización es conocida como malla
uniforme.
1
𝑥𝑖 = 𝑖ℎ, 𝑖 = 0,1,2, … . , 𝑛, ℎ= (1)
𝑛
En cada punto de la malla se encuentra su respectiva aproximación de la solución de
la ecuación diferencial.
3.2. Aproximación de la Ecuación Diferencial (Paso 2):
Se aproxima a la ecuación diferencial dada por su equivalente ecuación de diferencia, lo
cual significa que en lugar de resolver una ecuación diferencial para todo el problema se
tiene una ecuación de diferencia para cada punto de la malla.

3.2.1. Deducción de fórmulas de diferencias finitas

Las fórmulas de diferencias finitas pueden ser reducidas geométricamente usando series de
Taylor la derivada en un punto está dado por

𝐹 (𝑥 + ℎ) − 𝐹(𝑥) (2)
𝐹′ (𝑥 ) = lim
ℎ→0 ℎ

Si se asume que h es pequeña, da como resultado que

 𝐹 (𝑥 + ℎ ) − 𝐹 (𝑥 ) (3)
𝐹 ′ (𝑥 ) ≈ ,
ℎ

Esta fórmula es una aproximación a la primera derivada. Es importante notar que la Ec.(3)
establece que ambas partes son “casi” iguales, lo que implica que existe un error.

𝐹 (𝑥 + ℎ ) − 𝐹 (𝑥 )
𝐹 ′ (𝑥 ) = + 𝑒𝑟𝑟𝑜𝑟. (4)
ℎ

Figure 2: Comparación entre la derivada analítica (recta punteada) y la derivada aproximada

̅̅̅̅ recta continua
por diferencia hacia adelante 𝑃𝐵

A la Ec.(4) se le conoce como diferencia hacia delante para la primera derivada.

Geométricamente, como se observa en Fig.(2), esta fórmula representa la pendiente de la

línea secante que conecta los dos puntos (𝑥, 𝐹(𝑥)) y (𝑥 + ℎ, 𝐹(𝑥 + ℎ)).

Este método para obtener las aproximaciones de diferencia finita es intuitivo, un método más
general es la utilización de la Serie de Taylor

3.2.2. Aplicación de la Serie de Taylor.

La serie de Taylor es la herramienta más importante en el análisis del método de diferencias

finitas. Esta puede escribirse de la siguiente forma:

1 1
𝐹(𝑥 + ℎ) = 𝐹 (𝑥 ) + ℎ𝐹′ (𝑥 ) + (ℎ)2 𝐹′′ (𝑥 ) + (ℎ)3 𝐹′′′ (𝑥 ) + ⋯ (5)
2! 3!

Para obtener la primera derivada se despeja 𝐹’(𝑥),

 1 1
𝐹 (𝑥 + ℎ) − [𝐹 (𝑥 ) + (ℎ)2 𝐹′′ (𝑥 ) + (ℎ)3 𝐹′′′ (𝑥) + ⋯ ]
𝐹′(𝑥) = 2! 3! (6)
ℎ

reordenando,

𝐹 (𝑥 + ℎ ) − 𝐹 (𝑥 ) 1 1
𝐹′(𝑥) = − ( (ℎ)𝐹′′ (𝑥 ) + (ℎ)2 𝐹′′′(𝑥 ) + ⋯ ) (7)
ℎ 2! 3!

Finalmente, se trunca la serie, dando como resultado una fórmula que contiene un error cuyo
orden es el del valor de ℎ𝑛 de la parte no tomada en cuenta

𝐹 (𝑥 + ℎ ) − 𝐹 (𝑥 )
𝐹′(𝑥) = + 𝑂 (ℎ ) (8)
ℎ

La Ec.(8) corresponde a la fórmula de diferencia hacia adelante para la primera derivada con
orden de error 𝑂(ℎ1 ), lo que en este caso, significa que el error de la aproximación es
proporcional al tamaño del incremento ℎ (ver la Fig.(2)). Se lee como “error de orden ℎ𝑛 ”.
Es importante recordar que ℎ es un valor pequeño (menor a 1), por lo que para 𝑛1 <
𝑛2 , ℎ𝑛1 ≫ ℎ𝑛2 , de ahí que el orden del error obedezca solo al termino más significativo.

Por otro lado, al tomar el incremento ℎ hacia atrás en la serie de Taylor se obtiene,

1 1
𝐹(𝑥 − ℎ) = 𝐹 (𝑥 ) − ℎ𝐹′ (𝑥 ) + (ℎ)2 𝐹′′ (𝑥 ) − (ℎ)3 𝐹′′′ (𝑥 ) + ⋯ (9)
2! 3!

Figura 3: Comparación entre la derivada analítica (recta punteado) y la derivada aproximada

por diferencia hacia atrás ̅̅̅̅
𝐴𝑃 (recta continua).
Despejando F’ (x),
1 1
−ℎ𝐹′ (𝑥 ) = 𝐹(𝑥 − ℎ) − [𝐹(𝑥 ) + (ℎ)2 𝐹′′ (𝑥 ) − (ℎ)3 𝐹′′′ (𝑥 ) + ⋯ ] (10)
2! 3!

Al reagrupar,

𝐹 (𝑥 ) − 𝐹 (𝑥 − ℎ ) 1 1
𝐹 ′ (𝑥 ) = + [ (ℎ)𝐹′′ (𝑥 ) − (ℎ)2 𝐹′′′(𝑥 ) + ⋯ ] (11)
ℎ 2! 3!

Finalmente, se obtiene la fórmula de diferencia hacia atrás para la primera derivada con orden
de error h,

𝐹 (𝑥 ) − 𝐹 (𝑥 − ℎ ) (12)
𝐹 ′ (𝑥 ) = + 𝑂(ℎ)
ℎ

La cual se representa en la Fig.(3).

Sumando las series de Taylor hacia adelante y hacia atrás la aproximación a la primera
derivada por diferencia centrada es

𝐹 (𝑥 + ℎ ) − 𝐹 (𝑥 − ℎ )
𝐹 ′ (𝑥 ) = + 𝑂 (ℎ 2 ). (13)
ℎ

Cuya representación se muestra en la Fig. (4).

Al comparar las fórmulas para la primera derivada de diferencia hacia delante (Ec.(8)) hacia
atrás (Ec.(12)), y centrada (Ec.(13)), se aprecia que cada una de ellas requiere la misma
cantidad de puntos de la malla para aproximar la derivada en un punto (en este caso dos
puntos), la diferencia consiste en que la fórmula de diferencia centrada tiene un error de
𝑂(ℎ2 ) mientras que la diferencia hacia adelante y hacia atrás tiene un error de 𝑂(ℎ).

𝑂 (ℎ2 ) < 𝑂(ℎ)

Esto se aprecia en la Fig. (5), donde la línea que representa la solución para la diferencia
̅̅̅̅), es más cercana al resultado real que las líneas que representan el resultado a
centrada (𝐴𝐵
̅̅̅̅ 𝑦 ̅̅̅̅
partir de las otras opciones (𝐴𝑃 𝑃𝐵 ).

De lo anterior, se puede concluir que la exactitud de la aproximación puede incrementarse

básicamente de 2 formas:
(1) Incrementando el número de puntos tomados en la construcción de la fórmula. Se
obtiene una fórmula que involucre un mayor número de puntos de la malla y con esto
incremente su exactitud. Esto se logra al desarrollar y combinar cierto número de series de
Taylor tomando en cuenta un número mayor de puntos, por ejemplo 𝐹 (𝑥 + 2ℎ) , 𝐹 (𝑥 + 3ℎ),
etc.

Figure 3: Comparación entre la derivada analítica (recta punteada) y la derivada

̅̅̅̅ (𝑟𝑒𝑐𝑡𝑎 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑎 𝑖𝑛𝑓𝑒𝑟𝑖𝑜𝑟)
aproximada por diferencia centrada 𝐴𝐵

(1)
Figure 4: Comparación entre la derivada analítica (recta punteada) y las derivadas
̅̅̅̅, 𝐴𝑃
aproximadas (rectas 𝑃𝐵 ̅̅̅̅ Y 𝐴𝐵
̅̅̅̅).

(2) Incrementando el número de puntos en la malla. Dada una fórmula para cualquiera de
las aproximaciones, entre mayor sea el número de puntos en la malla, menor será el valor del
error.

Entre mayor sea el número de puntos que utilice una ecuación para aproximar la solución, y
en los que esté discretizada la región de solución, mayor será el trabajo computacional
requerido. Por esta razón, es importante que el grado de exactitud de la solución mínimo
requerido sea conocido, para evitar así cálculos innecesarios.

2.3. Método de Coeficientes Indeterminados

Supongamos que se desea utilizar diferencia finita hacia atrás para aproximar 𝐹′(𝑥) para un
error de segundo orden. Por el método de coeficientes indeterminados utilizando
𝐹 (𝑥 ), 𝐹 (𝑥 − ℎ), 𝐹(𝑥 − 2ℎ), escribimos,

𝐹′(𝑥) ≈ 𝛾1 𝐹(𝑥 ) + 𝛾2 𝐹 (𝑥 − ℎ) + 𝛾3 𝐹 (𝑥 − 2ℎ).

Podemos utilizar la Serie de Taylor en 𝑥 para obtener un sistema de ecuaciones,

𝛾1 𝐹(𝑥 ) + 𝛾2 𝐹 (𝑥 − ℎ) + 𝛾3 𝐹 (𝑥 − 2ℎ)
ℎ2 ′′ ℎ3
= 𝛾1 𝐹 (𝑥 ) + 𝛾2 (𝐹 (𝑥 ) − ℎ𝐹′ (𝑥 ) + 𝐹 (𝑥 ) − 𝐹′′′(𝑥))
2 6
4ℎ2 ′′ 8ℎ3
+ 𝛾3 (𝐹 (𝑥 ) − 2ℎ𝐹′ (𝑥 ) + (
𝐹 𝑥 − ) 𝐹′′′(𝑥))
2 6
(14)
+ 𝑂 (𝑚á𝑥 |𝛾𝑘 |ℎ4 )

La combinación lineal debe aproximar 𝐹′ (𝑥). Por lo que escribimos

𝛾1 + 𝛾2 + 𝛾3 = 0

−ℎ𝛾2 − 2ℎ𝛾3 = 1

ℎ2 𝛾2 + 2ℎ2 𝛾3 = 0

Debido a que se está buscando 𝐹′(𝑥), el renglón correspondiente a la primera derivada se

iguala a 1.

3 2 1
Resolviendo el sistema obtenemos, 𝛾1 = − , 𝛾2 = , 𝛾3 = − . Por lo tanto, obtenemos la
2ℎ ℎ 2ℎ

fórmula para diferencia hacia atrás con orden de 𝑂(ℎ2 ),

3 2 1
𝐹 ′ (𝑥 ) = 𝐹 (𝑥 ) − 𝐹 (𝑥 − ℎ ) + 𝐹 (𝑥 − 2ℎ) + 𝑂(ℎ2 ) (15)
2ℎ ℎ 2ℎ
Escrito de otra forma,

𝐹 (𝑥 − 2ℎ) − 4𝐹 (𝑥 − ℎ) + 3𝐹 (𝑥 )
𝐹 ′ (𝑥 ) = + 𝑂(ℎ2 ) (16)
2ℎ

3.2.4. Resumen del uso de la Serie de Taylor:

- La Serie de Taylor y el Método de Diferencias Finitas asumen que la función a ser resuelta
varía de forma continua (no existen discontinuidades).

- Es posible construir fórmulas para diferenciar hacia adelante, hacia atrás, y centrada a partir
de la Serie de Taylor.

3.3. Solución de las Ecuaciones de Diferencia (Paso 3)

𝑑𝑦
Se desea resolver la ecuación diferencial = sin(𝑥 2 ) en el intervalo [0,2𝜋 ], ver Fig.(6).
𝑑𝑥

Siguiendo el Paso 1 del procedimiento, la región de solución es discretizada, el intervalo

continuo [0,2𝜋 ] se convierte en un vector {𝑥𝑖 } de {𝑛 + 1} valores (o puntos en la malla),
desde 0 hasta 2𝜋.

A su vez, la ecuación diferencial en forma continua,

𝑑𝑦
= sin (𝑥 2 )
𝑑𝑥

pasa a la forma discreta,

[𝐴𝑖𝑗 ]{𝑦𝑖 } = {𝑠𝑖𝑛 (𝑥𝑖 )2 } (17)

donde,

[𝐴𝑖𝑗 ]: Matriz de diferenciación de [{𝑛 + 1} × {𝑛 + 1}] valores, forma discreta del operador
𝑑𝑦
continuo .
𝑑𝑥

{𝑦𝑖 }: Solución de la ecuación diferencial en forma discreta, vector de {𝑛 + 1} valores.

{𝑠𝑖𝑛 (𝑥𝑖 )2 }: Información del problema, vector de {𝑛 + 1} valores. e 𝑖, 𝑗 = 1, 2, 3 … 𝑛 + 1.

El segundo paso consiste en la Aproximación de la Ecuación Diferencial, por lo tanto, se

selecciona el tipo de fórmula de diferencia finita que se utilizara para resolver el problema.
En este caso se selecciona la fórmula de diferencia hacia adelante para la primera derivada
con 𝑂(ℎ2 ),

−𝐹 (𝑥 − 2ℎ) + 4𝐹 (𝑥 − ℎ) − 3𝐹 (𝑥 )
𝐹 ′ (𝑥 ) = + 𝑂(ℎ2 )
2ℎ

esta se aplica a cada uno de los puntos de la malla, formando así un sistema de 𝑛 + 1
ecuaciones simultaneas (una ecuación por cada punto),

1
(−3𝑦1 + 4𝑦2 − 𝑦3 + 0 + 0 + ⋯ + 0) = sin(𝑥1 )2
2ℎ
1
(0 − 3𝑦2 + 4𝑦3 − 𝑦4 + 0 + ⋯ + 0) = sin(𝑥2 )2
2ℎ
1
(0 + 0 − 3𝑦3 + 4𝑦4 − 𝑦5 + ⋯ + 0) = sin(𝑥3 )2
2ℎ
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . = . . . . . . . . ..
1
(0 + ⋯ + 0 + 0 − 3𝑦𝑛−1 + 4𝑦𝑛 − 𝑦𝑛+1 ) = sin(𝑥𝑛−1 )2
2ℎ

1
(0 + ⋯ + 0 + 𝑦𝑛−2 − 4𝑦𝑛−1 + 3𝑦𝑛 + 0) = sin(𝑥𝑛 )2
2ℎ

1
(0 + ⋯ + 0 + 0 + 𝑦𝑛−1 − 4𝑦𝑛 + 3𝑦𝑛+1 ) = sin(𝑥𝑛+1 )2
2ℎ

Este sistema puede escribirse en la forma [𝐴𝑖,𝑗 ]{𝑦𝑖 } = {sin(𝑥𝑖 )2 },

𝑦1 sin(𝑥1 )2
−3 + 4 − 1 + 0 + 0 + ⋯ + 0
0−3+4−1+0+⋯+0 𝑦2 sin(𝑥2)2
0+0−3+4−1+⋯+0 𝑦3 sin(𝑥3)2
1/2ℎ ............................ … = … (18)
0+ ⋯+0 +0 −3 + 4 −1 𝑦𝑛−1 sin(𝑥𝑛−1 )2
0+ ⋯+0 +1 −4 + 3 +0 𝑦𝑛
sin(𝑥𝑛)2
[
[ 0 + ⋯ + 0 + 0 + 1 − 4 + 3 ] 𝑛+1 𝑦 ]
[sin(𝑥 )2 ] 𝑛+1

Por lo que para encontrar el vector solución {𝑦𝑖 }, se despeja la ecuación, y obtenemos,

−1
{𝑦𝑖 } = [𝐴𝑖,𝑗 ] {sin(𝑥𝑖 )2 } (19)
 sin(𝑥1 )2
𝑦1 −3 + 4 − 1 + 0 + 0 + ⋯ + 0 −1
𝑦2 0 −3 +4 −1 +0 + ⋯+0 sin(𝑥2 )2
𝑦3 0 +0 −3 +4 −1 + ⋯+0 sin(𝑥3 )2
… = 1/2ℎ ............................ … (20)
𝑦𝑛−1 0 +⋯+0 + 0− 3 +4 −1 sin(𝑥𝑛−1 )2
𝑦𝑛 0 +⋯+0 + 1− 4 +3 +0 sin(𝑥𝑛)2
[𝑦𝑛+1 ] [ 0 +⋯+0 + 0+ 1 −4 +3 ]
[sin(𝑥𝑛+1 )2 ]

La matriz de diferenciación [𝐴𝑖,𝑗 ] es una matriz singular, al aplicar las condiciones de

frontera,

𝑦 (2𝜋 ) = 0 (21)

Figura 5:Gráfica de la ecuación [A i,j]{yi } = {sin(x i )2 } en el intervalo [0, 2π] .

−1
Figure 6: Gráfica de {yi } = [A i,j] {sin(x i )2 } en el intervalo [0, 2𝜋 ].

la singularidad se elimina de la forma

−1
{𝑦𝑖 } = [𝐴𝑖,𝑗 ] {sin(𝑥𝑖 )2 } − {𝐴𝑖,𝑁+1 𝑦𝑁+1 } 𝑖 = 1 … 𝑁, 𝑗 = 1 … 𝑁 (22)

donde,
{𝑦𝑁+1 } = {0}
Una vez realizada la operación, se obtiene un vector desde 𝑦1 hasta 𝑦𝑁+1 el cual es la solución
numérica de la ecuación diferencial (ver Fig. (7)). Y su solución numérica, la Ec. (19),
3.3.1. Técnica para la construcción de Matrices de Diferenciación.

Es importante recalcar que la matriz de diferenciación [𝐴𝑖𝑗 ] es independiente de la función

que se desea resolver.

La matriz de diferenciación [𝐴𝑖𝑗 ], puede utilizarse tanto para obtener la derivada, {𝑏𝑖 }, de
una función cualquiera {𝑢𝑖 },

[𝐴𝑖𝑗 ]{𝑢𝑖 } = {𝑏𝑖 } (23)

o como para resolver una ecuación diferencial, y obtener la función {𝑢 },

{𝑢𝑖 } = [𝐴_𝑖𝑗 ]−1 {𝑏𝑖 } (24)

A continuación, se muestran un par de ejemplos de matrices de diferenciación,

a)𝑛 + 1 = 6,

−𝐹 (𝑥 − 2ℎ) + 4𝐹 (𝑥 − ℎ) − 3𝐹 (𝑥 )
𝐹′(𝑥 ) = + 𝑂(ℎ2 )
2ℎ

−3 4 −1 0 0 0
0 −3 4 −1 0 0
1 0 0 −3 4 −1 0
→𝐴=
2ℎ 0 0 0 −3 4 −1
0 0 1 −4 3 0
[ 0 0 0 1 −4 3]

b)𝑛 + 1 = 8,

2𝐹 (𝑥 ) − 5𝐹 (𝑥 − ℎ) + 4𝐹 (𝑥 − 2ℎ) − 𝐹 (𝑥 − 3ℎ)
𝐹′′ (𝑥 ) = + 𝑂 (ℎ 2 )
ℎ3

2 −5 4 −1 0 0 0 0
0 2 −5 4 −1 0 0 0
0 0 2 −5 4 −1 0 0
1 −1 4 −5 2 0 0 0 0
→𝐴= 3
ℎ 0 −1 4 −5 2 0 0 0
0 0 −1 4 −5 2 0 0
0 0 0 −1 4 −5 2 0
[ 0 0 0 0 −1 4 −5 2]

Las siguientes tablas pueden ser utilizadas como guía en la construcción de algunas matrices
de diferenciación.
 Tipo Adelante

Derivada Error 𝐹(𝑥) 𝐹(𝑥+1) 𝐹(𝑥+2) 𝐹(𝑥+3) 𝐹(𝑥+4) 𝐹(𝑥+5) 𝐹(𝑥+6) 𝐹(𝑥+7)
Primera 𝑂(ℎ) (1/ℎ) × −1 1 0 0 0 0 0 0
2 (1/2ℎ) ×
Primera 𝑂(ℎ ) −3 4 −1 0 0 0 0 0
Primera 𝑂(ℎ3 ) (1/6ℎ) × −11 18 −9 2 0 0 0 0
4 (1/12ℎ) ×
Primera 𝑂(ℎ ) −25 38 −36 16 −3 0 0 0
Segunda 𝑂(ℎ2 ) (1/ℎ2 ) × 2 −5 4 −1 0 0 0 0
3 2
Segunda 𝑂(ℎ ) (1/12ℎ ) × 35 −104 114 −56 11 0 0 0
Segunda 𝑂(ℎ4 ) (1/12ℎ2 ) × 45 −154 214 −156 61 −10 0 0
3 3
Tercera 𝑂(ℎ ) (1/4ℎ ) × −17 71 −118 98 −41 7 0 0
Tercera 𝑂(ℎ4 ) (1/8ℎ3 ) × −49 232 −461 496 −307 104 −15 0
4 4
Cuarta 𝑂(ℎ ) (1/6ℎ ) × 56 −333 952 −1219 1056 −555 164 −21

Cuando se construye el sistema de ecuaciones, y consecuentemente, al obtener la matriz de

diferenciación, en el caso de diferencia hacia adelante,

1a) Se aplica la fórmula de diferencia hacia adelante desde el primer punto, 𝑦1 , hasta el nodo
𝑦(𝑛+1)−|𝑘 𝑚á𝑥|, donde 𝑛 + 1 es igual al numero de la malla, y donde |𝑘𝑚á𝑥 | es el número de
incremento máximo que contiene la fórmula de diferencia finita que está utilizando, por
ejemplo, para el caso de diferencia finita hacia adelante con orden de error 𝑂(ℎ2 ),

−𝐹 (𝑥 − 2ℎ) + 4𝐹 (𝑥 − ℎ) − 3𝐹 (𝑥 )
𝐹′(𝑥 ) = + 𝑂(ℎ2 )
2ℎ

el valor correspondiente a |𝑘𝑚á𝑥 | seria 𝑘 = 2, tomado del termino −𝐹(𝑥 + 2ℎ).

2a) Para completar la matriz, desde el punto 𝑦(𝑛+1)−(|𝑘 𝑚á𝑥|)+1 hasta el punto 𝑦𝑛+1 se utiliza
la fórmula de diferencia hacia atrás, con el mismo orden de error.
 Tipo Adelante

Derivada Error 𝐹(𝑥−7) 𝐹(𝑥−6) 𝐹(𝑥−5) 𝐹(𝑥−4) 𝐹(𝑥−3) 𝐹(𝑥−2) 𝐹(𝑥−1) 𝐹(𝑥)
Primera 𝑂(ℎ) (1/ℎ) × 0 0 0 0 0 0 -1 1
2 (1/2ℎ) ×
Primera 𝑂(ℎ ) 0 0 0 0 0 1 -4 3
Primera 𝑂(ℎ3 ) (1/6ℎ) × 0 0 0 0 -2 9 -18 11
4 (1/12ℎ) ×
Primera 𝑂(ℎ ) 0 0 0 3 -16 36 -48 25
Segunda 𝑂(ℎ2 ) (1/ℎ2 ) × 0 0 0 0 -1 4 -5 2
2
Segunda 𝑂(ℎ3 ) (1/12ℎ ) × 0 0 0 11 -56 114 -104 35
Segunda 𝑂(ℎ4 ) (1/12ℎ2 ) × 0 0 -10 61 -156 214 -154 45
3
Tercera 𝑂(ℎ3 ) (1/4ℎ ) × 0 0 -7 41 -98 118 -71 17
Tercera 𝑂(ℎ4 ) (1/8ℎ3 ) × 0 15 -104 307 -496 461 -232 49
4
Cuarta 𝑂(ℎ4 ) (1/6ℎ ) × -21 164 -555 1056 -1219 852 -333 56

Si se desea hacer diferencia hacia atrás,

1b) Se aplica diferencia hacia adelante al inicio de la malla, desde el punto 𝑦1 hasta el punto
𝑦(|𝐾𝑚á𝑥|), y

2b) Se utiliza diferencia hacia atrás desde modo 𝑦(|𝑖𝑚á𝑥|+1) hasta el modo 𝑦𝑛+1 .

Tipo Adelante

Derivada Error 𝐹(𝑥−3) 𝐹(𝑥−2) 𝐹(𝑥−1) 𝐹(𝑥) 𝐹(𝑥+1) 𝐹(𝑥+2) 𝐹(𝑥+3)

Primera 𝑂(ℎ2 ) (1/2ℎ) × 0 0 -1 0 1 0 0
Segunda 𝑂(ℎ2 ) (1/ℎ) × 0 0 -1 2 -1 0 0
4 2
Segunda 𝑂(ℎ ) (1/12ℎ ) × 0 -1 16 -30 16 -1 0

Finalmente, para el caso de diferencia centrada, se utilizan las dos técnicas anteriores para
completar la matriz de diferenciación.

1c) Se utiliza diferencia hacia adelante desde el punto 𝑦1 hasta el punto 𝑦(|𝐾𝑚á𝑥|)

2c) Desde el modo 𝑦(|𝑖𝑚á𝑥|+1) hasta el modo 𝑦(𝑛+1)−|𝐾 𝑚á𝑥| se utiliza diferencia centrada.
3c) Finalmente desde el punto 𝑦(𝑛+1)−(|𝐾 𝑚á𝑥|)+1 hasta el punto 𝑦𝑛+1 se usa la fórmula de
diferencia hacia atrás.

Es importante recordar que, para obtener el mayor grado de aproximación posible, todas las
ecuaciones involucradas en la generación de la matriz de diferenciación deben ser del mismo
orden de error.

A continuación, se presenta un resumen de varias fórmulas de diferencias finitas para

aproximar derivadas con diferentes órdenes de error:

Aproximación de la primera derivada:

Error de Primer Orden:

Hacia adelante:

𝐹 (𝑥 + ℎ ) − 𝐹 (𝑥 )
𝐹′ (𝑥) = + 𝑂 (ℎ ) (25)
ℎ

Hacia atrás:

𝐹 (𝑥 ) − 𝐹 (𝑥 − ℎ )
𝐹′ (𝑥) = + 𝑂 (ℎ )
ℎ (26)

Error de segundo orden:

Hacia adelante:

−𝐹(𝑥 + 2ℎ) + 4𝐹 (𝑥 + ℎ) − 3𝐹(𝑥) (27)

𝐹′ (𝑥) = + 𝑂 (ℎ 2 )
2ℎ

Hacia atrás:

𝐹 (𝑥 − 2ℎ) − 4𝐹 (𝑥 − ℎ) + 3𝐹(𝑥)
𝐹′ (𝑥) = + 𝑂 (ℎ 2 ) (28)
2ℎ

Centrada:

𝐹 (𝑥 + ℎ ) − 𝐹 (𝑥 − ℎ )
𝐹′ (𝑥) = + 𝑂 (ℎ 2 ) (29)
2ℎ

Error de tercer orden:

 Hacia adelante:

2𝐹 (𝑥 + 3ℎ) − 9𝐹 (𝑥 + 2ℎ) + 18𝐹 (𝑥 + ℎ) − 11𝐹(𝑥)

𝐹′ (𝑥) = + 𝑂 (ℎ 3 ) (30)
6ℎ

Hacia atrás:

11𝐹 (𝑥 ) − 18𝐹 (𝑥 − ℎ) + 9𝐹 (𝑥 − 2ℎ) − 2𝐹(𝑥 − 3ℎ)

𝐹′ (𝑥) = + 𝑂 (ℎ 3 ) (31)
6ℎ

Error de cuarto orden:

Hacia adelante:

−3𝐹(𝑥 + 4ℎ) + 16𝐹(𝑥 + 3ℎ) − 36𝐹(𝑥 + 2ℎ) + 48𝐹(𝑥 + ℎ) − 25𝐹(𝑥)

𝐹 ′ (𝑥) = (32)
12ℎ
+ 𝑂(ℎ4 )

Hacia atrás:

25𝐹(𝑥) − 48𝐹(𝑥 − ℎ) + 36𝐹(𝑥 − 2ℎ) − 16𝐹(𝑥 − 3ℎ) + 3𝐹(𝑥 − 4ℎ) (33)

𝐹′(𝑥) = + 𝑂(ℎ4 )
12ℎ

Aproximación de la Segunda Derivada:

Error de Segundo Orden:

Hacia adelante:

−𝐹 (𝑥 + 3ℎ) + 4𝐹 (𝑥 + 2ℎ) − 5𝐹 (𝑥 + ℎ) + 2𝐹 (𝑥 )
𝐹′′(𝑥 ) = + 𝑂(ℎ2 ) (34)
ℎ2

Hacia atras:

2𝐹(𝑥) − 5𝐹(𝑥 − ℎ) + 4𝐹(𝑥 − 2ℎ) − 𝐹(𝑥 − 3ℎ)

𝐹′′(𝑥) = + 𝑂(ℎ2 ) (35)
2
ℎ

Centrada:

𝐹(𝑥 − ℎ) − 2𝐹(𝑥) + 𝐹(𝑥 + ℎ) (36)

𝐹′′(𝑥) = 2 + 𝑂(ℎ2 )
ℎ

Error de Tercer Orden:

 Hacia adelante:

11𝐹(𝑥 + 4ℎ) − 56𝐹(𝑥 + 3ℎ) + 114𝐹(𝑥 + 2ℎ) − 104𝐹(𝑥 + ℎ) + 35𝐹(𝑥)

𝐹′′(𝑥) = + 𝑂(ℎ3 ) (37)
12ℎ2

Hacia atrás:

35𝐹(𝑥) − 104𝐹(𝑥 + ℎ) + 114𝐹(𝑥 + 2ℎ) − 56𝐹(𝑥 + 3ℎ) + 11𝐹(𝑥 + 4ℎ) (38)

𝐹′′(𝑥) = + 𝑂(ℎ3 )
12ℎ2

Error de Cuarto Orden:

Hacia adelante:

−10𝐹(𝑥 + 5ℎ) + 61𝐹(𝑥 + 4ℎ) − 156𝐹(𝑥 + 3ℎ) + 214𝐹(𝑥 + 2ℎ) − 154𝐹(𝑥 + ℎ) + 45𝐹(𝑥)
𝐹′′(𝑥) =
12ℎ2
+ 𝑂(ℎ4 )

(39)
Hacia atras:

45𝐹 (𝑥 ) − 154𝐹 (𝑥 − ℎ) + 214𝐹 (𝑥 − 2ℎ) − 156𝐹 (𝑥 − 3ℎ) + 61𝐹 (𝑥 − 4ℎ) − 10𝐹 (𝑥 − 5ℎ)
𝐹′′(𝑥 ) =
12ℎ2
+ 𝑂(ℎ4 ) (40)

Centrada:

−𝐹 (𝑥 − 2ℎ) + 16𝐹 (𝑥 − ℎ) − 30𝐹 (𝑥 ) + 16𝐹 (𝑥 + ℎ) − 𝐹(𝑥 + 2ℎ) (41)

𝐹′′(𝑥 ) = + 𝑂(ℎ4 )
12ℎ2

Aproximación de la Tercera Derivada:

Error de Tercer Orden:

Hacia adelante:

7𝐹(𝑥 + 5ℎ) − 41𝐹(𝑥 + 4ℎ) + 98𝐹(𝑥 + 3ℎ) − 118𝐹(𝑥 + 2ℎ) + 71𝐹(𝑥 + ℎ) − 17𝐹(𝑥)
𝐹 (3) (𝑥) = + 𝑂(ℎ3 )
4ℎ3

(42)
Hacia atrás:

17𝐹(𝑥) − 71𝐹(𝑥 − ℎ) + 118𝐹(𝑥 − 2ℎ) − 98𝐹(𝑥 − 3ℎ) + 41𝐹(𝑥 − 4ℎ) − 7𝐹(𝑥 − 5ℎ)
𝐹 (3) (𝑥) = + 𝑂(ℎ3 )
4ℎ3
(43)
Error de Cuarto Orden:

Hacia adelante:

1
𝐹 (3) (𝑥 ) = ( ) [−15𝐹 (𝑥 + 6ℎ) + 104𝐹 (𝑥 + 5ℎ) − 307𝐹 (𝑥 + 4ℎ) + 496𝐹 (𝑥 + 3ℎ)
8ℎ3
− 461𝐹 (𝑥 + 2ℎ) + 232𝐹 (𝑥 + ℎ) − 49𝐹 (𝑥 )] + 𝑂(ℎ4 )
(44)

Hacia atras:

1
𝐹 (3) (𝑥 ) = ( 3 ) [49𝐹 (𝑥 ) − 232𝐹 (𝑥 − ℎ) + 461𝐹 (𝑥 − 2ℎ) − 496𝐹 (𝑥 − 3ℎ)
8ℎ
+ 307𝐹 (𝑥 − 4ℎ) − 104𝐹 (𝑥 − 5ℎ) + 15𝐹 (𝑥 − 6ℎ)] + 𝑂(ℎ4 )

(45)
c
Aproximación de la Cuarta Derivada:

Error de Cuarto Orden:

Hacia adelante:

1
𝐹 (4) (𝑥 ) = ( 4 ) [−21𝐹 (𝑥 + 7ℎ) + 164𝐹 (𝑥 + 6ℎ) − 555𝐹 (𝑥 + 5ℎ)
6ℎ
+ 1056𝐹 (𝑥 + 4ℎ) − 1219𝐹 (𝑥 + 3ℎ) + 852𝐹 (𝑥 + 2ℎ)
− 333𝐹 (𝑥 + ℎ) + 56𝐹 (𝑥 )] + 𝑂(ℎ4 )
(46)
Hacia atras: c

1
𝐹 (4) (𝑥 ) = ( ) [56𝐹 (𝑥 ) − 333𝐹 (𝑥 − ℎ) + 852𝐹 (𝑥 − 2ℎ) − 1219𝐹 (𝑥 − 3ℎ)
6ℎ3
+ 1056𝐹 (𝑥 − 4ℎ) − 555𝐹 (𝑥 − 5ℎ) + 164𝐹 (𝑥 − 6ℎ) − 21𝐹(𝑥
− 7ℎ)] + 𝑂(ℎ4 )
(47)
c
1
𝐹 (4) (𝑥 ) = ( ) [56𝐹 (𝑥 ) − 333𝐹 (𝑥 − ℎ) + 852𝐹 (𝑥 − 2ℎ) − 1219𝐹 (𝑥 − 3ℎ)
6ℎ3
+ 1056𝐹 (𝑥 − 4ℎ) − 555𝐹 (𝑥 − 5ℎ) + 164𝐹 (𝑥 − 6ℎ) − 21𝐹(𝑥
− 7ℎ)] + 𝑂(ℎ4 )
(48)
c