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