get_legendre_par function

use legendre polynomials to fit a given data