parametric plot (-csc(a)^2 (-cos(a)+a csc(a)), cos(a)) from a=0.01*pi to a=0.99*pi