def roots_sh_chebyu(n, mu=False):
 (x, w, m) = roots_chebyu(n, True)
 x = ((x + 1) / 2)
 m_us = cephes.beta(1.5, 1.5)
 w *= (m_us / m)
 if mu:
  return (x, w, m_us)
 else:
  return (x, w)

def roots_sh_chebyu ( n , mu = False ) : ( x , w , m ) = roots_chebyu ( n , True ) x = ( ( x + 1 ) / 2 ) m_us = cephes . beta ( 1.5 , 1.5 ) w *= ( m_us / m ) if mu : return ( x , w , m_us ) else : return ( x , w )
gauss-chebyshev quadrature .