Uses Legendre Polynomials to determine the support points for the
design:
If \(a=-1\), \(b=1\): places \(k +1\) support points in
\([-1,1]\), located at the roots of \((1-x^{2})\frac{dP_{k}(x)}{dx}\)
where \(P_{k}(x)\) is the Legendre polynomial of degree \(k\)).
Distributes the n measurements almost equally over the
support points.