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

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