def chebyu(n, monic=False):
base = jacobi(n, <NUM_LIT:0.5>, <NUM_LIT:0.5>, monic=monic)<EOL>if monic:<EOL><INDENT>return base<EOL><DEDENT>factor = sqrt(pi) / <NUM_LIT> * _gam(n + <NUM_LIT:2>) / _gam(n + <NUM_LIT>)<EOL>base._scale(factor)<EOL>return base<EOL>
Return nth order Chebyshev polynomial of second kind, Un(x).  Orthogonal
    over [-1,1] with weight function (1-x**2)**(1/2).