Skip to main content
3 of 4
added 708 characters in body

Nasty integral advice

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