内联函数中怎样处理字母系数

sunyok

比如：对附件1中双重积分的处理——
[double_int.m]
function SS=double_int(fun,innlow,innhi,outlow,outhi)
y1=outlow;
y2=outhi;
x1=innlow;
x2=innhi;
f_p=fun;
SS=quad(@G_yi,y1,y2,[],[],x1,x2,f_p);

[G_yi.m]
function f=G_yi(y,x1,x2,f_p)
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(:);

[x_hi.m]
function f=x_hi(y)
f=y;

[Untitled.m]
ff=inline('x.^2+y.^2','x','y');      
SS=double_int(ff,3,@x_hi,3,2)

那对附件2中积分怎样处理呢？？？？？？

PS:我也不想重复发帖，只是刚才发第一次时有错误，抱歉。

[ 本帖最后由 sunyok 于 2008-3-12 10:59 编辑 ]

eight

原帖由 sunyok 于 2008-3-12 10:49 发表

                               
登录/注册后可看大图

比如：对附件1中双重积分的处理——
[double_int.m]
function SS=double_int(fun,innlow,innhi,outlow,outhi)
y1=outlow;
y2=outhi;
x1=innlow;
x2=innhi;
f_p=fun;
SS=quad(@G_yi,y1,y2,[],[],x1,x2,f_p);
...

发帖前，先看看置顶区所有帖子，特别是会员守则（对比一下你的标题，还有重复发帖）。如果不是见你的帖子内容较多、重写一次比较麻烦的话，我肯定删除。

[ 本帖最后由 eight 于 2008-3-12 10:55 编辑 ]

sigma665

不清楚，不过可以加个循环啊