Question

数值方法中要求用Matlab复合Simpson公式做此题（要求完整过程，急用，求高手！！！~~~）

Matlab复合Simpson公式(取n=4）计算积分∫0~1（(sinx)/(x))dx.

Answer 1

clc
format long
clear
a=0;b=1;p=0;
n=4;m=0;
h=(b-a)/n;
for i=1:n-1
x=0;
x=a+h*i;
m=m+sin(x)/x;
end
h=(b-a)/(2*n);
for k=1:2:2*n
p=p+sin(a+k*h)/(a+k*h);
end
Sn=h/3*[1+4*p+2*m+sin(b)/b] %辛普森公式结果
result=sinint(1) %计算精确值
format short

结果：
Sn =
0.94608331088847
result =
0.94608307036718
>>