code: def roots(cls, degree, domain, kind):
          
          basis_coefs = cls._basis_monomial_coefs(degree)
          basis_poly = cls.functions_factory(basis_coefs, domain, kind)
          return basis_poly.roots()
docs: Return optimal collocation nodes for some orthogonal polynomial.