legendre polynomials code snippets in matlab