Здорово вернуться к Fortran (95). Я закодировал орбиту на расстоянии 300 а.е. Ньютон предсказывает синюю кривую. QI предсказывает оранжевую. При QI ускорение к Солнцу немного выше, поэтому тело приближается и ускоряется: