polar plot r=1+cos(theta)