matlab 两种自适应数值积分算法

% Integral function f from a0 to b0 % tol0: Total Error limit % flag:0-->Trapezoid ; 1-->Simpson

function [intPointNumber] = AdaptiveIntegral（ fa0b0tol0flag ）
% Integral function f from a0 to b0
% tol0: Total Error limit
% flag:0-->Trapezoid ; 1-->Simpson
int=0;n=1;a（1）=a0;b（1）=b0;tol（1）=tol0;app（1）=trapezoid_or_Simpson（fa（1）b（1）flag）;PointNumber=0;
while n>0                                   %n是当前列表结束的位置
    c=（a（n）+b（n））/2;oldapp=app（n）;
    app（n）=trapezoid_or_Simpson（fa（n）cflag）;app（n+1）=trapezoid_or_Simpson（fcb（n）flag）;
    if abs（oldapp-（app（n）+app（n+1）））<3*tol（n）
        int=int+app（n）+ap