ÍNDICE

1. Unidad 3: Raíces de Ecuaciones e Integración Numérica .....................................3

Tema 2: Integración Numérica ..................................................................................... 3

Objetivo: ........................................................................................................................ 3

Introducción: ................................................................................................................. 3

2. Información de los subtemas .............................................................................6

2.1 Subtema 1: Método del Trapecio ....................................................................... 6

2.2 Subtema 2: Método de Simpson 1/3................................................................ 11

2.3 Subtema 3: Método de Simpson 3/8................................................................ 14

3. Preguntas de Comprensión de la Unidad ..........................................................21

4. Material Complementario ................................................................................23

5. Bibliografía ......................................................................................................24

2
 Integración Numérica

1. Unidad 3: Raíces de Ecuaciones e

Integracion Numerica
Tema 2: Integración Numérica
Objetivo:
Obtener la integral de una función continua en un intervalo [a, b] o de puntos x y f(x)
provenientes de mediciones o experimentos, aplicando fórmulas de integración
numérica como Trapecio, Simpson 1/3 y Simpson 3/8.

Introducción:
Las integrales son muy útiles en el área de ingeniería, por ejemplo, se pueden calcular
longitudes, áreas y volúmenes de distintas superficies, aunque sean irregulares. Otras
aplicaciones son el cálculo de centro de gravedad y valores eficaz rms de voltaje o
corriente. Hasta ahora se han visto la solución de integrales definidas y no definidas de
forma analítica. Sin embargo, existen problemas donde la función f(x) resulta
complicada integrar analíticamente o en otros casos cuando los datos son obtenidos
mediante mediciones y experimentos, en estos casos, no se podría obtener una
solución analítica. Para resolver los casos mencionados anteriormente, se utilizan
fórmulas de integración numérica, usualmente denominadas fórmulas de Newton
Cotes como las de Trapecio, Simpson 1/3 y Simpson 3/8, las cuales se basan en la
división en partes iguales de un intervalo cerrado [a, b], se va sumando el área que se
forma en cada una de las divisiones y se obtiene una valor, estas fórmulas no se
utilizan para obtener soluciones analíticas, puesto que convergen a un valor y
requieren necesariamente de un intervalo [a, b].
La función a integrar usualmente es una de las siguientes formas: (Chapra & Canale,
2006)
© Universidad Estatal de Milagro – UNEMI

1. Una función simple y continua en el intervalo [a, b], puede ser polinomial,
exponencial o trigonométrica.
2. Otro caso es una función trascendental continua, que de forma analítica es
difícil o imposible de integrar porque no existe una función antiderivada.
3. Una función con valores tabulados de t y f(t), proporcionados como un
conjunto discreto de puntos, esto cuando se trata de datos experimentales o
de campo.
La integración de Romberg y cuadratura de Gauss se utilizan para resolver integrales
en intervalos abiertos en los que se incluyen extremos de + o - infinito, o en casos

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 3

 Integración Numérica

donde exista una indeterminación, otra aplicación es en problemas con datos

irregularmente espaciados.
También en necesario analizar problemas de integración donde los datos están
irregularmente espaciados, esto es muy común en aplicaciones reales.

Integración Numérica

La forma general de una integral definida se muestra en la ecuación 1

b
I =  f ( x)dx (1)
a

donde f ( x) es una función continua en [a, b] que se desea integrar. La estrategia de

las fórmulas de Newton Cotes es reemplazar f(x) que puede ser una función
complicada o mediante datos x y f(x) por un polinomio, el cual es más sencillo de
integrar (Chapra & Canale, 2006).

b xn

I   pn ( x)dx =  p ( x)dx
n (2)
a x0

Donde pn(x) es un polinomio de la forma:

pn ( x) = a0 + a1 x + + an −1 x n −1 + an x n

El polinomio pn puede aproximarse mediante interpolación o mínimos cuadrados

(Sauer, 2012). La integral también se puede aproximar mediante un conjunto de
polinomios hallados con trazadores spline.
© Universidad Estatal de Milagro – UNEMI

Figura 1: Integración polinomio de interpolación

Fuente: (Nieves & Dominguez, 2014)

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 4

 Integración Numérica

Si la aproximación polinomial es buena, la integral de dicho polinomio dará una

aproximación excelente de la integral de la función f(x).

Las fórmulas de Newton Cotes se aplican cuando los valores de x son equidistantes, en
cambio para datos de x espaciados de forma irregular se utilizan las fórmulas de
cuadratura Guassiana.

Fórmulas de Newton Cotes

Para aplicar este tipo de fórmulas se divide el intervalo [a, b] en n partes iguales, los
valores extremos de cada intervalo se calculan con:

xi = a + ih, i = 0,1, 2,..., n (3)

donde x0=a y xn=b, n es el número de divisiones del intervalo, el valor de h se obtiene

con la siguiente expresión:

b−a
h=
n (4)

La función f(x) se aproxima por un polinomio pn(x) de grado n, el cual se integra para
obtener la aproximación de I (Nieves & Dominguez, 2014).
La integral puede ser reescrita como la suma de los sub-intervalos generados entre los
dos puntos generales de integración. De esta manera, la integral 1 puede ser
expresada como se muestra a continuación: (Lara, 2019)
b x1 x2 xn
I ( f ) =  f ( x)dx =  f ( x)dx +  f ( x)dx + ... +  f ( x)dx (5)
a x0 x1 xn−1

Para las fórmulas de Newton Cotes se puede utilizar el polinomio de diferencias finitas
o Newton, para determinar el error de la estimación se añade el error del polinomio de
interpolación, a continuación, se definen las formulas (Rodríguez, 2016):

x − x0
f ( x) = pn ( s ) + en ( s ), s=
© Universidad Estatal de Milagro – UNEMI

h
1 2 1 1 n
pn ( s ) = f 0 + f 0 s +  f 0 s ( s − 1) +  3 f 0 s ( s − 1)( s − 2) + +  f 0 s ( s − 1)( s − 2) ( s − n + 1)
2! 3! n!
 s  n +1 ( n +1)
en ( s ) =  h f ( z ), x0  x  xn
 n + 1

A partir de estas fórmulas se obtienen las de Trapecio, Simpson 1/3 y Simpson 3/8

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 5

 Integración Numérica

2. Informacion de los subtemas

2.1 Subtema 1: Método del Trapecio
Fórmula del Trapecio Simple

Se llama método del trapecio porque cuando f es una función con valores positivos de
b

 f ( x)dx se aproxima mediante el área de un trapecio:

a
(Burden &
Faires, 2011)
h h3
f ( x)dx =  f ( x0 ) + f ( x1 )  −
b
a 2 12
f ''( ) (6)

Figura 2: Método del trapecio para n=1

Fuente: (Burden & Faires, 2011)

Fórmula general del trapecio compuesto

© Universidad Estatal de Milagro – UNEMI

b
Si queremos aproximar la integral a
f ( x)dx primero dividimos el intervalo [a, b] en n
de segmentos trapezoides, luego aproximamos cada sub-intervalo de integración a un
segmento trapezoidal:

b x1 x2 xn
 a
f ( x)dx =  f ( x)dx +  f ( x)dx + ... + 
x0 x1 xn−1
f ( x)dx
(7)

f ( x0 ) + f ( x1 ) f ( x1 ) + f ( x2 ) f ( xn−1 ) + f ( xn )
b

 f ( x)dx  h
a
2
+h
2
+ ... + h
2 (8)

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 6

 Integración Numérica

Figura 3: Método del trapecio compuesto n intervalos

Fuente: (Burden & Faires, 2011)

Como h es igual para todos los sub-intervalos, la ecuación se puede escribir de forma
compacta como:

h
 f ( x0 ) + 2 f ( x1 ) + 2 f ( x1 ) + 2 f ( x2 ) + ... + 2 f ( xn−1 ) + f ( xn )
b
a
f ( x)dx 
2 (9)

La fórmula para el trapecio compuesta, incluyendo el término del error es:

h n −1
 b−a 2

b
a
f ( x)dx =
2 
f ( x0 ) + 2
i =1
f ( xi ) + f ( xn )  − 12 h f ''( )

(10)

Ejercicio 1

Integre numéricamente la función f(t) en el intervalo indicado, aplicando el método del

Trapecio con n=10, compare con el valor real de la solución analítica que es 9.2536.

Fuente: (Gutiérrez, Olmos, & Casillas, 2010)

3
I =  ( t 4 − 2t 3 + t 2e−t ) dt
© Universidad Estatal de Milagro – UNEMI

Solución

a=0 b=3 n = 10
b − a 3−0
h= = = 0.3
n 10
Colocamos los valores de t y f(t) en una tabla para luego aplicar la fórmula del Trapecio
compuesta.

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 7

 Integración Numérica

n ti f(ti)
0 0 0
1 0.3 0.02077364
2 0.6 -0.104827811
3 0.9 -0.472578576
4 1.2 -0.948680335
5 1.5 -1.18545714
6 1.8 -0.630831602
7 2.1 1.466132849
8 2.4 6.052135411
9 2.7 14.26802819
10 3 27.44808362
Tabla 1: Puntos ti y f(ti) de la función ejercicio 1
Fuente: Elaboración Propia

h n −1

f ( x)dx   f ( x0 ) + 2 f ( xi ) + f ( xn ) 
b
a 2 i =1 
h
 f ( x0 ) + 2 ( f ( x1 ) + f ( x2 ) + + f ( x9 ) ) + f ( x10 ) 
b

a f ( x)dx 
2
0.3   0.0207 − 0.10482 − 0.4725 − 0.9486 − 1.18545 − 0.6308  
I= 0 + 2   + 27.448
2   +1.4666 + 6.0521 + 14.2680  
I = 9.6566
Ahora se calcula el error relativo
I R − I Ap
eI = *100
IR
9.2536 − 9.6566
eI = *100
9.2536
eI = 4.35%
Ejercicio 2
Secrist y Hornbeck con el propósito de determinar las características térmicas de los
© Universidad Estatal de Milagro – UNEMI

frenos tuvieron que aproximar numéricamente la temperatura exterior promediada

del área T en el cojín del freno, basándose en la ecuación:

 T ( r ) r dr
r0
p
T=
re
r0
 r dr
re p

Donde r es el radio de inicio de contacto entre almohadilla-disco. El parámetro r0 es el

radio exterior entre el contacto almohadilla-disco,  p representa el ángulo de las
pastillas de freno y T(r) representa la temperatura para cada punto de la almohadilla,

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 8

 Integración Numérica

obtenida por medición. Si re = 0.308 ft, r0= 0.478 ft,  p = 0.7051 rad. Aproximado el
valor de la temperatura promedio.
Fuente: (Burden & Faires, 2011)

r ( ft ) T ( r )( F ) r ( ft ) T ( r )( F ) r ( ft ) T ( r )( F )
0.308 640 0.376 1034 0.444 1204
0.325 794 0.393 1064 0.461 1222
0.342 885 0.410 1114 0.478 1239
0.359 943 0.427 1152
Tabla 2: Mediciones de la temperatura en función del radio
Fuente: (Burden & Faires, 2011)

Figura 4: Freno de Disco

Fuente: (Burden & Faires, 2011)

Solución:
a = 0.308 b = 0.478
n = 10
h = 0.325 − 0.308 = 0.017
Se tabulan los valores de las funciones del numerador y denominador para resolver
individualmente cada integral.
n ri T r*teta T*r*teta
© Universidad Estatal de Milagro – UNEMI

0 0.308 640 0.2172 138.9893

1 0.325 794 0.2292 181.9511
2 0.342 885 0.2411 213.4126
3 0.359 943 0.2531 238.7024
4 0.376 1034 0.2651 274.1316
5 0.393 1064 0.2771 294.8390
6 0.41 1114 0.2891 322.0474
7 0.427 1152 0.3011 346.8415
8 0.444 1204 0.3131 376.9295
9 0.461 1222 0.3251 397.2124

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 9

 Integración Numérica

10 0.478 1239 0.3370 417.5898

Tabla 3: Calculo de valores de las funciones del numerador y denominador ejercicio 2.
Fuente: Elaboración propia.

r0

I1 =  T ( r ) r p dr
re

h n −1

f ( x)dx   f ( x0 ) + 2 f ( xi ) + f ( xn ) 
b
a 2 i =1 
 181.9551 + 213.4126 + 238.4126 + 238.7024 +  
0.017    
I1 = 138.9893 + 2  274.1316 + 294.8390 + 322.0474 + 346.8415 +  + 417.5898
2   376.9295 + 397.2124 
   

I1 = 49.7141

r0

I 2 =  r p dr
re

0.017   0.2292 + 0.2411 + 0.2531 + 0.2651 + 0.2771 +  

I2 = 0.2172 + 2   + 0.3370
2   0.2891 + 0.3011 + 0.3131 + 0.3251  
I 2 = 0.0471
I1
T=
I2
49.7141
T= F
0.0471
T = 1055.3272F
T = 568.515C
© Universidad Estatal de Milagro – UNEMI

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 10

 Integración Numérica

2.2 Subtema 2: Método de Simpson 1/3

Método de Simpson 1/3 para n=2

Una forma de obtener una estimación más exacta de una integral es utilizar un
polinomio de grado n para unir los puntos. El método de Simpson 1/3 se basa en
calcular el área de un polinomio cuadrático que pasa por los tres puntos.

Figura 4: Simpson 1/3 con n=2

Fuente: (Chapra & Canale, 2006)

Para obtener este polinomio se utiliza interpolación cuadrática mediante Lagrange o

Newton. Los coeficientes del polinomio cuadrático se obtienen a partir de tres puntos
igualmente espaciados en x. Si se cuenta con un dominio [a, b] los puntos que deben
a+b
x0 = a, x2 = b y x1 = .
ser usados son 2 . El polinomio puede ser escrito como
aparece en la expresión:

p( x) = a1 + a2 ( x − x0 ) + a3 ( x − x0 )( x − x1 )

Donde a0, a1 y a2 son constantes del polinomio y pueden ser calculadas como se
© Universidad Estatal de Milagro – UNEMI

muestra a continuación.

a1 = f ( x0 )

f ( x1 ) − f ( x0 )
a2 =
x1 − x0

f ( x2 ) − 2 f ( x1 ) + f ( x0 )
a3 =
2h 2

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 11

 Integración Numérica

b−a
Sea h = 2 . Sustituyendo el valor de las constantes se obtiene el siguiente arreglo
matemático.

h
 f ( x0 ) + 4 f ( x1 ) + f ( x2 )
x2 x2
x0
f ( x)dx  p( x)dx =
x0 3 (12)

La fórmula de Simpson 1/3 incluyendo el término del error para n=2 es:

h h5
 f ( x0 ) + 4 f ( x1 ) + f ( x2 ) − f (4) ( )
x2

x0
f ( x)dx =
3 90
(13)

Método Simpson 1/3 Compuesto

Para la formula compuesta de Simpson de divide la integral de f(x) en el intervalo [a, b]

en sub-intervalos iguales, cada uno con tres puntos regularmente espaciados, esta
regla se aplica en cada par consecutivo de puntos como se observa en figura 5, para
aplicar esta regla n debe ser par.
b x2 x4 xi+1 xn =b
 a
f ( x)dx =
x0 =a
f ( x)dx + f ( x)dx +... + 
x2 xi−1
f ( x)dx +... + 
xn−2
f ( x)dx (14)

Figura 5: Regla de Simpson 1/3 compuesta

Fuente: (Burden & Faires, 2011)
En forma general el método compuesto de Simpson 1/3 se expresa con la ecuación 15

h n −1 n−2  b − a 4 (4)
f ( x)dx =  f (a) + 4  f ( xi ) + 2  f ( x j ) + f (b)  −
b
 h f ( )
© Universidad Estatal de Milagro – UNEMI

(15)
a 3 i =1,3,5 j = 2,4,6  180
Ejercicio 3
Aplique la regla de Simpson 1/3 para integrar numéricamente la función tabular dada
por:
Fuente: (Chapra & Canale, 2006)

t 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4
f(t) 1 3 6 8 9 12 65 76 13 56 75 76 12 67 92
Tabla 4: Tabla de datos para ejercicio 3
Fuente: (Chapra & Canale, 2006)

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 12

 Integración Numérica

h  f ( t1 ) + f ( t3 ) + f ( t5 ) +   f ( t2 ) + f ( t4 ) + f (t6 ) +  
+ f ( t14 ) 
1.4
0
f (t )dt =  f (t0 ) + 4 
3 
 + 2
 f ( t ) + f ( t ) + f ( t ) + f (t )   f (t ) + f (t ) + f (t )  

 7 9 11 13   8 10 12 

0.1
I= 1 + 4 ( 3 + 8 + 12 + 76 + 56 + 76 + 67 ) + 2 ( 6 + 9 + 65 + 13 + 75 + 12 ) + 92
3 
I = 54.83
Ejercicio 4
Aplique la regla de Simpson 1/3 para resolver la siguiente integral con n=12
4
I =  log(t 3 ) coth(t )dt
1
Fuente: (Sauer, 2012)
Solución:
Primero se calcula el valor de h y luego se tabula la función
a =1 b=4
n = 12
4 −1
h= = 0.25
12

n t f(t)
0 1 0.0000
1 1.25 0.3427
2 1.5 0.5836
3 1.75 0.7745
4 2 0.9368
5 2.25 1.0803
6 2.5 1.2100
7 2.75 1.3288
8 3 1.4385
9 3.25 1.5403
10 3.5 1.6352
11 3.75 1.7240
© Universidad Estatal de Milagro – UNEMI

12 4 1.8074
Tabla 5: Tabulación de la función para aplicar Simpson 1/3.
Fuente: Elaboración propia.

h  f ( t1 ) + f ( t3 ) + f ( t5 ) +   f ( t2 ) + f ( t4 ) + f ( t6 ) +  
I=  f (t0 ) + 4  +
 f ( t ) + f ( t ) + f ( t )   f ( t ) + f ( t )
2  + f ( 12 )
t
3   7 9 11   8 10  
0.25   0.3427 + 0.7745 + 1.0803   0.5836 + 0.9368 + 1.21 +  
I= 0 + 4   + 2  + 1.8074
3   +1.3288 + 1.5403 + 1.7240  1.4385 + 1.6352  
I = 3.3815

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 13

 Integración Numérica

2.3 Subtema 3: Método de Simpson 3/8

Método de Simpson 3/8 para n=3

Este método utiliza un polinomio cúbico para aproximar la función a integrar. Para ello
se requiere cuatro puntos. En el intervalo [a, b] los puntos son x0=a, x3=b, los otros dos
puntos son intermedios x1 y x2 (Lara, 2019).

Figura 6: Simpson 3/8 con n=3

Fuente: (Chapra & Canale, 2006)
El polinomio cubico es:
p3 ( x) = b0 + b1 x + b2 x 2 + b3 x 3
Donde b0, b1, b2 y b3 son las constantes del polinomio y pasa por los puntos:
p( x0 ) = f ( x0 ), p( x1 ) = f ( x1 ), p( x2 ) = f ( x2 ) y p( x3 ) = f ( x3 )
Luego que se obtienen los valores de las constantes se tiene la integral.

3h
 f ( x0 ) + 3 f ( x1 ) + 3 f ( x2 ) + f ( x3 )
b b
 a
f ( x)dx  p( x)dx =
a 8
(16)

El método de Simpson 3/8 mejora la precisión de los descritos anteriormente. La

fórmula con el término del error es:

3h 3h5 (4)
f ( x)dx =  f ( x0 ) + 3 f ( x1 ) + 3 f ( x2 ) + f ( x3 )  −
b
 a 8 80
f ( ) (17)
© Universidad Estatal de Milagro – UNEMI

Método compuesto de Simpson 3/8

Este método divide el intervalo [a,b] en n sub-intervalos. Todos los sub-intervalos se

dividen en partes iguales. Se requieren cuatro puntos para construir un polinomio
cúbico, por lo que al ser aplicado a n sub-intervalos adyacentes, el método de Simpson
3/8 requiere los tres primeros intervalos. Luego al cuarto, quinto y sexto, y así
sucesivamente. Por tanto, el número total de sub-intervalo n debe ser divisible entre
3. La fórmula general usada en el método compuesto de Simpson es:
(Lara, 2019)

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 14

 Integración Numérica

3h  n −1 n −1 
b
f ( x ) dx   p ( x ) dx =  f ( a ) + 3   f ( xi ) + f ( xi +1 )  + 2  f ( x j ) + f (b )
b

a
a
8  i =1,4,7 j =3,6,9 
(18)

Ejercicio 6
Un automóvil recorre una pista en 90 segundos. La velocidad del automóvil se mide
cada 6 segundos mediante radar en pies / segundo, los datos de las mediciones se
registran en la siguiente tabla.
Fuente: (Burden & Faires, 2011)

n 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Tiempo 0 6 12 18 24 30 36 42 48 54 60 66 72 78 84 90
Velocidad 124 134 148 156 147 133 121 109 99 85 78 89 104 116 123 114
Tabla 6: Mediciones de velocidad en un instante de tiempo.
Fuente: (Burden & Faires, 2011).

Aplique el método de Simpson 3/8 para determinar la longitud de la pista

Solución:
La condición para aplicar Simpson 3/8 es que n debe ser divisible para 3, en este caso
n=15, por tanto, cumple la condición.
h = 6 , se reemplazan cada uno de los valores f(x) de la tabla 6, en la fórmula 18 y el
resultado queda como se muestra a continuación:
 f ( x1 ) + f ( x2 ) +  
   
 f ( x4 ) + f ( x5 ) +  
3h  n −1 n −1
 f ( x3 ) + f ( x6 ) +  
f ( x ) dx   f ( x0 ) + 3   f 
( x7 ) + f ( x8 ) +  + 2    + f ( x15 ) 
90
0 8 j =3,6,9,12  f ( x9 ) + f ( x12 ) + 
 i =1,4,7,10,13
f ( x10 ) + f ( x11 ) +  
  

  f ( x13 ) + f ( x14 )  


Al evaluar el resultado queda:

I = 10570.5
Ejercicio 7
Resuelva la siguiente integral para n=12 aplicando Simpson 3/8

12 x cos x
 ( 5 + senx ) dx
© Universidad Estatal de Milagro – UNEMI

2
0

Solución:

a=0 b= 
 −a
h= = 0.1477
12

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 15

 Integración Numérica

n x f(x)
0 0 0
1 0.1477045 0.0661734
2 0.295409 0.1211368
3 0.4431135 0.1629994
4 0.590818 0.1906689
5 0.7385224 0.2036134
6 0.8862269 0.2016596
7 1.0339314 0.1848322
8 1.1816359 0.1532377
9 1.3293404 0.1069878
10 1.4770449 0.0461583
11 1.6247494 -0.0292201
12 1.7724539 -0.1191404
Tabla 7: Tabulación de valores x y f(x) ejercicio 7.
Fuente: Elaboración propia.

 f ( x1 ) + f ( x2 ) +  
   
3h  n −1
f ( x4 ) + f ( x5 ) +  n −1
 f ( x3 ) + f ( x6 ) +  
I =  f ( x0 ) + 3    +2    + f ( x12 ) 
8  i =1,4,7,10,13 f
 ( x7 ) + f ( x8 ) +  j =3,6,9,12  f ( x9 )  
 f 
( x10 ) + f ( x11 ) +  
  
I = 0.2012
Integración con segmentos desiguales

Cuando se trabaja con segmentos desiguales se va aplicando las diferentes fórmulas de

integración numérica que se han visto en el compendio. Se divide el problema en
distintas integrales, se aplica la formula respectiva dependiendo del valor de n, luego
de que se tiene los valores de cada sub-intervalos, se suman todos los valores de la
integral y se obtiene el valor final.

Ejercicio 8: Tomado de (Chapra & Canale, 2006)

Para determinar el área transversal de un canal de agua se utiliza la siguiente integral:
B
Ac = Z ( y ) dy
© Universidad Estatal de Milagro – UNEMI

Siendo y= distancia desde uno de los márgenes (m), 𝐻= profundidad (m).

Análogamente, el flujo promedio Q ( m3 / s ) se calcula con:
B
Q = F ( y ) Z ( y ) dy
0

donde 𝐹= velocidad del agua ( m / s ). Aplique las diferentes fórmulas de integración

numérica para determinar Ac , y Q en base a los datos de la tabla 8:

y, m 0 2 4 5 6 9

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 16

 Integración Numérica

Z,m 0.5 1.3 1.25 1.7 1 0.25

F,m / s 0.03 0.06 0.05 0.12 0.11 0.02
Tabla 8: Mediciones de distancia, profundidad y velocidad del agua en un canal
Fuente: (Chapra & Canale, 2006)

Solución
La variable independiente es y (distancia desde uno de los márgenes)
y Z F Z*F
0 0.5 0.03 0.015
2 1.3 0.06 0.078
4 1.25 0.05 0.0625
5 1.7 0.12 0.204
6 1 0.11 0.11
9 0.25 0.02 0.005
Tabla 9: Calculo de U*H
Fuente: Elaboración propia

Los puntos 0, 2 y 4 están igualmente espaciados, se aplica Simpson 1/3, para n=2 y h=2
h
I1 =  f ( x0 ) + 4 f ( x1 ) + f ( x2 )
3
Calculo de Ac, el valor de H es f(x)
2
I A1 = 0.5 + 4(1.3) + 1.25
3
I A1 = 4.6333m 2
Calculo de Q, el valor de H*U es f(x)
2
I Q1 = 0.015 + 4(0.078) + 0.0625
3
m3
I Q1 = 0.2596
seg
Los puntos 4, 5 y 6 están igualmente espaciados, se aplica Simpson 1/3, para n=2 y h=1
h
I 2 =  f ( x0 ) + 4 f ( x1 ) + f ( x2 )
3
Calculo de Ac, el valor de H es f(x)
© Universidad Estatal de Milagro – UNEMI

1
I A2 = 1.25 + 4(1.7) + 1
3
I A 2 = 3.0166m 2
Calculo de Q, el valor de H*U es f(x)
1
IQ 2 = 0.0625 + 4(0.204) + 0.11
3
m3
IQ 2 = 0.3295
seg

Para los puntos 6 y 9 se aplica Trapecio simple, donde h=3

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 17

 Integración Numérica

h
I3 =  f ( x0 ) + f ( x1 )
2
Calculo de Ac, el valor de H es f(x)
3
I A3 = 1 + 0.25
2
I A1 = 1.875m 2
Calculo de Q, el valor de H*U es f(x)
3
IQ3 = 0.11 + 0.005
2
m3
IQ3 = 0.1725
seg
Ahora sumamos las respectivas integrales para obtener los valores de Ac y Q
Ac  I A1 + I A 2 + I A3
Ac  4.6333m 2 + 3.0166m 2 + 1.875m 2
Ac  9.5249m 2
Q  I Q1 + I Q 2 + I Q 3
Q  0.2596m3 / seg + 0.3295m3 / seg + 0.1725m3 / seg
m3
Q  0.7616
seg

Ejercicio 9
La siguiente tabla muestra las mediciones de concentración química de un reactor:
tf

M = Q  cdt
ti
t , min 0 1 4 6 8 12 16 20
c, mg / m 3
12 22 32 45 58 75 70 48
Tabla 10: Mediciones de concentración en función del tiempo.
Fuente: (Chapra & Canale, 2006)

Siendo Q=0.3m3/seg el flujo de salida, para el intervalo de tiempo [0, 20] min,
© Universidad Estatal de Milagro – UNEMI

determine en gramos la masa del producto.

Fuente: (Chapra & Canale, 2006)

Solución
Para que las unidades concuerden se convierte Q de 0.3 m3/seg a 18 m3/min
Al final se multiplicará el valor de la Integral por Q
tf

I =  cdt
ti
Entre los puntos para t=0 y t=1 se utiliza Trapecio simple con h=1

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 18

 Integración Numérica

h
I1 =  f ( x0 ) + f ( x1 )
2
1
I1 = 12 + 22
2
mg
I1 = 17 min
m3
Entre los puntos para t=1 y t=4 se utiliza Trapecio simple con h=3
h
I 2 =  f ( x0 ) + f ( x1 )
2
3
I2 =  22 + 32
2
mg
I 2 = 81 min
m3
Entre los puntos para t=4, t=6 y t=8 se utiliza Simpson 1/3 con h=2
h
I 3 =  f ( x0 ) + 4 f ( x1 ) + f ( x2 )
3
2
I 3 = 32 + 4(45) + 58
3
mg
I 3 = 180 min
m3
Entre los puntos para t=8, t=12, t=16 y t=20 se utiliza Simpson 3/8 con h=4
3h
I 4 =  f ( x0 ) + 3 f ( x1 ) + 3 f ( x2 ) + f ( x3 ) 
8
3(4)
I4 = 58 + 3(75) + 3(70) + 48
8
mg
I 4 = 811.5 min
m3
M  Q ( I1 + I 2 + I3 + I 4 )
m3 mg
M  18 (17 + 81 + 180 + 811.5 ) 3 min
min m
M  19611mg
M  19.61g
© Universidad Estatal de Milagro – UNEMI

Ejercicio 10
Se desea determinar el trabajo en un proceso termodinámico mediante:

T =  Pd v
Siendo T el trabajo, v=volumen y P=presión. Aplique las diferentes fórmulas de
integración para aproximar el trabajo en base a los datos de la tabla 11:
Fuente: (Chapra & Canale, 2006)

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 19

 Integración Numérica

(
P kJ / m3 ) 336 294.4 266.4 260.8 260.5 249.6 193.6 165.6

( )
v m3 0.5 2 3 4 6 8 10 11

Tabla 11: Mediciones de presión y volumen.

Fuente: (Chapra & Canale, 2006)

La variable con respecto a la que se integra es el volumen

Para los puntos donde v=0.5 y v=2, se aplica la regla del trapecio simple con h=1.5
h
I1 =  f ( x0 ) + f ( x1 )
2
1.5
I1 = 336 + 294.4
2
I1 = 472.8kJ
Para los puntos donde v=2, V=3 y v=4, se aplica la regla de Simpson 1/3 con h=1
h
I 2 =  f ( x0 ) + 4 f ( x1 ) + f ( x2 )
3
1
I 2 =  294.4 + 4(266.4) + 260.8
3
I 2 = 540.27kJ
Entre los puntos para v=4, v=6, v=8 y v=10 se utiliza Simpson 3/8 con h=2
3h
I 3 =  f ( x0 ) + 3 f ( x1 ) + 3 f ( x2 ) + f ( x3 ) 
8
3(2)
I3 =  260.8 + 3(260.5) + 3(249.6) + 193.6
8
I 3 = 1488,525kJ
Para los puntos donde v=10 y v=11, se aplica la regla del trapecio simple con h=1
h
I 4 =  f ( x0 ) + f ( x1 )
2
1
I 4 = 193.6 + 165.6
2
I 4 = 179.6kJ
T  I1 + I 2 + I 3 + I 4
© Universidad Estatal de Milagro – UNEMI

T  472.8kJ + 540.27kJ +1488.525kJ +179.6kJ

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 20

 Integración Numérica

3. Preguntas de Comprension de la
Unidad
¿Cuándo se recomienda aplicar integración numérica?
En los casos donde la función a integrar sea demasiado compleja o no exista una
antiderivada de dicha función, también se aplica integración numérica para resolver
problemas que tienen un conjunto de valores tabulados x, f(x) producto de mediciones
o experimentos. La integración numérica se aplica también para datos irregularmente
espaciados.

¿Explique la diferencia entre el método de Simpson 1/3 y Simpson 3/8?

El método de Simpson 1/3 parte de un polinomio cuadrático, por ello se necesita 3
puntos para obtener dicho polinomio que aproxima la función a integrar, llevando al
caso general para aplicar Simpson 1/3 n debe ser par. Mientras que el método de
Simpson 3/8 se basa en aproximar la función a integrar mediante un polinomio cubico,
se requieren de 4 puntos, para el caso general de Simpson 3/8 el valor de n debe ser
divisible para 3.

¿Cuál es el error al aplicar el método del trapecio compuesto con h=0.5 de la integral
2

 e dx ?
x

La fórmula del trapecio compuesta es:

h n −1
 b−a 2

b
a
f ( x)dx = 
2
f ( x0 ) + 2
i =1
f ( xi ) + f ( xn )  −
 12
h f ''( )

Donde el término del error es:

b−a 2
h f ''( )
12
El valor de zita es: max ( f ( ) ) ,0    2
f ''( x) = e x
© Universidad Estatal de Milagro – UNEMI

2−0
error = (0.5) 2 e2
12
error = 0.1539

¿Cuál es el error al aplicar el método Simpson 1/3 compuesto con n=12 de la integral
4

 sin xdx ?
1
La fórmula de Simpson 1/3 compuesta es:

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 21

 Integración Numérica

h n −1 n−2  b − a 4 (4)
 
b
a
f ( x)dx = 
3
f ( a ) + 4
i =1,3,5
f ( xi ) + 2
j = 2,4,6
f ( x j ) + f (b)  −
 180
h f ( )

Donde el término del error es:

b − a 4 (4)
h f ( )
180
b − a 4 −1
h= = = 0.25
n 12
f iv ( x) = sin( x) 1   4 f iv ( ) = sin(1.5707) = 1
4 −1
error = (0.25)4 (1)
180
error = 6.51*10−5

¿Cuál es el error al aplicar el método Simpson 3/8 con n=3 de la integral

3

 (1 + x )
−1
dx ?
1.5

3h 3h5 (4)
 
b
a 8
f ( x)dx =
f ( x0 ) + 3 f ( x1 ) + 3 f ( x2 ) + f ( x3 ) −
80
f ( )

Donde el término del error es:

3h5 (4)
f ( )
80
b − a 3 − 1.5
h= = = 0.5
n 3
f iv ( x) = 24 (1 + x ) f iv ( ) = 24 (1 + 1.5) = 0.24576
−5 −5
1.5    3
3(0.5)5
error = (0.24576)
80
error = 2.879*10−4
© Universidad Estatal de Milagro – UNEMI

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 22

 Integración Numérica

4. Material Complementario
Los siguientes recursos complementarios son sugerencias para que se pueda ampliar la
información sobre el tema trabajado, como parte de su proceso de aprendizaje
autónomo:

Videos de apoyo:
Integración Numérica: https://bit.ly/3kZ2qNw
Método del Trapecio: https://bit.ly/36fbN5H
Método de Simpson: https://bit.ly/2IgOZuq
https://bit.ly/3n3dYQt
Método de Simpson 3/8: https://bit.ly/3ke3WKy
Método de Simpson 1/3: https://bit.ly/2UcYBIN

Bibliografía de apoyo:
▪ Burden, R., & Faires, D. (2011). Numerical Analysis. (BOOKS/COLE, Ed.) (Ninth).
Boston.
▪ Chapra, S., & Canale, R. (2006). Métodos numéricos para ingenieros. (McGrawHill,
Ed.) (Quinta). México
Links de apoyo:
Link de libros: https://bit.ly/34WB0Td
© Universidad Estatal de Milagro – UNEMI

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 23

 Integración Numérica

5. Bibliografía

» Burden, R., & Faires, D. (2011). Numerical Analysis. (BOOKS/COLE, Ed.) (Ninth).
Boston.
» Chapra, S., & Canale, R. (2006). Métodos numéricos para ingenieros.
(McGrawHill, Ed.) (Quinta). México.
» Gutiérrez, J., Olmos, M., & Casillas, J. (2010). Análisis Numérico. (McGrawHill,
Ed.). México. Retrieved from http://librosysolucionarios.net/
» Lara, O. (2019). Métodos Numéricos: Teoría y Aplicaciones. 2018.
» Nieves, A., & Dominguez, F. (2014). Métodos Numéricos aplicados a la
ingeniería. (Patria, Ed.) (Primera). México.
» Rodríguez, L. (2016). Análisis Numérico Básico Un enfoque algorítmico con el
soporte de Python. Guayaquil. Retrieved from
http://calicreativa.com/natunova-flores-comestibles/
» Sauer, T. (2012). Análisis Numérico. (PEARSON, Ed.) (Segunda). México.
Retrieved from http://librosysolucionarios.net/
© Universidad Estatal de Milagro – UNEMI

FORMATO CONTROLADO: FR0018/ v1.0 / 18-11-2019 24