% 1.核心程序%
%designed by CONST from BUAA%
% 可根据电脑性能更改参数
th=linspace(-pi,pi,100);
ph=linspace(0,2*pi,100);
r=linspace(0,1.5,1000);
[t,p]=meshgrid(th,ph);
f=zeros(1,1000);
sf=zeros(100,100);
for i=1:100
for j=1:100
x=sin(ph(i))*cos(th(j));
y=cos(ph(i))*cos(th(j));
z=sin(th(j));
f=((x*r).^2 + (9/4)*(y*r).^2 + (z*r).^2 - 1).^3 - (x*r).^2.*(z*r).^3 - (9/80)*(y*r).^2.*(z*r).^3; [minf,ind]=min(abs(f));
sf(i,j)=r(ind);
end
end
[x,y,z]=sph2cart(p,t,sf);
% 2.mesh红心%
mesh(x,y,z)
axis equal
Colormap([1,0,0])
% 3.surf红心%
surf(x,y,z)