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


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

где 

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

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

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

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

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