sphericalbesselj(l, x) = NaN for large l and small nonzero x

[hersle]

using Bessels
sphericalbesselj(1000, 0.0) # 0.0, good
sphericalbesselj(1000, 1e-20) # NaN, bad