function trapezoid
  
n=20;
  
a=0;
b=pi;

h = (b-a)/n;
sum = .5*h*(sin(a)+sin(b));

for i=1:(n-1)
  sum = sum + h*sin(a+i*h);
end

sum
err=abs(2-sum)