У massa есть формула для скорости кругового движения вокру тяготеющего центра под действием радиапьной силы

V=sqrt( (1-2GM / r) (GM / r^2 - F1 / m) / ( ( (1-2GM/ r ) / r ) - F1 / m ) ).

На самом деле так выражается не сама скорость, а отношение V/c. Просто у massa c=1. Мне кажется, что удобнее эту формулу представить вот так

V/c = sqrt{[(1+A)(1-B)B/2] / [1-(1+A/2)B]}

где B=Rg/R, a Rg - гравитационный радиус
A -коэффициент пропорционалности между силой тяжести и радиально силой, который может (в зависимости от R) меняться в пределах

-1=< A < (1-B)^2 / [(3B/2)-B^2].

Видно,что при R-> Rg верхний предел А стремится к 0.