自适应Simpson法

function [qsrmerr]=naadapt（fnameabe）
srm=zeros（306）;
it=0;done=1;m=1;
sr=simpson（fnameabe）;
srm（11:6）=sr;
state=it;
while（state==it）
    n=m;
    for j=n:-1:1
        p=j;
        sr0=srm（p:）;
        err=sr0（5）;e=sr0（6）;
        if e<=err
            state=done;
            a=sr0（1）;;b=sr0（2）;err=sr0（5）;e=sr0（6）;
            c=0.5*（a+b）;e2=0.5*e;
            sr1=simpson（fnameace2）;
            sr2=simpson（fnamecbe2）;
            err=abs（sr0（3）-sr1（3）-sr2（3））/10;
            if err