scipy.special.besselpoly

scipy.special.besselpoly(x1, x2, x3[, out])

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