UNIVERSIDAD CENTRAL DEL ECUADOR

FACULTAD DE INGENIERÍA QUÍMICA

CARRERA DE INGENIERÍA QUÍMICA

Tema: Regla de Simpson 1/3

NOMBRE:
• Ximena Salazar
• Melany Cruz
• Vanessa Pucuji
• Klaus Huilca
• Benancio Tipantaxi

GRUPO: 3

MATERIA: Análisis Numérico

DOCENTE: Ing. Fernando Vaca

Semestre: 2020-2021

Quito- Ecuador
 UNIVERSIDAD CENTRAL DEL ECUADOR
FACULTAD DE INGENIERIA QUIMICA
ANÁLISIS NUMÉRICO

Regla de Simpson 1/3

1. Marco Teórico
1.1. ¿Qué es la regla de Simpson?
Además de aplicar la regla del trapecio con una segmentación más fina, otra
forma de obtener una estimación más exacta de una integral consiste en usar
polinomios de grado superior para unir los puntos. Por ejemplo, si hay otro
punto a la mitad entre f(a) y f(b), los tres puntos se pueden unir con una parábola.
Si hay dos puntos igualmente espaciados entre f(a) y f(b), los cuatro puntos se
pueden unir mediante un polinomio de tercer grado. Las fórmulas que resultan
de tomar las integrales bajo esos polinomios se conocen como reglas de
Simpson. (Chapra, 2006)

1.2. Regla de Simpson 1/3

Este método consiste en la aproximación del cálculo del área plana bajo una
curva utilizando trapecios curvilíneos a partir de una interpolación con una
función cuadrática.
Esta aproximación es denominada “simple” debido a que utiliza tan sólo un
polinomio. Requiere el conocimiento de tres puntos equiespaciados: los
extremos y un punto central. Aplicando esta expresión utilizando mayor
cantidad de puntos intermedios (es decir, realizando más de un Simpson 1/3
dentro del intervalo) puede definirse la variante “compuesta” del método para
el cual se utilizan N puntos que corresponden a n = N − 1 subintervalos. Este
caso requiere que la cantidad de subintervalos sean pares (el caso simple utiliza
dos, por lo tanto, este debe ser un múltiplo).
A modo de interpretación, las sumatorias corresponden a los valores de la
función en los puntos impares y pares respectivamente sin contar los extremos.
Esto es debido a que, en los puntos impares, se encuentra el factor de 4 que
introduce la fórmula de Simpson 1/3 simple para los puntos intermedios y, en
los puntos pares, el factor de 2 se debe a que ese punto es compartido por los
trapecios curvilíneos adyacentes. (Nieves, 2014)
 UNIVERSIDAD CENTRAL DEL ECUADOR
FACULTAD DE INGENIERIA QUIMICA
ANÁLISIS NUMÉRICO

1.3. Descripción Gráfica de la regla de Simpson 1/3

En la gráfica a) Consiste en tomar el área bajo una parábola que une 3 puntos,
en cambio en b) consiste en tomar el aérea bajo una ecuación cubica que une
cuatro puntos. (Chapra, 2006)

Figura 1. Regla de Simpson 1/3

Fuente: Chapra, S. C. (2006). Métodos Numéricos para Ingenieros. Mexico:
The McGraw-Hill Companies, Inc.

1.4. Descripción grafica de la regla de Simpson 1/3 de aplicación múltiple

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.

Figura 2. Regla de Simpson 1/3 de aplicación múltiple

Fuente: Chapra, S. C. (2006). Métodos Numéricos para Ingenieros. Mexico:
The McGraw-Hill Companies, Inc.
 UNIVERSIDAD CENTRAL DEL ECUADOR
FACULTAD DE INGENIERIA QUIMICA
ANÁLISIS NUMÉRICO

2. Modelo Matemático
2.1. Regla de Simpson 1/3
La regla de Simpson 1/3 resulta cuando un polinomio de interpolación de segundo
grado se sustituye en la ecuación
𝑏 𝑏
𝐼 = ∫ 𝑓 (𝑥 )𝑑𝑥 ≅ ∫ 𝑓2 (𝑥 )𝑑𝑥 (1)
𝑎 𝑎

Si se designan a y b como 𝑥0 𝑦 𝑥2 , 𝑦 𝑓2 (𝑥 )se representa por un polinomio de

Lagrange de segundo grado, la integral de transforma en:
𝑥2 (
𝑥 − 𝑥1 )(𝑥 − 𝑥2 ) (𝑥 − 𝑥0 )(𝑥 − 𝑥2 )
𝐼=∫ [ 𝑓(𝑥0 ) + 𝑓(𝑥1 )
𝑥0 (𝑥0 − 𝑥1 )(𝑥0 − 𝑥2 ) (𝑥1 − 𝑥0 )(𝑥1 − 𝑥2 )
(𝑥 − 𝑥0 )(𝑥 − 𝑥1 )
+ 𝑓(𝑥2 )] 𝑑𝑥 (2)
(𝑥2 − 𝑥0 )(𝑥2 − 𝑥1 )
Después de la integración y de las manipulaciones algebraicas, se obtiene la siguiente
formula:
ℎ
𝐼≅ [𝑓(𝑥0 ) + 4𝑓(𝑥1 ) + 𝑓(𝑥2 )] (3)
3
Donde, en este caso, h = (b – a)/2. Esta ecuación se conoce como regla de
Simpson 1/3, y es la segunda fórmula de integración cerrada de Newton-Cotes. La
especificación “1/3” se origina del hecho de que h está dividida entre 3 en la ecuación
donde se integra el polinomio de Newton Gregory para llegar a la misma fórmula. La
regla de Simpson 1/3 también se puede expresar usando el formato de la ecuación

Donde 𝑎 = 𝑥0 , 𝑏 = 𝑥2 𝑦 𝑥1 = 𝑒𝑙 𝑝𝑢𝑛𝑡𝑜 𝑎 𝑙𝑎 𝑚𝑖𝑡𝑎𝑑 𝑒𝑛𝑡𝑟𝑒 𝑎 𝑦 𝑏, que esta dado

por (b+a)/2. Observe que, de acuerdo con la ecuación, el punto medio está ponderado
por dos tercios; y los dos puntos extremos, por un sexto. Se puede demostrar que la
aplicación a un solo segmento de la regla de Simpson 1/3 tiene un error de truncamiento
de:
1 5 ( 4)
𝐸𝑡 = − ℎ 𝑓 (𝜉 ) (4)
90
O como h=(b-a)/2,
 UNIVERSIDAD CENTRAL DEL ECUADOR
FACULTAD DE INGENIERIA QUIMICA
ANÁLISIS NUMÉRICO

(𝑏 − 𝑎 ) 5 ( 4)
𝐸𝑡 = − 𝑓 (𝜉 ) (5)
2880

Donde ξ está en algún lugar en el intervalo de a - b. Así, la regla de Simpson 1/3 es más
exacta que la regla del trapecio. No obstante, una comparación con la ecuación (21.6)
indica que es más exacta de lo esperado. En lugar de ser proporcional a la tercera
derivada, el error es proporcional a la cuarta derivada. En consecuencia, la regla de
Simpson 1/3 alcanza una precisión de tercer orden aun cuando se base en sólo tres
puntos. En otras palabras, ¡da resultados exactos para polinomios cúbicos aun cuando
se obtenga de una parábola!

2.2. Regla de Simpson 1/3 de aplicación múltiple

Así como en la regla del trapecio, la regla de Simpson se mejora al dividir el intervalo
de integración en varios segmentos de un mismo tamaño

𝑏−𝑎
ℎ= (6)
𝑛

La integral total se puede representar como

𝑥2 𝑥4 𝑥𝑛
𝐼 = ∫ 𝑓 (𝑥 )𝑑𝑥 + ∫ 𝑓 (𝑥 )𝑑𝑥 + ⋯ + ∫ 𝑓(𝑥 )𝑑𝑥 (7)
𝑥0 𝑥2 𝑥𝑛−2

Al sustituir la regla de Simpson 1/3 en cada integral se obtiene

𝑓(𝑥0 ) + 4𝑓 (𝑥1 ) + 𝑓(𝑥2 ) 𝑓(𝑥2 ) + 4𝑓 (𝑥3 ) + 𝑓(𝑥4 )

𝐼 ≅ 2ℎ + 2ℎ +⋯
6 6
𝑓(𝑥𝑛−2 ) + 4𝑓 (𝑥𝑛−1 ) + 𝑓(𝑥𝑛 )
+ 2ℎ (8)
6

o, combinando términos y usando la ecuación

 UNIVERSIDAD CENTRAL DEL ECUADOR
FACULTAD DE INGENIERIA QUIMICA
ANÁLISIS NUMÉRICO

Se debe utilizar un número par de segmentos para implementar el método.

Además, los coeficientes “4” y “2”) a primera vista parecerían peculiares. No obstante,
siguen en forma natural la regla de Simpson 1/3. Los puntos impares representan el
término medio en cada aplicación y, por lo tanto, llevan el peso de 4 Los puntos pares
son comunes a aplicaciones adyacentes y, por lo tanto, se cuentan dos veces. Un error
estimado en la regla de Simpson de aplicación múltiple se obtiene de la misma forma
que en la regla del trapecio: sumando los errores individuales de los segmentos y
sacando el promedio de la derivada para llegar a

3. Obtención y estimación del error de la regla de Simpson 1/3

la regla de Simpson 1/3 se obtiene al integrar el polinomio de interpolación de Newton-
Gregory hacia adelante
𝑥2
Δ2 𝑓(𝑥0 ) Δ3 𝑓(𝑥0 )
𝑰 = ∫ [𝑓(𝑥0 ) + Δ𝑓(𝑥0 )𝛼 + ( )
𝛼 𝛼−1 + 𝛼 (𝛼 − 1)(𝛼 − 2)
𝑥0 2 6
(f)4 𝜉
+ 𝛼(𝛼 − 1)(𝛼 − 2)(𝛼 − 3)ℎ4 ] 𝑑𝑥 (9)
24
Observe que se escribió el polinomio hasta el término de cuarto grado, en lugar
de hasta el de tercer grado como se esperaría. La razón de esto se verá un poco después.
Advierta también que los límites de integración van de 𝑥0 𝑎 𝑥2 , Por lo tanto, cuando se
realizan las sustituciones para simplificarla integral es de α=0 a 2:
𝑥2
Δ2 𝑓(𝑥0 ) Δ3 𝑓(𝑥0 )
𝐼 = ℎ ∫ [𝑓(𝑥0 ) + Δ𝑓(𝑥0 )𝛼 + 𝛼(𝛼 − 1) + 𝛼 (𝛼 − 1)(𝛼 − 2)
𝑥0 2 6
(f)4 𝜉
+ 𝛼(𝛼 − 1)(𝛼 − 2)(𝛼 − 3)ℎ4 ] 𝑑𝑥 (10)
24
Que al integrarse tiene:
 UNIVERSIDAD CENTRAL DEL ECUADOR
FACULTAD DE INGENIERIA QUIMICA
ANÁLISIS NUMÉRICO

𝛼2 𝛼2 𝛼2 𝛼4 𝛼3 𝛼2 3
𝐼 = ℎ [𝛼𝑓 (𝑥0 ) + Δ𝑓(𝑥0 ) + ( − ) Δ2 𝑓(𝑥0 ) + ( − + ) Δ 𝑓 (𝑥0 )
2 6 4 24 6 6
2
𝛼5 𝛼 4 11𝛼 3 𝛼 2 (4)
+( − + − ) 𝑓 ( 𝜉 )ℎ 4 ] (11)
120 16 72 8 0

Y evaluando en los límites se obtiene

Δ2 𝑓(𝑥0 ) 1 ( 4)
𝐼 = ℎ [2𝑓(𝑥0 ) + 2Δ𝑓(𝑥0 )𝛼 + + (0)Δ3 𝑓 (𝑥0 ) − 𝑓 (𝜉)ℎ4 ] (12)
3 90
Observe el resultado significativo de que el coeficiente de la tercera diferencia dividida
es cero. Debido a que Δ𝑓 (𝑥0 ) = 𝑓 (𝑥1 ) − 𝑓 (𝑥0 ) 𝑦 Δ2 𝑓 (𝑥0 ) = 𝑓 (𝑥2 ) − 2𝑓(𝑥1 ) + 𝑓(𝑥0 )

Así, el primer término es la regla de Simpson 1/3 y el segundo es el error de

truncamiento. Puesto que se suprime la tercera diferencia dividida, se obtiene el
resultado significativo de que la fórmula tiene una precisión de tercer orden. (Chapra,
2006)

4. Ejercicio
Ejemplo 1
Método de Simpson 1/3
Para ello tenemos en cuenta la formula cuando se quiere integrar [x0, xi]
2
𝐹(𝑥) = con 𝑋0 = 0 𝑋2 = 0.35
𝑥 2 −4
UNIVERSIDAD CENTRAL DEL ECUADOR
FACULTAD DE INGENIERIA QUIMICA
ANÁLISIS NUMÉRICO

Primero calculamos h
b−a
h=
2
X2 − X0
h=
2
0.35 − 0
h=
2
ℎ = 0.175
Entonces también conocemos que
h = X1 − X 0 = X 2 − X1
h − X 0 = X1
X1 = 0.175
Entonces resolviendo F(x)
2 2 1
F(x0 ) = = =−
x02 −4 0−4 2
2 2
𝐹(𝑥1 ) = = = −0.5039
𝑥12 − 4 0.1752 − 4
2 2
𝐹(𝑥2 ) = = = −0.5158
𝑥22 − 4 0.352 − 4
Aplicando el Método de Simpson
𝑋2
ℎ
∫ 𝐹(𝑥) 𝑑𝑥 = ∗ [𝐹(𝑥0 ) + 4𝐹( 𝑥1 ) + 2𝐹(𝑥2 )]
𝑋0 3
 UNIVERSIDAD CENTRAL DEL ECUADOR
FACULTAD DE INGENIERIA QUIMICA
ANÁLISIS NUMÉRICO

0.35
2
∫ 𝑑𝑥 =
0 𝑥2 −4
ℎ 0.175 1
∗ [𝐹(𝑥0 ) + 4𝐹(𝑥1 ) + 𝐹(𝑥2 )] = ∗ [− 2 + 4 ∗ (−0.5039) + (−0.5158)]
3 3
0.35
2
∫ 𝑑𝑥 = −0.1768
0 𝑥2 −4

Ejemplo 2
Aplicar La regla de Simpson para aplicación múltiple para:
𝑓(𝑥) = 0.2 + 25𝑥 − 200𝑥 2 + 675𝑥 3 − 900𝑥 4 + 400𝑥 5
Con a=0 hasta b=0.8 teniendo en cuenta que la integral exacta es 1.640533 y con (h=0.2:
n=4)
𝑓(0) = 0.2 + 25(0) − 200(0)2 + 675(0)3 − 900( 0)4 + 400(0)5
𝑓 (0) = 0.2

x F(x)
0 0.2
0.2 1.288
0.4 2.456
0.6 3.464
0.8 0.232

Aplicando la formula del método de Simpson

𝑓(𝑥0 ) + 4 ∑𝑛−1 𝑛−2
𝑖=1.3.5 𝑓(𝑥𝑖 ) + 2 ∑𝑖=2.4.6 𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑛 )
𝐼 = (𝑏 − 𝑎) ∗
3∗𝑛
0.2 + 4 ∗ (1.288 + 3.464) + 2 ∗ (2.456) + 0.232
𝐼 = (0.8 − 0) ∗
12
𝐼 = 1.623467
Cálculo del error estimado
(𝑏 − 𝑎)5
𝐸𝑎 = − ∗ 𝑓 ̅(4)
180 ∗ 𝑛4
𝑓 ̅(4) = 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑜 𝑑𝑒 𝑙𝑎 𝑐𝑢𝑎𝑟𝑡𝑎 𝑑𝑒𝑟𝑖𝑣𝑎𝑑𝑎
Entonces calculamos la 4 derivada de la función
𝑓 𝑖 = 2000𝑥 4 − 3600𝑥 3 + 2025𝑥 2 − 400𝑥 + 25
𝑓 𝑖𝑖 = 8000𝑥 3 − 10800𝑥 2 + 4050𝑥 − 400
 UNIVERSIDAD CENTRAL DEL ECUADOR
FACULTAD DE INGENIERIA QUIMICA
ANÁLISIS NUMÉRICO

𝑓 𝑖𝑖𝑖 = 24000𝑥 2 − 21600𝑥 + 4050

𝑓 𝑖𝑣 = 48000𝑥 − 21600
𝑓(0.4) = −2400
(0.8)5
𝐸𝑎 = − ∗ −2400
180 ∗ 44
𝐸𝑎 = 0.017067
𝐸𝑖 = 1.640533 − 1.640533 = 0.017067

5. APLICACIONES
La regla de Simpson 1/3 es muy útil utilizar parábolas para ajustar una función
curva obtendrás un resultado más exacto (con el mismo número de subintervalos). Con
algunos ejemplos, se puede observar que las parábolas solapadas se ajustan mejor a la
curva de la función que el segmento lineal y, por lo tanto, el área bajo la curva de una
mejor manera. Las parábolas hacen un mejor trabajo que los trapezoides a la hora de la
concavidad de la curva de la función.

6. CONSIDERACIONES FINALES

• El método de Simpson 1/3 solo se puede aplicar a segmentos equidistantes

• Es posible aplicar Simpson 1/3 a grupos de segmentos equidistantes y el aporte del resto
se calcula con la regla del Trapecio

• La aproximación se mejora al dividir el intervalo en varios segmentos

 UNIVERSIDAD CENTRAL DEL ECUADOR
FACULTAD DE INGENIERIA QUIMICA
ANÁLISIS NUMÉRICO

Bibliografía
• Chapra, S. C. (2006). Métodos Numéricos para Ingenieros. Mexico: The McGraw-
Hill Companies, Inc. .
• Nieves, A. (2014). Métodos Numéricos. Mexico: EBOOK.