trapezoid.m

algotithm to calculate pi

function y=trapezoid(n)
for k=1: n-1
    a(k)=2*fun(k/n);
end;
vpa(1/(2*n)*(sum(a)+fun(0)+fun(1)))