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);