目录
clc;
clear all;
close all;
clear;
[x,y,z] = meshgrid(linspace(-1.3,1.3));
val=(x.^2 + (9/4)*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (1/9)*y.^2.*z.^3;
p = patch(isosurface(x, y, z, val, 0));
isonormals(x, y, z, val, p);
set(p, 'FaceColor', 'red', 'EdgeColor', 'none');
axis equal;view(-10, 24);
colormap(hot);axis off;
camlight;rotate3d on;