Hello,
I have a system of two differential equations. The system is derived through Hamilton's equations from a specific Hamiltonian. I have been trying to integrate the system for 10^8 years. However, after 10^6 years the integration is failing and I am getting an error saying:
"at t = 2.472503629*10^6, the step size is effectively zero. Singularity or stiff system suspected."
When I integrate up to 10^6 years, the integration goes fine. However, I am trying to graph the full phase space of the system. The entire phase space requires integration up to ~ 10^8 years or more.
I have attached the notebook that I am using to integrate the system.
I would appreciate any help with that.
Thanks!