code: def roots_chebyt(n, mu=False):
      r"""Gauss-Chebyshev (first kind) quadrature.
  
      Computes the sample points and weights for Gauss-Chebyshev quadrature.
      The sample points are the roots of the n-th degree Chebyshev polynomial of
      the first kind, :math:`T_n(x)`.  These sample points and weights correctly
      integrate polynomials of degree :math:`2n - 1` or less over the interval
      :math:`[-1, 1]` with weight function :math:`f(x) = 1/\sqrt{1 - x^2}`.
  
      Parameters
      ----------
      n : int
          quadrature order
      mu : bool, optional
          If True, return the sum of the weights, optional.
  
      Returns
      -------
      x : ndarray
          Sample points
      w : ndarray
          Weights
      mu : float
          Sum of the weights
  
      See Also
      --------
      scipy.integrate.quadrature
      scipy.integrate.fixed_quad
      numpy.polynomial.chebyshev.chebgauss
      """
      m = int(n)
      if n < 1 or n != m:
          raise ValueError("n must be a positive integer.")
      x = _ufuncs._sinpi(np.arange(-m + 1, m, 2) / (2 * m))
      w = np.full_like(x, pi / m)
      if mu:
          return x, w, pi
      else:
          return x, w
  
imports: from __future__ import division, print_function, absolute_import
  import numpy
  from numpy import exp, inf, pi, sqrt, floor, sin, cos, around, int, hstack, arccos, arange
  from scipy import linalg
  from scipy.special import airy
  from .None import _ufuncs
  from .None import _ufuncs
  from .None import specfun
  from ._ufuncs import binom, eval_jacobi, eval_sh_jacobi, eval_gegenbauer, eval_chebyt, eval_chebyu, eval_chebys, eval_chebyc, eval_sh_chebyt, eval_sh_chebyu, eval_legendre, eval_sh_legendre, eval_genlaguerre, eval_laguerre, eval_hermite, eval_hermitenorm
label_desc: Use image node_id 23 for calling a global function with example usage: roots_chebyt(n, mu) and returns: x,  w,  pi, x,  w