第18章 chapter 18.2（1 / 1）

为了表达我的爱意，特别粘贴一段Matlab代码，有条件的童鞋跑跑看。

[x,y,z]=meshgrid(linspace(-1.5,1.5));

f=@(x,y,z)(x.^2 + (9/4)*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9/80)*y.^2.*z.^3;

p=patch(isosurface(x,y,z,f(x,y,z),0)); view(3);axis equal

set(p,'FaceVertexCData',jet(size(get(p,'faces'),1)) ,'FaceColor', 'flat', 'EdgeColor', 'none');

好了，终于看我墨迹完了，看完洗洗睡吧，晚安