www.pudn.com > matlab8.zip > 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=quad8('G_yi',y1,y2,[],[],x1,x2,f_p);