https://www.bragitoff.com/2017/08/spherical-bessel-functions-c-program/