simpson.m

该程序实现自适应的复化simpson积分。 含有三个M文件： main naadapt simpson

function z=simpson(func,a,b,e)
%%本函数将实现复化simpson数值积分：
%%输入：
%%func:积分函数
%%a:积分下限
%%b:积分上限
%%e:积分精度
%%输出：
%%z:积分过程
h=b-a;
c=0.5*(a+b);
f=feval(func,[a c b]); 
s=h*(f(1)+4*f(2)+f(3))/6; %%计算simpson积分
s2=s;
err=e;
z=[ a b s s2 err e];