UNIVERSIDAD TÉCNICA DE MACHALA
Calidad, Pertinencia y Calidez
UNIDAD ACADEMICA DE CIENCIAS AGROPECUARIAS
ESCUELA DE INGENIERÍA AGRONÓMICA
LABORATORIO COMPUTACIONAL MATEMÁTICO
NÚMERO DE PRACTICA: AC.02.04.9
TEMA: INTEGRACION NUMERICA
1. DATOS INFORMATIVOS
CARRERA: INGENIERÍA AGRONOMICA CICLO/NIVEL: TERCER CICLO “B”
FECHA: ENERO DEL 2019 ALUMNO: Wagner Ramos
DOCENTE RESPONZABLE: ING. BLADIMIR SERRANO RUGEL
TÉCNICO DOCENTE RESPONSABLE: ING. HENRRY AGUILAR AGUILAR
2. FUNDAMENTACION TEÓRICA.
FÓRMULAS DE INTEGRACIÓN DE NEWTON-COTES
Estas fórmulas se basan en la idea de integrar una función polinomial en vez de
Donde es un polinomio de aproximación de grado n para
ciertos valores de f(x) que se escogen apropiadamente (se suele conocer también como
polinomio de interpolación, ya que la condición es que tome los mismos valores que la
función original en los puntos elegidos). Estas fórmulas se pueden aplicar también a una
tabla de datos, siendo éstos los puntos a considerar.
Dentro de las fórmulas de Newton-Cotes, existen las formas cerradas y abiertas. En las
formas cerradas se conocen los valores de f(a)y f(b), en caso contrario, se llaman formas
abiertas. Nos remitiremos a estudiar únicamente las formas cerradas, y por lo tanto,
siempre supondremos que conocemos los valores de los extremos, f(a) y f(b).
METODO DE SIMPSON
Este método consiste en la aproximación del cálculo del área plana bajo una
curva utilizando trapecios curvilíneos a partir una interpolación con una función
cuadrática:
Esta aproximación es denominada “simple” debido a que utiliza tan solo un
polinomio. Requiere el conocimiento de tres puntos equiespaciados: los
extremos y un punto central. Aplicando esta expresión utilizando mayor cantidad
de puntos intermedios (es decir, realizando más de un Simpson 1/3 dentro del
intervalo) puede definirse la variante “compuesta” del método para el cual se
utilizan N puntos que corresponden a n = N − 1 subintervalos. Este caso requiere
que la cantidad de subintervalos sean pares (el caso simple utiliza dos, por lo
tanto ´este debe ser un múltiplo). Por lo tanto, se define un valor h que
corresponde a el ancho del subintervalo o el paso que hay entre puntos. Se
calcula como h = b−a/n. Finalmente, la aproximación del ´área se puede calcular
como:
A modo de interpretación, las sumatorias corresponden a los valores de la
función en los puntos impares y pares respectivamente sin contar los extremos.
Esto es debido a que, en los puntos impares, se encuentra el factor de 4 que
introduce la fórmula de Simpson 1/3 simple para los puntos intermedios y, en los
puntos pares, el factor de 2 se debe a que ese punto es compartido por los
trapecios curvilíneos adyacentes.
[Link]:
Desarrollar un programa para aproximar integrales definidas aplicando la fórmula
de Simpson, implementado Phyton.
4. METRIALES
COMPUTADOR
SOFTWARE DE CÓDIGO ABIERTO
MODULOS SymPy, NumPy y Matplotlib
5. PROCEDIMIENTO
Vamos a Aproximar la siguiente integral definida:
, Para M= 5;10;15;20
Abrimos Python, e importamos los módulos SymPy, NumPy y Matplotlib, de la
siguiente manera from sympy import*, import numpy as np, import
[Link] as plt, Luego declaramos las variables de forma simbólica, con
el comando Symbol de esta forma x=Symbol('x') para la variable x, con el mismo
procedimiento vamos a declarar a la variable y. Para nuestro problema la variable
“Y”, definirá la función, la vamos a solicitar en nuestro programa por medio del
teclado, de esta forma f=input('Ingrese la función: '), ahora pedimos por teclado
el valor de x, el cual nos servirá para evaluar la aproximacionde la función en
dicho punto.
Abrimos Python, e importamos los módulos SymPy, NumPy y Matplotlib, de la
siguiente manera from sympy import*, import numpy as np, import
[Link] as plt, Luego declaramos las variables de forma simbólica, con
el comando Symbol de esta forma x=Symbol('x') para la variable x, con el mismo
procedimiento vamos a declarar a la variable y. Para nuestro problema la variable
“Y”, definirá la función, la vamos a solicitar en nuestro programa por medio del
teclado, de esta forma f=input('Ingrese la función: '), ahora pedimos por teclado
el valor de x, el cual nos servirá para evaluar la aproximacion de la función en
dicho punto.
Ahora pedimos por teclado los extremos superior e inferior de la integral definida
de la siguiente manera
a=float(input("ingrese el extremo inferior: "))
b=float(input("ingrese el extremo superior: "))
Se tomara “M” como el numero de intervalos de la integral,
M=int(input("ingrese el numero de intervalos: "))
F=sympify(f)
Para el calculo del valor de h se utilizara la siguiente formula
h=(b-a)/(2*M)
s1=0
s2=0
En este punto se establecere doble derivada por lo que estableceremos siguiente
con el programa.
for k in range(1,M):
x1=a+h*(2*k-1)
s1=s1+[Link](x,x1)
for k in range(1,M-1):
x1=a+h*2*k
s2=s2+[Link](x,x1)
Ahora se establece la formula para el calculo
s=h*([Link](x,a)+[Link](x,b)+4*s1+2*s2)/3
print("el valor aproximado de la integral es:",' ', k) , se utiliza print para que el
programa pueda imprimir la respuesta.
5. RESULTADOS
Este programa nos permitió calcular una aproximación a una función usando la
Diferenciación Numérica aplicada en Phyton.
CÓDIGO
Figura 1; Aquí se muestra el código a utilizarse para desarrollar el calculo de las
Aproximaciones de la siguiente integral.
CÁLCULOS
Figura 2
En la figura 2 se muestra los resultados de evaluar la integral
para M=5
Figura 3
En la figura 3 se muestra los resultados de evaluar la integral
para M=10
Figura 4
En la figura 4 se muestra los resultados de evaluar la integral
para M=15
Figura 5
En la figura 5 se muestra los resultados de evaluar la integral
para M=20
6. CONCLUSIONES
La integración con métodos numéricos es una herramienta útil cuando se trata
de integrar una función muy complicada o datos tabulados.
Con el método de Simpson se puede aproximar una integral compleja a la
integral de un polinomio.
Con el método de Simpson se logra convertir matemáticas superiores en
aritméticas básicas.
Mediante el uso de la tecnología podemos realizar operaciones más complejas
en menos tiempo.
7. BIBLIOGRAFÍA
[Link]
ma5cTeoria(09-10).pdf
[Link]
comp/[Link]
