Ok, я понял. А что не так с распределением плотности тёмной материи по ?
Именно при таком распределении невозможно по наблюдаемому вращению отличить сферическое распределение от дискового. эффект будет тот же самый.
А у меня тут еще мысли: если гравитационное поле имеет собственную гравитирующую массу пропорциональную квадрату напряженности (как в электромагнитном), то тоже получится именно такое распределение. Я тут подобрал гидростатическую модель гравитации для этого случая и все вроде бы сходится.
Для сферически симметричного случая получается:
где
- полная гравитирующая плотность массы (она же - гидростатическое давление)
- ускорение (напряженность гравитационного поля)
- плотность обычной массы
- константа гравитации (вероятно что-то типа
)
Решаю численно в маткаде - картинка вполне правдоподобная.
Можно дальше попробовать расширить систему для динамики.