double_int.m

这是我最近学习的程序

function SS=double_int(fun,innlow,innhi,outlow,outhi)
%double_int		采用8阶多项式近似的二重数值积分
%						内积分变量为x，其积分区间上下限可以是变量y的表达式。
%fun				被积函数，以x,y为积分变量。
%innlow,innhi		内积分区间上下限，可以是标量，或函数文件名字符串
%outlow,outhi		外积分区间上下限，是标量。
y1=outlow;y2=outhi;x1=innlow;x2=innhi;f_p=fun;
SS=quad(@G_yi,y1,y2,[],[],x1,x2,f_p);