scipy.special.jacobi
scipy.special.genlaguerre
Return the nth order Laguerre polynoimal, L_n(x), orthogonal over [0,inf) with weighting function exp(-x)