I am trying to solve this integral which has come up as part of some other work, but it is proving to be much harder than I had originally thought. For $0 < |a| \le 1$ being some constant, I am trying to solve the following
$$f_Z(z)=\lim_{z\to \infty} \int_0^z{\dfrac{\exp{\left(-\left(\frac{z-w}{2} + \frac{w}{2a^2}\right)\right)}}{\sqrt{\left(z - w\right)w}}}dw$$
which I reduced to $$f_Z(z)=\lim_{z\to \infty} \exp{\left(\frac{-z}{2}\right)}\int_0^z{\dfrac{\exp{\left(\frac{w\left(a^2 - 1\right)}{2a^2}\right)}}{\sqrt{\left(z - w\right)w}}}dw$$
Using the substitution $t=\sqrt{\left(z - w\right)}$ which gives $w=z - t^2$ and $\frac{dw}{dt}=-2t$. I was able to restate it as $$f_Z(z)=\lim_{z\to \infty}\exp{\left(\frac{-z}{2}\right)}\int_\sqrt{z}^0{\dfrac{-2\exp{\left(\frac{\left(z - t^2\right)\left(a^2 - 1\right)}{2a^2}\right)}}{\sqrt{\left(z - t^2\right)}}}dt$$ Using the substitution $t=\sqrt{z}\sin \theta$ and $\frac{dt}{d\theta}=\sqrt{z}\cos \theta$ and the trigonometric identity $\sin^2 \theta + \cos^2 \theta = 1$. I was able to restate it as. $$f_Z(z)=\lim_{z\to \infty}\exp{\left(\frac{-z}{2}\right)}\int_\frac{\pi}{2}^\pi{-2\exp{\left(z\cos^2 \theta \frac{\left(a^2 - 1\right)}{2a^2}\right)}}d\theta$$.
I am not quite sure how to proceed here. For all I know, trying to get rid of those square root parts, I may have made this worse, not only I am unsure about how to deal with the cosine in the power part of the $e$, I also have this negative sign to deal with now.
Any help would be much appreciated.
Edit: Using Laplace transform and then the inverse I was able to compute this to $$f_Z(z)=\dfrac{exp\left(\frac{-z}{4}\left(1 + \frac{1}{a^2}\right)\right)}{4a^2}J_0\left(\dfrac{z\sqrt{1-a}\sqrt{a+1}\sqrt{a^2 -1}}{2a}\right)$$ where $J_0$ is the Bessel function of first kind. The link for representation Wolfram Alpha Inverse Laplace Transform. I would love to see a derivation without the Laplace transform usage. :-)
For the integral $$I:=\int_0^z {\frac {\exp {(w\cdot \frac{a^2-1}{2a^2})}}{\sqrt {(z-w) w}} dw} $$ by completing the square in the denumenator and simplify we get $$\sqrt{\frac{1}{4} z-\frac{1}{4} z+zw-w^2}=\frac{1}{2} z\sqrt {1-(1-\frac{2w}{z})^2}$$ Now using the substitution $$\sin\theta = 1-\frac {2w}{z} $$ Therefore $ I $ reduces to $$ I=\int_{\frac {-\pi}{2}}^{\frac{\pi}{2}}{ e^{\alpha_z (1-\sin\theta)}d\theta}$$ where $\alpha_z=\frac{z}{2}\frac{a^2-1}{2a^2} $. Now, expand $\alpha_z (1-\sin\theta) $ in exponential Taylor series we get $$ I=\int_{\frac {-\pi}{2}}^{\frac {\pi}{2}}{\sum_n { \frac {\alpha_z^n}{n!}(1-\sin\theta)^n }d\theta}.$$ assuming the uniform convergence of the series then term by term integration yields that $$ I=\sum_n {\frac{\alpha_z^n}{n!}\cdot \int_{\frac{-\pi}{2}}^{\frac{\pi}{2}} {(1-\sin\theta)^n d\theta}} $$ Using binomial series we may write $$(1-\sin\theta)^n=\sum_{k=0}^{n}{(-1)^k \left ({\begin{array}{*{20} c} n \\ k \\ \end{array}}\right) \int_{\frac {-\pi}{2}}^{ \frac {\pi}{2 }}{\sin^{n-k}\theta d\theta}}.$$
It's well known that $\int\sin^m x$ given by:
If $ m$ even i.e., $ m =2n $ then $$ \int{\sin^{2n} xdx}= 2^{-2n} \left ({\begin {array}{*{20} c} 2n \\ n \\ \end {array}}\right) x+(-1)^{n} 2^{1-2n}\sum_{ k=0 }^{n-1}(-1)^k\left ({\begin {array}{*{20} c} 2n\\ k\\ \end{array}}\right) \frac {\sin (2n-2k) x}{2n-2k} $$ and for $ m $ odd; $ m=2n+1$ we have $$\int { \sin^{2n+1} xdx }=2^{-2n}(-1)^n \sum_{k=0}^n \left ({\begin {array}{*{20} c} 2n+1\\ k\end {array}}\right) \frac { \cos (2n+1-2k) x }{2n+1-2k} $$
For the both cases substitute the result in the original sum and simplify to get the result. It remains to calculate the limit as $ z→\infty $.