simpson.m

algotithm to calculate pi

function y=simpson(m)
for j=1: m
    a(j)=4*fun((2*j-1)/(2*m));
end
for i=1: m-1
    b(i)=2*fun((2*i)/(2*m));
end
vpa(1/(6*m)*(sum(a)+sum(b)+fun(0)+fun(1)))