Returns the nth order Legendre polynomial, P_n(x), orthogonal over [-1,1] with weight function 1.
scipy.special.eval_sh_jacobi
scipy.special.chebyt