y=y1(x) returns the Bessel function of the second kind of order 1 at x.
scipy.special.y0
scipy.special.i0