急！急！如何画程序中三维隐函数的图像

[lv19860312]

我编了一个程序，要画出其中隐函数的图像，但总是出错，请高手帮我看看！程序如下:
syms x y z w0 lad m e o %e为真空中介电常数，o为光的角频率
x=-0.005:0.0001:0.005;
y=-0.005:0.0001:0.005;
z=-100:2:100;
[x,y,z]=meshgrid(x,y,z);
m=1;
w0=2*10^-3;
lad=632*10^-9;
e=8.854*10^-12;
k=2*pi/lad;
o=3*10^8*k;
w=sqrt(w0^2+4*z.^2/k^2*w0^2);
R=z+k^2*w0^4./(4*z);
u=(w0./w).*(sqrt(x.^2+y.^2)./w).^abs(m).*exp(-(x.^2+y.^2)./w.^2).*exp(i*m*atan(y./x))...
.*exp(0.5*i*k*(x.^2+y.^2)./R).*exp(-i*(abs(m)+1)*atan(2*z/(k*w0^2)));
Mz=m/(pi^2+z.^2*lad^2).*exp(-2*(x.^2+y.^2)*pi^2/w0^2./(pi^2+z.^2*lad^2)).*pi^(2*abs(m)+2)...
*((x.^2+y.^2).^(1/2)/(w0^2*(pi^2+z.^2*lad^2)).^(1/2)).^(2*abs(m))*e*o;
p=patch(isosurface(x,y,z,Mz,0));
set(p, 'FaceColor', 'red', 'EdgeColor', 'none');
daspect([1 1 1])
view(3)
camlight; lighting phong
程序能运行但是就是出来的图却只有坐标轴。请高手帮我分析分析是什么问题，怎么改？