Use new subroutine to calculate the modified Bessel functions of
the first kind with integer order and real argument. The subroutine
is based on http://www.netlib.org/toms/644.
Add the .gitignore file to ignore certain files.