So wie Sie sich bisher mit Formeln angestellt haben, wage ich es ja zu bezweifeln, daß Sie über ein idealisiertes g=G·M/r² bei einem infinitesmal kleinem Gravitationszentrum hinauskommen werden. In der Tat müssen Sie dabei aber mindestens die modifizierte Neutonnendynamik anwenden, für g den Wert g=(G·M·a)/r  setzen und in die Formel für g'=(g/a)/[1+(g/a)²]*g einsetzen. Die Variable a hat die Einheit einer Beschleunigung und kann jeden Wert zwischen 1 und cosh(g/x)*g annehmen. Das genaue Verhalten des x-Parameters ist jedoch noch Gegenstand aktueller Forschung, was die Sache ein wenig verkompliziert. Nichtsdestotrotz, sollten Sie so weit kommen, müssen Sie als nächstes nach infinitesmal kleinen Teilstücken und deren Feldlinien, auch hier je mehr desto besser, integrieren. Hier eine grobe Skizze mit drei Teilstücken [0:0, 4:0, 8:0] in zwei Dimensionen: