g_yi.m

来自「这是本人收集的一些程序源代码」· M 代码 · 共 14 行

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