4744.sagews

var('x,y,z')

\newcommand{\Bold}[1]{\mathbf{#1}}\left(x, y, z\right)

p = (x^2+9/4*y^2+z^2-1)^3-x^2*z^3-9/(80)*y^2*z^3
r = 1.5

implicit_plot3d(p, (x, -r, r), (y, -r, r), (z, -r, r), plot_points=200, color='red', smooth=True, axes=True).show(viewer='tachyon')