Обобщенная задача кеплера

Добрый день, друзья! Такая проблема. Сказали написать программу для обобщенной задача Кеплера в геометрии Шварцшильда (некоторое тело вращается вокруг черной дыры). Заданы уравнения:
dr/dt = квадратный корень из (E^2 - (1 - 2m/r)(1 + L^2/r^2));
df/dt = L/r^2.
Суть программы в том, что значения t0 - собственное время и f0 - угол поворота должны менять непрерывно, а r0 - радиус траектории вращения должен менять свое значение на определенном отрезке (в расчетах получается отрезок, если округлить, [4;8]). Сначала он увеличивается, а потом уменьшается. Увеличение и уменьшение связанно с условием: E^2 - (1 - 2m/r)(1 + L^2/r^2) > 0 (т.к. dr/dt = квадратный корень из (E^2 - (1 - 2m/r)(1 + L^2/r^2)) ). Это реализовано с помощью цикла for - как только условие не выполняется, цикл с увеличением(уменьшением) r останавливается , и начинается другой, с уменьшением(увеличением) r. Значения r и f ищутся с помощью метода Адамса-Башфорта.
Проблема заключается в том, что значения в одну сторону и в обратную (с увеличением и уменьшением r, соответственно) по отдельности ищутся, но как только пытаешься соединить эту задачи вместе и зациклить (да хотя бы совместить) программа перестает работать (пишет ошибку: "Необработанное исключение в "0x000b24a0" в "BH++.exe": 0xC00001A5: Обнаружен недопустимый обработчик исключений."). Как мне кажется, все связанно с переобозначениями при переходе во второй цикл do (именно на этом моменте программа стопориться). По отладки ничего не смог увидеть.
Друзья, прошу помощи! Помогите исправить ошибку! Заранее благодарю!
P.S. Так как получается "длиннопост", прикладываю к посту файл с программой.

Вложения

BH++.rar (445 байт, 7 просмотров)

0

Прошу прощения за неполный файл с лишними записями в коде.

Вложения

BH++.rar (2.09 Мб, 16 просмотров)

0