d/dx -1/2.cot(x)^2