clc,clear
[x,y,z]=meshgrid(linspace(-3,3,120));
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(gca,'Color','g')
set(p,'FaceColor','b')
grid on
daspect([1 1 1])
view(3)
camlight('right')
camlight('left')
camlight('headlight')
lighting phong
xlabel('MATLAB万岁')
ylabel('数学建模')
zlabel('相信自己')
title('亲爱的,证明自己')
text(1.2,0,0,'Just for you!','FontSize',14)