clc 
%	Regla de Simpson de 1/3
% 	Ing. Angel Hernández Corvera
%	Posgrado en Electrónica
% 	MC J. Gilberto Mateos Suárez
%	Noviembre de 1997
format long
echo on
% La función a evaluar esta contenida dentro del 
% archivo 'funcion.m' y se requiere introducir los valores
% de los límites, como se indica a continuación
pause % Pulsa una tecla

a= input('límite inferior:');
b= input('límite superior:');
c= (b-a)/2;

% Formato general de Simpson de 1/3 que da el resultado
% de la integración de la función en el archivo 'funcion.m'.

I=(b-a)*((funcion(a)+4*funcion(c)+funcion(b))/6)