Я первый раз когда увидел тоже сильно удивился, что система так себя ведет. В крайнем случае рассчитывал на то, что планеты стабилизируются в устойчивых точках Лагранжа L4 и L5, но не тут то было.

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