Question

用matlab 画(x^2 + (9/4)y^2 + z^2 - 1)^3 - x^2z^3 - (9/80)y^2z^3 = 0的图应该怎么编代码？

Answer 1

[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');

Answer 2

[x,y,z]=meshgrid(linspace(-1.5,1.5,50),linspace(-.8,.8,50),linspace(-1.5,1.5,50));
isosurface(x,y,z,(x.^2 + (9/4)*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (9/80)*y.^2.*z.^3,0)
axis equal