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


Для сферически симметричного случая получается:

где

- полная гравитирующая плотность массы (она же - гидростатическое давление)

- ускорение (напряженность гравитационного поля)

- плотность обычной массы

- константа гравитации (вероятно что-то типа

)
Решаю численно в маткаде - картинка вполне правдоподобная.
Можно дальше попробовать расширить систему для динамики.