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?????????????