title: Integration of Bessel functions in C++/Fortran
body: <p>How can I integrate an equation including bessel functions numerically from "0" to "infinity" in Fortran or/and C? 
  I did in matlab, but it's not true for larger inputs and after a specific values , the bessel functions give completely wrong results(there is a restriction in Matlab)</p>