function f=G_yi(y,x1,x2,f_p)
%G_yi			實現"內重"積分的函數文件，專供double_int.m調用。
%y				"外重"積分的積分變量
%x1,x2			內積分區間的下、上限。
%				它們由double_int.m中quad8指令以參數方式傳遞過來。
%f_p			被積原函數，它由double_int.m中quad8指令以參數方式傳遞過來。
y=y(:);n=length(y);
if isstr(x1)==1;xx1=feval(x1,y);else xx1=x1*ones(size(y));end
if isstr(x2)==1;xx2=feval(x2,y);else xx2=x2*ones(size(y));end
for i=1:n
	f(i)=quad8(f_p,xx1(i),xx2(i),[],[],y(i));
end
f=f(:);