DSolve[Derivative[2][y][x] == Tan[x]^2*y[x] + Csc[x]*Sec[x]*(1 + 3*Sin[x]^2)*Derivative[1][y][x],y[x],x]