Олег Игоревич, попробовал я получить информацию о прецессии Земли и Нептуна.

За основу взял эфемериды DE441 (от JPL лаборатории НАСА, покрывают временной интервал -13200--17191 гг.). В данном случае взял интервал времени от -10000 года до 1900 года.

С помощью пакета SkyField (для пайтона) получил таблицу, где на каждое 1 января каждого года из интервала, указанного выше, просто запросил вычислить эклиптические координаты Земли и Нептуна (также дополнил информацией по Солнцу, для возможной отладки). Угол между Землёй и Солнцем вычислял только по долготе как разность значения координат Земли и Нептуна (игнорировал при этом широту). Помимо таблицы, построил по этим данным также график изменения угла между ними ("Earth_Neptune_-10000..1900").

Чтобы попробовать извлечь из этого массива данных информацию о прецессии, с помощью пакета astropy сделал перевод в гелиоцентрическую систему координат с усреднённой эклиптикой, попробовал построить аналогичные графики, но пропуская по 165 лет между соседними точками и стартуя с разных сдвигов в пределах этих 165 лет. Результат изображён на втором графике ("Earth_Neptune_сдвиги_-10000..1900"). Возможно, точки пересечения кривых смогут дать какие-то сведения.

Все файлы выложил на гугл-диске в виде архива, ссылка [обновил - добавил прокоментированный файл с вычислениями, сохранённый в html формате, а также лист в xls файле с данными, по которым строил второй график]:
https://drive.google.com/file/d/1L0jUrqWn-ORtl8jvEJbkCjEKy4-Xbb7B

Картинки в двух форматах - pdf и jpg.