Unidad 6.1

Cátedra de Métodos Numéricos
CÁTEDRA DE MÉTODOS NUMÉRICOS

DEPARTAMENTO DE COMPUTACIÓN
DERIVACIÓN E INTEGRACIÓN NUMÉRICA
DICIEMBRE 2004
(Revisado Octubre 2009)
Ing. Héctor Bridera

Ing. Beatriz Pedrotti
ÍNDICE DE CONTENIDO
Introducción ............................................................................................................ 1
Derivación Numérica .............................................................................................. 1
Introducción .......................................................................................................... 1
Obtención de Derivada Primera ............................................................................. 3
Obtención de Derivada Segunda ............................................................................ 3
Generalización de Conceptos................................................................................. 4
La Tabla de Derivación Numérica ......................................................................... 5
Derivadas de Orden Superior................................................................................. 7
Extrapolación de Richardson para Derivadas Numéricas .................................. 10
Integración Numérica ........................................................................................... 13
Introducción ........................................................................................................ 13
La Tabla de Integración Numérica....................................................................... 17
Extrapolación de Richardson para Integrales Numéricas .................................. 20
ANEXO A: TABLA DE DERIVACIÓN NUMÉRICA ....................................... 24
ANEXO B: TABLA DE INTEGRACIÓN NUMÉRICA .................................... 25
Introducción
Este capítulo versa sobre el cálculo de las derivadas e integrales de funciones f(x) en
forma numérica.
La función f(x) es una función del tipo analítica, es decir conocida por una fórmula que
la satisface o una función discreta, conocida solo por un conjunto de pares de puntos de
la misma.
Si la función f(x) es del tipo analítica, deberá ser discretizada para aplicar la derivación
y la integración numérica.
Ampliando el concepto de las funciones discretas, podemos decir que:
Una función discreta esta formada por un conjunto finito de argumentos xk, cada uno de
los cuales posee su correspondiente yk.
La distribución de argumentos xk puede ser igualmente o uniformemente espaciada o no

uniformemente espaciado.
En el caso de argumentos uniformemente espaciados, se cumple con la siguiente

relación:
h= xk-1 - xk
Las funciones discretas se utilizan:
• cuando no se conoce la función analítica, por ejemplo una

medición de laboratorio.
• cuando la función analítica es muy compleja.
• cuando se desea desarrollar un algoritmo computacional para
resolver operaciones, por ejemplo derivación, integración,
interpolación, aproximaciones polinómicas.
Discretizar un función significa expresarla a través de un conjunto finito de puntos.
Derivación Numérica
Introducción
Las derivadas numéricas puede obtenerse a través de diversos métodos, a saber:
1
• por Expansión en Serie de Taylor, cuya ventaja es que el error se obtiene en

forma explícita y se puede aplicar a argumentos uniformemente distribuidos.
• por Diferencias Finitas u Operadores de Diferencias, la desventaja es que para

analizar el error que genera el método se debe expandir en Serie de Taylor.
• por Polinomios de Interpolación, derivando los polinomios, su desventaja es la

dificultad de aplicación con argumentos no uniformemente espaciados.
Se opta entonces por utilizar el método de derivación numérica a través del Polinomio
de Taylor.
Conociendo la función en un punto y sus n-derivadas se puede conocer en forma

aproximada el valor de la función en otro punto.
Si xi+m=xi+ m h
fi+m=f( xi+ m h)
expandiendo en Serie de Taylor a la función f en el punto x i+m,
fi ± m = fi ± m h f’i + (mh)2 f’’i ± (mh)3 f’’’i + (mh)4 fIVi ± ...

2! 3! 4!
Haciendo m= ± 1 y m= ± 2
a) fi+1 = fi + h f’i + (h)2 f’’i +(h)3 f’’’i + (h)4 fIVi + ...

2 6 24
b) fi+2 = fi + 2 h f’i + 2 (h)2 f’’i +4(h)3 f’’’i + 2(h)4 fIVi + ...

3 3
c) fi-1 = fi - h f’i + (h)2 f’’i - (h)3 f’’’i + (h)4 fIVi + ...

2 6 24
d) fi-2 = fi - 2 h f’i + 2 (h)2 f’’i - 4(h)3 f’’’i + 2(h)4 fIVi + ...

3 3
Combinando linealmente las ecuaciones y truncando la serie hasta un orden de error

preestablecido se obtienen fórmulas para distintos ordenes de derivación y dentro de un
mismo orden, fórmulas de diferentes tipos, por ejemplo:
2
Obtención de Derivada Primera
De la ecuación a) se despeja la derivada primera, obteniéndose
f’i= (fi+1 – fi) 1/h ei= - (h/2)fi’’ + (h2 /6)f i’’’ +...
Esta fórmula es del tipo Lateral Inicial o Derecha con un orden de error h.
El orden de error es el paso elevado a la potencia correspondiente al primer término de

la serie de Taylor que se desprecia, siendo éste además el término más significativo del
error ya que
Si h → 0, entonces el término de error ei → 0

2
Si h < 1, entonces el término h << h, entonces ei ≈ -(h/2) fi’’
De la ecuación c) se despeja también la derivada primera obteniéndose en este caso una
fórmula del tipo Lateral Final o Izquierda con orden de error h, a saber:
f’i= (fi - fi-1 ) 1/h ei ≈ (h/2)fi’’
Combinando las ecuaciones a) y c) a través de la operación resta se obtiene una fórmula

2
del tipo Central con un orden de error h .
De a) – c)
(fi+1 - fi-1 ) = 2h fi’ + (2h3/6) fi’’’
fi’ = (fi+1 - fi-1 ) ei ≈ - (h2/6) fi’’’

2h
Obtención de Derivada Segunda
Combinando las ecuaciones a) y c) a través de la operación suma se obtiene una fórmula

2
del tipo Central con un orden de error h .
De a) + c)
fi’’= (fi+1 -2 fi + fi-1 ) ei ≈ (h2/12) fiIV

h2
Sumando las ecuaciones a, b, c y d se obtiene una fórmula del tipo Central con orden
2
de error h .
3
De a) + b) + c) + d)
fi’’= (-fi-2 - fi-1 + 4 fi - fi+1 - fi+2)

(5h2)
ei ≈ (2/24 + 4/3) h4 fiIV = -17/60 h2 fiIV

(5h2)
Generalización de Conceptos
Combinando linealmente el conjunto de cuatro ecuaciones planteadas anteriormente por

desarrollo en serie de Taylor se pueden obtener derivadas de mayor orden o nuevas
fórmulas de las derivadas primera o segunda.
Por ejemplo la siguiente combinación lineal de ecuaciones, da una derivada tercera:
-2 a +b +2 c - d
Un concepto general es que los tipos de fórmulas de derivadas numéricas pueden ser
del tipo Central, Lateral Izquierda o Inicial y Lateral Derecha o Final.
El orden de error se puede calcular con la siguiente fórmula:
sea m el número de puntos y n el orden de derivada
• En derivadas laterales, el exponente del h del término de error

es : m-n
• En derivadas centrales, el exponente del h del término de error

es : m-n+1
Ejemplo:
Sea f(x)= x3-2x2+2
Evaluar la función en forma discreta (discretizar) entre x=0 y x=1 con h=0.1 y calcular
la derivada 2º con la fórmula obtenida anteriormente:
fi’’ = fi-1 – 2 fi + fi+1

h2
Comparar con los valores exactos de la derivada 2º.
La derivada 2º de la función es:
4
derivando una vez f’(x)= 3x2 –4x
derivando nuevamente f’’(x)= 6x-4
Se realiza la resolución utilizando la siguiente tabla:
x 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

f(x) 2 1.981 1.928 1.847 1.744 1.625 1.496 1.363 1.232 1.109 1
f''(x) exacta -4 -3.40 -2.80 -2.20 -1.60 -1.00 -0.40 0.20 0.80 1.40 2
der num 2º NA -3.40 -2.80 -2.20 -1.60 -1.00 -0.40 0.20 0.80 1.40 NA
(NA: no aplicable)
2 IV
El orden de error de la fórmula es (h /12) fi , como la derivada cuarta es nula, el
error de esta fórmula aplicada a esta función es cero.
La Tabla de Derivación Numérica
La tabla de derivación numérica que se utilizará (ver tabla completa en Anexo A)

contiene todos los factores necesarios para poder calcular derivadas de orden 1º al 4º,
cada una con fórmulas de distinto tipo (centrales o laterales derecha) con diversos
ordenes de error.
La tabla se origina realizando, al igual que en la introducción, combinaciones lineales

entre ecuaciones representativas de las funciones expresadas como desarrollos en serie
de Taylor.
En forma esquemática se puede representar de la siguiente forma:
Derivada x factor Factor y Coeficientes de Error

(no olvidar traspasar cada fi+/-1 aproximado
al segundo miembro (el coeficiente marcado con (Primer término
dividiendo) círculo corresponde a fi) de Taylor que se
trunca)
hn f(n) Factor Coeficientes de la Error en la derivada f(n) Nº de

fórmula formula
hn f(n) 1 -1 1 -(1/2) h f’’ 1
h2f’’
h3f’’’ 1 -1 3 -3 1 -(3/2) h fIV 16

1/2 -5 18 -24 14 -3 (7/4) h2 fV 17
h4fIV
5
La interpretación de la tabla conlleva a los siguientes conceptos:
• La suma de todos lo coeficientes es nula ya que si la función f(x) es una

constante, todas sus derivadas deben ser nulas.
• Los tipos de fórmulas que se presentan son iniciales o inclinadas a la derecha

y centradas.
• Se pueden deducir las fórmulas finales o inclinadas a la izquierda realizando

los siguientes cambios en las fórmulas laterales derecha:
o invertir el signo de m en fi +/- mh

o invertir el signo del coeficiente de cada fi +/-m h
Ejemplo:
Transformar la fórmula Nº1 de la Tabla de derivación numérica del tipo inclinada a

derecha al tipo opuesto, inclinado a izquierda.
De la tabla, la fórmula Nº1 es fi’= (- fi + fi+1)/ h
Para obtener la fórmula lateral izquierda se cambian los signos de los coeficientes y los
de las +/- i, como sigue:
fi’= (+ fi - fi-1)/ h
Reordenando los términos, la fórmula es: fi’= (- fi-1+ fi)/ h
Ejemplo:
Transformar a lateral izquierda las fórmulas Nº 17 y 18 de la tabla de derivadas

numéricas.
La fórmula Nº 17 es para calcular una derivada tercera con 5 puntos y es
fi’’’= (-5 fi + 18 fi+1 –24 fi+2 + 14 fi+3 – 3 fi+4 ) / (2 h3)

La fórmula lateral izquierda es
fi’’’= (3 fi-4 - 14 fi-3+24 fi-2 - 18 fi-1 + 5 fi) / (2 h3)
La fórmula Nº 18 también es para calcular una derivada tercera pero con 4 puntos y su
ecuación es
fi’’’= (- fi-1 + 3 fi– 3fi+1 + fi+2) / h3

6
Transformada en lateral derecha la fórmula es

fi’’’= (- fi-2 + 3 fi-1– 3fi+ fi+1) / h3
Derivadas de Orden Superior
Las derivadas de orden superior se obtienen aplicando sucesivamente dos operadores de

orden inferior.
Es decir si por ejemplo se desea calcular una derivada de orden 6º se puede aplicar una
derivada 4º a una 2º o una 2º a una 4º o tres veces una 2º.
Ejemplo:
Desarrollar la fórmula numérica para derivada 6º aplicando una derivada 2º a una

derivada 4º.
Se eligen las fórmulas Nº23 y Nº 12 de la tabla:
fiIV= (fi-2 - 4 fi-1 + 6 fi – 4 fi+1 + fi+2) / h4 orden de error h2
fi’’= (fi-1 - 2 fi + fi+1) / h2 orden de error h2
Si se llama gi = fiIV entonces gi’’ = fiVI
Desarrollando la derivada 2º de gi según la fórmula elegida
gi’’= (g i-1 – 2 gi + gi+1) / h2
Reemplazando cada gi por las derivadas 4º, reordenando y simplificando términos se

llegará a la fórmula buscada.
Se puede evitar un desarrollo engorroso y largo realizando lo mismo en una tabla, como
sigue:
h4 factor fi-3 fi-2 fi-1 fi fi+1 fi+2 fi+3

2
h
gi-1= fi-1IV 1 1 -4 6 -4 1
gi-1= fiIV -2 1 -4 6 -4 1
gi-1= fi+1IV 1 1 -4 6 -4 1
h2 h4 1 -6 15 -20 15 -6 1
La fórmula buscada es
7
gi’’= fiIV = (fi-3 – 6 fi-2 +15 fi-1 –20 fi + 15 fi+1 –6 fi+2 + fi+3) / (h2 h4)
El exponente orden de error calculado por la fórmula m-n+1 ( para fórmula central)
si m (orden de derivación) = 6
si n (número de puntos )= 7
entonces m (exponente del orden de error) = 7 – 6 +1 = 2
Se ha mantenido el orden de error de las fórmulas generadoras.
Ejemplo:
Dada la siguiente tabla de valores de una función f(x) discreta, evaluar lo solicitado
i 0 1 2 3 4 5 6
x 0 0.1 0.2 0.3 0.4 0.5 0.6
f(x) 1 0.90484 0.81873 0.74082 0.67032 0.60653 0.54881
a) Encontrar f’(0) aplicando las fórmulas de 2, 3 y 4 puntos de la tabla, es decir:
i i+1 i+2 i+3 factor orden error

-1 1 1/h 1 fórmula 1
-3 4 -1 1/2h 2 fórmula 2
-11 18 -9 2 1/6h 3 fórmula 3
según fórmula 1 f'(0) es = ( -1 + 0.9048 + 0 + 0 ) *1/ 0.1 = -0.9516
según fórmula 2 f'(0) es = ( -3 + 3.6194 + -0.8187 + 0 ) *1/ 0.2 = -0.9969
según fórmula 3 f'(0) es = ( -11 + 16.287 + -7.3686 + 1.5 ) *1/ 0.6 = -0.9997
b) Encontrar f’’’(0) aplicando las fórmulas de 5 puntos de la tabla, es decir:
i i+1 i+2 i+3 i+4 factor orden error

-5 18 -24 14 -3 1/2(h*h*h) 2 fórmula 17
No son aplicables las formulas 19 y 20 que son también de 5 puntos.
según fórmula 17 f'''(0) es = ( -5 + 16.287 + -19.65 + 10 + -2 ) *1/ 0.002 = -0.9400

8
iv
c) Encontrar f (0) aplicando las fórmulas de 5 puntos de la tabla, es decir:
i i+1 i+2 i+3 i+4 factor orden error

1 -4 6 -4 1 1/(h*h*h*h)) 1 fórmula 21
No son aplicables las fórmulas 22 y 23.
según fórmula 21 der4(0) es = ( 1 + -3.6194 + 4.9124 + -3 + 1 ) *1/ 0.00010 = 0.6000
Ejemplo:
Dada la siguiente función discreta,
i 0 1 2 3 4 5 6
x 0.8 0.9 1 1.1 1.2 1.3 1.4
f(x) 0.71736 0.78333 0.84147 0.89121 0.93204 0.96356 0.98545
a) Encontrar f’’(1.0) o sea i=2, aplicando las fórmulas centrales de 3 y 5 puntos, es

decir:
i-2 i-1 i i+1 i+2 factor orden error

1 -2 1 1/(h*h) 2 fórmula 12
1 -2 1 0 1/(h*h) 2 fórmula13
-1 16 -30 16 -1 1/12(h*) 4 fórmula 15
La fórmula 14 no es central.
según fórmula 12 f''(1.0) es = ( 0 + 0,7833 + -1,6829 + 0,89121 + 0 ) *1/ 0,010 = -0,8400
según fórmula13 f''(1.0) es = ( 0 + 0,7833 + -1,6829 + 0,89121 + 0 ) *1/ 0,010 = -0,8400
según fórmula 15 f''(1.0) es = ( -1 + 12,533 + -25,244 + 14,2594 + -0,9320 ) *1/ 0,120 = -0,8405
b) Encontrar fIII(1.4) o sea i=6, aplicando la fórmula 17.
Para aplicar la fórmula 17 que es lateral derecha o inicial en el punto 1.4 y según los
datos que poseo, debo transformarla en fórmula lateral izquierda o final
9
formula 17 lateral derecha o inicial

i i+1 i+2 i+3 i+4
-5 18 -24 14 -3
cambio signo de indices y signo de factores
i i-1 i-2 i-3 i-4
5 -18 24 -14 3
reordeno y obtengo
formula 17 lateral izquierda o final
i-4 i-3 i-2 i-1 i factor orden error
3 -14 24 -18 5 1/2(h*h*h) 2
f'''(1.4) es = ( 2.524 + -12.477 + 22.369 + -17.344 + 4.927 ) *1/ 0.002 = -0.2000
Extrapolación de Richardson para Derivadas Numéricas
Una de las ventajas de utilizar las derivadas numéricas a partir del desarrollo en serie de
Taylor es que se conoce el orden del error que se comete al aplicar el método.
Conociendo las variaciones del error se pueden realizar mejores estimaciones, en este
concepto se basa la extrapolación de Richardson, la cual permite mejorar la precisión
obtenida por dos cálculos aproximados.
Para demostrar esto se expresa una función derivación como la aproximación obtenida
por serie de Taylor mas el error
Q = Q1 + e1
n
para un h suficientemente pequeño el error es proporcional a h a través de una
constante C
Q = Q1 + C h1n ∴ C = (Q - Q1) / h1n (1)
Q = Q2 + C h2n ∴ C = (Q – Q2) / h2n (2)
igualando (1) con (2)
(Q - Q1) = (Q – Q 2)
h1n h2n
10
Q - Q1 = (Q – Q 2) (h1 / h2)n
Q – Q (h1 / h2)n = Q1 – Q 2 (h1 / h2)n
Q [ 1– (h1 / h2)n ] = Q 1 – Q2 (h1 / h2)n
Q= Q1 – Q2 (h1 / h2)n
[ 1– (h1 / h2)n ]
o bien
Q= Q 2 (h1 / h2)n- Q 1
[ (h1 / h2)n-1 ]
Ejemplo:
Dada la función discretizada f(x) calcular la derivada primera en x=0.2 extrapolando el

cálculo por Richardson.
x 0 0.1 0.2 0.3 0.4 0.5

f(x) 1 1.105171 1.221403 1.349859 1.491825 1.648721
Para extrapolar se necesitan realizar dos cálculos con diferente paso h, se puede utilizar
h=0.1 y h=0.2 y la fórmula Nº5 de la tabla de derivación de acuerdo a la información
disponible.
i-1 i i+1 factor
Fórmula Nº 5
-1 0 1 1/2h
Error ≈ (-1/6) h2 f’’’ y orden de error es h2
a) La derivada en x=2 con una paso h1= 0.2 es
f1’(0.2) = (-1 + 1.491825) / ( 2 * 0.2) = 1.2295625
b) La derivada en x=2 con una paso h2= 0.1 es
f2’(0.2) = (-1.105171 +1.349859) / (2 * 0.1) = 1.223440
11
c) La derivada en x=2 mejorada por Richardson es
Q= Q2 (h1 / h2)n- Q1 = 22 * 1.223440 – 1.2295625 =

[ (h1 / h2)n-1 ] 22 - 1
= 1.2213992
Ejemplo:
Calcular la derivada primera de f(x)=seno(x)/x en los puntos x=1.2 y x=1.3 aplicando

la formula “8” con pasos de cálculo h1=0.1 y h2=0.2 y Extrapolando por Richardson.
La fórmula “8” es:
i-2 i-1 i i+1 i+2 factor

1 -8 0 8 -1 1/12h
Los datos de la función discretizada que necesito para el cálculo son:
x 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7

f(x) 0.896695 0.870363 0.841471 0.810189 0.776699 0.741199 0.703893 0.664997 0.624734 0.583332
a) La derivada en x=1.2 con h1=0.1 es:

(0.841470985 -6.481508073 + 0 + 5.929588833 - 0.703892664)/ (12* 0.1)= -0.345284099
La derivada en x=1.2 con h2=0.2 es:

(0.896695114 - 6.731767878 + 0 + 5.631141314 - 0.624733502)/(12*0.2)= -0.345277064
Extrapolando,
h2/h1= 2; orden de error= 4
(h2/h1)4= 16
f’(1.2)ex=(16 *( -0.345284099) – (-0.345277064))/(16-1)= -0.34452845691
b) La derivada en x=1.3 con h1=0.1 es:
(0.810188509 - 6.213593906 + 0 + 5.631141314 - 0.6649967)/(12*0.1)= -0.364383951
12
La derivada en x=1.1 con h2=0.2 es:
(0.870363233 - 6.481508073 +0+ +5.319973262 - 0.583332241)/(12*0.2)= -0.36437659
Extrapolando,
h2/h1= 2; orden de error = 4
(h2/h1)4= 16
f’(1.3)ex=(16 *( -0.364383951) – (-0.36437659))/(16-1)= -0.3643844452
Integración Numérica
Introducción
Para desarrollar las fórmulas de integración numérica se utilizará también el desarrollo

en serie de Taylor.
Una integral definida entre un valor xi y xi+nh se representa por la siguiente fórmula
xi +nh
I= ∫ f ( x)dx
xi
x z xi+ n
13
Realizando un cambio de variables y considerando un eje Z que marque el inicio del

intervalo de integración, un x genérico es
x= xi + z entonces
dx=dz
nh
I= ∫ f (xi + z)dz
0
desarrollando por Taylor la función f(xi + z)
= f(xi) + z f’(x i) + z2 f’’(x i ) + z3 f’’’(xi ) + z4 fIV(xi ) + ...=

2! 3! 4!
entonces la integral se calcula
nh nh nh
I= ∫
0
f(x i) dz + ∫
0
z f’(xi) dz + ∫
0
(z2 f’’(xi )/2) dz +
nh nh
+ ∫
0
(z3 f’’’(xi )/ 6)dz + ∫
0
4 IV
(z f (xi )/ 24) dz
se debe tener cuidado ya que se está integrando respecto a z por lo que las f(x) son
constantes para este caso, operando
I= z f(xi) + (z2/2) f’(xi) + (z3/6) f’’(xi) +(zIV/24) f’’’ (x i) +

+ (z5/120) fIV (x i) ] nh 0
la expresión anterior valuada en z=0 y z= nh da la siguiente fórmula para integrales

numéricas
I = nh fi + ((nh)2/2) f’i + ((nh)3/6) f’’i + ((nh)4/24) f’’’i +

+ ((nh)5/120) fiIV+ ...
14
Esta expresión es una expresión generadora de integrales numéricas ya que

reemplazando las derivadas por diversas fórmulas de derivadas numéricas se obtienen
también diversas fórmulas de integración, existiendo una variedad de posibilidades.
El error que se produce en esta metodología depende del primer término de Taylor que
se trunca y del error de las fórmulas de derivación numérica utilizadas.
Ejemplo:
Se toma la fórmula generadora y se reemplaza solo derivada primera utilizando la

fórmula Nº1 de la tabla derivación numérica incluido su término de error
fi’= (- fi + fi+1)/ h - h fi´´ /2 (fórmula Nº1)
I = nh fi + (n2h2/2) ((- fi + fi+1)/ h - h fi´´ /2 ) – n3 h3fi´´/6
Reagrupando lo que se puede evaluar y lo que no se puede evaluar se obtiene la

siguiente fórmula de integración
I= h/2 ( fi (2n - n2) + n2 fi+1) – n2 h3 fi´´/4 + n3 h3 fi´´´/6
I ≈ h/2 ( fi (2n - n2) + n2 fi+1) ( integral aproximada)
e ≈ n2 h3 fi´´(3-2n)/12 ( error aproximado)
Continuando con el análisis y aplicación de la fórmula generadora, si n es el número de

franjas y m es la cantidad de puntos a evaluar para calcular la fórmula, por ejemplo
dando a n el valor 1, se calcula la integral en una franja y su fórmula es
I = h ( fi + fi+1) / 2
ésta representa el área de un trapecio y su error es
e= - h3 fi´´ / 12
En la tabla de integrales numéricas esta fórmula se encontrará como tipo de fórmula I12.
En general las fórmulas de la tabla se dicen que son de tipo Inm, donde n es la cantidad
de franjas de integración que resuelve la fórmula y m es la cantidad de puntos a
evaluar.
No existe una relación directa entre n y m ya que depende de la fórmula de derivación

que se utilice en la fórmula generadora y si hay términos que se reducen.
15
Una interpretación gráfica de la fórmula anterior nos muestra que la integral se calcula
como el área de un trapecio definido entre los puntos xi y x i+1.
fi+1
fi
xi xi+1
Si en la fórmula generadora se toma ahora n=2 (dos franjas) se observa que se obtiene
una fórmula del tipo 2 franjas y 1 punto pivote, como sigue:
I21= 2 h fi+1
Realizando el análisis gráfico de la misma esta representa el área de un rectángulo de

altura constante fi+1.
fi+1
fi
xi xi+1
Hay que tener siempre en cuenta las siguientes consideraciones:
• la derivada se calcula en un punto y la integral en un intervalo

• el intervalo de integración puede estar compuesto por una o mas
franjas.
• no hay una relación directa entre el número de franjas y la cantidad de
puntos pivotes.
• el intervalo de integración siempre comienza en xi
16
La Tabla de Integración Numérica
La tabla de integración numérica que se utilizará (ver tabla completa en Anexo B)

contiene todos los factores necesarios para poder calcular integrales numéricas , con
distintos tipos de fórmulas que se diferencian por la cantidad de franjas, cantidad de
puntos pivotes y orden de error.
La tabla se origina reemplazando en la fórmula generadora diferentes fórmulas de
derivación numérica.
En forma esquemática se puede representar de la siguiente forma:

Factor y Coeficientes de Error aproximado
cada fi+/-1 (depende del
Tipo de Fórmula
(los coeficientes con truncamiento en la
(n indica la cantidad de
circulo marcan inicio y fin formula generadora y
franjas y m la cantidad de
del intervalo de de la derrivada
puntos pivotes)
integración) numérica aplicada)
Inm(i) Factor Coeficientes de la fórmula Error en la Nº de formula

fi-1 fi fi+1 fi+5 derivada f(n)
I11(i) h 1 0 (h2/2)f´ 1
.......... ...........
I14(i) h/24 -1 13 13 -1 (11/720)h5fiv 6
Hay que tener cuidado con las fórmulas de la tabla de integración utilizada ya que al
tener esta un fin didáctico hay fórmulas en ella que no son muy precisas.
Las fórmulas se obtienen armándolas con los factores y coeficientes presentados en ella,
por ejemplo la fórmula I12, vista ya anteriormente, es la fórmula Nº 2 de la tabla y
calcula la integral en una franja utilizando dos puntos pivotes, su ecuación es
I12 = h (fi + fi+1) / 2 y un término de error de h3 f’’/12
La fórmula Nº9 de la tabla es del tipo I21 y calcula la integral en dos franjas con 1 punto
pivote, observar aquí que solo se cuentas como puntos pivotes los que tienen
coeficientes distinto de cero, la ecuación en
3
I21 = 2 h fi+1 y un término de error de h f’’/3
También se encontrarán fórmulas del tipo de demasiados puntos pivotes respecto a las
franjas calculadas como la fórmula Nº3 del tipo I13 la cual calcula solo una franja con 3
puntos pivotes, su gráfica es
17
fi+ 2
fi+1
fi
xi xi+ h xi+ 2 h
y su fórmula es
I13= h( 5 fi + 8 fi+1 – fi+2 ) / 12 y su término de error de h4 f’’’/24
Si se compara los términos de error, la que tiene menor error es la fórmula Nº 3 y es

correcto ya que se trabaja con muchos mas puntos que los de la franja de integración o
sea se trabaja con información de otros puntos fuera del intervalo de integración.
La fórmula Nº 9 es la que tiene el mayor error entre las tres fórmulas analizadas y es
correcta la apreciación dado que calcula la integral en forma muy aproximada solo con
un punto de ordenada.
Analizando las tres fórmulas como el cálculo del área de un rectángulo hipotético
podemos inferir que:
En la fórmula I12 = h (fi + fi+1) / 2 la altura del rectángulo es la media entre las
funciones en los dos puntos pivotes.
En la fórmula I21 = 2 h fi+1 la altura del rectángulo es la ordenada media.
En la fórmula I13= h( 5 fi + 8 fi+1 – fi+2 ) / 12 la altura es un promedio ponderado

entre las tres ordenadas de los puntos pivotes.
La fórmula Nº6 de la tabla calcula 1 franja con 4 puntos pivotes, la fórmula Nº 15

calcula 2 franjas con 3 puntos pivotes y la fórmula Nº19 calcula 5 franjas con 2 puntos
pivotes, siendo esta la meno precisa de las tres.
Ejemplo:
18
x
Calcular la integral entre x= 0.1 y x= 0.3 de la función f(x) = e – 1 utilizando la
fórmula Nº4 con un paso de cálculo h=0.1.
Se debe tener en cuenta que no todo problema se puede resolver con cualquier fórmula
o con cualquier paso de cálculo.
La fórmula Nº4 es I14 = (9 fi + 19 fi+1 – 5fi+2 + fi+3 ) h/24
El tipo de fórmula no permite obtener la integral con una sola aplicación de la misma
sino que se deben realizar dos cálculos como sigue:
0 .3 0 .2 0 .3
I= ∫ f (x) dx = ∫ f (x) dx + ∫ f (x) dx =

0 .1 0 .1 0 .2
= (9 f0.1 + 19 f0.2 – 5f0.3 + f0.4 ) 0.1/24 + (9 f0.2 + 19 f0.3 – 5f0.4 + f0.5 ) 0.1/24
La tabla de valores a utilizar de la función es:
X 0.1 0.2 0.3 0.4 0.5

F(X) 0.1051709 0.2214028 0.3498588 0.4918247 0.6487213
La integral numérica valuada es
I = 0.016232172 + 0.028456416 = 0.044688589
Si se pudiera elegir h y para realizar un cálculo más rápido se puede utilizar un paso
h=0.2 . En este caso la integral da como resultado
I = 0.044700016
Ya con un paso h= 0.3 la fórmula no es aplicable.
Ejemplo:
Imaginemos una función f(x) cualquiera, se desea integrar entre x=0 y x=1, utilizar la
fórmula Nº18 de la tabla de integrales, proponer por lo menos dos pasos de cálculo con
los que se pueda resolver la integral.
La fórmula Nº18 es del tipo I42 es decir de cuatro franjas con dos puntos pivotes y su
ecuación es
I42= (-8 fi + 16 fi+1) h/2
19
Su representación gráfica es la siguiente
fi+1
fi
xi xi+ 4 h
El caso planteado se puede resolver con un paso máximo de h=0.25 y cualquier otro
submúltiplo, por ejemplo 0.125 aplicando dos veces la fórmula.
Extrapolación de Richardson para Integrales Numéricas
Repasando los conceptos vistos para derivación numérica, se puede decir que una de las
ventajas de utilizar las derivadas numéricas o integrales a partir del desarrollo en serie
de Taylor es que se conoce el orden del error que se comete al aplicarlas.
Conociendo las variaciones del error se pueden realizar mejores estimaciones, en este
concepto se basa la extrapolación de Richardson, la cual permite mejorar la precisión
obtenida por dos cálculos aproximados.
Para demostrar esto se expresa una función derivación o integración como la

aproximación obtenida por serie de Taylor mas el error
Q = Q1 + e1
n
para un h suficientemente pequeño el error es proporcional a h a través de una
constante C
Q = Q1 + C h1n ∴ C = (Q - Q1) / h1n
si se calcula la integral con un segundo paso de cálculo h2 k-veces mas pequeño que h1
entonces la fórmula de integración se aplicará k-veces con este paso h2 apareciendo k
términos de error,
20
si h1= k h2
Q = Q2 + C h2n +....+ C h2n = Q 2 + k C h2n
k-veces
∴ C = (Q – Q2) /k h2n
(Q - Q1) / h1n = (Q – Q 2) /k h2n
(Q - Q1) = (Q – Q2) h1n /k h2n
pero h1= k h2
(Q - Q1) = (Q – Q2) h1n-1 / h2n-1
Q= Q1 – Q2 (h1 / h2)n-1
[ 1– (h1 / h2)n-1 ]
o bien
Q= Q 2 (h1 / h2)n-1- Q 1
[ (h1 / h2)n-1-1 ]
Ejemplo:
Dado el siguiente conjunto de datos x y f(x), calcular la integral numérica en el

intervalo [0,1] usando la fórmula Nº19 de la Tabla de Integrales, hacer el cálculo con 2
asos diferentes y extrapolar luego por Richardson para mejorar el resultado obtenidos.
21
x f(x) x f(x)
0 2 0.8 1.93
0.1 1.8 0.9 2.38
0.2 1.67 1 3
0.3 1.98 1.1 3.79
0.4 1.39 1.2 4.78
0.5 1.37 1.3 5.99
0.6 1.44 1.4 7.43
0.7 1.62 1.5 9.12
La fórmula Nº19 está representada por la ecuación siguiente:
I52= (- 15 fi + 25 fi+1) * h/2 con orden de error n=3
esta fórmula no es de las mas precisas de la tabla ya que calcula un conjunto grande de
franjas utilizando pocos puntos que dan información sobre la función pero vale su
utilización con el fin de promover una correcta interpretación de la utilización de la
tabla.
resuelve la integral en 5 franjas utilizando dos puntos pivotes.
Gráficamente la fórmula se puede esquematizar de la siguiente forma:
Para este caso se utilizarán los pasos de cálculo

h1=0.2 aplicando una sola vez la fórmula y h2=0.1 h
aplicando dos veces la fórmula.
Con h1= 0.2
I1 = (-15 f0 + 25 f0.2) * 0.2/ 2
= ( -15 * 2 + 25 * 1.67) * 0.2/2 = 1.175
x xi+ 5
Con h2= 0.1
I2 = [(-15 f0 + 25 f0.1) * 0.1/ 2] + [(-15 f0.5 + 25 f0.6) * 0.1/ 2]
= [(-15 * 2 + 25 * 1.80) * 0.1/ 2] + [(-15 * 1.37 + 25 * 1.44) * 0.1/ 2]

= 1.5225
Aplicando la fórmula de Extrapolación de Richardson resulta una integral mejorada
Iex = (h1/h2)n-1 I2 – I1
(h1/h2)n-1 –1
22
Iex = (0.2/0.1)3-1 I2 – I1 = (4 I2 – I1)/3 = ( 4 * 1.5225 – 1.175)/3= 1.6383

(0.2/0.1)3-1 –1
23
ANEXO A: TABLA DE DERIVACIÓN NUMÉRICA
hn f(n) Factor Coeficientes de la fórmula Error en Nº de

la derivada f(n) fórmula
1 -1 1 -(1/2) h f ´´ 1
1/2 -3 4 -1 (1/3) h2 f ´´´ 2
1/6 -11 18 -9 2 -(1/4) h3 f iv 3
hf´ 1/12 -25 48 -36 16 -3 (1/5) h4 f v 4
1/2 -1 0 1 -(1/6) h2 f ´´´ 5
1/6 -2 -3 6 -1 (1/12) h3 f iv 6
1/12 -3 -10 18 -6 1 -(1/20) h4 f v 7
1/12 1 -8 0 8 -1 (1/30) h4 f v 8
1 1 -2 1 -h f ´´´ 9
1 2 -5 4 -1 (11/12) h2 f iv 10
1/12 35 -104 114 -56 11 -(5/6) h3 f v 11
h2 f ´´
1 1 -2 1 -(1/12) h2 f iv 12
1 1 -2 1 0 -(1/12) h2 f iv 13
1/12 11 -20 6 4 -1 (1/12) h3 f v 14
1/12 -1 16 -30 16 -1 (1/90) h4 f v 15
1 -1 3 -3 1 -(3/2) h f iv 16
1/2 -5 18 -24 14 -3 (7/4) h2 f v 17
h3 f ´´´
1 -1 3 -3 1 -(1/2) h f iv 18
1/2 -3 10 -12 6 -1 (1/4) h2 f v 19
1/2 -1 2 0 -2 1 -(1/4) h2 f v 20
1 1 -4 6 -4 1 -2 h f v 21
1 1 -4 6 -4 1 -hfv 22
h4 f IV
1 1 -4 6 -4 1 -(1/6) h2 f v 23
24
ANEXO B: TABLA DE INTEGRACIÓN NUMÉRICA

Coeficientes de la fórmula
Inm(i) Factor fi-1 fi fi+1 fi+2 fi+3 fi+4 fi+5 Error en la Nº de
integral formula
I11(i) h 1 0 (h2/2)f´ 1
I12(i) h/2 1 1 -(h3/12)f´´ 2
4
I13(i) h/12 5 8 -1 (h /24)f´´´ 3
I14(i) h/24 9 19 -5 1 -(19/720)h5f iv 4
I15(i) h/720 251 616 -264 106 -19 +(27/1440)h6f v 5
I14(i) h/24 -1 13 13 -1 (11/720) h5f iv 6
I15(i) h/720 -19 346 456 -74 11 -(11/1440) h6f v 7
I13(i) h/12 -1 8 5 -(h4/24)f´´´ 8
I21(i) 2h 0 1 0 (h3/3)f´´ 9
I21(i) 2h 1 0 0 2 h2 f´ 10
5 iv
I23(i) h/3 1 4 1 -(h /90)f 11
I24(i) h/3 1 4 1 0 -(h5/90)fiv 12
I34(i) 3 h/8 1 3 3 1 -(3/80)h5f iv 13
I43(i) 4 h/3 0 2 -1 2 0 +(14/15)h5f iv 14
I23(i) h/12 4 -8 28 0 (h4/3)f´´´ 15
I22(i) h/2 0 4 0 (h3/3)f´´ 16
I32(i) h/2 -3 9 0 0 (27/12) h3 f´´ 17
I42(i) h/2 -8 16 0 0 0 (80/12) h3 f´´ 18
I52(i) h/2 -15 25 0 0 0 0 (175/12) h3 f´´ 19
2
I31(i) h 3 0 0 0 (9/2) h f´ 20
I41(i) h 4 0 0 0 0 8 h2 f´ 21
I51(i) h 5 0 0 0 0 0 (25/2) h2 f´ 22
I32(i) h/12 9 0 27 0 (9/24) h4 f´´´ 23
I43(i) h/12 32 -64 80 0 0 (64/24) h4 f´´´ 24
I33(i) h/12 27 -72 81 0 0 (63/24) h4 f´´´ 25
I43(i) h/12 80 -208 176 0 0 0 (28/3) h4 f´´´ 26
25
26

Unidad 6.1

Cargado por

Copyright:

Formatos disponibles

Unidad 6.1

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Unidad 6.1

Cargado por

Copyright:

Formatos disponibles

Cátedra de Métodos Numéricos

CÁTEDRA DE MÉTODOS NUMÉRICOS

DERIVACIÓN E INTEGRACIÓN NUMÉRICA

Ing. Héctor Bridera

Ampliando el concepto de las funciones discretas, podemos decir que:

La distribución de argumentos xk puede ser igualmente o uniformemente espaciada o no

En el caso de argumentos uniformemente espaciados, se cumple con la siguiente

Las funciones discretas se utilizan:

• cuando no se conoce la función analítica, por ejemplo una

Discretizar un función significa expresarla a través de un conjunto finito de puntos.

Las derivadas numéricas puede obtenerse a través de diversos métodos, a saber:

• por Expansión en Serie de Taylor, cuya ventaja es que el error se obtiene en

• por Diferencias Finitas u Operadores de Diferencias, la desventaja es que para

• por Polinomios de Interpolación, derivando los polinomios, su desventaja es la

Conociendo la función en un punto y sus n-derivadas se puede conocer en forma

expandiendo en Serie de Taylor a la función f en el punto x i+m,

fi ± m = fi ± m h f’i + (mh)2 f’’i ± (mh)3 f’’’i + (mh)4 fIVi ± ...

a) fi+1 = fi + h f’i + (h)2 f’’i +(h)3 f’’’i + (h)4 fIVi + ...

b) fi+2 = fi + 2 h f’i + 2 (h)2 f’’i +4(h)3 f’’’i + 2(h)4 fIVi + ...

c) fi-1 = fi - h f’i + (h)2 f’’i - (h)3 f’’’i + (h)4 fIVi + ...

d) fi-2 = fi - 2 h f’i + 2 (h)2 f’’i - 4(h)3 f’’’i + 2(h)4 fIVi + ...

Combinando linealmente las ecuaciones y truncando la serie hasta un orden de error

Obtención de Derivada Primera

De la ecuación a) se despeja la derivada primera, obteniéndose

El orden de error es el paso elevado a la potencia correspondiente al primer término de

Si h → 0, entonces el término de error ei → 0

f’i= (fi - fi-1 ) 1/h ei ≈ (h/2)fi’’

Combinando las ecuaciones a) y c) a través de la operación resta se obtiene una fórmula

(fi+1 - fi-1 ) = 2h fi’ + (2h3/6) fi’’’

fi’ = (fi+1 - fi-1 ) ei ≈ - (h2/6) fi’’’

Obtención de Derivada Segunda

Combinando las ecuaciones a) y c) a través de la operación suma se obtiene una fórmula

fi’’= (fi+1 -2 fi + fi-1 ) ei ≈ (h2/12) fiIV

fi’’= (-fi-2 - fi-1 + 4 fi - fi+1 - fi+2)

ei ≈ (2/24 + 4/3) h4 fiIV = -17/60 h2 fiIV

Combinando linealmente el conjunto de cuatro ecuaciones planteadas anteriormente por

Por ejemplo la siguiente combinación lineal de ecuaciones, da una derivada tercera:

El orden de error se puede calcular con la siguiente fórmula:

sea m el número de puntos y n el orden de derivada

• En derivadas laterales, el exponente del h del término de error

• En derivadas centrales, el exponente del h del término de error

Sea f(x)= x3-2x2+2

fi’’ = fi-1 – 2 fi + fi+1

Comparar con los valores exactos de la derivada 2º.

La derivada 2º de la función es:

derivando una vez f’(x)= 3x2 –4x

derivando nuevamente f’’(x)= 6x-4

Se realiza la resolución utilizando la siguiente tabla:

x 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

La Tabla de Derivación Numérica

La tabla de derivación numérica que se utilizará (ver tabla completa en Anexo A)

La tabla se origina realizando, al igual que en la introducción, combinaciones lineales

En forma esquemática se puede representar de la siguiente forma:

Derivada x factor Factor y Coeficientes de Error

hn f(n) Factor Coeficientes de la Error en la derivada f(n) Nº de

h3f’’’ 1 -1 3 -3 1 -(3/2) h fIV 16

La interpretación de la tabla conlleva a los siguientes conceptos:

• La suma de todos lo coeficientes es nula ya que si la función f(x) es una

• Los tipos de fórmulas que se presentan son iniciales o inclinadas a la derecha

• Se pueden deducir las fórmulas finales o inclinadas a la izquierda realizando