Что-то созвездия не на местах... Я к примеру родился в августе... под созвездием Льва, а не Водолея. А вообще-то можно посмотреть неплохое видео...
Надо только не забывать что наше Солнце рядовая звезда в галактике под названием Млечный путь.

Валерий, мне нравится направление вашей деятельности. Мне кажется согласовать все взаимодействия гравитационное, электрическое, магнитное возможно только в графическо-векторном варианте. Вы и ваш сын поставили задачу невероятной сложности, упрощать её конечно возможно для показа конкретных взаимодействий и частных случаев, так в основном и делают в современной науке. Но упрощать лучше когда вы сами знаете о том, что пограмма решает частный случай понимания и это оговаривается. Накапливайте максимальное количество знаний в этом направлении, результат будет положительный. Предложите сыну учавствовать в астрономическом форуме, там немало грамотных людей и есть раздел астрономических программистов. Анатолий Волович прав движение солнечной системы в сторону северного полюса астрономы не отрицают. Значит траектория движения планет спираль

Рассматривать движение планет в декартовой системе координат и составлять математический алгоритм, я думаю менее эффективно, чем радиально векторная система координат. Посмотрите в астрономии там есть и другие системы координат пространства и времени тоже. Пример сутки Это оборот Земли вокруг своей оси, продвижение по орбите, Движение в системе солнечной системы и тд. Описать эти движения в векторном выражении лучше в радиальной системе координат. Это только мои мысли не более.

Валера, извини но это не правдоподобно:

Более правдоподобно будет выглядеть система если у Солнца будет один заряд а у планет противоположный, но это всё ровно только приближённо. Нужно не забывать что и отдельно взятые планеты тоже должны позиционировать себя также как и Солнце к другим участникам системы. То-есть, только силы притяжения между всеми и ни каких сил отталкивания.
Не пойму тебя, почему в свою программу не вобьёшь закон всемирного тяготения вместо кулоновского? Зачем эти излишества с отталкиванием?

(А у зарядов - там бывают оба варианта, в зависимости от знаков основных зарядов. Там пробным является стандартный - положительный.) Поясните вашу мысль, что значит в вашем понимании " стандартный - положительный"

Валерий посмотрите как на картах выражается давление это эквопотенциальные поверхности на сфере. Вы даете изображение в плоскости, но понять гравитацию в плоскости, очень трудно, вы стараетесь привязать к декартовой системе координат . В объеме расчитывать эквопотенциальные сферы, это конечно круто, но это реально или нет?, особенно с декартовой системе координат