А как в SE определяется область стабильных планетных орбит для двойных звезд? 

Для орбит S-типа (вокруг одной из звёзд) - по предельным возмущениям от второй звезды (этот случай аналогичен спутнику у планеты: орбита стабильна, если её радиус меньше 0.5 радиуса сферы Хилла).

Для орбит P-типа (вокруг обоих звёзд) - по эмпирическому правилу: радиус стабильных орбит спутников больше утроенного радиуса взаимной орбиты звёзд.

Ага, вот эти цифры мне и были нужны. Спасибо

А орбиты самих звезд в парке как вычисляются?

А как их вычислить? Рандомно генерируются. Ну почти, учитывается наблюдательная статистика, но она есть только в виде % тесных пар - % широких пар.

В смысле в паре а не в парке. Как это рандомно? То есть вот есть барицентр и тупо от балды ставится расстояние и скорости этих звезд?

QR_BBPOST  но она есть только в виде % тесных пар - % широких пар.

Этого я не нашел. Нашел упоминания о максимуме орбит с периодом взаимного обращения на уровне 30-100 лет и все.

Скорость не от балды конечно, SE рассчитывает движение по орбитам честно. Но взаимная большая полуось - с некоторым псевдослучайным гауссовым распределением.

Как так-то? Вот Земля тоже немного вращается против Луны. Если бы мы не знали точные цифры, то SE смог бы их посчитать или нет? Я имею в виду что вот есть два тела, там же есть какая то формула чтобы посчитать на каком они расстоянии от барицентра и все такое?

Нет, размер орбиты может быть любым. На какую орбиту вывел спутник, на такой он и останется.

Вы это серьезно? Тогда бы Земля могла вращаться на таком же расстоянии от барицентра как и Луна. А барицентр это центр масс, то есть там формула на каком расстоянии от двух тел он получится.

Не, масса тел учитывается конечно. Я о том, что радиус взаимной орбиты ничем не ограничивается, разве что внешней планетой или второй звездой.

Вот, я в начале как раз спрашивал, если есть 2 звезды то как вычислить их взаимное обращение? Я попробовал такие создать без указания орбит, просто обе вокруг барицентра. Так они просто свалились в одну точку и никак не вращаются. Что я им прописать должен?

Хотя бы период. И массы обязательно.

А еще такой вопрос. Как в SE решается уравнение Кеплера? Разложением в ряд Фурье или итерационным методом, или еще как? 

Итерационным. Для любых орбит не более 50 итераций, для близких к круговым - вообще 1-2 итерации.