Расчёт орбиты движения звёзд

В архиве находится код программы для расчёта орбиты движения звёзд. В файле STPIC10.FOR находится основной код программы с графикой и считыванием данных с клавиатуры. В файле torb1.for код программы без графики и считыванием данных из файла. Файлы ORBIT.FOR RKF46.FOR являются подпрограммами. Физика задачи хорошо описана в этих 2х статьях

http://www.astronet.ru/db/msg/1245721/lec.16.1.htm...

http://www.astronet.ru/db/msg/1245721/lec.16.2.htm...

P.S.

STPIC10.mp3 это STPIC10.FOR

Необходимо переписать программу с fortran'а на с или c# или любой другой современный язык программирования с хорошими комментариями. Программа использует 2D графику, В переработанном варианте должна быть 3D графика или хотя бы 2D. Также должна быть возможность в процессе вывода траектории звезды на экран останавливать работу программы и получать текущие координаты.Программу можно и не переписывать если напишите свою отвечающюю всем требованиям