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