polar plot -ln(tan((phi)/(2)))