The besj1(x) function returns the J1st Bessel function of its argument. besj1 expects its argument to be in radians.