MATLAB绘画三维心形

今天是520，单身狗只能用code来表达对matlab的钟爱了= =

clc;
const=0;
x=-5:0.05:5;y=-5:0.05:5;z=-5:0.05:5;
[x,y,z]=meshgrid(x,y,z);%生成采样网格
%笛卡尔心形函数
f=(x.^2 + (9/4)*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9/80)*y.^2.*z.^3-const;
p=patch(isosurface(x,y,z,f,0));   %创建补片图形对象,隐函数实现
set(p, 'FaceColor', 'red', 'EdgeColor', 'none');%设置颜色
daspect([1 1 1]);; %设置坐标轴的长宽高之比等于1
view(3);  %视角的控制
camlight;%光源位置
lighting phong;%光照模式