变步长复化辛普森公式计算积分

变步长复化辛普森公式计算积分变步长复化辛普森公式计算积分变步长复化辛普森公式计算积分V:1.0精细整理，仅供参考变步长复化辛普森公式计算积分日期：20xx年X月2.编写用变步长复化辛普森公式计算积分的程序。用上面编写的程序计算下列积分并分析计算结果（1）（2）（3）程序：functionS=bianfuhuasimpson(fx,a,b,eps,M)%变步长复合simpson求积公式%调用方式：S=fuhuasimpson(@fx,a,b,epsilon)%fx--求积函数（函数文件）%a,b--求积区间%eps--计算精度%M--最大允许输出划分数n=1;h=(b-a)/n;T1=h*(feval(fx,a)-feval(fx,b))/2;Hn=h*feval(fx,(a+b)/2);S1=(T1+2*Hn)/3;n=2*n;%最好与倒数第三行保持一致（变步长）whilen<=MT2=(T1+Hn)/2;Hn=0;h=(b-a)/n;forj=1:nx(j)=a+(j-1/2)*h;y(j)=feval(fx,x(j));Hn=Hn+y(j);endHn=h*Hn;S2=(T2+2*Hn)/3;fprintf('n=%2dS2=%S2-S1=%\n',n,S2,abs(S2-S1));ifabs(S2-S1)<epsbreak;elseT1=T2;S1=S2;n=2*n;endendS=S2;%达到下列条件之一，则运算终止:%(1).abs(S2-S1)<eps%(2).下一次的n>M%输入1：S=bianfuhuasimpson(inline('sqrt(x)*cos(x)'),0,pi,10e-6,2000)%输入2：S=bianfuhuasimpson(inline('2*x^2*cos(x^2)'),0,sqrt(pi),10e-6,2000)%输入3：S=bianfuhuasimpson(inline('sqrt(x)'),0,1,10e-6,2000)输出结果：（1）S=bianfuhuasimpson(inline('sqrt(x)*cos(x)'),0,pi,10e-6,2000)n=2S2=S2-S1=n=4S2=S2-S1=n=8S2=S2-S1=n=16S2=S2-S1=n=32S2=S2-S1=n=64S2=S2-S1=n=128S2=S2-S1=n=256S2=S2-S1=n=512S2=S2-S1=n=1024S2=S2-S1=S=（2）S=bianfuhuasimpson(inline('2*x^2*cos(x^2)'),0,sqrt(pi),10e-6,2000)n=2S2=S2-S1=n=4S2=S2-S1=n=8S2=S2-S1=n=16S2=S2-S1=n=32S2=S2-S1=n=64S2=S2-S1=n=128S2=S2-S1=n=256S2=S2-S1=n=512S2=S2-S1=n=1024S2=S2-S1=S=（3）S=bianfuhuasimpson(inline('sqrt(x)'),0,1,10e-6,2000)n=2S2=S2-S1=n=4S2=S2-S1=n=8S2=S2-S1=n=16S2=S2-S1=n=32S2=S2-S1=n=64S2=S2-S1=n=128