clear;clc[x,y,z]=meshgrid(linspace(-1.3,1.3));v=(x.^2+(9*y.^2)./4+z.^2-1).^3-(x.^2).*(z.^3)-((9*y.^2).*z.^3)/80;isosurface(x,y,z,v,0);axis equal
