clc;
clear all;
close all;

tini=0;
tinc=pi/1000;
tfinal=pi;
x=tini:tinc:tfinal;
f=sin(x);
I_trapz= trapz(x,f);
I_a= 1-cos(pi);
disp('Resultado')
%despliega los resultados comparativos
%del  método trapezoidal con el analítico. 
disp([I_trapz I_a])