with(plots):
> implicitplot3d(x^2+(y-abs(x)^(2/3))^2+5*(abs(z*(1-y/5+x^2/10+y^2/10)))^3=1, x=-1..1, y=-1..5/3, z=-1..1, scaling=constrained, numpoints=10000, axes=none, thickness=3, color=red, style=patchnogrid, light=[30,-30,0.75,0,0], orientation=[-100,20]);