一、二維圖形

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
         
 

實現效果: