untitledfuhua2.m

数值计算实践是我们必备的一门学科！这是我在学习之余利用MATLAB编写的程序

%复化辛普森;
Iexact=0.78539816325;   %取精确值;
a=0;b=1;
fprintf('\n Simpson 法计算积分\n')
fprintf('\n 精确值Iexact=0.78539816325;\n');
fprintf('  n          1        Error\n');
n=1;
for k=1:9;n=2*n;
    h=(b-a)/n;i=1:n+1;
    x=a+(i-1)*h;f=1./(1+x.^2);
    I=(h/3)*(f(1)+4*sum(f(2:2:n))+2*sum(f(3:2:n))+f(n+1));
    fprintf('%3.0f %12.9f %12.9f\n',n,I,Iexact-I);
end