code: def chebgauss(deg):
      """
      Gauss-Chebyshev quadrature.
  
      Computes the sample points and weights for Gauss-Chebyshev quadrature.
      These sample points and weights will correctly integrate polynomials of
      degree :math:`2*deg - 1` or less over the interval :math:`[-1, 1]` with
      the weight function :math:`f(x) = 1/\\sqrt{1 - x^2}`.
  
      Parameters
      ----------
      deg : int
          Number of sample points and weights. It must be >= 1.
  
      Returns
      -------
      x : ndarray
          1-D ndarray containing the sample points.
      y : ndarray
          1-D ndarray containing the weights.
  
      Notes
      -----
  
      .. versionadded:: 1.7.0
  
      The results have only been tested up to degree 100, higher degrees may
      be problematic. For Gauss-Chebyshev there are closed form solutions for
      the sample points and weights. If n = `deg`, then
  
      .. math:: x_i = \\cos(\\pi (2 i - 1) / (2 n))
  
      .. math:: w_i = \\pi / n
  
      """
      ideg = pu._deprecate_as_int(deg, "deg")
      if ideg <= 0:
          raise ValueError("deg must be a positive integer")
  
      x = np.cos(np.pi * np.arange(1, 2 * ideg, 2) / (2.0 * ideg))
      w = np.ones(ideg) * (np.pi / ideg)
  
      return x, w
  
imports: import numpy
  import numpy.linalg
  from numpy.core.multiarray import normalize_axis_index
  from .None import polyutils
  from ._polybase import ABCPolyBase
label_desc: Use image node_id 31 for calling a global function with example usage: chebgauss(deg) and returns: x,  w