polar plot r=cot(phi) * tan(phi) + phi / pi, -pi <= phi <= pi