I am getting weird results when I try to compute the bessel_k function to 40 or 50 places. I think that the bessel_k function is just returning a float and the big float is just random digits from the 17th digit on. Is this the correct analysis? Rich