Análisis Numérico

para Ingeniería

Clase Nro. 11
 Integración Numérica

Temas a tratar:
Métodos de Newton-Cotes.
Método de los Trapecios.
Método de 1/3 de Simpson.
Método de 3/8 de Simpson.
Método de Romberg.

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 2

Problema a resolver

Se desea calcular numéricamente:

b

I f  =∫ f  x⋅dx
a

Cuando la función a integrar no tiene anti-

derivada explícita.
La función antiderivada es muy complicada de
calcular.
La función está expresada en forma discreta.
Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 3
Método de Resolución
El método básico utilizado para obtener el valor de una
Integral definida, se denomina cuadratura numérica y
está basado en la siguiente fórmula:

∑ H i⋅f  x i 
i=0

Entonces tendremos:
b n
I f  =∫ f  x⋅dx =∑ H i⋅f  x i Error
a i=0

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 4

Estrategias para Integrar Numéricamente

Si se conocen los valores de las abscisas x i con

i=0,1,...,n y éstos son equiespaciados, se
utilizarán las fórmulas de Newton Cotes para hallar
los valores de H i
Otras forma de calcular los H i se logran utilizando las
fórmulas Gauss-Legendre, o de Cuadratura
Gaussiana.

Otra forma consiste en hallar la Serie de Taylor e

integrar el polinomio resultante al truncar la misma.

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 5

Métodos de Newton-Cotes
La estrategia consiste en aproximar los puntos de la
función f(x) y luego integrar el polinomio resultante:

b b

I f  =∫ f  x⋅dx ≈∫ pn  x⋅dx
a a

Y considerando el error de interpolación, tendremos:

b b b

I f  =∫ f  x⋅dx =∫ p n  x⋅dx∫ E  x⋅dx

a a a

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 6

Métodos de Newton-Cotes

Teniendo en cuenta el error de interpolación

para puntos equiespaciados, obtenemos la
siguiente expresión para el Error de
Integración:

a a
s⋅(s−1)⋅(s−2)⋯( s−n) n+1 (n+1)
Error =∫ E (x)⋅dx= ∫ ⋅h ⋅f (ξ)⋅dx
b b (n+1)!

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 7

 Métodos de Newton-Cotes
Esta fórmula puede utilizarse de diferentes formas, ya
que el rango de ajuste del polinomio no
necesariamente tiene que coincidir con el intervalo de
integración.
b b b

I f  =∫ f  x⋅dx =∫ p n  x⋅dx∫ E  x⋅dx

a a a

En el caso en que los extremos de integración

coincidan con el intervalo de interpolación se
llaman fórmulas cerradas de Newton-Cotes, de
lo contrario se denominan fórmulas abiertas.
Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 8
Métodos de Newton-Cotes
Las fórmulas de Newton-Cotes más utilizadas son
aquellas que poseen mejor ajuste y menor error de
redondeo, es decir, aquellas que utilizan un polinomio
de interpolación de grado bajo.

n=1 Trapecios
NEWTON −COTES
CERRADAS
{n=2

n=3
1
3
3
8
de Simpson

de Simpson

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 9

Método de Trapecios
Dados dos puntos es posible aproximar por un
polinomio interpolante de Newton de grado 1.

y
p1(x)

f(x)

x0 x1 x

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 10

 Método de Trapecios
Aproximando con un polinomio de 1er. Grado,
obtenemos:

x1 x1

∫ f  x⋅dx ≈∫  y 0s⋅ y 0 ⋅dx =

x0 x0
2 1
s  y0
2 0
1
= h⋅[ y 0⋅s∣0 y 0⋅ ] = h⋅[ y 0
2∣ ] =
y 1− y 0 h
= h⋅[ y 0 ] = ⋅[ y 1 y 0 ]
2 2

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 11

Teorema del Valor Medio para Integrales
Si f es una función contínua en el
intervalo [a, b] , entonces existe un
número z en [a, b] tal que:

∫ f  x⋅dx = f  z⋅b−a
a

Veamos una interpretación geométrica del teorema:

Supongamos que f(x) ≥ 0 para todo x que pertenece a [a, b] ,

en este caso, la integral de f(x) se toma como el área de la
región encerrada bajo la curva f(x) y las rectas x=a y x=b,
entonces se demuestra que existe al menos un valor z que
pertenece a [a, b] , que verifica que f(z)*(b-a) es igual al área
antes mencionada.
Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 12
Error en el Método de Trapecios
Por lo tanto, el error cometido en la integración es:

b x1
s⋅s−1 2
Error = ∫ E  x⋅dx = ∫ ⋅h f ' ' ⋅dx =
a x 2 0
1
s⋅s−1
3
(por TVM) = h ⋅f ' ' 1 ∫ ⋅ds =
0 2
3 2 1 3
s s
3
= h ⋅f ' ' 1 ⋅ −  =
6 4 0
−h
12 ∣
⋅f ' ' 1 

Donde x0 ≤ ξ1 ≤ x1

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 13

Método de Trapecios Compuestos
Si subdividimos el intervalo [a, b] en sub-intervalos,
equiespaciados o no, la integral aproximada, será
igual a la suma de las áreas de los trapecios
correspondientes a cada uno de los sub-intervalos.

Por ejemplo, para dos puntos cualesquiera dentro

del intervalo [a, b] , tenemos:

x i1
h
∫ f  x⋅dx ≈ 2⋅[ y i y i1 ]
x i

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 14

Método de Trapecios Compuestos
Por lo tanto, si extendemos esta expresión para calcular
la integral aproximada dentro del intervalo [a, b] y
considerando n puntos equiespaciados, tenemos :

b n−1
h
∫ f ( x)⋅dx ≈ ∑ 2
⋅[ y i + y i +1 ] =
a i=0
h
= ⋅[ y 0 + y 1 + y 1 + y 2 + y 2 +⋯+ y n−1 + y n−1 + y n ] =
2
h
= ⋅[ y 0 +2⋅y 1 +2⋅y 2 +⋯+2⋅y n−1 + y n ]
2

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 15

Error en Trapecios Compuestos

El error local correspondiente a un trapecio es :

−h3
Error local = ⋅f ' ' 1 
12

Por lo tanto, podemos calcular el error global del

método, sumando todos los errores locales :

3
−h
Error global = ⋅[ f ' ' (ξ 1 )+ f ' ' (ξ2 )+⋯+f ' ' (ξ n−1 )+ f ' ' (ξ n )]
12

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 16

Error en Trapecios Compuestos
Los ξi corresponden a cada uno de los subintervalos. Si
asumimos que f ''(x) es contínua en (a, b), entonces habrá
algún valor ξ perteneciente al intervalo (a, b), para el cual:
3
−h
Error global = ⋅[ f ' ' (ξ1 )+ f ' ' (ξ2 )+⋯+f ' ' (ξ n−1 )+ f ' ' (ξn )]=
12
−h3 −h3 (b−a)
= ⋅n⋅f ' ' (ξ) = ⋅ ⋅f ' ' (ξ)
12 12 h

Por lo tanto, nos queda:

(a−b) 2
Error global = ⋅h ⋅f ' ' (ξ)
12
Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 17
Ejemplo de Trapecios Compuestos

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 18

Método de los Trapecios Compuestos

FUNCTION trapecios(n, h, v)
! Función que calcula la integral numérica
! por el método de los Trapecios
INTEGER n
REAL(8) v(n)
REAL(8) trapecios, h
trapecios = h*(v(1)+2*SUM(v(2:n-1))+v(n))/2.0
END FUNCTION

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 19

Método de 1/3 de Simpson

Aproximando con un polinomio de 2do. Grado,

obtenemos:

x2 x2
s⋅s−1 2
∫ f  x⋅dx ≈∫  y 0s⋅ y 0 2 ⋅ y 0 ⋅dx =
x 0 x 0
2
s⋅s−1 2
= h⋅∫  y 0s⋅ y 0 ⋅ y 0 ⋅ds =
0 2
h
= ⋅ y 04⋅y 1 y 2 
3

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 20

Error en 1/3 de Simpson
Si denominamos :
x2

f ( x) = F ' ( x) ⇒ ∫ f ( x)⋅dx = F ( x 2)−F ( x 0 )

x0

Desarrollando la integral de f(x) en Serie de Taylor,

obtenemos:

x2 2 3
(2h ) (2h)
∫ f (x )⋅dx =F ' (x 0 )⋅(2h)+ 2 ⋅F ' ' (x 0 )+ 6 ⋅F ' ' ' ( x 0 ) +
x 0

(2h)4 iv (2h )5 v
+ ⋅F (x 0 )+ ⋅F ( x 0 )+⋯
24 120

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 21

Error en 1/3 de Simpson

Reemplazando F'(x0) = f(x0) = y0 , finalmente

nos queda :

x2
8 32
∫ f  x⋅dx =2h⋅y 02h ⋅y ' 0 6⋅h ⋅y 0 ' ' +
x 0

16 4 32 5 iv
+ ⋅h ⋅y 0 ' ' ' ⋅h ⋅y 0 ⋯
24 120

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 22

Error en 1/3 de Simpson
Para estimar el error de truncamiento del método,
utilizaremos el desarrollo en Serie de Taylor del valor
aproximado de la integral, obtenida por 1/3 de Simpson.

h
⋅( y 0 +4⋅y 1 + y 2 ) =
3
2 3 4
h h h h iv
= ⋅[ y 0 +4⋅( y 0 +h⋅y 0 ' + ⋅y 0 ' ' + ⋅y 0 ' ' ' + ⋅y 0 +⋯)+
3 2! 3! 4!
(2 h)2 (2 h)3 (2 h)4 iv
+( y 0 +(2 h)⋅y 0 ' + ⋅y 0 ' ' + ⋅y 0 ' ' ' + ⋅y 0 +⋯)=
2! 3! 4!

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 23

Error en 1/3 de Simpson

Realizando las simplificaciones correspondientes en la

expresión anterior, tenemos:

h 4 4 2
= ⋅[ 6⋅y 0 +6⋅h⋅y 0 ' +( + )⋅h ⋅y 0 ' ' +
3 2 2
4 8 3 4 16 4 iv
+( + )⋅h ⋅y 0 ' ' ' +( + )⋅h ⋅y 0 +⋯] =
6 6 24 24
h 2 3 5 4 iv
= ⋅[ 6⋅y 0 +6⋅h⋅y 0 ' +4⋅h ⋅y 0 ' ' +2⋅h ⋅y 0 ' ' ' + ⋅h ⋅y 0 +⋯]
3 6

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 24

Error en 1/3 de Simpson

Ahora, para obtener una expresión del error hacemos la

diferencia entre la integral exacta y la integral aproximada.

x2
h
∫ f ( x)⋅dx− 3⋅( y 0 +4⋅y 1 + y2 ) =
x 0

6 2 6 2
=(2⋅h⋅y 0− ⋅h⋅y 0 )+(2⋅h ⋅y 0 ' − ⋅h ⋅y 0 ' )+
3 3
16 2 32 5
+( ⋅h 4⋅y 0 ' ' '− ⋅h 4⋅y 0 ' ' ' )+( ⋅h5⋅y 0vi− ⋅h5⋅y iv0 )=
24 3 120 18
−1 5 iv
= ⋅h ⋅y 0 +⋯
90

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 25

Error en 1/3 de Simpson
Este es el error de truncamiento de la serie.

−1 5 iv
Error de Truncamiento= ⋅h ⋅y 0 +⋯
90

Para estimar el error local, recordemos el término del

error de interpolación para puntos equiespaciados.
Al integrarlo, nos queda:

b n (n+1)
n+2 f (ξ)
Error=∫ E (x )⋅dx=h ⋅∫ s⋅(s−1)⋅(s−2)⋯(s−n)⋅ ⋅ds
a 0 (n+1)!

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 26

Error en 1/3 de Simpson
Uno esperaría que la integración de 1/3 de Simpson
sea exacta si f(x) es un polinomio de grado menor o
igual a 2, sin embargo, resulta exacta, también para un
polinomio de grado 3, ya que :

ERROR DE INTERPOLACIÓN
PARA UN POLINOMIO DE GRADO 3

2
f ' ' '  4
Error = h ⋅ ⋅∫ s⋅s−1⋅s−2⋅ds = 0
3! 0

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 27

Error en 1/3 de Simpson
Por lo tanto, si calculamos el error de interpolación
para un polinomio de grado 4, obtenemos :

(iv ) 2
f 5 (ξ)
Error = h ⋅ ⋅∫ s⋅(s−1)⋅(s−2)⋅(s−3)⋅ds =
4! 0
(iv )
(ξ) 5 f
= −h ⋅
90

Donde x0 ≤ ξ ≤ x2
Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 28
Método de 1/3 de Simpson Compuesto
Si tenemos un número impar de puntos
equiespaciados dentro de nuestro intervalo de
integración, podemos aplicar la regla de 1/3 de Simpson
en forma repetida, y si luego sumamos las áreas,
obtenemos :

∫ f ( x)⋅dx =
a
h
= ⋅( y 0 + 4⋅y 1 +2⋅y 2 + 4⋅y 3 +⋯+2⋅y n−2 +4⋅y n−1 + y n )+
3
+ Error Global

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 29

Error en 1/3 de Simpson Compuesto
Sumando los errores locales de cada uno de los arcos
de parábola, obtenemos una expresión para el error
global.
5
−h (iv) (iv ) (iv)
Error Global= ⋅[f (ξ1 )+f (ξ 2)+⋯+f (ξn )]=
90 2
b−a
5 5
−h n (iv) −h h (iv )
= ⋅ ⋅f (ξ)= ⋅ ⋅f (ξ)=
90 2 90 2
5
−h b−a (iv ) −(b−a) 4 (iv )
= ⋅ ⋅f (ξ) = ⋅h ⋅f (ξ)
90 2⋅h 180

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 30

Ejemplo de 1/3 de Simpson Compuesto

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 31

Método de 1/3 de Simpson Compuesto

FUNCTION simpson(n, h, v)
! Función que calcula la integral numérica
! por el método de 1/3 de Simpson
INTEGER n
REAL(8) v(n)
REAL(8) simpson, h

simpson = h*(v(1)+4*SUM(v(2:n-1:2))+2*SUM(v(3:n-2:2))+v(n))/3.0

END FUNCTION

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 32

Método de 3/8 de Simpson

Está basado en la integral de una aproximación

polinomica de 3er. Grado:

x3 x3
3⋅h
∫ f ( x)⋅dx ≈∫ p3 ( x)⋅dx= 8 ⋅( y 0+ 3⋅y 1+ 3⋅y 2+ y 3 )
x 0 x 0

−3 5 (v )
Error Local= ⋅h ⋅f (ξ)
80

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 33

Método de 3/8 de Simpson Compuesto

Sumando los errores locales de cada una de las

aproximaciones, obtenemos una expresión para el
error global.
x3 x3
3⋅h
∫ f  x⋅dx ≈∫ p3  x⋅dx= 8 ⋅( y 03⋅y 13⋅y 2 +
x 0 x 0

+ 2⋅y 3 3⋅y 4 3⋅y 5⋯2⋅y n−23⋅y n−1 y n ]

−b−a 4  v 
Error Global= ⋅h ⋅f 
80

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 34

Comparación de Errores

Si comparamos los errores globales de los diferentes

métodos, tenemos :
(a−b) 2
Trapecios Compuestos = ⋅h ⋅f ' ' (ξ)
12

1 −(b−a) 4 (iv )
Simpson Compuestos= ⋅h ⋅f (ξ)
3 180

3 −(b−a) 4 (v )
Simpson Compuestos= ⋅h ⋅f (ξ)
8 80

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 35

Método de Romberg
Aunque el método de los trapecios es muy sencillo de aplicar,
carece de la exactitud requerida habitualmente. La integración de
Romberg es un método que utiliza inicialmente los valores
obtenidos por el método de los trapecios y posteriormente
aplica el proceso conocido como ”extrapolación de
Richardson” para obtener correcciones a las aproximaciones
anteriores.

Sabemos que el error de Trapecios Compuestos es:

3
a−b 2 a−b
Error global = ⋅h ⋅f ' '  = ⋅f ' ' 
12 12⋅n 2

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 36

Método de Romberg

Sean n1 y n2 dos cantidades de trapecios distintas.

Y sea I* el verdadero valor de la integral y In el valor
de la integral aproximada por n trapecios.

Entonces:

*
I = I n1E n1 = I n2 E n2

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 37

Método de Romberg

Entonces:

3
b−a
2
⋅f ' ' 2 
En2
12⋅n2
=
En b−a 3
1
2
⋅f ' ' 1 
12⋅n1

Donde ξ1 y ξ2 se encuentran en (a, b)

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 38

Método de Romberg
Asumiendo que f ''(ξ1) = f ''(ξ2) nos queda :

2
n1
En 2
=
n2 [ ]
⋅E n 1

I n1E n1 −I n2 −E n2 = 0
2
n1
I n1E n1 −I n2 −
n2 [ ]
⋅E n = 0 1

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 39

Método de Romberg

2
n1
 I n −I n E n ⋅ 1−
1 2
n2
1
{ [ ]} = 0

Por lo tanto :

I n −I n
En = 2 1

2
n1
1

1−
n2 [ ]
Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 40
Método de Romberg

*
I n −I n
I =I n  2 1

2
n1
1

1−
n2 [ ]
Por lo tanto, si hacemos n2 = 2*n1 :

* 4 1
I = ⋅I n − ⋅I n
3 3 2 1

Es de esperar que I* sea mejor que In1 e In2

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 41
Método de Romberg

Sistematizando este proceso denominado

Extrapolación de Richardson, obtenemos el
método de integración de Romberg.

Si tenemos 1 Trapecio h = (b - a), por lo tanto:

b−a 1
T 0, 1 = ⋅ ⋅[ f af b]
1 2

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 42

Método de Romberg

Si tenemos 2 Trapecios h = (b - a)/2, por lo tanto:

b−a 1 b−a
T 1, 1 =
2 2 {
⋅ ⋅[ f af b]f a
2
 = }
1 b−a
=
2 {
⋅ T 0,1b−a⋅f a
2
 }
Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 43
Método de Romberg

Si tenemos 2N Trapecios h = (b - a)/(2N), por lo tanto:

N
2 −1
T N ,1
1
2 {
b−a
2 i=1
b−a
= ⋅ T N −1, 1 N −1 ⋅ ∑ f a N ⋅i
2
 i=2
}
Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 44
Método de Romberg
Si tomamos dos valores consecutivos de estas
aproximaciones, podemos aplicar la fórmula ya vista
para encontrar una mejor aproximación.

4⋅T N 1,1−T N , 1
T N ,2 =
3

Por ejemplo, si tomamos N=0, nos queda:

4⋅T 1,1−T 0, 1 b−a b−a

T 0, 2 =
3
=
6 [
⋅ f a4⋅f a
2
f b ]
Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 45
Método de Romberg

Ahora podemos generalizar el esquema visto, para j ≥ 2.

j−1
4 ⋅T N 1, j−1 −T N , j−1
TN,j = j−1
4 −1

Donde j es el grado del polinomio aproximante y 2N es

la cantidad de Trapecios utilizados para obtener dicha
aproximación.
Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 46
 Método de Romberg paso 1
Cálculo del área de 1 trapecio, es decir, T0,1

a b

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 47

 Método de Romberg paso 2
Cálculo del área de 2 trapecios, es decir, T1,1

a b

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 48

 Método de Romberg paso 3
Cálculo del área de 4 trapecios, es decir, T2,1

a b

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 49

 Método de Romberg paso 4
Cálculo del área de 8 trapecios, es decir, T3,1

a b

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 50

Método de Romberg

VALORES CALCULADOS VALORES CALCULADOS

POR TRAPECIOS POR ROMBERG

T0,1
T0,2
T1,1 T0,3
T1,2 T0,4
T2,1 T1,3
T2,2
VALOR FINAL
T3,1 DE LA INTEGRAL

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 51

Método de Romberg
Los valores de la segunda columna en adelante se
calculan por medio de esta fórmula:

Fórmula de
T0,1
cálculo
T0,2
T1,1 T0,3 j−1
4 ⋅T N 1, j−1 −T N , j−1
T1,2 T0,4 TN,j =
4 j−1 −1
T2,1 T1,3
T2,2
T3,1 VALOR FINAL
DE LA INTEGRAL

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 52

 Ejemplo de Romberg

-56.5195
-50.8040
-52.2328 -11.9705
-14.3976 -0.0828
-23.8564 -0.2685 -1.33069
-1.1516 -1.3258 -1.42647
-6.8278 -1.3093 -1.42638 -1.42604
-1.2994 -1.42598 -1.42604 -1.42603
-2.6815 -1.42416 -1.42604 -1.42603 -1.42603
-1.4164 -1.42604 -1.42603 -1.42603
-1.7327 -1.42600 -1.42603 -1.42603
-1.4254 -1.42603 -1.42603
-1.5022 -1.42602 -1.42603
-1.4226 -1.42603
-1.4450 -1.42602
-1.4260
-1.4308

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 53

Ejemplo de Romberg (II)
En este caso vemos que aún tomando muchos menos
puntos, 64 en lugar de los 256 de la tabla anterior,
obtenemos casi el mismo resultado final.

-56.5195
-50.8040
-52.2328 -11.9705
-14.3976 -0.0828
-23.8564 -0.2685 -1.33069
-1.1516 -1.3258 -1.42647
-6.8278 -1.3093 -1.42638 -1.42604
-1.2994 -1.42598 -1.42604
-2.6815 -1.42416 -1.42604
-1.4164 -1.42604
-1.7327 -1.42600
-1.4254
-1.5022

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 54

 Método de Romberg

FUNCTION romberg(n, h, v)
! Función que calcula la integral numérica
! por el método de Romberg
INTEGER n, i, j, nt, nv, step
REAL(8) h, romberg, v(n)
REAL(8), ALLOCATABLE :: t(:)

step = n-1
nt = LOG(REAL(step))/LOG(2.)
ALLOCATE(t(0:nt))

DO i=0, nt
t(i) = trapecios(INT(2.**i)+1, step*h, v(::step))
step = step/2
ENDDO Sique en la
próxima página
Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 55
Método de Romberg

nv = nt-1

DO j=2, nt+1
DO i=0, nv
t(i) = (4**(j-1)*t(i+1) - t(i))/(4**(j-1)-1)
ENDDO
nv = nv - 1
ENDDO

romberg = t(0)

END FUNCTION

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 56

Comparación entre métodos
TRAPECIOS
Es sencillo de aplicar.
Puede utilizarse con cualquier cantidad de puntos.
No integra con demasiada exactitud, a menos que la
distancia entre puntos sea muy pequeña.
1/3 de SIMPSON
Es sencillo de aplicar.
Solo puede utilizarse con una cantidad impar de
puntos.
El error de integración es muy pequeño.

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 57

Comparación entre métodos

3/8 de SIMPSON
Su fórmula es apenas un poco más complicada que la de
1/3 de Simpson.
Solo puede utilizarse con una cantidad de puntos igual a
3*n+1, siendo n el número de arcos.
Tiene un error ligeramente mayor que 1/3 de Simpson.
ROMBERG
Se inicia con valores obtenidos por Trapecios.
Solo puede utilizarse con una cantidad de puntos igual a
2N+1, siendo 2N el número máximo de trapecios.
El error de integración es muy pequeño.

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 58

PREGUNTAS ...

Mg. Ing. Francisco A. Lizarralde Facultad de Ingeniería - UNMDP - 2021 59