使用mesh函数画复杂公式的三维图问题

[小小书生]

仿真矩形载流线圈空间磁场分布时，由于公式复杂一些，出现错误，哪位大神看看这个公式该怎么仿真？

[小小书生]

仿真的时候总是出现矩阵不匹配的问题
代码：
x=linspace(-2, 2, 25); % 在x轴上取25点

y=linspace(-2, 2, 25); % 在y轴上取25点

[xx,yy]=meshgrid(x, y); % xx和yy都是21x21的矩阵

%[x,y]=meshgrid(0:0.25:4); %在[0,4pi]×[0,4pi]区域生成网格坐标
zz=((0.5+yy)/((0.5+yy).^2+0.04))*(((0.5+xx)/sqrt((0.5+xx).^2+(0.5+yy).^2+0.04))+((0.5-xx)/sqrt((0.5-xx).^2+(0.5+yy).^2+0.04)))+((0.5-yy)/((0.5-yy).^2+0.04))*(((0.5-xx)/sqrt((0.5-xx).^2+(0.5-yy).^2+0.04))+((0.5+xx)/sqrt((0.5+xx).^2+(0.5-yy).^2+0.04)))+((0.5-xx)/((0.5-xx).^2+0.04))*(((0.5+yy)/sqrt((0.5-xx).^2+(0.5+yy).^2+0.04))+((0.5-y)/sqrt((0.5-xx).^2+(0.5-yy).^2+0.04)))+((0.5+xx)/((0.5+xx).^2+0.04))*(((0.5-yy)/sqrt((0.5+xx).^2+(0.5-yy).^2+0.04))+((0.5+yy)/sqrt((0.5+xx).^2+(0.5+yy).^2+0.04)));
%z=x^2+y^2;
%a1=(1+y)/((1+y)^2+0.04);
%{
a2=(1+x)/sqrt((1+x)^2+(1+y)^2+0.04);
a3=(1-x)/sqrt((1-x)^2+(1+y)^2+0.04);
b1=(1-y)/((1-y)^2+0.04);
b2=(1-x)/sqrt((1-x)^2+(1-y)^2+0.04);
b3=(1+x)/sqrt((1+x)^2+(1-y)^2+0.04);
c1=(1-x)/((1-x)^2+0.04);
c2=(1+y)/sqrt((1-x)^2+(1+y)^2+0.04);
c3=(1-y)/sqrt((1-x)^2+(1-y)^2+0.04);
d1=(1+x)/((1+x)^2+0.04);
d2=(1-y)/sqrt((1+x)^2+(1-y)^2+0.04);
d3=(1+y)/sqrt((1+x)^2+(1+y)^2+0.04);
%}
%z=a1*(a2+a3)+b1*(b2+b3)+c1*(c2+c3)+d1*(d2+d3);
mesh(xx,yy,zz);
axis([0 4*pi 0 4*pi -2.5 1]);

[小小书生]

mesh是不是不能够画画出类似z=(x+y)/((x+y)^2+a)这样的三维图呢？