37

Guía de estudio Matemática V

TEMA 4
INTEGRACION NUMERICA

En ingeniería se presenta con frecuencia la necesidad de integrar una función con

alguna de las siguientes características:
(a) Una función complicada y continua que es difícil o imposible de integrar
directamente.
(b) Una función tabulada en donde los valores de x y f(x) se dan en un conjunto
de puntos discretos, como es el caso a menudo, de datos experimentales.
En estos dos casos, se deben emplear métodos aproximados. Las fórmulas de
integración de Newton-Cotes son los esquemas más comunes dentro de la
integración numérica. Se basan en la estrategia de reemplazar una función
complicada o un conjunto de datos tabulares con alguna función aproximada que
sea más fácil de integrar. La integral se puede aproximar usando una serie de
polinomios aplicados por partes a la función o a los datos sobre intervalos de
longitud constante.
Se dispone de las formas abierta y cerrada de las fórmulas de Newton-Cotes. Las
formas cerradas son aquellas en donde los puntos al principio y al final de los
límites de integración se conocen. Las fórmulas abiertas tienen los límites de
integración extendidos más allá del rango de los datos. Las fórmulas abiertas de
Newton-Cotes, en general, no se usan en la integración definida. Sin embargo, se
usan extensamente para evaluar integrales impropias y en la solución de
ecuaciones diferenciales ordinarias.

4.1. REGLA DEL TRAPECIO (SIMPLE Y COMPUESTA)

La regla del trapecio o regla trapezoidal es la primera de las fórmulas cerradas de
Newton-Cotes.

4.1.1. Regla del trapecio simple:

Considérese la función f(x), cuya gráfica esta entre los extremos x=a y x=b como
se muestra en la figura. Si utilizamos un polinomio P(x) de primer grado como una
aproximación de f(x) tenemos:
xb xa
P( x)  f ( a )  f (b) , el cual es equivalente a:
a b ba
f (b)  f (a)
P( x)  f (a )  ( x  a)
ba
El área bajo esta línea recta será una aproximación del área bajo la curva entre los
límites a y b
 38
Guía de estudio Matemática V

Integrando este polinomio:

f(x) b b f (b)  f (a) 
f(a)
a f ( x)dx  a  f (a)  b  a ( x  a) dx
b
f (b)  f (a) ( x  a) 2
 f (a) x 
ba 2 a
f(b)
f (b)  f (a) (b  a) 2
 f (a)(b  a) 
ba 2
(b  a)
 f (a)(b  a)  ( f (b)  f (a))
a b 2
 f (b)  f (a) 
 (b  a) f (a)  
 2 
f (a)  f (b)
 (b  a)
2
Que es la conocida Regla del Trapecio Simple. Geométricamente, la Regla del
trapecio aproxima el área bajo una curva mediante el área del trapecio bajo la
línea recta que une f(a) y f(b)
Bb
b A h
2
f(a) f (a)  f (b)
h A (b  a)
2
y f(b)
B
b-a

Ejemplo:
1
Utilizar la regla del trapecio simple para aproximar la integral

1 / 2 arcsen( x)dx
Solución:
1

1 / 2 arcsen( x)dx  (1  1/ 2) f (1) 2f (1/ 2)
 /2 /6 
   0.5235988
x
    
4 
6  
La solución exacta de esta integral es:
 5  6 3
 0.4429715
12





 39
Guía de estudio Matemática V

El error relativo porcentual que se cometió al aplicar la regla del trapecio simple
esta dado por:
p  p* 0.4429715  0.5235988
Er =  100   100  18.2%
p 0.4429715
El error de la estimación es muy alto.

4.1.2. Regla del trapecio compuesta:

Una aproximación suficiente al área bajo la curva se obtiene dividiéndola en n
segmentos de ancho h  (b  a) / n y aproximando el área de cada segmento
mediante un trapecio, como se indica en la figura:
Sea P  {x0 , x1 , xn } la partición que se
f(x) forma al hacer dicha subdivisión. Usando
f(a)
propiedades de la integral tenemos que:

b x1 x2
a f ( x)dx  x 0
f ( x)dx  
x1
f ( x)dx  
f(b) xn
 f ( x)dx
xn 1
Aplicando la regla del trapecio en cada una
de las integrales, obtenemos:
a b f ( x0 )  f ( x1 ) f ( x1 )  f ( x2 )
h h 
2 2
f ( xn1 )  f ( xn )
h
2
Agrupando términos:
n1
h h 
( f ( x0 )  2 f ( x1 )  2 f ( x2 )    f ( xn ))   f ( x0 )  2 f ( xi )  f ( xn ) 
b
a f ( x)dx 
2 2 i 1 
n 1
f ( x 0 )  2  f ( xi )  f ( x n )
b i 1
a f ( x)dx  (b  a) 2n
Que es la conocida como la Regla del Trapecio múltiple o compuesta.

Ejemplo:
Utilizar la regla del trapecio compuesta con n=5 subintervalos para aproximar la
1
integral 1 / 2 arcsen( x)dx
 y
40
Guía de estudio Matemática V

Solución:
b  a 1  1/ 2
 h   0.1
n 5
P= {0.5, 0.6, 0.7, 0.8, 0.9, 1}

1
1 / 2 arcsen( x)dx 
0.1
2
( f (0.5)  2 f (0.6)  2 f (0.7) 
x
     2 f (0.8)  2 f (0.9)  f (1)) =0.4513161

 La solución exacta de esta integral es:

5  6 3
 0.4429715
12

El error relativo porcentual que se cometió al aplicar la regla del trapecio
compuesta esta dado por:
p  
p* 0.4429715  0.4513161
Er =  100   100  1.884%
p 0.4429715
Vemos que  efectivamente se ha obtenido una mejor aproximación con este
método, reduciendo el error relativo verdadero de un 18.2% hasta un 1.884%.
Si aumentamos el valor de n obtendremos los siguientes resultados:
n Snumérica Er%
10 0.4460196420 0.688%
50 0.4432559383 0.0642%
100 0.4430730772 0.0229%
200 0.4430076838 0.00816$
250 0.4429974465 0.00585%
1000 0.4429747968 0.000736%

Se ha obtenido una mejor aproximación al aumentar el número de subdivisiones.

1. Aplique la regla del trapecio compuesta con los valores indicados de n para
aproximar las siguientes integrales:
2 2
 x2
a)  e dx , n=5 b)  x ln xdx , n=4
2 1
 41
Guía de estudio Matemática V

2 
  x 2 cos( x)dx , n=6
3 x
c) x e dx , n=4 d)
2 0
2 5
2 1
e)  x2  4
dx , n=8 f)  dx , n=8
0 3 x 4
2

2. Utilizar la regla del trapecio compuesta para aproximar el trabajo W realizado

por una fuerza f ( x)  e  x  x 2  3x  1 , en el intervalo 0  x  3 . Utilice n=1, 2, 3
y 5 subintervalos. Resuelva analíticamente y determine el error relativo
porcentual de la aproximación en cada caso.

3 La velocidad de un paracaidista que cae esta dada por:

gm   t
c
v 1 e m 
c  
 
a) Aplicando la regla del trapecio simple, aproxime la distancia recorrida al
cabo de 9s del lanzamiento si la masa m es de 68,1 kg y su velocidad v es de
35 m/s. Tome g=9,8 m/s2 y el coeficiente de arrastre c=15 kg/s. [Sugerencia:
dx
v ]
dt
b) Utilice la regla del trapecio compuesta con n=2, 3 y 5 subintervalos.

4 En estadística, la probabilidad de que un valor aleatoriamente seleccionado

descrito por la Distribución Normal se encuentre en [a, b] está dada por:
2
1 x 
b   
1
 e 2   con media   0 y desviación estándar  .
a  2
Utilice la regla del trapecio compuesta con n=2, 3 y 5 subintervalos para
aproximar la probabilidad de que un valor aleatoriamente seleccionado descrito
por la Distribución Normal se encuentre en:
a)   ,   b)  2 , 2  c)  3 , 3 

5. Utilice la regla del trapecio para aproximar el desplazamiento de un móvil si la

velocidad en cada intervalo de 6 seg se muestra en la siguiente tabla:

t (s) 0 6 12 18 24 30
v (m/s) 124 134 148 156 147 133

4.2. REGLAS DE SIMPSON

 42
Guía de estudio Matemática V

A través de la Regla de Simpson (nombrada así en honor de Thomas Simpson) se

puede obtener una estimación más exacta de la integral. El método consiste en
usar polinomios de grado superior para aproximar la curva de la función y tomar
las integrales bajo tales polinomios.

4.2.1. Regla de Simpson 1/3 simple:

Considérese la función integrando f(x), cuya gráfica está entre los extremos x0 = a
y x2 = b, si hay otro punto a la mitad x1 = (x0+ x2)/2 como se muestra en la figura.
Si utilizamos un polinomio de lagrange P2(x) de segundo grado como una
aproximación de f(x):
( x  x1 ) ( x  x2 ) ( x  x0 ) ( x  x 2 ) ( x  x0 ) ( x  x1 )
P2 ( x)  f ( x0 )   f ( x1 )   f ( x2 ) 
( x0  x1 ) ( x0  x2 ) ( x1  x0 ) ( x1  x2 ) ( x2  x0 ) ( x2  x1 )

El área bajo este polinomio será una aproximación del área bajo la curva entre los
límites a y b
f(x) Integrando este polinomio:
b x2  ( x  x1 ) ( x  x2 )
a f ( x)dx  x0  f ( x0 ) ( x0  x1 )  ( x0  x2 )  f ( x1 )
f(x2) ( x  x0 ) ( x  x 2 ) ( x  x0 ) ( x  x1 ) 
  f ( x2 ) 
( x2  x0 ) ( x2  x1 ) 
f(x0) dx
P(x) ( x1  x0 ) ( x1  x 2 )

Después de la integración y manipulación

algebraicas, se obtiene la siguiente formula:
x0 x1 x2 b f ( x0 )  4 f ( x1 )  f ( x2 )
a f ( x)dx  ( x2  x1 ) 6
b h
a f ( x)dx  3 ( f ( x0 )  4 f ( x1 )  f ( x2 ))
Que es la conocida Regla de Simpson 1/3 Simple donde h = (b − a)/2.
Geométricamente, la Regla de Simpson 1/3 simple aproxima el área bajo una
curva mediante el área bajo una parábola que une tres puntos.

Ejemplo:
Utilizar la regla de Simpson 1/3 simple para aproximar la integral
1
1 / 2 arcsen( x)dx
Solución:
 y
43
Guía de estudio Matemática V


5  6 3
La solución exacta de esta integral es:  0.4429715
12

1
1 / 2 arcsen( x)dx  (1 61/ 2) ( f (1)  4 f (3 / 4)  f (1/ 2))
 1
 (1.5707963  0.8480621  0.5235988)  0.4572203
12
El error relativo porcentual
x
que se cometió al aplicar
    laregla del
 trapecio
 simple esta dado por:

p  p*

Er =  100
p

0.4429715  0.4572203
 100  3.22%
 0.4429715
El error de la estimación es menor que el obtenido
con la Regla del Trapecio simple.


4.2.2. Regla de Simpson 1/3 compuesta:


En el caso de que el intervalo [a, b] no sea lo suficientemente pequeño, el error al
calcular la integral puede ser muy grande. Para ello, se recurre a la fórmula
compuesta de Simpson 1/3. Dividiremos el intervalo [a, b] en n subintervalos
iguales, de manera que xi = a + ih, donde h = (b − a) / n para i = 0,1,...,n.
f(x) b x2 x4
a f ( x)dx  x 0
f ( x)dx   f ( x)dx  
x2
xn
 f ( x)dx
xn  2
f(b)
f(a) Aplicando la regla de Simpson 1/3 en cada
una de las integrales, obtenemos:
h h
 ( f ( x0 )  4 f ( x1 )  f ( x2 ))  h ( f ( x2 )  4 f ( x3 )
3 3
h
 f ( x4 ))    ( f ( xn2 )  4 f ( xn1 )  f ( xn ))
a b 3

Agrupando términos:
n / 21
h  n/2 
 
b

a f ( x)dx 
3 
f ( x 0 )  2 f ( x 2i )  4 f ( x 2i 1 )  f ( x n ) 
i 1 i 1 
 44
Guía de estudio Matemática V

n / 21 n/2

b
f ( x0 )  2  f ( x2i )  4  f ( x2i 1 )  f ( xn )
i 1 i 1
a f ( x)dx  (b  a) 3n
Que es la conocida como la Regla del Trapecio múltiple o compuesta. Se debe
utilizar un número par de divisiones para implementar el método.

Ejemplo:
y
Utilizar la regla de Simpson 1/3 compuesta con n=4 subintervalos para aproximar
1
la integral 1 / 2 arcsen( x)dx


Solución:
b  a 1  1/ 2
 h   0.125
n 4
P= {0.5, 0.625, 0.75, 0.875, 1}

1
1 / 2 arcsen( x)dx
0.125 x
     ( f (0.5) 4 f (0.625 )  2 f (0.75) 
3
 4 f (0.875)  f (1)) =0.4480329


La solución exacta de esta integral es:


5  6 3
 0.4429715
12
El error relativo
 porcentual que se cometió al aplicar la regla de Simpson 1/3
compuesta esta dado por:
p  p* 0.4429715  0.4480329
Er =   100   100  1.14%
p 0.4429715
Vemos que efectivamente se ha obtenido una mejor aproximación con este
método, reduciendo el error relativo verdadero de un 3.22% hasta un 1.14%.

Si aumentamos el valor de n obtendremos los siguientes resultados:

n Snumérica Er%
10 0.4442541593 0.2895500498
50 0.4430863307 0.02591491573
 45
Guía de estudio Matemática V

100 0.4430121240 0.009162891245

200 0.4429858860 0.003239711554
250 0.4429818040 0.002318207647
1000 .4429728180 0.0002896348633

Se ha obtenido una mejor aproximación al aumentar el número de subdivisiones.

1. Resuelva aplicando la regla de Simpson los ejercicios de la ACTIVIDAD No.

11.

2. El trabajo producido por un proceso termodinámico a temperatura constante

esta dado por dW = P(v)dv. A continuación están tabulados los resultados
obtenidos experimentalmente:
v (m3) 2 3 4 5 6 7
P (kPa) 294.4 266.4 260.8 249.6 193.6 165.6
Utilice la regla de Simpson para estimar el trabajo en kJ.

3. El estudio de la difracción de la luz en una apertura rectangular implica el uso

de las integrales de Fresnel
t t
 
c(t )   cos w dw y s(t )   sen
2
w 2 dw
0
2 0
2
Aplique la regla de Simpson compuesta para aproximar c(1) y s(1).

4. La fuerza total del viento ejercida sobre un mástil de un bote de

vela de carreras se expresa como la integral de una función
contínua:
L 2
z Lz
F 200 e dw
0
5 z
Si la longitud del mástil es de 30 pies Aplique la regla de
Simpson compuesta para aproximar F, utilizando n=1 y 5.

5. La raíz media cuadrática de la intensidad de una corriente alterna esta dada

por:
T
1 2 t
I RMC  
T 0
i 2 (t )dt , donde i(t )  sen
T
, T=1 s

Aplique la regla de Simpson compuesta para aproximar IRMC, utilizando n=1 y 5.