clc
%	UNIVERSIDAD DE GUADALAJARA
%	Posgrado en Electrónica
%	Método de Simpson 1/3 de segmentos múltiples
%	Noviembre de 1997
%	Ing. Angel Hernández Corvera
%	M.C. J.Gilberto Mateos Súarez
format long
echo on
a=input('limite inferior:');
b=input('limite superior:');
n=input('segmentos pares:');
h=(b-a)/n;
clc
%	Subrutina para la sumatoria impar
% 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)
%	Que es el valor de la integración