y=besselpoly(a,lam,nu) returns the value of the integral: integral(x**lam * jv(nu,2*a*x),x=0..1).
scipy.special.it2i0k0
scipy.special.jvp