step=0.2;
[x y z]=meshgrid(-2:step:2, -2:step:2, -2:step:2);
F=(x.^2+(9*y.^2)./4+z.^2-1).^3-((9*y.^2).*(z.^3))./80-(x.^2).*(z.^3);
isosurface(x,y,z,F,0);
lighting phong; axis off; axis equal; colormap('flag'); view([30 30]);
clear all
[x,y,z]=meshgrid(linspace(-2,2,80));
f=(x.^2+(9*y.^2)./4+z.^2-1).^3-((9*y.^2).*(z.^3))./80-(x.^2).*(z.^3);
p=patch(isosurface(x,y,z,f,0));
set(p,'FaceColor','r')
grid on
daspect([1 1 1])
view(3)
camlight('right')
camlight('left')
camlight('headlight')
lighting phong
title('Heart of Math')