timflow.bessel.besselnumba.bessells_gauss_ho_qxqy_d1d2#
- timflow.bessel.besselnumba.bessells_gauss_ho_qxqy_d1d2(x, y, z1, z2, lab, order, d1, d2)#
Returns integral from d1 to d2 along real axis.
While strength is still Delta^order from -1 to +1.
implicit none integer, intent(in) :: order real(kind=8), intent(in) :: x,y,d1,d2 complex(kind=8), intent(in) :: z1,z2,lab complex(kind=8), dimension(0:2*order+1) :: qxqy, qxqyc integer :: n, m real(kind=8) :: xp, yp, dc, fac complex(kind=8) :: z1p,z2p,bigz1,bigz2