文章来自于MATLAB论坛，见此链接：http://www.ilovematlab.cn/thread-264471-1-1.html，感谢原作者winner245的辛勤总结！

实现原理：先定义3元显函数v=f(x, y, z), 则 v = 0 定义的等值面就是z= g(x,y)的3D曲面。利用isosurface函数获取v=0 的等值面，将得到的等值面直接输入给patch函数，得出patch句柄p，并画出patch曲面的平面视角图形。对p用isonormals函数设置曲面顶点数据的法线，最后设置颜色、亮度、3D视角，得到3D曲面。

代码如下：

f = @(x,y,z) x.*y.*z.*log(1+x.^2+y.^2+z.^2)-10;      % 函数表达式
[x,y,z] = meshgrid(-10:.2:10,-10:.2:10,-10:.2:10);       % 画图范围
v = f(x,y,z);
h = patch(isosurface(x,y,z,v,0)); 
isonormals(x,y,z,v,h)              
set(h,'FaceColor','r','EdgeColor','none');
xlabel('x');ylabel('y');zlabel('z'); 
alpha(1)   
grid on; view([1,1,1]); axis equal; camlight; lighting gouraud

代码说明：

• alpha函数用于设置patch曲面的透明度（可以是0~1任意数值），1 表示不透明，0 表示最大透明度。如果想设置透明度为0.7，可以修改alpha(1)为alpha(0.7)。
  

• 使用此代码解决特定问题时，只需将第1行的函数表达式替换为特定问题的函数表达式，将第2行数据（x、y、z）范围换成合适的范围，后续代码无需任何变动。
  

得到图形：

plot(plot::Implicit3d(x*y*z*ln(1+x^2+y^2+z^2)-10,
x = -10..10, y = -10..10, z = -10..10), Scaling = Constrained)

得到如下图形：

解决办法三：第三方工具包ezimplot3

http://www.mathworks.com/matlabcentral/fileexchange/23623-ezimplot3-implicit-3d-functions-plotter