MÉTODOS NUMÉRICOS

INTEGRACIÓN NUMÉRICA
PARTE TEÓRICA

INTRODUCCIÓN

La Integración Numérica es uno de los métodos más citados dentro del campo de la Ingeniería, en general. De
gran aplicación en todas las áreas de ciencia y tecnología, pues consiste en encontrar el área bajo la curva de
una función tabular, usando diversas técnicas que permitan aproximar de la mejor manera el valor correcto.

Toda Integración Numérica parte de la fórmula de Interpolación de Newton, para el intervalo X0 ≤ X ≤ Xn:

( ) ( ) ( ) ( ) ( )

Se integra la anterior función de la forma:

∫ ( ) ∫[ ( ) ( ) ( ) ( ) ]

Pero dx = hdk, entonces de tiene:

∫ ( ) [ ( ) ( ) ]

REGLA DEL TRAPECIO:

En la anterior ecuación, si se considera sólo las primeras diferencias, se tiene:

∫ ( ) [ ]

Y la integral sólo se calcula entre los dos primeros valores, es decir, entre x 0 y x1, obteniéndose como resultado
final:

∫ ( ) [ ]

Que es el área bajo la curva entre los puntos de coordenadas (x0,y0) y (x1,y1).

En forma general, se tiene la Fórmula del Trapecio Compuesta igual a:

Luis Cabezas Tito Métodos Numéricos Página 1

 ∫ ( ) [ ∑ ]

REGLA DE SIMPSON 1/3:

Si se desea obtener una mejor aproximación de la integral definida entre dos puntos, se considera el
polinomio Interpolante hasta las segundas diferencias, es decir:

∫ ( ) [ ( ) ]

Realizando las integraciones respectivas, se llega a la fórmula general de la Regla de Simpson 1/3, que es igual
a:

∫ ( ) [ ∑ ∑ ]

A diferencia de la fórmula de Integración Trapezoidal, la ecuación anterior se aplica solamente para valores
pares de n, ya que de lo contrario, no se podrá obtener el valor de la última integral. Si el número de puntos
fuese impar, entonces el último tramo puede integrarse mediante la Regla del Trapecio.

REGLA DE SIMPSON 3/8:

Se aplica para tener una mejor aproximación de la integral, para lo cual el polinomio Interpolante respectivo,
se aproxima hasta el tercer grado.

∫ ( ) [ ( ) ( ) ]

Obteniéndose como fórmula general de esta Regla, lo siguiente:

∫ ( ) [ ∑ ∑ ]

Si n = 4:

∫ ( ) [ ]

Si n = 5:

Luis Cabezas Tito Métodos Numéricos Página 2

 ∫ ( ) [ ]

Si n = 6:

∫ ( ) [ ]

EJEMPLO:

Integrar la función tabular mostrada a continuación, desee x = 0 hasta x = 6, usando:

a. La Regla del Trapecio.

b. La Regla de Simpson 1/3.
c. La Regla de Simpson 3/8.
d. La fórmula de 6 puntos.

X 0 1 2 3 4 5 6
Y -5 1 9 25 55 105 181

SOLUCIÓN:

a. Regla del Trapecio:

∫ ( ) [ ∑ ]

∫ ( ) [ ( )]

b. Regla de Simpson 1/3:

∫ ( ) [ ∑ ∑ ]

∫ ( ) [ ( ) ( )]

c. Regla de Simpson 3/8:

∫ ( ) [ ∑ ∑ ]

Luis Cabezas Tito Métodos Numéricos Página 3

 ∫ ( ) [ ( ) ( )]

d. Fórmula de 6 Puntos:

∫ ( ) [ ]

NOTA IMPORTANTE:

Si sólo existirían valores desde x = 0 hasta x = 5, entonces se puede integrar por partes, de la siguiente manera:

 Desde x = 0 hasta x = 3, aplicar la Regla de Simpson 3/8.

 Desde x = 3 hasta x = 5, aplicar la Regla de Simpson 1/3.

O también se puede hacer:

 Desde x = 0 hasta x = 4, aplicar la Regla de Simpson 1/3.

 Desde x = 4 hasta x = 5, aplicar la Regla del Trapecio.

PRÁCTICA 1:

Integrar la función tabular mostrada a continuación, desee x = 0 hasta x = 8, usando:

a. La Regla del Trapecio.

b. La Regla de Simpson 1/3.
c. La Regla de Simpson 3/8.
d. La fórmula de 6 puntos.

X 0 1 2 3 4 5 6
y -3 2 5 23 35 51 79

PRÁCTICA 2:

Para la siguiente función tabular:

X 0 1 2 3 4 5
y -3 2 5 23 35 51

Calcular la integral, usando las dos técnicas sugeridas en la NOTA IMPORTANTE.

PRÁCTICA 3:

Con los datos del Ejemplo, calcular las integrales desde x = 0 hasta x = 5, tomando en cuenta las dos técnicas
sugeridas en la NOTA IMPORTANTE y compare los resultados finales obtenidos con los encontrados en dicho
Ejemplo.

Luis Cabezas Tito Métodos Numéricos Página 4

PARTE PRÁCTICA

INTEGRACIÓN NUMÉRICA CON MATLAB

Se desea integrar la siguiente expresión mediante MATLAB:

Introducir en el workspace el siguiente código:

>> quad('1./(x.^3-2*x-5)',0,2)

ans =
-0.4605

O bien este otro:

>> f='1./(x.^3-2*x-5)'

f=
1./(x.^3-2*x-5)

>> quad(f,0,2)

ans =
-0.4605

REGLA DEL TRAPECIO:

EJEMPLO 1: aproxime el área bajo la curva de la siguiente función tabular, entre x = -1 y x = 4 utilizando la
regla trapezoidal compuesto.

x -1 0 1 2 3 4
f(x) 8 10 10 20 76 238

SOLUCIÓN:

Copiar previamente a una carpeta la función Trapezoida_L.p, luego, escribir en el workspace lo siguiente:

>> x=[-1 0 1 2 3 4];

>> y=[8 10 10 20 76 238];
>> Trapezoidal_L(x,y)

SOLUCION
I = 239.0000000
ans =
239

Luis Cabezas Tito Métodos Numéricos Página 5

O bien el formato avanzado:

>> sol=Trapezoidal_L(x,y)

sol =
239

EJEMPLO 2: encuentre la integral aproximada, por medio de este método, de la función:

( )
√

Que da lugar a la curva normal tipificada entre los límites -1 y 1.

SOLUCIÓN:

Copie en un archivo script el siguiente código y guárdelo con el nombre de Integral.

function y = Integral(x)
y=(1/sqrt(2*pi))*exp((-x.^2)/2);

En el workspace escriba:

>> Trapezoidal_L(-1,1)

METODO TRAPEZOIDAL
Arch. Ecuación ? Integral

SOLUCION
Nro Franjas = 1024
I = 0.6826893
ans =
0.6827

Luis Cabezas Tito Métodos Numéricos Página 6

Usando el formato avanzado:

>> sol=Trapezoidal_L(-1,1,'Integral')

sol =
0.6827

REGLA DE SIMPSON 1/3:

EJEMPLO 3: mediante el método de Simpson 1/3 aproxime el área bajo la curva de la siguiente función dada
en formato tabular, entre x = -1 y x = 3:

x -1 0 1 2 3
f(x) 8 10 10 20 76

SOLUCIÓN:

Copiar previamente a una carpeta la función Simpson13_L.p, luego, escribir en el workspace lo siguiente:

>> x=[-1 0 1 2 3];

>> y=[8 10 10 20 76];
>> Simpson13_L(x,y)

SOLUCION
I = 74.6666667

ans =
74.6667

Luis Cabezas Tito Métodos Numéricos Página 7

Usando el formato avanzado:

>> sol=Simpson13_L(x,y)

sol =
74.6667

EJEMPLO 4: encuentre la integral aproximada, por medio de este método, de la función:

( )
√

Que da lugar a la curva normal tipificada entre los límites -1 y 1.

SOLUCIÓN:

Copie en un archivo script el siguiente código y guárdelo con el nombre de Integral.

function y = Integral(x)
y=(1/sqrt(2*pi))*exp((-x.^2)/2);

En el workspace escriba:

>> Simpson13_L(-1,1)

METODO SIMPSON 1/3

Arch. Ecuación ? Integral

SOLUCION
Nro Franjas = 64
I = 0.6826895

ans =
0.6827

Luis Cabezas Tito Métodos Numéricos Página 8

Usando el formato avanzado:

>> sol=Simpson13_L(-1,1,'Integral')

sol =
0.6827

REGLA DE SIMPSON 3/8:

EJEMPLO 5: mediante el método de Simpson 3/8 aproxime el área bajo la curva de la siguiente función dada
en formato tabular, entre x = 0.7 y x = 0.4:

x 0.70 0.65 0.60 0.55 0.50 0.45 0.40

f(x) 2.20 2.17 2.13 2.09 2.04 1.99 1.94

SOLUCIÓN:

Copiar previamente a una carpeta la función Simpson38_L.p, luego, escribir en el workspace lo siguiente:

>> x=[0.70,0.65,0.60,0.55,0.50,0.45,0.40];
>> y=[2.20,2.17,2.13,2.09,2.04,1.99,1.94];
>> Simpson38_L(x,y)

SOLUCION
I = 0.6245625

ans =
0.6246

Luis Cabezas Tito Métodos Numéricos Página 9

Evaluando con el formato avanzado:

>> sol=Simpson38_L(x,y)

sol =
0.6246

EJEMPLO 6: encuentre la integral aproximada, por medio de este método, de la función:

( )

Desde a = 0 hasta b = 0.8.

SOLUCIÓN:

Copie en un archivo script el siguiente código y guárdelo con el nombre de Integracion.

function y = Integracion(x)
y=0.2+25*x-200*x.^2+675*x.^3-900*x.^4+400*x.^5;

En el workspace escriba:

>> Simpson38_L(0,0.8)

METODO SIMPSON 3/8

Arch. Ecuación ? Integracion

SOLUCION
Nro Franjas = 192
I = 1.6405333

ans =
1.6405

Luis Cabezas Tito Métodos Numéricos Página 10

Resolviendo en formato avanzado:

>> sol=Simpson38_L(0,0.8,'Integracion')

sol =
1.6405

CUADRATURA DE GAUSS - LEGENDRE:

Es una clase de técnica que permite utilizar valores de f(x) en abscisas desigualmente espaciadas, y que se
encuentran determinadas por ciertas propiedades de familias de polinomios ortogonales.

El objetivo de la cuadratura de Gauss es determinar los coeficientes de una ecuación de la forma:

( ) ( )

Donde w son coeficientes a ser determinados.

La técnica utilizada generada la fórmula:

∫ ( ) ( ) ( )

Donde z1 = -0.57735 y z2 = 0.57735, respectivamente. El integrando f(x)dx en términos de la nueva variable

queda definida como:

∫ ( ) [ ( ( ) ) ( )
( ))]

Este método es ideal para funciones cúbicas.

Luis Cabezas Tito Métodos Numéricos Página 11

EJEMPLO 7: integre la función f(x) en el intervalo (-0.8, 1.5) por Cuadratura de Gauss – Legendre, utilizndo do
2 puntos.

( )
√

SOLUCIÓN:

Copie en un archivo script el siguiente código y guárdelo con el nombre de Integral.

function y = Integral(x)
y=(1/sqrt(2*pi))*exp((-x.^2)/2);

En el workspace escriba:

>> CuadraturaGauss_L(-0.8,1.5,2)

METODO CUADRATURA DE GAUSS

Arch. Ecuación ? Integral
Orden Polinomio ? 2

SOLUCION
I = 0.7218250

ans =
0.7218

INTEGRALES MÚLTIPLES:

Cualquiera de los métodos de integración existentes se puede utilizar en la aproximación de integrales dobles o
triples, realizando una integración repetida y manteniendo, en cada caso, una determinada variable como si se
tratase de una constante.

EJEMPLO 8: encuentre la integral aproximada de la función:

∫∫

Utilizando el método de Simpson 1/3 para ambos ejes. La integral interna con 4 franjas y la integral externa
con 6 franjas.

SOLUCIÓN:

Copie en un archivo script el siguiente código y guárdelo con el nombre de Multiple.

function valor = Multiple(x,y)

valor=exp(x+y);

Luis Cabezas Tito Métodos Numéricos Página 12

En el workspace escriba:

>> IntegralMultiple_L

METODO INTEGRAL MULTIPLE

Arch. Ecuación ? Multiple
DATOS INTEGRAL INTERNA
Limite Inferior ? 0
Limite Superior ? 4
Nro. de Franjas ? 4
METODO DE TRABAJO
-----------------------
Trapezoidal 1
Simpson 1/3 2
-----------------------
Método ? 2
DATOS INTEGRAL EXTERNA
Limite Inferior ? 1
Limite Superior ? 3
Nro. Franjas ? 6
METODO DE TRABAJO
-----------------------
Trapezoidal 1
Simpson 1/3 2
-----------------------
Método ? 2

SOLUCION
-----------------------
Solución I = 935.5304715

ans =
935.5305

En formato avanzado se escribe:

>> sol=IntegralMultiple_L('Multiple',0,4,4,2,1,3,6,2)

sol =
935.5305

INTEGRALES DOBLES CON REGIONES DELIMITADAS POR ECUACIONES:

Hasta ahora solo se han visto integraciones dobles sobre regiones R rectangulares. No obstante, también
pueden resolverse integrales del tipo:

( ) ( )

∫ ∫ ( ) ∫ ∫ ( )
( ) ( )

Luis Cabezas Tito Métodos Numéricos Página 13

EJEMPLO 9: encuentre la integral aproximada de la función:

∫ ∫( )

Utilizando el método de Simpson 1/3 para ambos ejes. La integral interna con 4 franjas y la integral externa
con 6 franjas.

SOLUCIÓN:

Copie en un archivo script el siguiente código y guárdelo con el nombre de Ecuaciones.

function [L,fxy] = Ecuaciones(x,y)

L(1)=x.^2;
L(2)=2*x;
fxy=x.^3+4*y;

En el workspace escriba:

>> IntegralMultipleE_L

METODO INTEGRAL MULTIPLE

Arch. Ecuación ? Ecuaciones
DATOS INTEGRAL INTERNA
Nro. de Franjas ? 2
METODO DE TRABAJO
-----------------------
Trapezoidal 1
Simpson 1/3 2
-----------------------
Método ? 2
DATOS INTEGRAL EXTERNA
Limite Inferior ? 0
Limite Superior ? 2
Nro. Franjas ? 4
METODO DE TRABAJO
-----------------------
Trapezoidal 1
Simpson 1/3 2
-----------------------
Método ? 2

SOLUCION
-----------------------
Solución I = 10.5833333

ans =
10.5833

Luis Cabezas Tito Métodos Numéricos Página 14

En formato avanzado:

>> sol=IntegralMultipleE_L('Ecuaciones',2,2,0,2,4,2)

sol =
10.5833

PROBLEMAS DE APLICACIÓN:

PROBLEMA 1:

Calcule las integrales de todos los ejemplos dados en la parte teórica, usando los diferentes métodos
estudiados para resolver mediante MATLAB.

PROBLEMA 2:

Haga cuatro capturas de datos mediante la conexión ARDUINO – MATLAB, utilizando sensores: un
potenciómetro, un LM35, una fotoresistencia y un sensor de efecto Hall. Para cada uno de ellos obtendrá las
variables (t,y). Con esos datos obtenga un polinomio p e intégrelo por todos los métodos de resolución de
integrales usando MATLAB.

PROBLEMA 3:

Usando los datos capturados del Problema 2, realice Arte Digital al enviar los mismos mediante la conexión
MATLAB – ARDUINO a un cubo de LEDS 3x3x3 y 8x8x8. Añada música electrónica al evento y grabe en un video
la secuencia enviada.

Luis Cabezas Tito Métodos Numéricos Página 15

 ----- oOo -----

Luis Cabezas Tito Métodos Numéricos Página 16