> retstart:

> plot((1+cos(3*t)+sin(3*t)^2),t=0..2*Pi,coords=polar,axes=none,scaling=constrained);

>