函數y = traint(a，b，n，f)
h =(b-a)/n；
x = linspace(a，b，n+1)；
y1 = h * feval(f，x)；
y 1(1)= y 1(1)/2；
y 1(n+1)= y 1(n+1)/2；
y = sum(y 1)；
%復合辛普森公式復合辛普森公式
函數y = sraint(a，b，n，f)
h =(b-a)/n；
x = linspace(a，b，2 * n+1)；
y1 = feval(f，x)；
y 1(2:2:2 * n)= 4 * y 1(2:2:2 * n)；
y 1(3:2:2 * n-1)= 2 * y 1(3:2:2 * n-1)；
y = h/6 * sum(y 1)；
& gt& gtf =內聯(' 4。/(1+x . * x)')；
& gt& gttraint(-1，1，8，f)
ans =
6.2624
& gt& gtsraint(-1，1，8，f)
ans =
6.2832