如何绘制普林格尔形状的函数

Question

我想在 Octave 中重现品客形状。 请查看附件 - 这是在 Mathematica 中使用“区域”函数完成的。 在 Octave 中也可以切割圆形吗？到目前为止，我已经得到了这个（下面的代码），但它具有矩形形状/平台，而不是像 Mathematica 中那样的圆形（如下图）。

a = 10
b = 10
x = [-10:.1:10];
y = [-10:.1:10];
[xx, yy] = meshgrid (x, y);
z = xx.^2/a.^2 - yy.^2/b.^2
h = surf(xx,yy,z);
colormap hsv;
set(h,'linestyle','none');