def chebyc(n, monic=False):
if n < <NUM_LIT:0>:<EOL><INDENT>raise ValueError("<STR_LIT>")<EOL><DEDENT>if n == <NUM_LIT:0>:<EOL><INDENT>n1 = n + <NUM_LIT:1><EOL><DEDENT>else:<EOL><INDENT>n1 = n<EOL><DEDENT>x, w, mu0 = c_roots(n1, mu=True)<EOL>if n == <NUM_LIT:0>:<EOL><INDENT>x, w = [], []<EOL><DEDENT>hn = <NUM_LIT:4> * pi * ((n == <NUM_LIT:0>) + <NUM_LIT:1>)<EOL>kn = <NUM_LIT:1.0><EOL>p = orthopoly1d(x, w, hn, kn,<EOL>wfunc=lambda x: <NUM_LIT:1.0> / sqrt(<NUM_LIT:1> - x * x / <NUM_LIT>),<EOL>limits=(-<NUM_LIT:2>, <NUM_LIT:2>), monic=monic)<EOL>if not monic:<EOL><INDENT>p._scale(<NUM_LIT> / p(<NUM_LIT:2>))<EOL>p.__dict__['<STR_LIT>'] = lambda x: eval_chebyc(n, x)<EOL><DEDENT>return p<EOL>
Return nth order Chebyshev polynomial of first kind, Cn(x).  Orthogonal
    over [-2,2] with weight function (1-(x/2)**2)**(-1/2).