說明,用下面程式時,必須要先確定m的值.這兩個程式都有誤差估計.設n=2m+1(這是因為,辛普森公式要求的)梯形公式:h=π/2/2m=π/4m;a=0;b=pi/2;x=a:h:b;y=sin(x);z1=(y(1)+y(n))*h/2;z2=sum(y(2:n-1))*h;z=z1+z2,symstf=sin(t);intf=int(f,t,a,b)(計算精確值),Fs=double(intf),Juewucha8000=abs(z8000-Fs)(計算誤差)辛普森公式a=0;b=π/2;;h=(b-a)/(2*m);x=a:h:b;y=sin(x);z1=y(1)+y(2*m+1);z2=2*sum(y(2:2:2*m));z3=4*sum(y(3:2:2*m));z=(z1+z2+z3)*h/3,symst,f=exp((-t^2)/2)/(sqrt(2*pi));intf=int(f,t,a,b),Fs=double(intf);Juewucha=abs(z-Fs).
說明,用下面程式時,必須要先確定m的值.這兩個程式都有誤差估計.設n=2m+1(這是因為,辛普森公式要求的)梯形公式:h=π/2/2m=π/4m;a=0;b=pi/2;x=a:h:b;y=sin(x);z1=(y(1)+y(n))*h/2;z2=sum(y(2:n-1))*h;z=z1+z2,symstf=sin(t);intf=int(f,t,a,b)(計算精確值),Fs=double(intf),Juewucha8000=abs(z8000-Fs)(計算誤差)辛普森公式a=0;b=π/2;;h=(b-a)/(2*m);x=a:h:b;y=sin(x);z1=y(1)+y(2*m+1);z2=2*sum(y(2:2:2*m));z3=4*sum(y(3:2:2*m));z=(z1+z2+z3)*h/3,symst,f=exp((-t^2)/2)/(sqrt(2*pi));intf=int(f,t,a,b),Fs=double(intf);Juewucha=abs(z-Fs).