4.
Integración numérica
4.1 Regla del rectángulo
4.2 Regla del trapecio
4.3. Regla del trapecio mejorado
4.2 Método de Simpson 1/3
4.3 Método de Simpson 3/8
4.4 Método de Romberg
4.5 Cuadratura Gaussiana
Introducción
Con los métodos, podremos integrar funciones dadas en
• forma tabular o bien
• funciones analíticas muy complejas e, incluso,
• integrar aquéllas cuya integral “no existe”* como es el caso de
"#$
sen 𝑥 cos 𝑥
𝑒 , , , 1 − 𝑘 0 sen0 𝑥, etc
𝑥 𝑥
*Se dice que la integral de una función 𝑓(𝑥) no existe, cuando no hay una
función elemental (polinomial, racional, trascendente y sus combinaciones
finitas) cuya derivada sea 𝑓(𝑥).
Los métodos de integración comúnmente usados pueden clasificarse en dos
grupos:
• los que emplean valores dados de la función 𝑓(𝑥) en abscisas equidistantes y
que se conocen como fómulas de Newton-Cotes,
• y aquellos que utilizan valores de 𝑓(𝑥) en abscisas desigualmente espaciadas,
determinadas por ciertas propiedades de familias de polinomios ortogonales,
conocidas como fórmulas de cuadratura gaussiana.
Regla del rectángulo
:
5 𝑓 𝑥 𝑑𝑥 ≈ 𝑏 − 𝑎 𝑓(𝑎)
;
(:";)$
Error: 𝑓 < 𝜉 , 𝑎 < 𝜉 < 𝑏
0
(Ver: Samuel Daniel Conte, Elementary Numerical Analysis An Algorithmic
Approach, McGraw-Hill (1980), Third Edition).
Regla del punto medio
:
𝑎+𝑏
5 𝑓 𝑥 𝑑𝑥 ≈ 𝑏 − 𝑎 𝑓
; 2
(:";)A
Error: 𝑓 << 𝜉 , 𝑎 < 𝜉 < 𝑏 ●
0B
(Ver: Samuel Daniel Conte, Elementary Numerical Analysis An Algorithmic
Approach, McGraw-Hill (1980), Third Edition).
Métodos de Newton-Cotes
:
Para estimar 𝐼 = ∫; 𝑓 𝑥 𝑑𝑥 , los métodos de Newton-Cotes funcionan, en
general, en dos pasos:
1. Se divide el intervalo [𝑎, 𝑏] en 𝑛 intervalos de igual amplitud, cuyos valores
extremos son sucesivamente
Para quedar en la nueva notación 𝑥I = 𝑎 y 𝑥J = 𝑏.
2. Se aproxima 𝑓(𝑥) por un polinomio 𝑝J (𝑥) de grado 𝑛, y se integra para
obtener la aproximación de 𝐼.
Obtención de la Regla del trapecio (𝒏 = 𝟏)
Usaremos el polinomio lineal de Lagrange
𝑝N (𝑥)
Así,
: #O
5 𝑓(𝑥) 𝑑𝑥 ≈ 5 𝑝N 𝑥 𝑑𝑥
; #P
ℎ
= [𝑓 𝑥I + 𝑓(𝑥N )] donde ℎ = 𝑥N − 𝑥I .
2
Regla del trapecio (𝒏 = 𝟏)
:
ℎ área del trapecio con
5 𝑓(𝑥) 𝑑𝑥 ≈ 𝑓 𝑥I + 𝑓 𝑥N =
; 2 vértices 𝑥I , 𝑥N , 𝑓 𝑥I , 𝑓(𝑥N )
donde ℎ = 𝑏 − 𝑎 = 𝑥N − 𝑥I
𝑦 = 𝑝N 𝑥
Regla del trapecio corregido (o mejorado)
(:";)S
Error: 𝑓 (B) 𝜉 , 𝑎 < 𝜉 < 𝑏
T0I
(Ver: Samuel Daniel Conte, Elementary Numerical Analysis An Algorithmic
Approach, McGraw-Hill (1980), Third Edition).
Regla de Simpson 1/3 (𝒏 = 𝟐)
𝑏−𝑎
𝑥N = 𝑥I + = 𝑥I + ℎ 𝑦 = 𝑝0 𝑥
2
𝑏−𝑎
≈ donde ℎ=
2
Fórmulas de Newton-Cotes
Fórmulas de Newton-Cotes
Ejemplo 1 Uso de la regla del trapecio.
a) Aproxime el área 𝐴N bajo la curva de la función dada por la tabla siguiente,
en el intervalo 𝑎 = 500, 𝑏 = 1800.
Solución
Ejemplo 2 Uso de la regla de Simpson 1/3.
a) Aproxime el área 𝐴N bajo la curva de la función dada por la tabla siguiente,
en el intervalo 𝑎 = 500, 𝑏 = 1800.
Solución
[ 𝑓(𝑥N ) se obtiene interpolando con un polinomio de Newton de segundo
grado ]
El grado de exactitud o precisión de una fórmula de cuadratura es el
entero positivo más grande 𝑛, tal que la fórmula sea exacta para 𝑥 [ ,
cuando 𝑘 = 0,1, … , 𝑛.
• La regla del trapecio tiene un grado de precisión de uno.
• La regla de Simpson 1/3 tiene un grado de precisión de tres.
Esto es, produce el resultado exacto con cada polinomio de grado tres o
menor.
Ilustración de la
Regla compuesta
del trapecio
Ilustración de la
Regla compuesta
del trapecio (con
su error)
𝑏 −𝑎 0
ℎ 𝑓′′(𝜉)
12
Regla compuesta del trapecio
En notación de sumatoria
𝑏 −𝑎
donde ℎ = = 𝑥^_N − 𝑥^
𝑛
para 𝑖 = 0, 1, … , (𝑛 − 1)
Ejemplo 1
Mediante la regla compuesta del trapecio, aproxime el área bajo la curva de la
siguiente función dada en forma tabular, entre 𝑥 = −1 y 𝑥 = 4.
Solución
Ejemplo 1
Comparamos este resultado con la solución analítica (los datos corresponden a
la función 𝑓 𝑥 = 𝑥 B − 2𝑥 0 + 𝑥 + 10).
B
5 (𝑥 B −2𝑥 0 + 𝑥 + 10) 𝑑𝑥 = 219.167
"N
Ejemplo 3
Encuentre la integral aproximada de la función
que da lugar a la curva normal estándar, entre los límites −1 y 1.
Utilice la regla trapezoidal con varios trapezoides y compare con el resultado
(0.682) obtenido de tablas.
Solución
Existen dos categorías de fórmulas de Newton-Cotes, abiertas y cerradas.
Fórmulas cerradas comunes de Newton Cotes (con su error)
Error de truncamiento
Error de truncamiento
Fórmulas cerradas comunes de Newton Cotes (con su error)
Fórmulas abiertas comunes de Newton Cotes (con su error)
Fórmulas abiertas comunes de Newton Cotes (con su error)
Extrapolación de Richardson
• Con el nombre de extrapolación de Richardson se conoce a un conjunto de
técnicas que generan mejores aproximaciones a los resultados buscados o
aproximaciones equivalentes a métodos de alto orden, a partir de las
aproximaciones obtenidas por medio de algún método de bajo orden y pocos
cálculos.
• Dichas técnicas están basadas en el análisis del error de truncamiento, cuya
aplicación a la integración numérica se presenta a continuación.
Integración de Romberg
Supóngase que el error de truncamiento de cierto algoritmo de aproximación
de
se expresa como 𝐸 = 𝑐ℎh 𝑓 h
(𝜉), donde
𝑐 es una constante, 𝑟 es un entero positivo, ℎ es el tamaño de paso y 𝑎 < 𝜉 < 𝑏.
Por ejemplo, observe el error en la regla compuesta del trapecio,
𝑏 −𝑎 0
ℎ 𝑓′′(𝜉)
12
Luego de obtener dos aproximaciones de 𝐼, con tamaños de paso distintos: ℎN
y ℎ0 , de llamar a dichas aproximaciones 𝐼N y 𝐼0 , respectivamente, y despreciar
errores de redondeo, se puede escribir
h h
Dividiendo y suponiendo que 𝑓 (𝜉N ) y 𝑓 (𝜉0 ) son casi iguales, se tiene
Despejando 𝐼:
• Si en particular ℎ0 = ℎN /2, la ecuación
• se simplifica a
• Este proceso, conocido como integración de Romberg, es efectivo cuando
𝑓 h (𝑥) no varía bruscamente en (𝑎, 𝑏), y no cambia de signo en dicho
intervalo.
• En estos casos, las ecuaciones anteriores permiten obtener una mejor
aproximación a 𝐼 a partir de 𝐼N y 𝐼0 , sin repetir el proceso de integración y con
cálculos breves.
• En la regla compuesta del trapecio, por ejemplo, 𝑟 = 2 y la ecuación
toma la forma
Para sistematizar, usamos la notación:
(I)
• 𝐼[ = las aproximaciones de 𝐼 obtenidas empleando 2[ trapezoides,
𝑘 = 0, 1, 2, …
estos resultados sirven para producir la primera extrapolación ( ver la tabla).
(P) (P)
(N) 0$ jklO "jk
• 1a extrapolación: 𝐼[ = , 𝑘 = 0, 1, 2, …
0$ "N
estos resultados sirven para producir la segunda extrapolación (ver la tabla).
(O) (O)
(0) 0m jklO "jk
• 2da extrapolación: 𝐼[ = , 𝑘 = 0, 1, 2, …
0m "N
estos resultados sirven para producir la 3ra extrapolación.
Así se continua. Ver la siguiente fórmula general:
(I)
𝐼[ = las aproximaciones de 𝐼 obtenidas empleando 2[ trapezoides, 𝑘 = 0, 1, 2, …
(I) (I) (I) (I)
(N) 20 𝐼[_N − 𝐼[ 4𝐼[_N − 𝐼[
𝐼[ = =
20 − 1 4−1
B (N) (N) 0 (N) (N)
(0) 2 𝐼[_N − 𝐼[ 4 𝐼[_N − 𝐼[
𝐼[ = =
2B − 1 40 − 1
• En general, tenemos la siguiente fórmula del método de Romberg
𝑘 = 0, 1, 2, …
Tabla. Aplicación del método de Romberg
(I)
Cuando los valores de 𝐼[ → 𝐼 al crecer 𝑘, los valores de la diagonal superior
de la tabla convergen a 𝐼.
Ejemplo 1 Uso del método de Romberg.
Encuentre una aproximación de la integral
N
5 sen 𝜋𝑥 𝑑𝑥
I
empleando 1, 2, 4, 8 y 16 trapecios.
Con los resultados obtenidos y la fórmula del método de Romberg, obtenga
mejores aproximaciones. Compare los valores obtenidos con el valor, calculado
analíticamente: 0.6366197.
Solución
• Función, usando mathematica, para integrar con la regla compuesta del
trapecio.
• Con ella se obtienen lo valores de la 3er columna en la siguiente
diapositiva.
(B)
El valor 𝐼B es el valor analítico de la integral.
Ilustración de la cuadratura gaussiana
(o cuadratura de Gauss)
a) Representación gráfica de la regla del
trapecio como el área bajo la recta
que une los puntos extremos fijos.
b) Se obtiene una mejor estimación de la
integral tomando el área bajo la recta
que pasa por dos puntos intermedios.
Estos puntos se ubican en una forma
adecuada, de tal manera que se equilibran
los errores positivo y negativo.
a) Representación gráfica de la regla del trapecio como el área bajo la recta
que une los puntos extremos fijos.
b) Se obtiene una mejor estimación de la integral tomando el área bajo la
recta que pasa por dos puntos intermedios.
Estos puntos se ubican en una forma adecuada, de tal manera que se
equilibran los errores positivo y negativo.
a) Regla del trapecio b) Cuadratura gaussiana con dos puntos
Observación. Cualquiera de las fórmulas de integración desarrolladas
anteriormente (las de Newton-Cotes) puede escribirse en la forma
J
:
5 𝑓 𝑥 𝑑𝑥 ≈ r 𝑤^ 𝑓(𝑥^ )
; ^ sN
• Por ejemplo, en la regla del trapecio
ℎ ℎ ℎ
𝐼 ≈ 𝑓 𝑥I + 𝑓 𝑥N = 𝑓 𝑥I + 𝑓 𝑥N = 𝑤I 𝑓 𝑥I + 𝑤N 𝑓 𝑥N
2 2 2
ℎ
donde 𝑤I = 𝑤N =
2
ℎ 4ℎ
• En la regla de Simpson 1/3: 𝑤I = 𝑤0 = y 𝑤N =
3 3
Desarrollo de la fórmula de cuadratura gaussiana de dos puntos
Considérese primero que se
desea integrar la función
𝐹 𝑧 entre los límites –1 y
+1.
−1
Motivado por la fórmula de la regla del trapecio, Gauss se propuso desarrollar
una fórmula del tipo
N
5 𝑓 𝑥 𝑑𝑥 ≈ 𝑤N 𝐹 𝑧N + 𝑤0 𝐹(𝑧0 )
"N
donde se tienen cuatro parámetros por determinar y, en consecuencia, se
requieren cuatro condiciones para determinarlas.
Étas se eligen de manera que el método dé resultados exactos
cuando la función por integrar sea alguna de las cuatro siguientes,
o combinaciones lineales de ellas:
Los valores exactos de integrar estas cuatro funciones entre –1 y +1 son
2/3
Si la ecuación 𝑤N 𝐹 𝑧N + 𝑤0 𝐹(𝑧0 ) da resultados exactos para las funciones
1, 𝑧, 𝑧 0 , 𝑧 v , se tiene el siguiente sistema de ecuaciones:
Solución del sistema
1 1
𝑤N = 𝑤0 = 1 𝑧N = − = −0.57735 … 𝑧0 = = 0.57735 …
3 3
Solución usando Mathematica
De forma numérica
Así,
1 1
𝑤N = 𝑤0 = 1, 𝑧N = − = −0.57735 … , 𝑧0 = = 0.57735 …
3 3
Sustituyendo, obtenemos la fórmula de cuadratura gaussiana de dos puntos:
N
1 1
5 𝐹 𝑧 𝑑𝑧 = 𝑤N 𝐹 𝑧N + 𝑤0 𝐹 𝑧0 = 𝐹 − +𝐹
3 3
"N
= 𝐹 −0.57735 … + 𝐹(0.57735 … )
O simplemente
N • Grado de precisión tres
1 1 (produce el resultado exacto
5 𝐹 𝑧 𝑑𝑧 ≈ 𝐹 − +𝐹
3 3 con cada polinomio de grado
"N
tres o menor).
• Observe que los límites de integración son de −1 a 1.
:
• Una integral ∫; 𝑓 𝑥 𝑑𝑥 en un intervalo arbitrario [𝑎, 𝑏] se puede
transformar en otra en [−1, 1] usando el cambio de variable
ya que si 𝑥 = 𝑎, 𝑧 = −1, y si 𝑥 = 𝑏, 𝑧 = 1.
Si despejamos 𝑥, obtenemos
𝑏−𝑎 𝑎+𝑏 𝑏−𝑎 𝑎+𝑏 𝑏−𝑎
𝑥= 𝑧+ y 𝑑𝑥 = 𝑑 𝑧+ = 𝑑𝑧
2 2 2 2 2
Esto nos permite aplicar la cuadratura gaussiana a cualquier intervalo [𝑎, 𝑏]:
El método de Gauss puede extenderse a tres o más puntos
• Por ejemplo, si se escogen tres puntos no equidistantes en el segmento de
la curva 𝑓(𝑧), comprendida entre –1 y 1, se podría pasar una parábola por
los tres como en la regla de Simpson, excepto que dichos puntos se
escogerían de modo que minimicen o anulen el error.
• Similarmente es factible elegir cuatro puntos y una curva cúbica,
• cinco puntos y una curva de grado cuatro, etc.
• En general, el algoritmo tiene la forma
• La siguiente tabla presenta valores hasta para seis puntos.
• Con dos puntos se tiene un grado de exactitud de tres.
• Con tres puntos se tiene un grado de exactitud de cinco.
• Con cuatro puntos se tiene un grado de exactitud de siete y así
sucesivamente.
• Los coeficientes y las abscisas dadas en la tabla sirven para integrar sobre
todo el intervalo de interés, o bien puede dividirse el intervalo en varios
subintervalos (como en los métodos compuestos de integración) y aplicar el
método de Gauss a cada uno de ellos.
Coeficientes y abscisas en el método de la cuadratura de Gauss.
N $ /0
Ejemplo 1 Integre la función 𝑒 "# en el intervalo (−0.8, 1.5) por
0w
cuadratura de Gauss.
Solución
Con el cambio de la función en términos de la nueva variable 𝑧, queda:
N.x
1 $ /0
𝐼= 5 𝑒 "# 𝑑𝑥
2𝜋 "I.y
$
N N.x" "I.y "I.y_N.x
1 1.5 − (−0.8) "
0
z_
0
/0
= 5 𝑒 𝑑𝑧
2𝜋 "N 2
N
2.3 $ /y
= 5 𝑒 " 0.vz_I.T 𝑑𝑧
2 2𝜋 "N
De la tabla,
Al evaluar la función del integrando en 𝑧N , 𝑧0 :
" 0.v(I.xTTvxI0{|0 )_I.T $ /y
𝐹 0.5773502692 = 𝑒 = 0.5980684
" 0.v("I.xTTvxI0{|0 )_I.T $ /y
𝐹 −0.5773502692 = 𝑒 = 0.9519115
Se aplica la fórmula
y se obtiene
se tiene
Ejemplo 2
Solución
𝑥 −𝜋
𝜋
Con el empleo de la fórmula con 𝑛 = 3 y los valores de la tabla, queda:
Recordar la fórmula
Se puede poner en forma más general como
