Integración Numérica

MSc. Paulina Morillo

 INTRODUCCION

La integración numérica es una herramienta esencial que se usa

en la ciencia y en la ingeniería para obtener valores aproximados de
integrales definidas que no pueden calcularse analíticamente.

Pero… QUÉ ES INTEGRAR?

De acuerdo con la definición del diccionario,

integrar significa “llevar junto, como partes, en un
todo, unir, indicar la cantidad total…”
 INTRODUCCION

Matemáticamente la integración se representa por:

Ec 1

que se tiene para la integral de la función 𝑓(𝑥) con

respecto a la variable independiente 𝑥 , evaluada en los
límites 𝑥 = 𝑎 y 𝑥 = 𝑏
Como lo sugiere la definición del diccionario, el
“significado” de la Ec 1 es el valor total o sumatoria de f(x)
sobre el rango 𝑥 = 𝑎 a 𝑥 = 𝑏.
De hecho, el símbolo es una letra S estilizada que intenta
representar la conexión cercana entre la integración y la
sumatoria.
 INTRODUCCION

Observe que el proceso representado en la Ec 1 y en la Fig 1 es

llamado integración definida
 INTRODUCCION

La integra por tanto se puede representar como el límite de las sumas

de Riemann
INTRODUCCION
 INTRODUCCION

METODOS de INTEGRACIÓN NUMÉRICA

Regla
Trapezoidal
Fórmulas de
Integración Regla 1/3 de
De Simpson
Newton-Cotes Regla de
Simpson
Métodos Regla 3/8 de
de Simpson
Integración
Numérica

Integración Método de
De Extrapolación
Romberg De Richadson
 FORMULAS DE NEWTON-COTES

Las fórmulas de integración de Newton-Cotes son los esquemas de integración

numérica más comunes.

Se basan en la estrategia de reemplazar una función complicada o datos tabulados

con una función aproximada que sea fácil de integrar:

b b
I   f ( x)dx   f n ( x)dx Ec 2
a a

donde fn(x) es igual a un polinomio de la forma:

Ec 3

donde n es el orden del polinomio.

• Se llama grado de precisión de la fórmula de integración al máximo grado de los
polinomios que son integrados exactamente por dicha fórmula.
• La integración numérica es un proceso bastante más estable y preciso que la
derivación numérica vista previamente.
FORMULAS DE NEWTON-COTES
 FORMULAS DE NEWTON-COTES

Por ejemplo, en la Fig. 2 se usa el polinomio de primer orden (una línea recta) como
una aproximación. Mientras que en la Fig. 3 se emplea una parábola para el mismo
propósito.

Fig 2 Fig 3
 FORMULAS DE NEWTON-COTES

Por ejemplo, en la Fig. 4 se usan tres segmentos de línea recta para aproximar la
integral. Pueden utilizarse polinomios de orden superior para los mismos propósitos.

Fig 4
 LA REGLA
Base legal DEL TRAPECIO O TRAPEZOIDAL

La Regla Trapezoidal es la primera de las fórmulas de integración cerrada de Newton-

Cotes. Corresponde al caso donde el polinomio de la Ec 2. es de primer orden.
ℎ

 b  f (a )  f (b) 
)dx  I   f ( x)dx  (b  a )   Ec 2
 a
 2

Pero…
QUÉ SIGNIFICA LA REGLA TRAPEZOIDAL?

Geométricamente, la
regla del trapecio es
equivalente a aproximar
el área del trapecio bajo
la línea recta que conecta
a f(a) y f(b) como se
muestra en Fig. 7.

Para ver el trapecio

gira la cabera 90°
Fig. 7
 LA REGLA
Base legal DEL TRAPECIO O TRAPEZOIDAL

Recuerde, que la fórmula para calcular el área de un trapezoide es la altura por el

promedio de las bases, tal y como se muestra en la Fig. 8.

Fig. 8 Fig. 9

En la Fig. 8 se muestra la fórmula para calcular el área de un trapezoide (altura por el

promedio de las bases).
En la Fig. 9 para la regla trapezoidal el concepto es el mismo pero ahora el trapezoide
está sobre su lado
 EJERCICIOS DE APLICACIÓN

Use la Regla del Trapecio para aproximar los valores de las siguientes integrales:

a) b)
 APLICACIÓN MULTIPLE
Base
DE LAlegal
REGLA DEL TRAPECIO o COMPUESTA

La Regla del Trapecio se puede ampliar si subdividimos el intervalo [a,b] en n

subintervalos, todos de la misma longitud

Sea la partición que se forma al hacer dicha subdivisión.

Usando las propiedades de la integral, tenemos que:

Aplicando la Regla del Trapecio a cada una de las integrales, obtenemos:

 APLICACIÓN MULTIPLE
Base
DE LAlegal
REGLA DEL TRAPECIO

Ahora bien, ya que los subintervalos tienen la misma longitud h, tenemos que:

Sustituyendo el valor de h y haciendo uso de la notación sigma (sumatoria), tenemos

finalmente:

Ec. 6

Esta es la regla del trapecio para n subintervalos.

Obviamente, esperamos que entre más subintervalos usemos, mejor sea la

aproximación a la integral.
 APLICACIÓN MULTIPLE
Base
DE LAlegal
REGLA DEL TRAPECIO

Ilustración de la Regla Trapezoidal de aplicación múltiple: a) dos segmentos, b) tres

segmentos, c) cuatro segmentos y d) cinco segmentos

Fig. 10 Fig. 11
 APLICACIÓN MULTIPLE
Base
DE LAlegal
REGLA DEL TRAPECIO

Función en Matlab

Resolviendo la siguiente integral en Matlab

 EJERCICIOS DE APLICACIÓN

Use la Regla del Trapecio para aproximar el valor de la siguiente integral:

Si subdividimos en 5 intervalos
 PARA LA PRÓXIMA CLASE

ESTUDIAR LA REGLA DE SIMPSON

Regla 1/3 de
Simpson
Regla de
Simpson
Regla 3/8 de
Simpson
 Base
REGLAlegal
DE SIMPSON

Además de aplicar la Regla Trapezoidal con segmentación más fina, otra forma de
obtener una estimación más exacta de la integral es con el uso de polinomios de
orden superior para conectar los puntos.

Fig. 12 Fig. 13

Por ejemplo, si hay un punto extra a la Si hay dos puntos igualmente

mitad del camino f(a) y f(b), los tres espaciados entre f(a) y f(b), los cuatro
puntos se pueden conectar en una puntos se pueden conectar con un
parábola, tal y como se muestra en la Fig. polinomio de tercer orden, tal y como
12. se muestra en la Fig. 13.

Las fórmulas que resultan al tomar las integrales bajo estos polinomios son conocidas
como Regla de Simpson.
 Base
REGLAlegal
DE SIMPSON

La Regla de Simpson 1/3 resuelta cuando una interpolación polinomial de segundo

orden es sustituida en la ecuación:

Ec. 7

Si a y b se designan como xo y x2 y f2(x) es representada por un polinomio de

Lagrange de segundo orden y la integral se transforma en:

 ( x  x1 )( x  x2 ) ( x  x0 )( x  x2 ) ( x  x0 )( x  x1 ) 
x2

I   f ( x0 )  f ( x1 )  f ( x2 )dx
x0 
( x0  x1 )( x0  x2 ) ( x1  x0 )( x1  x2 ) ( x2  x0 )( x2  x1 ) 

Después de la integración y manejo algebraico, resulta la siguiente fórmula:

(b  a )
I   f ( x0 )  4 f ( x1 )  f ( x2 )  2  f ( x0 )  4 f ( x1 )  f ( x2 )
h
3 3
 Base
REGLAlegal
DE SIMPSON

Recuerde que x1 es el
punto medio entre a y b.
f ( xo )  4 f ( x1 )  f ( x2 )
I  (b  a ) Ec. 8
6

Esta ecuación es conocida como Regla de Simpson de 1/3. Es la segunda fórmula de

integración cerrada de Newton-Cotes.

La especificación “1/3” surge del hecho de que h está dividida entre 3 en la ecuación
anterior.
 EJERCICIOS DE APLICACIÓN

Use la Regla de Simpson de 1/3 para aproximar el valor de las siguientes integrales:

a) b)
 REGLA DE SIMPSON 1/3
Base legal
DE APLICACIÓN MÚLTIPLE

La Regla de Simpson se puede ampliar si subdividimos el intervalo [a,b] en n

subintervalos, todos de la misma longitud

Sea la partición que se forma al hacer dicha subdivisión y

sea Pm  xi 1 , xi  el conjunto de los untos medios de los subintervalos.

Usando las propiedades de la integral, tenemos que:

Al sustituir la Regla de Simpson de 1/3 a cada una de las integrales, obtenemos:

 REGLA DE SIMPSON
Base legal
DE APLICACIÓN MÚLTIPLE

Combinando términos y sustituyendo nos queda:

n n 1
f ( xo )  4 f ( xm )  2 f ( xi )  f ( x2 )
I  (b  a) i 1 i 1 Ec. 9
6n

Fig. 14: Representación

gráfica de la Regla de
Simpson 1/3 de aplicación
múltiple. Observe que el
método se puede emplear
sólo si el número de
segmentos es par
 REGLA DE SIMPSON
Base legal
DE APLICACIÓN MÚLTIPLE

Función en Matlab

Resolviendo la siguiente integral en Matlab

 EJERCICIOS DE APLICACIÓN

Use la Regla de Simpson de 1/3 para aproximar el valor de la siguiente integral y

sibdividiendo en 5 intervalos

a)

Use la Regla de Simpson de 1/3 para aproximar el valor de la siguiente integral y

sibdividiendo en 4 intervalos

b)
 Base
REGLAlegal
DE SIMPSON de 3/8

En una manera similar a la derivación de la Regla Trapezoidal y Regla de Simpson

1/3, un polinomio de Lagrange de tercer orden se puede ajustar a cuatro puntos e
integrarse:

Ec. 10

Para obtener:

Donde . Esta ecuación se llama Regla e Simpson de 3/8 debido a

que h se multiplica por 3/8.

NOTE QUE x1 Y x2 SON LOS PUNTOS QUE DIVIDEN EN TRES PARTES IGUALES
EL INTERVALO [a,b]
 Base
REGLAlegal
DE SIMPSON de 3/8

Ésta es la tercera fórmula de integración cerrada de Newton-Cotes. La Regla de

Simpson 3/8 se puede expresar también de la forma:

Ec. 11

Fig. 15: Ilustración de cómo se puede usar en conjunto

las Reglas de Simpson de 1/3 y 3/8 para menejar
aplicaciones múltiples con números nones de
intervalos.
 EJERCICIOS DE APLICACIÓN

Aproximar la siguiente integral usando la Regla de Simpson de 3/8:

a)
 Base
REGLAlegal
DE SIMPSON de 3/8 MÚLTIPLE

Al igual que en los casos anteriores, la Regla de Simpson de 3/8 se puede extender si
subdividimos el intervalo [a.b] en n intervalos de la misma longitud h.

Sea la partición determinada de esta forma. Cada subintervalo

lo dividimos en tres partes iguales, y sean y

los puntos determinados así:

Aplicado la Regla de Simpson de 3/8 en cada uno de los intervalos, tenemos:

ba  n  n 1

b

 f ( x)dx   f ( x o )  3   f ( y i )  f ( z i )   2  f ( x i )  f ( x n 
)
a
8n   i 1  i 1 
 EJERCICIOS DE APLICACIÓN

Aproximar la siguiente integral usando la Regla de Simpson de 3/8, subdividiendo en

3 intervalos:

a)
 RESUMEN DE FÓRMULAS

REGLA DEL TRAPECIO SIMPLE REGLA DEL TRAPECIO COMPUESTA

b f (a )  f (b)
dx  I   f ( x)dx  (b  a )
a 2

REGLA DE SIMPSON DE 1/3 SIMPLE REGLA DE SIMPSON DE 1/3 COMPUESTA

n n 1

I  (b  a)
f ( xo )  4 f ( x1 )  f ( x2 ) f ( xo )  4 f ( xm )  2 f ( xi )  f ( x2 )
6 I  (b  a) i 1 i 1
6n

REGLA DE SIMPSON DE 3/8 SIMPLE

REGLA DE SIMPSON DE 3/8 COMPUESTA

ba  n  n 1

b

 f ( x)dx   f ( x o )  3  f ( y i )  f ( z i 
)   2 f ( x i )  f ( x n 
)
a
8n   i 1  i 1 
 TALLER

CALCULAR EL VALOR DE LA INTEGRAL:

 8  3Senx  dx
0

HACIENDO USO DE:

1. REGLA DEL TRAPECIO SIMPLE
2. REGLA DEL TRAPACIO COMPUESTO EN n=3
3. REGLA DE SIMPSON DE 1/3 SIMPLE
4. REGLA DE SIMPSON DE 1/3 COMPUESTO CON n=3
5. REGLA DE SIMPSON DE 3/8 SIMPLE
Métodos adaptativos de
cuadratura
Métodos adaptativos • Las fórmulas conocidas como ecuaciones de
de cuadratura Newton-Cotes, se caracterizan porque la
estimación de la integral, se basa en valores
igualmente espaciados (h) para la función. Por lo
tanto, dichos valores son predeterminados o fijos.
(Chapra, S. y Canale, R. 2015)

Sin embargo, esta técnica no es adecuada cuando se

integra una función en un intervalo que contiene
regiones con variación funcional muy grande o muy
pequeña. (Burden, 2002)

Para este tipo de problemas, un método eficiente

a b podría estimar el grado de variación funcional y
adaptar el tamaño de paso a las diferentes
necesidades, a éstas técnicas se les conoce como
métodos adaptativos de cuadratura. (Burden, 2002).
Cuadratura Gaussiana

• Suponga que se elimina la restricción de los puntos fijos y que tiene la libertad
de evaluar el área bajo una línea recta que une dos puntos cualesquiera de la
curva. Si se eligen esos puntos adecuadamente, se obtendría una línea recta que
equilibraría los errores positivos y negativos, obteniendo una mejor
aproximación a la integral (Figura 3).

• A una de las técnicas que emplean éste principio, se le conoce como Cuadratura
de Gauss y se basa en las fórmulas de Gauss-Legendre.
 f(x f(x
) )

x x

a) Representación gráfica de la Regla del b) Estimación de la integral tomando el

Trapecio como el área bajo la línea recta área bajo la recta que pasa por dos puntos
que une los puntos extremos fijos intermedios. Estos puntos se ubican de tal
forma que se equilibra el error.
Figura 3. Representación gráfica de la integración de una función
• La cuadratura gaussiana selecciona los puntos de la evaluación de manera
óptima y no de forma igualmente espaciada, eligiendo los nodos
𝑥1 , 𝑥2 , ⋯ , 𝑥𝑛 en el intervalo 𝑎, 𝑏 y los coeficientes 𝑐1 , 𝑐2 , ⋯ , 𝑐𝑛 , para
reducir en lo posible el error esperado que se obtiene al efectuar la
aproximación:

𝑏 𝑛

𝐼= 𝑓 𝑥 𝑑𝑥 ≈ 𝑐𝑖 𝑓 𝑥𝑖
𝑎 𝑖=1
• El conjunto relacionado con este problema es el de los polinomios de
Legendre, es decir, un conjunto 𝑃0 𝑥 , 𝑃1 𝑥 , ⋯ , 𝑃𝑛 𝑥 , ⋯ , con las
siguientes propiedades:

1. Para cada 𝑛, 𝑃𝑛 𝑥 es un polinomio de grado 𝑛

1
2. −1
𝑃 𝑥 𝑃𝑛 𝑥 𝑑𝑥 = 0 siempre que 𝑃 𝑥 sea un polinomio de un grado
menor que 𝑛
Los primeros polinomios de Legendre son:
1
𝑃0 𝑥 = 1, 𝑃1 𝑥 = 𝑥, 𝑃2 𝑥 = 𝑥 2 −
3

3 6 3
𝑃3 𝑥 = 𝑥 3 − 𝑥 𝑦 𝑃4 𝑥 = 𝑥 4 − 𝑥 2 +
5 7 35

Las raíces de estos polinomios son diferentes: se encuentran en el intervalo

−1, 1 , tienen simetría con respecto del origen y lo más importante es que
son la opción correcta para determinar los parámetros que resuelven el
problema. (Burden, 2002). (Figura 4)
Representación t
gráfica de las
variables 𝑓(𝑥0 ) 𝑓(𝑥1 )
desconocidas 𝑥0 y
𝑥1 para la
integración por
medio de la
cuadratura de Gauss
x
-1 𝑥0 𝑥1 1

Figura 4. Representación de la Cuadratura de Gauss

Teorema

• Supongamos que 𝑥1 , 𝑥2 , ⋯ , 𝑥𝑛 son las raíces de los polinomios de Legendre

𝑃𝑛 𝑥 de 𝑛 − é𝑠𝑖𝑚𝑜 𝑔𝑟𝑎𝑑𝑜 y que para cada 𝑖 = 1, 2, ⋯ , 𝑛, los números 𝑐𝑖 están
definidos por:

1 𝑛 𝑥−𝑥𝑗
• 𝑐𝑖 = −1 𝑗=1 𝑥 −𝑥 𝑑𝑥 (3)
𝑖 𝑗
𝑗≠𝑖

• Si 𝑃 𝑥 es un polinomio cualquiera de un grado menor que 2𝑛, entonces:

1 𝑛
• −1
𝑃 𝑥 𝑑𝑥 ≈ 𝑖=1 𝑐𝑖 𝑃 𝑥𝑖
 • Las constantes 𝑐𝑖 necesarias para que la
cuadratura funcione, pueden ser generadas a
partir del teorema anterior. Sin embargo,
Raíces y coeficientes tanto las constantes como las raíces de los
de los polinomios de polinomios de Legendre han sido
Gauss-Legendre ampliamente tabuladas y, generalmente, se
encuentran para 𝑛 = 2, 3, 4 𝑦 5.

Tabla 1. Raíces y coeficientes de los

polinomios de Gauss-Legendre
 𝑏
• Una integral 𝑎
𝑓 𝑥 𝑑𝑥 en un intervalo arbitrario 𝑎, 𝑏 se puede transformar en
otra en −1, 1 (ver figura 5), usando un cambio de variables:

2𝑥−𝑎−𝑏 1
• 𝑡= ↔ 𝑥 =2 𝑏−𝑎 𝑡+𝑎+𝑏
𝑏−𝑎

𝑏−𝑎
• 𝑑𝑥 = 𝑑𝑡
2

• Esto permitirá aplicar cuadratura gaussiana a cualquier intervalo 𝑎 + 𝑏 ya que

𝑏 1 1 𝑏−𝑎 𝑡+ 𝑏+𝑎 𝑏−𝑎

• 𝑎
𝑓 𝑥 𝑑𝑥 = −1
𝑔 𝑡 𝑑𝑡 ≈ −1
𝑓 𝑑𝑡
2 2
 t 𝑏, 1
Cambio de
variables para 2𝑥 − 𝑎 − 𝑏
aplicar la 𝑡=
1 𝑏−𝑎
Cuadratura
Gaussiana

x
𝑎 𝑏

-1
𝑎, −1
Figura 5. Representación de la Cuadratura de Gauss
• Resumiendo, una integral definida puede aproximar su valor
aplicando el siguiente procedimiento:

𝑏 1 𝑛
• 𝑎
𝑓 𝑥 𝑑𝑥 = −1
𝑔 𝑡 𝑑𝑡 ≈ 𝑖=1 𝑐𝑖 ∙ 𝑔 𝑡𝑖

1 𝑏−𝑎
• Donde: 𝑥 = 𝑏 − 𝑎 𝑡 + 𝑎 + 𝑏 y 𝑑𝑥 = 𝑑𝑡
2 2

•𝑐𝑖 y 𝑡𝑖 son los coeficientes y raíces de los polinomios de Gauss-

Legendre
Ejemplo:
𝜋 2
Aproxime la integral 0
cos 𝑥 𝑑𝑥 usando cuadratura gaussiana con 𝑛 = 2

Solución:
1
1°. Se obtiene: 𝑥 = 2 𝑏 − 𝑎 𝑡 + 𝑎 + 𝑏
1
𝑥= 𝜋−0 𝑡+0+𝜋
2
1
𝑥= 𝜋𝑡 + 𝜋
2

𝑏−𝑎
2°. Se obtiene: 𝑑𝑥 = 𝑑𝑡
2

𝜋
∴ 𝑑𝑥 = 2 𝑑𝑡
3°. Los valores obtenidos se sustituyen en la función que se quiere integrar
para realizar el cambio de variable:

2 2
𝜋 2 1 1 𝜋 2 1 𝜋
0
cos 𝑥 = −1
𝑐𝑜𝑠 𝜋𝑡 + 𝜋 𝑑𝑡 ≈ 𝑖=1 𝑐𝑖 ∙ 𝑐𝑜𝑠 𝜋𝑡𝑖 + 𝜋
2 2 2 2

4°. Se desarrolla la sumatoria para 𝑛 = 2

2 2
1 𝜋
𝑐𝑖 ∙ 𝑐𝑜𝑠 𝜋𝑡𝑖 + 𝜋
2 2
𝑖=1
2 2
1 1 𝜋
= 𝑐1 𝑐𝑜𝑠 𝜋𝑡 + 𝜋 + 𝑐2 𝑐𝑜𝑠 𝜋𝑡2 + 𝜋
2 1 2 2
5°. Se sustituyen los valores de la tabla 1 en la sumatoria desarrollada en 𝑐𝑖 (coeficientes) y 𝑡𝑖
(raíces)

2 2 2 2
1 𝜋 1 1 𝜋
𝑐𝑖 ∙ 𝑐𝑜𝑠 𝜋𝑡𝑖 + 𝜋 = 𝑐1 𝑐𝑜𝑠 𝜋𝑡 + 𝜋 + 𝑐2 𝑐𝑜𝑠 𝜋𝑡2 + 𝜋
2 2 2 1 2 2
𝑖=1

2 2
1 1 𝜋
∴ 1 𝑐𝑜𝑠 𝜋 0.57735 + 𝜋 + 1 𝑐𝑜𝑠 𝜋 −0.57735 + 𝜋
2 2 2

𝜋
2
cos 𝑥 𝑑𝑥 = 0.6423172
0