I have a pretty ferocious integral to solve, and since it doesn't seem I'll be able to do much analytically, I've taken to Mathematica to get some information. Mainly, I want to see if there are any constant terms in the following integral:
$$ I = \int_{r}^{\infty} r_0^{-5/2} W_{-i\alpha'/2, \nu}\left(\frac{-ir_0^2\omega'}{L}\right) W_{-i\alpha/2, \nu}\left(\frac{-ir_0^2\omega}{L}\right) J_{\nu}\left(-i\hat{k}r_0\right) dr_0$$
Where $W$ is WhittakerW
, $J$ is BesselJ
, and everything besides $r/r_0$ is a constant.
So far, I've simply been doing NIntegrate
and breaking the integral up over several pieces. I plotted the function on a log plot, and it seems after about $r \approx 10$ it drops off significantly; as such, the integral from 10 to infinity seems to be zero. After this, I've just plotted the integral and tried to look at dependence using LogPlots. Does anyone have a better method to deal with Numerical Integrals like this to get a bit more information?
Here is my mathematica code:
g[r_, \[Alpha]_, v_, w_] := WhittakerW[-I*\[Alpha]/2, v, -I*r^2*w]
Plot[Log[(Re/Im)[
r^(-5/2)*g[r, 10, 2.1, .01 ]*g[r, 15, 2.1, .01]*
BesselJ[3.2, 5*r]]], {r, 1, 100}, PlotLabel -> "Log(I2) vs r0"]
% see that the function being integrated goes to zero for both Real and Imaginary parts
I2bh[r_? NumericQ, v_] :=
Re[NIntegrate[
r0^(-5/2)*g[r0, 10, 2.1, .01 ]*g[r0, 15, 2.1, .01]*
BesselJ[3.2, 5*r0], {r0, r, 15}, WorkingPrecision -> 15]]
I've done some more plotting and probing of the features of the integral, but for the purposes of this question I suppose what I have here is all I'm asking about