polar plot rho=1+2*cos(phi)