DSolve[Derivative[2][y][x] == Csc[x]^2*y[x] - Cot[x]*Derivative[1][y][x],y[x],x]