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