UNIVERSIDAD TECNOLGICA

DE LOS ANDES
FILIAL ANDAHUAYLAS
Escuela profesional de ingeniera ambiental
Semestre 2017-1

MTODOS NUMRICOS II
Regla de Simpson
Ing. Max Fredy Quispe Aguilar
Presentado por:
Ivoshka k. Ramos Becerra
Jhordan Chuhi Lopez
Vladdy Gabriel Salazar
Arteaga
Wilfredo Altamirano
morales
 INTRODUCCIN
En integracin numrica, una forma de aproximar una integral definida en
un intervalo [a,b] es mediante la regla del trapecio, es decir, que sobre cada
subintervalo en el que se divide [a,b] se aproxima f por un polinomio de primer
grado, para luego calcular la integral como suma de las reas de los trapecios
formados en esos subintervalos . El mtodo utilizado para la regla de Simpson
sigue la misma filosofa, pero aproximando los subintervalos de f mediante
polinomios de segundo grado
Uno de los problemas matemticos ms frecuentes es el clculo del rea que se
forma al graficar una funcin. Por ejemplo, se necesita calcular el rea A que
aparece en la siguiente figura:

En donde la funcin f(x) y los valores a y b son conocidos.

En este tipo de problemas se pueden obtener dos tipos de soluciones:
Soluciones algebraicas: se obtiene una frmula precisa y exacta para el rea
solicitada.
Soluciones numricas: se calcula numricamente una estimacin del rea.
Desde luego, la soluciones algebraicas son mejores que las numricas, porque
son exactas. Pero a veces, la complejidad de las funciones hace imposible (o
difcil) obtener la solucin algebraica, por lo que una solucin numrica permite
ahorrar tiempo.
En integracin numrica, una forma de aproximar una integral definida en
un intervalo [a,b] es mediante la regla del trapecio, es decir, que sobre cada
subintervalo en el que se divide [a,b] se aproxima f por un polinomio de primer
grado, para luego calcular la integral como suma de las reas de los trapecios
formados en esos subintervalos . El mtodo utilizado para la regla de Simpson
sigue la misma filosofa, pero aproximando los subintervalos de f mediante
polinomios de segundo grado.
1.-OBJETIVOS
1. Comprender los rasgos generales de la integracin aproximada utilizando el mtodo

de Simpson.

2. Comprender la aproximacin del error por truncamiento de la integracin aproximada

utilizando el mtodo de Simpson, frente al valor exacto. 4. Resolver problemas de

integracin aproximada utilizando el mtodo de Simpson.

2.-METODO DE SIMPSON
Adems de aplicar la regla trapezoidal o Rectangular con segmentos o sub reas

cada vez ms pequeas, otra manera de obtener una estimacin an ms exacta

de una integral, es la de usar polinomios de orden superior para conectar los

puntos, en el caso particular del mtodo que usa orden 2, es decir de la forma

ax + bx + c 2 .
A las frmulas resultantes de calcular la integral bajo estos polinomios se les

conoce como reglas de Simpson.

En este procedimiento, se toma el intervalo de anchura 2h, comprendido entre

Xi y xi+2, y se sustituye la funcin f(x) por la parbola que pasa por tres puntos

(xi, yi), (xi+1, yi+1), y (xi+2, yi+2). El valor del rea aproximada, sombreada en la

figura, se calcula con un poco ms de trabajo y el resultado es

(yi+4yi+1+yi+2), que se demuestra en seguida.

DESARROLLO DEL MODELO DE SIMPSON:

Para efectos de la demostracin del mtodo de Simpson, se asume cada sub

rea como un pequeo arco de parbola de la forma ax2 + bx + c con lmites as:

Limite inferior en h, limite superior en h, por ende la mitad de la pequea sub

rea se encontrar en el Punto 0, tal como se ilustra en Fig. 2.

Se procede a integrar dicho arco de parbola entre los lmites descritos se

tendr:

Reemplazando cada uno de los

lmites, se tiene:

, ahora destruyendo parntesis se

tendr:

, simplificando un poco la solucin

se
obtendr la ecuacin 1 que se muestra a continuacin.
 .Ec 1

Observando la Fig 3, en lo que respecta a las notaciones, se puede decir que ,

, Entonces se podra

obtener el siguiente sistemas de ecuaciones, evaluando la ecuacin general de

la parbola ax2 + bx + c en cada uno de los puntos de la pequea sub rea

[h,0-h]:

De lo anterior se puede decir que:

y0+y2=2ah2+2c. Ec 2

y1=c. Ec 3

Retomando la Ec 1 se puede expresar igualmente de la siguiente manera:

.. Ec 4
Reemplazando las ecuaciones 2 y 3 en la Ec 4 se tiene que:

.Ec 5

Interpretando la ecuacin Ec 5 con base en la sub rea seleccionada A1 para

desarrollar el modelo de Simpson, se dira que el rea del segmento es igual a

la suma de la altura o funcin evaluada en el lado izquierdo mas cuatro veces la

funcin evaluada en la parte central de la sub rea mas la funcin evaluada en

el lado derecho de la sub rea, todo esto multiplicado por el ancho del sub rea

y dividido por 3.

La simple inspeccin visual de esta figura y la que describe el procedimiento de

los trapecios o los rectngulos, confirma que el mtodo de Simpson deber ser

mucho ms exacto que los procedimientos mencionados.

Si a y b se denominan como x0 y x2,y fi (xi) se representa mediante un polinomio

de Lagrange de segundo orden, entonces la integral es:

Despus de integrar y de reordenar los trminos, resulta la siguiente ecuacin:

.Ec 5a

Si se toma entonces se

tiene como solucin de la sub rea , que sera lo mismo

mostrado en la ecuacin 5.
Ahora, se sabe que el rea que se desea encontrar sera la sumatoria de todas

las sub reas que se calculen. Al igual que los mtodos de la regla trapezoidal y

de la regla rectangular, entre ms sub reas tenga la integral a calcular, ms

exacto ser el valor encontrado. El rea aproximada en el intervalo [a, b] es:

, ahora dejando esta ecuacin en trminos de

la ecuacin 5 se tendr:

Simplificando h /3 y sumando los trminos se tendr:

donde n sera el nmero de sub reas en el cual se ha dividido el rea que se

desea calcular.

A manera de ejemplo, si el rea a calcular se hubiera dividido en 4 Sub reas

entonces en trminos de y la solucin sera:

Bien, dependiendo como se agrupen los trminos se llegara a expresar la

solucin de dos maneras:

. Ec 6

. Ec 7

Los primeros trminos del parntesis, contienen los valores de la evaluacin de

la funcin en los extremos, el segundo, la suma de los trminos de ndice impar,

y el tercero la suma de los trminos de ndice par.

Las dos ecuaciones se pudieran representar en trminos de sumatorias de la

siguiente manera.

La Ec 6:

.. Ec 8

La Ec 7:

.. Ec 9

Para efectos de programacin y en lo que respecta a mi concepto personal, es

mejor la solucin representada como Ec 9 y con ella se contina el trabajo. Hay

que tener en cuenta que n es el nmero de sub reas en la que se divide el rea

total a calcular y h = dx / 2 .

Ahora lo que se conoce en un momento determinado, cuando se desea calcular

el valor de la integral definida, son los siguientes trminos:

a = Lmite inferior

b = Lmite Superior

n = Nmero de sub reas

f (x) La funcin sobre la cual se desea integrar.

Con los valores anteriores se pude calcular el valor de dx as: dx = (b a)/ n

y h = dx / 2 .

En necesario entonces dejar la ecuacin en trminos de f (x), a, b y dx h as:

Los primeros trminos: y0 = f(a) y y2n = f (b)

Analizando ahora los trminos impares y1 = f(a+1dx/2), y3 = f(a+3dx/2),

y5=f(a+5dx/2), por tanto se tendra de manera general:

Y2-i = f(a+(2i-1)dx/2) y2i-1 = f(a+(2i-1)h)..Ec 10

Analizando ahora los trminos pares: y2 =f(a+1dx), y4 = f(a+2dx), y5 = f(a+3dx) , por

tanto se tendra de manera general:

Y2i = f(a+idx) Y2i = f(a+2ih) ..Ec 11

Dejando la ecuacin Ec 9 en trminos de lo expresado en las ecuaciones Ec 10

y Ec 11 se tendr en forma definitiva la solucin as:

...Ec 12

EJERCICIOS:
1.- Calcula la integral , usando la siguiente tabla de datos:

Solucin combinando varios mtodos de integracin.

Para este caso, se puede aplicar la regla de Simpson de 1/3 en el intervalo [-

1,0], la regla del trapecio en el intervalo [0,1] y la regla de Simpson de 3/8 en el

intervalo [1,3.25]. As, se tiene las siguientes integrales:

Por lo tanto, la integral buscada es la suma de las tres integrales anteriores:

2.-Utilizar la regla de Simpson para aproximar la integral:

Tenga en cuenta que el valor real es 1.4626

Solucin: Usando la frmula directamente con los siguientes datos:

Si se asume el rea a calcular como un solo arco de parbola, se tendra

entonces que

dx = (b a)/1=1 h = (dx / 2) = 0.5

EJERCICIO APLICATIVO:

Un agricultor desea conocer la superficie aproximada de un prado limitado por

una carretera, dos caminos perpendiculares a ella y la ribera de un ro, de

manera que si colocamos unos ejes cartesianos sobre la carretera (eje OX) y

uno de los caminos (eje OY, abscisa x = 0), el segundo camino ser a la recta
vertical x = 2 (unidades en cientos de metros). Se toman varias medidas desde

la carretera hasta la ribera, obtenindose las siguientes coordenadas para los

puntos de la ribera: (0, 1.5), (0.5, 1.8), (1, 2.1), (1.5, 1.75), (2, 1.3).

Calcular aproximadamente el rea de dicho terreno utilizando las reglas de los

trapecios y de Simpson. Determinar el rea si extendemos el terreno hasta la

abscisa x = 2.5 sabiendo que el ro en tal caso pasa por el punto (2.5, 1.1).

En este caso desconocemos la funcin de forma explcita, teniendo en cuenta

tan solo los valores de la tabla que nos han sido facilitados. Se tiene:

de modo que usando el mtodo de los trapecios podemos escribir

mientras que si usamos el mtodo de Simpson se llega a

Si se aade un nuevo punto, la tabla queda dada por

Ahora la regla de los trapecios proporcionar:

mientras que si el mtodo de Simpson no es aplicable de forma directa dado que

estamos considerando un nmero impar de subintervalos en este caso. Lo que

podemos hacer es considerar la regla de Simpson para los 4 subintervalos

primeros y estimar el quinto subintervalo mediante la regla de los trapecios.

As queda:

BIBLIOGRAFIA:
http://www.ilustrados.com/documentos/eb-
simpson.pdf
https://www.easycalculation.com/es/integration/si
mpson-rule.php
http://campus.usal.es/~mpg/Personales/PersonalM
AGL/Docencia/CalculoCATema5cTeoria(09-10).pdf
 https://es.wikipedia.org/wiki/Regla_de_Simpson
https://www.youtube.com/watch?v=qWpulFYFGEc