Numerical Computation and FORTRAN Program on the Modified Bessel Functions I and K

This paper has given a computation method on the modified Bessel functions I_n(x) and K_n(x) for positive real x and integern, and a FORTRAN program used in microcomputer. In the range of the floating point figure of the microcomputer, I_n(x) and K_n(x) havc accuracy of 10 digits.