e = 1.0 problem

[Trifonoff]

Hi Stefano!

I just remember another thing…. When you are checking for stability with the implemented n-body engines I noticed the following:

Integration: 3652500.00 days [10000.00 years], survival time: 10000.00 years, int.method: 4 Planet 1 Pmin: 294.66, Pmax: 360.80, dev(P)/med(P): 5.61e-02 emin: 0.05, emax: 1.00, dev(e)/med(e): 5.69e-01… Planet 2…. …

well although the survival time is 10000.00 years, I think this system is unstable (in the first 500 years) because e1 = 1.0 (emax). When e1 = 1.0 then the peri distance q1 = (1.0-e1)*a1, which will be 0 and this is not physical. Everybody can judge for him self what would be a stable system, but I think you can set more strict criteria when the stability simulation to be interrupted (to save CPU time and give adequate stability answer).

All the best, Trifon