MAINT: to support older versions of numpy, use scipy.special.p_roots ins...

...tead of np.polynomial.legendre.leggauss