[x,y,z]=meshgrid(linspace(-3,3));
v=(x.^2+(9/4).*y.^2+z.^2 - 1).^3-x.^(2).*z.^3-(9/80).*y.^2.*z.^3;
isosurface(x,y,z,v,0);
set(patch(isosurface(x,y,z,v,0)),'facecolor','r','edgecolor','none','facelighting','gouraud')
axis equal;
v=(x.^2+(9/4).*y.^2+z.^2 - 1).^3-x.^(2).*z.^3-(9/80).*y.^2.*z.^3;
isosurface(x,y,z,v,0);
set(patch(isosurface(x,y,z,v,0)),'facecolor','r','edgecolor','none','facelighting','gouraud')
axis equal;
grid on;