clc;
clear all;
close all;
h=0.001; % periodo de muestreo
t=0:h:10; % intervalo de integración
f=sin(t); % función a integrar
I_e=0; % condición inicial de la integral por
        % método de Euler
% método comparativo trapezoidal 
I_trap=trapz(t,f);
% método comparativo regla de Simpson
I_simp=quad('sin',0,10);
[m, n]=size(t); % Orden del vector de tiempo\n
for k=1:n
     I_e=I_e+h*f(k);
end  
fprintf('Método Euler= %f \n Método Trapezoidal= %f \n Regla de Simpson =%f \n',I_e,I_trap, I_simp)