C++: Отладить решение системы нелинейных обыкновенных дифференциальных уравнений с граничными условиями

Есть программа на C++, в которой записана система нелинейных обыкновенных дифференциальных уравнений и граничные условия (значения параметров в начальный и конечный моменты времени). В программе не получается учесть условия на правой границе (в конечный момент времени). Нужно отладить программу, чтобы были учтены условия на обеих границах.

Программа доступна по ссылке: https://cloud.mail.ru/public/ggye/tQuMG3PVB

PS: Использовать библиотеку heyoka (применяется в программе) не советую ввиду энергозатратного процесса установки.