]> bessel_y1_small_z Y 1 z = z π log z 2 − 2 π z − z 2 π 1 − 2 γ ; z < ε