程式碼如下：

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

• ezimplot3(f)畫函式f(X,Y,Z)= 0 在-2*pi< X < 2* pi, -2* pi < Y < 2* pi, -2* pi < Z < 2* pi上的圖形
• ezimplot3(f, [A,B])畫函式f(X,Y,Z)= 0 在A< X < B, A < Y < B, A < Z < B上的圖形
• ezimplot3(f, [XMIN,XMAX,YMIN,YMAX,ZMIN,ZMAX])畫函式f(X,Y,Z)= 0 在XMIN< X < XMAX, YMIN < Y < YMAX, ZMIN < Z < ZMAX上的圖形
  
  
  

f = @(x,y,z) x*y*z*log(1+x^2+y^2+z^2)-10;
ezimplot3(f,[-10,10]);  % [-10, 10] 表示圖形範圍x、y、z都在區間[-10, 10]

• ezimplot3和方法一本質上完全相同。即ezimplot3實際上也是基於isosurface+ patch + isonormals的實現
• ezimplot3與方法一的圖形視覺效果相同，唯一的區別是，ezimplot3的使用了0.7的透明度：alpha(0.7)
• ezimplot3在方法一基礎上增加了一些外包功能，如：允許函式控制代碼f是非向量化的函式（即函式定義無需.*  ./  .^），這在ezimplot3內部會自動呼叫vectorize實現函式向量化。另外，ezimplot3可以在呼叫的時候方便的設定座標範圍。
  

• 常見問題：很多人在使用以上方法後，經常出現的問題是程式碼沒有任何錯誤，程式可以執行，就是出來的圖形只有一個空座標軸，看不到圖形。
  

• 問題分析：出現這種問題的原因是圖形的顯示區域沒設對。比如，我們上述三種方法都是在x為-10到10的範圍內，如果你設的範圍內本身就沒有圖形，那當然就看不到圖形了。
• 解決辦法：把圖形顯示範圍重新設定對即可，如果不知道圖形的大致範圍，就手工多改幾次，直到看到圖形為止
  

• 方法一，圖形範圍是在第2句的meshgrid函式決定的，meshgrid裡給出的x、y、z範圍就是最終畫圖範圍，修改meshgrid語句即可。
• 方法二（Mupad），x =-10..10, y = -10..10, z = -10..10是表示顯示範圍，修改這裡即可。
• 方法三，用ezimplot3(f,[A, B]) ezimplot3(f, [XMIN,XMAX,YMIN,YMAX,ZMIN,ZMAX])兩種方式控制圖形顯示範圍。
  

[x,y,z] = meshgrid(linspace(-2,2));
v = x.*y.*z.*exp(-(x.^2+y.^2+z.^2));
xslice = 0; yslice = 1; zslice = 1;
slice(x,y,z,v,xslice,yslice,zslice)
xlabel('x'); ylabel('y'); zlabel('z');
colormap hsv