试用matlab语言编写一个小程序，并说明该程序实现的功能。

Question

如题所述

Answer 1

function s = simpson(a,b,N);
if mod(N,2)==1
'error!'
return
end
h=(b-a)/N;
s0=fun(a)+fun(b);
s1=0;
s2=0;
for i=1:N-1
x=a+i*h;
if mod(i,2)==1
s1=s1+fun(x);
else
s2=s2+fun(x);
end
end
s=(h/3)*(s0+4*s1+2*s2)

simpson法求积分