matlab 7.0三维心形函数

% 二维心形图
figure(1)
ezplot('(x^2-1+y^2)^3-x^2*y^3 = 0')
axis([-2.5 2.5 -2 2])

% 三维心形图
figure(2)
n=100;
x=linspace(-3,3,n);
y=linspace(-3,3,n);
z=linspace(-3,3,n);
[X,Y,Z]=ndgrid(x,y,z);
F=((-(X.^2) .* (Z.^3) -(9/80).*(Y.^2).*(Z.^3)) + ((X.^2) + (9/4).* (Y.^2) + (Z.^2)-1).^3);
isosurface(F,0)
lighting phong
caxis
axis equal
colormap('flag');
view([55 34]);