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中可运行后可得结果：