r(t) = 4*cos(t/3)^3

plot( r, 0, 4*pi, 200, "polar", "red", "frame:#00cc00" )