代码如下
%%
clear;clc;
f = @(x,y,z) (x.^2+4.*y.^2+z.^2-1).^3-x.^2.*z.^3 - 0.4*y.^2.*z.^3;
fimplicit3(f);
axis equal;
%%
clear;clc;
figure
[x,y,z]=meshgrid(linspace(-3,3));
p=(x.^2+4*y.^2+z.^2-1).^3-x.^2.*z.^3-0.4*y.^2.*z.^3;
isosurface(x,y,z,p,0,x);
axis equal;
axis off;
view(0,45);
colormap([1 0 0]);
brighten(0.5);
camlight right;
lighting phong;
第一段代码效果如下:
第二段代码效果如下: