[SciPy-dev] Gauss-Lobatto points

[Ondrej Certik]

Hi,

does anyone know a way to get Gauss-Lobatto points out of scipy? I
know how to get Gauss integration points, but that's different:

In [1]: from scipy.special.orthogonal import p_roots

In [2]: p_roots(2)
Out[2]: [array([-0.57735027,  0.57735027]), array([ 1.,  1.])]

In [3]: p_roots(3)
Out[3]:
[array([ -7.74596669e-01,  -4.78946310e-17,   7.74596669e-01]),
 array([ 0.55555556,  0.88888889,  0.55555556])]

In [4]: p_roots(4)
Out[4]:
[array([-0.86113631, -0.33998104,  0.33998104,  0.86113631]),
 array([ 0.34785485,  0.65214515,  0.65214515,  0.34785485])]

In [5]: p_roots(5)
Out[5]:
[array([ -9.06179846e-01,  -5.38469310e-01,  -1.25197490e-16,
        5.38469310e-01,   9.06179846e-01]),
 array([ 0.23692689,  0.47862867,  0.56888889,  0.47862867,  0.23692689])]

the Gauss-Lobatto points are used for interpolation.

Thanks!
Ondrej