smartsimpson.m

来自「Matlab语言常用算法程序集.rar 里面有很多有用的Matlab数值运」· M 代码 · 共 18 行

function q=SmartSimpson(f,a,b,eps)
if(nargin==3)
    eps=1.0e-4;
end;
e=5*eps;
q=SubSmartSimpson(f,a,b,e);

function q=SubSmartSimpson(f,a,b,eps)
QA=IntSimpson(f,a,b,1,eps);
QLeft=IntSimpson(f,a,(a+b)/2,1,eps);
QRight=IntSimpson(f,(a+b)/2,b,1,eps);

if(abs(QLeft+QRight-QA)<=eps)
    q=QA;
else
    q=SubSmartSimpson(f,a,(a+b)/2,eps)+SubSmartSimpson(f,(a+b)/2,b,eps);   %递归公式
end