用MeteoInfoLab绘制三维心。

a = linspace(-3, 3, 100)
x,y,z = meshgrid(a, a, a)
p = (x**2+(9/4.)*y**2+z**2-1)**3-x**2*z**3-(9/80.)*y**2*z**3
lighting(position=[0,0,-1,0])
isosurface(a, a, a, p, 0, facecolor='r', edgecolor=None, nthread=4)