求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:Mx=a+(i-1)*h;y=a+i*h;to=abs(simpson_integral(x,y,

求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:Mx=a+(i-1)*h;y=a+i*h;to=abs(simpson_integral(x,y,
求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

求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:Mx=a+(i-1)*h;y=a+i*h;to=abs(simpson_integral(x,y,
while(to>=tol) %循环直到to