UNIVERSIDAD NACIONAL DE Jaén

Escuela Profesional de Ingeniería Civil

Ingeniería - UNJ
Ingeniería Civil
y más….
CURSO : MÉTODOS NUMÉRICOS

ALUMNOS:
 AGREDA CRUZ, Milcer Duberli
 CUNAIQUE ALBERCA, Alicia Stefanny
 DAVILA SARMINENTO, Julio Fernando
 GUERRERO NUÑEZ, Jhon Carlos.
 JIMENEZ TORO, Karen
 OJEDA SANCHEZ, Pedro Delmys
 ROSORIO BOCANEGRA, Stefano
 ROMERO RUIZ, Olandy
 VILLEGAS LIZANA, Eishten
DOCENTE :
Mg. Mat. NUÑEZ RIVAS, CARLOS

JAÉN, 24 septiembre 2018

ÍNDICE
1
I. INTRODUCCIÓN 3
II. FÓRMULAS DE NEWTON-COTES 4
2.1 REGLAS DE SIMPSON 7
2.1.1 REGLA DE SIMPSON 1/ 3 7
2.1.2 REGLA DE SIMPSON 3/8 11
2.1.3 EJEMPLO 13
2.2 REGLA DEL TRAPECIO 14
2.2.1 TRAPECIO SIMPLE (PARA UN TRAMO) 15
2.2.2 TRAPECIO COMPUESTA (FÓRMULA GENERAL PARA “n” TRAMOS) 15
2.2.3 ALGORITMO DE LA REGLA DEL TRAPECIO 16
2.2.4 EJEMPLO 17
2.3 REGLA DE BOOLE 18

2
I. INTRODUCCIÓN

En este trabajo trataremos sobre uno de los problemas más vastos de la

aproximación numérica la solución de ecuaciones no lineales analizado
de diferentes maneras desde la óptica analítica y su interpretación
geométrica.
En el campo de la tecnología principalmente en la ingeniería nos
encontramos generalmente con el siguiente problema determinar las
raíces de la ecuación f(x) = 0.
Como la teoría de la difracción de la luz se precisa de la siguiente
ecuación; x-tanx=0 . así como también, el método de bisección se puede
dividir el intervalo de x i a x u en mitades iguales, no se toman en
consideración las magnitudes de f (xi ) y f (xu ) . Por ejemplo, si
f ( xi ) está mucho más cercana a cero que f (xu ) , es lógico que la
raíz se encuentre más cerca de x i que de f (xu ) .

3
II. FÓRMULAS DE NEWTON-COTES

Las fórmulas de Newton-Cotes son los tipos de integración numérica más comunes.
Se basan en la estrategia de reemplazar una función complicada o datos tabulados por
un polinomio de aproximación que es fácil de integrar:

b b
I =∫ f ( x ) dx ≅ ∫ f n ( x ) dx ……… (1)
a a

donde f n ( x ) = un polinomio de la forma

f n ( x )=a0 +a1 x+ ⋯+a n−1 x n−1+ an x n

donde n es el grado del polinomio. Por ejemplo, en la figura 1.a, se utiliza un

polinomio de primer grado (una línea recta) como una aproximación. En la figura 1.b,
se emplea una parábola con el mismo propósito.

La integral también se puede aproximar usando un conjunto de polinomios aplicados

por pedazos a la función o datos, sobre segmentos de longitud constante. Por ejemplo,
en la figura 2, se usan tres segmentos de línea recta para aproximar la integral.

Imagen 1: La aproximación de una integral mediante el área bajo a) una sola línea recta y b) una parábola.

4
 Aunque pueden utilizarse polinomios de grado superior con los mismos propósitos.
Con este antecedente, reconocemos que el “método de barras” emplea un conjunto de
polinomios de grado cero (es decir, constantes) para aproximar la integral.

Existen formas cerradas y abiertas de las fórmulas de Newton-Cotes. Las formas

cerradas son aquellas donde se conocen los datos al inicio y al final de los límites de
integración (figura 3.a). Las formas abiertas tienen límites de integración que se
extienden más allá del intervalo de los datos (figura 3.b). en este sentido, son similares a
la extrapolación. Por lo general, las formas abiertas de Newton-Cotes no se usan para
integración definida. Sin embargo, se utilizan para evaluar integrales impropias y para
obtener la solución de ecuaciones diferenciales ordinarias. Este capítulo enfatiza las
formas cerradas. No obstante, al final del mismo se presenta brevemente una
introducción a las fórmulas abiertas de Newton-Cotes.

5
 Imagen 2: La aproximación de una integral mediante el área
bajo tres segmentos de línea recta.

Imagen 3: La diferencia entre las fórmulas de integración a) cerradas y b) abiertas.

II.1 REGLAS DE SIMPSON

6
 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 (figura 4.a). Si hay
dos puntos igualmente espaciados entre f (a) y f (b) , los cuatro puntos se pueden
unir mediante un polinomio de tercer grado (figura 4.b). Las fórmulas que resultan de
tomar las integrales bajo esos polinomios se conocen como reglas de Simpson.

Imagen 4: a) Descripción gráfica de la regla de Simpson 1/3, que consiste en tomar el área bajo una
parábola que une tres puntos. b) Descripción gráfica de la regla de Simpson 3/8, que consiste en tomar
el área bajo una ecuación cúbica que une cuatro puntos.

II.1.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 (1):

b b
I =∫ f ( x ) dx ≅ ∫ f 2 ( x ) dx
a a

Si se designan a y b como x0 y x2 , y f 2(x ) se representa por un

polinomio de Lagrange de segundo grado, la integral se transforma en

x2

I =∫
x0 [ ( x−x 1 )( x −x2 )
( x 0−x 1 )( x 0 −x2 )
f ( x0)+
( x−x 0 )( x −x2 )
( x1 −x 0 )( x 1−x 2 )
f ( x 1) +
( x−x 0 ) ( x−x 1 )
( x 2−x 0 ) ( x 2−x 1 ) ]
f ( x 2 ) dx

Después de la integración y de las manipulaciones algebraicas, se obtiene la

siguiente fórmula:

7
 h
3[ 0
I≅ f ( x ) +4 f ( x 1 ) +f ( x2 )] ……… (2)

donde, en este caso, h=(b−a)/2 . Esta ecuación se conoce como la 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 (2). Una
alternativa para obtenerla se muestra en el título 2.1.1.1, 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 I ≅ ( b−a ) × altura promedio

……… (3)

donde a=x 0 , b=x2 y x 1=¿ el punto a la mitad entre a y b , que está

dado por (b+ a)/2 . Observe que, de acuerdo con la ecuación (3), el punto medio está
ponderado por dos tercios; y los 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 (título 2.1.1.1):

−1 5 4
Et = h f (ξ)
90

o, como h=(b−a)/2 ,

−(b−a)5 4
Et = f (ξ ) ……… (4)
2880

donde ξ está en algún lugar en el intervalo de a 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 de error de la regla del trapecio 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 basa en sólo tres puntos. En otras palabras, da resultados exactos para
polinomios cúbicos aun cuando se obtenga de una parábola.

8
 II.1.1.1 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.

[ ]
x2
∆2 f ( x 0 ) ∆ 3 f ( x 0) f 4 (ξ )
I =∫ f ( x0 )+ ∆ f ( x0 ) α + α ( α−1 )+ α ( α−1 )( α −2 ) + α ( α−1 ) ( α −2 ) ( α −3 ) h4 dx
x0 2 6 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 x0 a x 2 . Por lo tanto,
cuando se realizan las sustituciones para simplificar, la integral es de α=0 a 2 :

[ ∆ 2 f ( x0 ) ∆ 3 f ( x0 )
]
2
f 4 (ξ )
I =h∫ f ( x0 )+ ∆ f ( x0 ) α + α ( α−1 )+ α ( α−1 ) ( α −2 ) + α ( α −1 ) ( α −2 ) ( α −3 ) h4
0 2 6 24

que al integrarse tiene

[ α2
I =h αf ( x0 ) + ∆ f ( x 0 )+
2
α2 α2 2
−
6 4 (
∆ f ( x0 )+
α 4 α3 α2 3
− +
24 6 6) ∆ f ( x0 )+( α 5 α 4 11 α 3 α 2 4
− +
120 16 72
−
8) f (ξ) h4 ( ) ]
0

y evaluando en los límites se obtiene

[
I =h 2 f ( x 0 ) +2 ∆ f ( x 0 ) +
∆ 2 f ( x0 )
3
1
+ ( 0 ) ∆3 f ( x 0 )− f 4 (ξ) h4
90 ]

Observe el resultado significativo de que el coeficiente de la tercera diferencia

dividida es cero. Debido a que ∆ f ( x0 ) =f ( x 1 )−f ( x 0 ) y

∆ 2 f ( x 0 ) =f ( x 2 )−2 f ( x 1 ) + f ( x 0) , la ecuación de reescribe como

9
 II.1.1.2 LA 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 (Imagen 5):

b−a
h= ……… (5)
n

La integral total se puede representar como

x2 x4 xn

I =∫ f ( x ) dx+∫ f ( x ) dx +…+ ∫ f ( x ) dx
x0 x2 x n−2

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

f ( x 0 ) +4 f ( x 1 ) +f ( x 2 ) f ( x 2 ) +4 f ( x 3 ) +f ( x 4 ) f ( x n−2 ) + 4 f ( x n−1 ) +f ( x n )
I ≅2h + 2h + …+ 2 h
6 6 6
o, combinando términos y usando la ecuación (5)

……… (6)

Observe que, como se ilustra en la figura (5), se debe utilizar un número par de
segmentos para implementar el método. Además, los coeficientes “4” y “2” en la
ecuación (6) 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 en la ecuación (3). 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

(b−a)5 (4 )
Ea = f́ ……… (7)
180 n4

10
donde f́ (4 ) es el promedio de la cuarta derivada en el intervalo.

Imagen 5: Representación gráfica de la regla de Simpson 1/3 de aplicación múltiple.

Observe que le método se puede emplear sólo si el número de segmentos es par.
II.1.2 REGLA DE SIMPSON 3/8
De manera similar a la obtención de la regla de Simpson 1/3, es posible ajustar un
polinomio de Lagrange de tercer grado a cuatro puntos e integrarlo:

b b
I =∫ f ( x ) dx ≅ ∫ f 3 ( x ) dx
a a

para obtener

3h
8 [
I≅ f ( x 0 ) +3 f ( x 1 ) +3 f ( x 2 ) +f ( x 3 ) ]

donde h=(b−a)/3 . Esta ecuación se llama regla de Simpson 3/8 debido a que h
se multiplica por 3/8. La regla 3/8 se expresa también en la forma:

……… (8)

Así los dos puntos interiores tienen pesos de tres octavos que los puntos extremos
tienen un peso de un octavo. La regla de Simpson 3/8 tiene un error de

11
 −3 5 (4)
Et = h f (ξ )
80

o, como h=(b−a)/3 ,

−(b−a)5 (4)
Et = f ( ξ ) ……… (9)
6480

Puesto que el denominador de la ecuación (9) es mayor que la ecuación (4), la regla
es más exacta que a regla 1/3.

Por lo común, se prefiere la regla de Simpson 1/3, ya que alcanza una exactitud de
tercer orden con tres puntos en lugar de los cuatro puntos requeridos en la versión 3/8.
No obstante, la regla de 3/8 es útil cuando el número de segmentos es impar.
Supongamos que se desea una estimación de 5 segmentos. Una opción podría ser
utilizar una versión de la regla del trapecio de aplicación múltiple, quizá esto no sea
recomendable, sim embargo, debido al gran error de truncamiento asociado con dicho
método. Una alternativa sería aplicar la regla de Simpson 1/3 a los dos primeros
segmentos y la regla de Simpson 3/8 a los últimos tres (Imagen 6). De esta forma,
podríamos obtener un estimado con una exactitud de tercer orden durante todo el
intervalo

Imagen 3: Ilustración de cómo se utilizan en

conjunto las reglas de Simpson 1/3 y 3/8
para manejar aplicaciones múltiples con
números impares de intervalos.

II.1.3 EJEMPLO

12
 1. calcular la siguiente integral por el método de simpsom
2

∫ √ 1+ x 2 dx donde n = 4 Nota: n siempre tiene

0
Solución que ser par para poder
resolver por este
a) Hallamos nuestro h método

b−a 2−0 2
h= = = =0.5
n 4 4

b) Hallamos imagen de ; xi

x 0=a ; x 2=x 0 +h ; x 3=x 2+2 x h; x 4 =x 3+ h …. ; x n=b

xi f ( xi ¿
x 0=¿ 0 1
x 1=¿ 0.5 1.118
x 2=¿ 1 1.4142
x 3=¿ 1.5 1.8027
x 4=¿ 2 2.2360

c) Reemplazamos los valores de la tabla en ecuación:

∫ f (x) x ≡ b−a
3n
[ f ( x 0 ) + 4 f ( x 1 ) +2 f ( x 2 ) + 4 f ( x 3 ) + f ( x 4 ) ]
a

x 0=a ; x 2=x 0 +h , x3 =x2 +2 x h ; x 4 =x3 + h …. ; x n=b

Entonces quedará

∫ √ 1+ x 2 dx= 2−0
3x4
[f ( 0 ) +4 f ( 0.5 ) +2 f ( 1 ) +4 f ( 1.5 ) + f (2 )]
0

2
¿ [1+ 4 ( 1.118 )+2 ( 1.4142 ) +4 ( 1.8027 ) + ( 2.2360 ) ]
12

= 3.004533

13
 II.2 REGLA DEL TRAPECIO

El método del trapecio es un método que consiste en reemplazar la gráfica de

una función f ( x) por la de un polinomio “polinomio aproximante” el cual va a
coincidir con la función original f(x) en un tramo [a,b]. Entonces se tendría que integrar
el polinomio aproximante desde el punto (a, f(a)) hasta (b, f(b)). (Ver figura).

La función f(x) (en azul) es aproximada por la función

lineal (en rojo)

II.2.1 TRAPECIO SIMPLE (PARA UN TRAMO)

La integral de ésta es igual al área del trapecio bajo la gráfica de la función lineal. Se
sigue que:
b
f ( a )+ f (b)
I= ∫ f (x)dx ≡( b−a) 2
a

Cuyo error se calcula por la siguiente fórmula:

14
 f}(ξ
3 ´
donde f } left (ξ right ) = {int from {a} to {b} {f (x)dx ¿
−( b−a ) ´ ;
¿ b−a
12

Siendo ξ un numero entre a y b

II.2.2 TRAPECIO COMPUESTA (FÓRMULA GENERAL PARA “n”

TRAMOS)

Representación gráfica de la
regla de trapecio compuesta

La regla del trapecio compuesta o regla de los trapecios es una forma de aproximar una
integral definida utilizando n trapecios. En la formulación de este método se supone que f es
b
continua y positiva en el intervalo [a, b]. De tal modo la integral definida ∫ f (x)dx
a
representa el área de la región delimitada por la gráfica de f y el eje x, desde x=a hasta
x=b . Primero se divide el intervalo [a,b] en n subintervalos, cada uno de ancho
b−a
Δx =h= ; n es el número de divisiones.
n
Después de realizar todo el proceso matemático se llega a la siguiente fórmula:
h
f ( x)dx ≡ (¿ f ( a+0 xh ) +2 f ( a+ 1 x h ) +2 f ( a+2 h ) +… …+ f ( b))
2
b

∫¿
a

x 0=f ( a+0 xh ) ; x2=a+1 x h , … … ; x i=a+ixh

Cuyo error se calcula por la siguiente formula:

15
 f } (ξ
3 ´
donde f } left (ξ right ) = {int from {a} to {b} {f (x) dx ¿
−(b−a) ´ ;
2
¿ b−a
12n

II.2.3 ALGORITMO DE LA REGLA DEL TRAPECIO

i. Defino mi integral ∫ f ( x)dx , declaro mi a , b , y el número de trapecios =n

a
b−a
ii.Calculo el H= ,
n

iii. Calculo , x i=a+ iH ; f ( x i ) ,

iv. Remplazo mis datos en la formula siguiente

[ ]
b n−1

∫ f ( x ) dx= h2 f ( x 0 ) +2 ∑ f ( x i ) + f ( x n )
a i=1

v. Finalmente obtengo mi integración.

II.2.4 EJEMPLO

1. Calcular al valor aproximado de la integral definida, por el método de los

trapecios para el valor de n indicado.
3

∫ √ 1+ x 2 dx ; n= 6
2

Solución:
d) Hallamos nuestro h

b−a 3−2 1
h= = = =0.1 6
n 6 6

16
 e) Hallamos imagen de ; xi
x i=a+ixh

xi f ( xi ¿
x 0=¿ 2 2.2360
x 1=¿ 2.16 2.380
x 2=¿ 2.32 2.526
x 3=¿ 2.48 2.674
x 4=¿ 2.64 2.823
x 5=¿ 2.80 2.9773
x2.96
6=¿ 3.0435

f) Reemplazamos los valores de la tabla en ecuación:

h
f ( x) x ≡ (¿ f ( a+0 xh ) +2 f ( a+1 x h )+ 2 f ( a+2 h ) +… …+ f (b))
2
b

∫¿
a

x 0=f ( a+0 xh ) ; x2=a+1 x h , x 3=a+2 x h

Entonces quedará

∫ √ 1+ x 2 dx= 0.16
2
[f ( x 0 ) +2 f ( x 1 ) +2 f ( x 2 ) +2 f ( x 3 ) +2 f ( x 4 ) +2 f ( x 5 ) + f ( x6 ) ]
2

0.16
¿ [2.236+2 ( 2.380 ) +2 ( 2.526 ) +2 ( 2.674 ) +2 ( 2.823 ) +2 ( 2.9773 ) + ( 3.0435 ) ]
2

= 2.81825

II.3 REGLA DE BOOLE

La regla de Boole (llamada así debido a George Boole) utiliza cinco puntos
consecutivos igualmente separados para calcular la integral aproximada de la función
utilizando un polinomio de cuarto grado, considerando cinco valores de la
misma equiespaciados en [a ,b] ..

17
 2h
f ( x)dx ≡ (¿ 7 f ( x 0 ) +3 2 f ( x 1 ) +12 f ( x 3 ) +7 f (x 4 ))
45
b

∫¿
a

Donde

b−a
h= ; f i =f (a+ixh)
4

Y el error es

−(b−a)7 (6)
f (ξ )
1935360

Siendo ξ un numero entre a y b

18