matlab正余弦画心形图案

for b=0:0.08:20 
    clf;
    x=[0:0.01:1.8];
    y=x.^(2/3)+0.9*sqrt(3.3-x.^2).*sin(b.*pi.*x);
    plot(x,y,'r');
    hold on
    plot(-x,y,'r')
    plot(x*2/3,y*2/3,'b')
    plot(-x*2/3,y*2/3,'b')
    axis equal
    pause(0.0005)
    hold off
end