def roots_sh_chebyt(n, mu=False):
 xw = roots_chebyt(n, mu)
 return ((((xw[0] + 1) / 2),) + xw[1:])

def roots_sh_chebyt ( n , mu = False ) : xw = roots_chebyt ( n , mu ) return ( ( ( ( xw [ 0 ] + 1 ) / 2 ) , ) + xw [ 1 : ] )
gauss-chebyshev quadrature .