def gauss_chebyshev_t(n, n_digits):
 xi = []
 w = []
 for i in range(1, (n + 1)):
  xi.append(cos(((((2 * i) - S.One) / (2 * n)) * S.Pi)).n(n_digits))
  w.append((S.Pi / n).n(n_digits))
 return (xi, w)

def gauss_chebyshev_t ( n , n_digits ) : xi = [ ] w = [ ] for i in range ( 1 , ( n + 1 ) ) : xi . append ( cos ( ( ( ( ( 2 * i ) - S . One ) / ( 2 * n ) ) * S . Pi ) ) . n ( n_digits ) ) w . append ( ( S . Pi / n ) . n ( n_digits ) ) return ( xi , w )
computes the gauss-chebyshev quadrature [1]_ points and weights of the first kind .