import numpy as np
import numpy.linalg as la
import scipy.special as sp
import matplotlib.pyplot as pt
n = 5
nodes = np.linspace(-1, 1, n)
#nodes = sp.legendre(n).weights[:, 0]
pt.plot(nodes, 0*nodes, "o")
exponents = np.arange(n)
exponents
exponents[:, np.newaxis]
vdm = nodes ** exponents[:, np. newaxis]
if n < 5:
print(vdm)
print(la.cond(vdm))
vdm = np.array([sp.eval_legendre(i, nodes) for i in range(n)]).T
print(la.cond(vdm))