x=linspace(-50, 50, 50); % 在x軸上取50點
y=linspace(-25, 25, 25); % 在y軸上取25點
[xx,yy]=meshgrid(x, y); % xx和yy都是矩陣
zz=8000-2.*xx.*xx-5.*yy.*yy; % 計算函數值，zz也是21x21的矩陣
surf(xx, yy, zz); % 畫出立體曲面圖
colorbar; %如下圖，右邊那個色卡
可按回車
contour(xx,yy,zz)
colorbar
可按回車
surfc(xx,yy,zz)
colorbar
可按回車
在octave-3.6.4-gcc4.6.2中可運行后可得結果：
