以前中二的时候，特别喜欢鼓捣各种乱七八糟的东西，这是翻出来的Matlab绘制心形的代码

f=inline('(x.^2+(9/4)*y.^2+z.^2-1).^3-(x.^2).*(z.^3)-(9/80)*(y.^2).*(z.^3)');
gd=100;
xlimit=[-1.5 1.5];
ylimit=xlimit;zlimit=xlimit;
x=linspace(xlimit(1),xlimit(2),gd);
y=linspace(ylimit(1),ylimit(2),gd);
z=linspace(zlimit(1),zlimit(2),gd);
[x,y,z]=meshgrid(x,y,z);val=f(x,y,z);
[f v]=isosurface(x,y,z,val,0);
 newplot;
p=patch('Faces',f,'Vertices',v,'CData',v(:,3),'facecolor','w','EdgeColor','flat');
isonormals(x,y,z,val,p);
  view([29 7]);
 grid on;
 colormap(autumn);