polar plot r=3/(1+cos(α))