qdODE_fit function

legendre polynomials fit to qdODE model