The besi1(x) function is the modified Bessel function or order 1. besi1 expects its argument to be in radians.