cubic interpolation. i use this to speed up evaluation of bessel functions