function y=trapezios(fun,xinf,xsup,m)
%
% Integracao numerica usando a regra dos trapezios
% 
%
h=(xsup-xinf)/m;
x=linspace(xinf,xsup,m+1);
f=feval(fun,x);
y=0.5*h*(f(1)+2*sum(f(2:end-1))+f(end));