scipy.special.eval_legendre
scipy.special.eval_chebyu
Evaluate Chebyshev T polynomial at a point.
This routine is numerically stable for x in [-1, 1] at least up to order 10000.