问题描述：
matlab 自适应辛普森
function s=Self_Adaptive_integral(a,b,tol,M)
%input:a--下限
% b--下限
% tol--the tolerance（容差）
% m--初始设置的步数
h=(b-a)/M;%步距
s=0;
for i=1:M
x=a+(i-1)*h;
y=a+i*h;
to=abs(simpson_integral(x,y,2)+simpson_integral(x,y,1))/10;
j=1;
while(to>=tol) %循环直到to
function s=Self_Adaptive_integral(a,b,tol,M)
%input:a--下限
% b--下限
% tol--the tolerance（容差）
% m--初始设置的步数
h=(b-a)/M;%步距
s=0;
for i=1:M
x=a+(i-1)*h;
y=a+i*h;
to=abs(simpson_integral(x,y,2)+simpson_integral(x,y,1))/10;
j=1;
while(to>=tol) %循环直到to
问题解答：
我来补答展开全文阅读