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)
在基础上作了一点修改:
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);