[求助]matlab求积分的问题？

tonyboy005

我是MATLAB新手，今天写了下面一个积分的测试程序，却总是出不来结果，请高手赐教


x = 0 : 0.1 : 3.2;
y = -1.6 : 0.1 : 1.6;
[X,Y] = meshgrid(x,y);
func = inline('L/((X-L).^2+Y.^2+eps)','L');
Z = quad(func,0,1.5);
mesh(Z);

zjliu

warning off
syms x y x1 x2;
F = int(1/(x^2+y^2),x,x1,x2);
ff=inline(char(F),'x1','x2','y');
x = 0:0.1:3.2;
y = -1.6:0.1:1.6;
m=length(x);
n=length(y);
for p=1:m;
for q=1:n;
Fd(p,q)=ff(-x(p),-x(p)+1.5,y(q));
end
end
Fd(17,17)=nan;
[X,Y]=meshgrid(x,y);
mesh(X,Y,Fd)

tonyboy005

在基础上作了一点修改：
warning off
syms x y r1 r2 l;
F = int(1/((x-l)^2+y^2)^1.5,l,r1,r2);
func=inline(char(F),'x','y','r1','r2');
x = 0:0.1:3.2;
y = -1.6:0.1:1.6;
m=length(x);
n=length(y);
for p=1:m;
for q=1:n;
Fd(p,q)=func(x(p),y(q),0,1.5);
end
end
[X,Y]=meshgrid(x,y);
mesh(X,Y,Fd);

但不知为什么Fd(0,17)（即x=y=0时）为无限大啊？这个值应当是L的-2次方从0到1.5的积分值啊？

tonyboy005

搞定，谢谢