base = sh_jacobi(n,2.0,1.5,monic=monic)
if monic: return base
factor = 4**n
return base * factor
// Legendre
def p_roots(n,mu=0):
[x,w] = p_roots(n)
After Change
Return nth order shifted Chebyshev polynomial of second kind, Un(x).
Orthogonal over [0,1] with weight function (x-x**2)**(1/2).
base = sh_jacobi(n,2.0,1.5,monic=monic)
if monic: return base
factor = 4**n
base._scale(factor)return base
// Legendre
def p_roots(n,mu=0):
[x,w] = p_roots(n)
