func5ok.m

数值分析最常用的四十种算法

n=500
a=linspace(0,1,n)
s=0;
for i=2:n-1       %i can't start from 1,else the denominator will be 0
s=s+func5(a(i))+4*func5(a(i)+1/n/2)+func5(a(i+1))
disp(s)
end
s=s*1/n/6