%
% Programa de chamada da function trapezios
%
clear;
clc;
% dados
xinf=-10;
xsup=10;
m=2000;
h=(xsup-xinf)/m;
% comparacao dos resultados trapezios vs Simpson
trap=trapezios('f',xinf,xsup,m);
simp=quad('f',xinf,xsup);
erro=abs(simp-trap);
% tabela
disp('----------------------------------------------');
disp('   Trapezios       Simpson          Erro      ');
disp('----------------------------------------------');
fprintf('\n%12.5f%15.5f%15.5f%15.5f\r',simp,trap,erro); 
disp(' ');
disp('----------------------------------------------');
% graficos
subplot(2,1,1), fplot('f',[xinf xsup],'r'), title('Grafico da funcao')
subplot(2,1,2), area(feval('f',[xinf:h:xsup])), 
	title(['Integral = ',num2str(trap),'; Erro = ',num2str(erro)]), set(gca,'XTickLabel','','YTickLabel','')

