x=linspace(0,5,100);y1=zeros(1,100);y2=zeros(1,100);for i = 1:100 y1(i)=2/sqrt(pi)*Simpson(g,0,x(i),1000); %幸普森公式 y2(i)=2/sqrt(pi)*trapezoidal(g,0,x(i),1000); %复合梯形公式 end plot(x,y1,'r-',x,y2,'b.');legend('辛普森公式','复合梯形公...
a=1;b=10;eps=0.00001;DDTraprl(fun,a,b,eps)ans = 14.0258 代码见附件。