I finally figure out the problem of my SBML ODE Solver programmer, I didn't make the right tolerance setting for Sundials cvodes function. I feel so good!