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