复合梯形公式与复合辛普森公式MATLAB

采用复合梯形公式与复合辛普森公式，计算 sin(x)/x 在[0, 1]范围内的积分。采样点数 目为 5、9、17、33。

function Sn=CompoundSimpson（fabn）
h=（b-a）/n;
for k=0:n   
    x（k+1）=a+k*h;     
    xk（k+1）=x（k+1）+1/2*h;   
    if （x（k+1）==0）    
        x（k+1）=10^（-10）;       
    end
    if （xk（k+1）==0）            
        xk（k+1）=10^（-10）;
    end
end
p1=h/6*（f（x（1））+f（x（n+1）））;
for i=2:n     
    fx1（i）=h/3*f（x（i））;
end
for j=1:n     
    fx2（j）=2/3*h*f（xk（j））;
end
p2=sum（fx1）;
p3=sum（fx2）;
Sn=p1+p2+p3;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-03-18 00:43  5\
     文件         271  2018-07-05 20:23  5\CompoundTrapezoid.m
     文件         422  2018-07-05 20:37  5\CompoundSimpson.m