In my case instead of J_n(ax)^2 I have J_n(ax)' *J_n(ax) , where ' is the complex conjugate.
Basically, my first type Bessel functions are complex and I don't know
if the formula I found holds anymore.