Программа расчета движения планет (я новичок) - помогите
 

Мне нужно решить задачу n тел на примере солнечной системмы методом численного интегрирования. Вроде бы все просто, но я не могу разделаться с ошибками. Помогите пожалуйста. Вот формула, по которой я считаю: (во воложении). Подскажите плз, чему равна в ней постоянная Гаусса, если расстояния в ае, время в эфемеридных сутках. Массы - обратные массы планет (m0/m). Координаты - гелеоцентрические.Также не могли бы вы подсказать как получить координаты и скорости больших планет на опр. промежуток времени в нормальном виде (скорости - ае/сутки координаты - ае или подобное (без склонений всяких, я в них утону) ). Пытался разобраться с DE200, там вроде они есть, но либо я совсем дурак, либо там не то, что мне нужно (скорее первое в плане астрономии). Я мало что смыслю в Астрономии, так что извините если что.