[交流] Matlab自适应数值积分，用了递归，最后怎样清理全局变量？

求解问题是：对f在区间[a,b]上实现自适应剖分，并在每个小区间上用三点高斯数值积分公式。自适应剖分的条件可以参看下面的if条件.要求函数输出最终数值积分值，以及所有剖分得到的小区间及每个小区间上的三点高斯积分值 以下是m文件源码： function [int, abt] = gadap(a, b, f, p, tol) % format long global p global abt xc=(a+b)/2;int=gauss3(f,a,b); t1=gauss3(f,a,xc);t2=gauss3(f,xc,b);  %%[a,xc]、[xc,b]上f的三点高斯积分值 if (abs(int-t1-t2))>tol*max(abs(int),abs(t1)+abs(t2))   %%判定条件，若满足则将[a,b]等分为两部分     int=gadap(a,xc,f,p,tol)+gadap(xc,b,f,p,tol); else     p=p+1;     abt(p,=[a,b,int];     if b==1         clear global p     %%清除全局变量p     end end function I = gauss3(f,a,b) %%三点高斯数值积分 ta = (b-a)/2; tb = (a+b)/2; I = ta*(0.55555556*subs(sym(f),findsym(sym(f)),ta*0.7745967+tb)+...             0.55555556*subs(sym(f),findsym(sym(f)),-ta*0.7745967+tb)+...             0.88888889*subs(sym(f),findsym(sym(f)),tb)); 现在问题是这代码只能清除全局变量p，若是还清除abt，则函数根本返回不了abt，只能每次在command window里用clear global来清除，各位有没有更好的办法，如果不用到全局变量可以吗

回复此楼