DSolve[Derivative[2][y][x] == -(((e + d*Cos[x] + c*Cos[x]^2)*Csc[x]^2*y[x])/a) - (b*Cot[x]*Derivative[1][y][x])/a,y[x],x]