整理的代码如下：function plot2heart
[x,y,z]=meshgrid(linspace(-3,3));
val=(x.^2 + (9/4)*y.^2 + z.^2 - 1).^3 - x.^2.*z.^3 - (1/9)*y.^2.*z.^3;
isosurface(x,y,z,val,0);
axis equal;
colormap([1 0 0]);
brighten(0.5);      
hold on
[x1,y1,z1]=meshgrid(linspace(-4,2));
val1=((x1+2).^2 + (9/4)*(y1).^2 + (z1+0.5).^2 - 1).^3 - (x1+2).^2.*(z1+0.5).^3 - (1/9)*(y1).^2.*(z1+0.5).^3;
isosurface(x1,y1,z1,val1,0);
axis equal;view(30,15);
colormap([1 0 0]);
hold on
t=linspace(-pi,pi,25);
y=linspace(-1,1,30);
[T,Y]=meshgrid(t,y);
X=3.5*sin(T);
Z=2*cos(T);
mesh(X,Y,Z);
axis equal
end