DSolve[Derivative[2][y][x] == -(Csc[x]^2*(c + a*Cos[x]^2 + b*Sin[x]^2)*y[x]),y[x],x]