I doubt this is still needed, but:
If you want to write your own script, you can use outward recursion for small arguments and then inward recursion for large arguments. However, this requires the knowledge of and in order to do the outward recursion and they are also helpful for the inward recursion for scaling (since they will be very large when unscaled). However, there are other ways to do the scaling because the amplitudes all add to one, which can be reasonably approximated if you take n high enough. The issue with the outward recursion is that it becomes unstable after approximately n=10, while the inward recursion doesn't, but it gets huge as you iterate in and you might run into precision issues if you don't work carefully - i.e. use doubles.
If you use the below algorithm and code it in matlab, you get excellent results for nearly any order spherical bessel for, at least a wide range, of the argument. I have tested it out to 50, but I can not confirm it past that. Though, they have scaling done so that it should work very well for large arguments.
If you make use of a parfor loop over the argument, you can interate the recursion through your desired argument range quite quickly, though perhaps not the most efficiently.