xt=@(t) cos(2*t).*cos(t);
yt=@(t) cos(2*t).*sin(t);
fplot(xt,yt)