Математическая модель трех небесных тел

Задача Коши:

Где:

Я нашел интересную формулу для вычисления ошибки метода Мерсона:

План для модернизации:

  1. Изучить вычисление ошибки в методе Мерсона (локальная , глобальная — )

    Есть теоретическая оценка, есть через ричардсон-экстраполяцию:

  2. Добавить отрисовку с matplotlibcpp.h
  3. Оптимизировать метод Мерсона с avx
  4. Переписать задачу в виде классов для того, чтобы можно было удобно работать с несколькими сущностями
  5. Сделать нормальный адаптивный шаг с учетом ошибки из пункта (1)
  6. Сделать возможность работать с несколькими объектами
  7. Разработать логику для расчета столкновений