g_yi.m

matlab mfiles实例分析

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