clc
%	UNIVERSIDAD DE GUADALAJARA
%	Posgrado en Electrónica
%	Método de Simpson 3/8 de segmentos múltiples impares
%	Noviembre de 1997
%	Ing. José Martín Villegas González
%	M.C. J.Gilberto Mateos Súarez

echo on
format long
a=input('limite inferior:');
b=input('limite superior:');
n=input('segmentos impares :');
h=(b-a)/n
%	Subrutina para la regla de Simpson de 3/8....
% Pulsa una tecla...
pause
if (-1)^n==-1, b=b-(3*h),n=(n-3),
e=b+h,
f=b+2*h,
g=b+3*h,
k=(e)*((funcion(b)+3*(funcion(e)+funcion(f))+funcion(g))/8);end
pause	
clc
%	Subrutina para la sumatoria impar de Simpson 1/3
% pulsa una tecla...
pause
p=0;
for i=1:n;
	if (-1)^i==1,p=p+funcion(a+(i-1)*h); end
	p
end
clc
%	Subrutina para la sumatoria par
% pulsa una tecla...
pause
s=0;
for i=2:n;
	if (-1)^i==-1,s=s+funcion(a+(i-1)*h); end
	s
end
clc
	
%	Y se tiene la fórmula de Simpson 1/3 de segmentos múltiples
I=(b-a)*(funcion(a)+4*p+2*s+funcion(b))/(3*n)
%	Finalmente la suma de las reglas de 1/3 y 3/8 es....
I=I+k
%	Que es el valor de la integración