Re:Heart

on 5 Oct 2022

t= 0:pi/2048:2*pi;

x = 16*sin(t).^3;

y = 14*cos(t) - 5*cos(2*t) - 2*cos(3*t) - cos(4*t);

m = min(y);

hold on

axis equal off

for i=0:1/100:1

p(x*i,y*i,i);

p(x*i,-y*i+2*m,i);

p(y*i-m,x*i+m,i);

p(-y*i+m,x*i+m,i);

% p(0.5*x*i,0.5*y*i-10,i);

% p(0.5*x*i,0.5*(-y*i+2*m)-10,i);

% p(0.5*(y*i-m),0.5*(x*i+m)-10,i);

% p(0.5*(-y*i+m),0.5*(x*i+m)-10,i);

end

hold off;

function p(x,y,i)

plot(x,y, 'Color', [i,i,0]);

end