var('x,y,z')
heart = implicit_plot3d((x^2 + 9*y^2/4 + z^2 - 1)^3 == x^2*z^3 + 9*y^2*z^3/80,(x,-1.5,1.5),(y,-1,1),(z,-1,1.5),color='red',plot_points=120,axes=False,bounding_box=False)
heart.show()
saveheart = file("heart.x3d", "w")
saveheart.write(heart.x3d())
saveheart.close()
