数值方法中要求用Matlab复合Simpson公式做此题(要求完整过程,急用,求高手!!!~~~)

Matlab复合Simpson公式(取n=4)计算积分∫0~1((sinx)/(x))dx.
2024-12-19 23:40:44
推荐回答(1个)
回答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
>>