maltab实现心形图形

一、二维图形

ezplot('x^2+(y-(x^2)^(1/3))^2-1',[-2,2])

实现效果 

二、三维图形

function fun1_2()
n=100;
x=linspace(-3,3,n); %nirv
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([60 30]); %视角
end
         

实现效果: