February 20, 2023 at 1:45 pm
sidjou
Subscriber

Thanks for your suggestion!
I reduce the time step and the issue for the position repetition has been corrected.
About the exact position issue at t=0, I'm already using the double precision for my calculation.