Using method of variation of parameters , solve the following differential equation $$\frac{d^2y}{dx^2}+y=\csc x.$$
I know that, if, the given differential equation is of the form, $\frac{d^2y}{dx^2}+P\frac{dy}{dx}+Qy=g(x),$ and the two fundamental solution of $\frac{d^2y}{dx^2}+P\frac{dy}{dx}+Qy=0,$ are $u$ and $v$ then, the particular integral of $$\frac{d^2y}{dx^2}+P\frac{dy}{dx}+Qy=g(x),$$ is given by $$y=-u\int\frac{vg(x)}{w(u,v)}dx +v\int\frac{ug(x)}{w(u,v)}dx\tag 1,$$ where $w(u,v)$ is the wronskian of the two functions $u,v$ i.e $w(u,v)=\begin{vmatrix} u & v\\ u'&v'\end{vmatrix}.$
However, the problem is, in this case, the complementary function, or the fundamental set of solutions of $\frac{d^2y}{dx^2}+y=0$ are $u=e^{ix}$ and $v=e^{-ix}.$
Here, I am experiencing a couple of issues while applying the above mentioned formula $(1)$ of particular integrals.
All in all, I have to calculate the wronskian firstly. But then, I have only worked with calculus problems in real functions. This is the first time, I am trying to differentiate a complex function as, $u=e^{ix}.$
I am unsure whether, $u'=ie^{ix}.$ I tried emulating the same thing that we do in case of real functions.
Then, comes the next hurdle similar to these. How do we integrate $$y=-u\int\frac{vg(x)}{w(u,v)}dx +v\int\frac{ug(x)}{w(u,v)}dx,$$ when, $u,v$ are complex functions?
Say, for example, $\int\frac{vg(x)}{w(u,v)}=\int\frac{e^{-ix}\csc(x)}{w(e^{ix},e^{-ix})}$ I have no whatsoever idea about this whole process, since never did I perform an integration on complex functions.
Any clarification regarding this apparent issues will be highly appreciated.