
n=6;
x=linspace(-1,1,10); x=x';
px = legendre_polynomials ( n, x);
