qdODE_ls function

least-square fit for qdODE model