% comet plot of the trifolium: (x^2+y^2)(x^2+y^2+x)=4xy^2
% F. Waleffe, UW-Madison, 09/20/01

figure(1), clf

t=[0:0.001:1]*2*pi;
r=(4*(sin(t)).^2-1).*cos(t);
x=r.*cos(t); y=r.*sin(t);

comet(x,y,0.1)
     
% plot(x,y), grid, axis image