Проблемы реализации интегратора чехарды в C++ ⇐ C++

[Anonymous]

Я пытаюсь реализовать интегратор чехарды как средство численного решения задачи двух тел. Мой код не прошел тест из-за того, что он слишком далек от аналитических значений. Судя по результатам теста, компонент y вектора положения отличается в 2 раза. Просмотрите код и помогите мне.
Это код:

Код: Выделить всё

#include "ODESolverLeapfrog.hpp"
#include 

ODESolverLeapfrog::ODESolverLeapfrog(const SetNorm &normalisation, std::vector &mass,
double t_final, double dt_initial, double dt_out) :
ODESolverSingleStepBase(normalisation, mass, t_final, dt_initial, dt_out) {
std::cout num_steps_++;
this->i_step_ = 0;
}

ODESolverLeapfrog::~ODESolverLeapfrog() {
std::cout 

Подробнее здесь: [url]https://stackoverflow.com/questions/79157409/issues-implementing-the-leapfrog-integrator-in-c[/url]