How can I handle BesselI(0,a)/BesselI(0,b) to prevent NaN in Fortran or MATLAB?
I expect a method of simplifying the ratio or a mathematical trick to handle the problem. I try a series representation of BesselI(0) but sometimes, it fails to converge.