[x,y,z]=meshgrid(-1:.1:1, -1:.1:1, -1:.1:1);
a=(1+sqrt(5))/2;
b1=0;
b2=1;
b=b1+b2;
f=4*(a^2*x.^2-y.^2).*(a^2*y.^2-z.^2).*(a^2*z.^2-x.^2)-b^2*(1+2*a)*(x.^2+y.^2+z.^2-b^2).^2;
const=0;
surfcv(x,y,z,f,const)
shading interp
axis off
colormap(winter)
light
lighting phong
