It seems there is no ::Interval support for besselk and besseli of SpecialFunctions.jl since IntervalRootFinding.roots() gives a StackOverflowError when applied to (a combination of) these functions. This issue may be relevant for both IntervalArithmetic and IntervalRootFinding.