\sqrt{r^2 + (vt)^2} = r \sqrt{1 + (vt/r)^2} = r + (1/2) (vt)^2 / r + o(t^2),
так что вторая производная равна

v^2/r = (vr)^2 / r^3.

Числитель — (vr)^2 — это квадрат углового момента. Так что он вдоль орбиты всегда один и тот же!
А знаменатель r^3 — как раз и соответствует закону всемирного тяготения: куб, как и раньше. Так что, если взять константу
l = (vr)^2/ (GM),
то для разницы (r-l) будет
(r-l)’’ = r’’ = hr + (vr)^2 / r^3
= - GM r / r^3 + (vr)^2 / r^3
= - GM / r^3 * (r- l)
= h (r-l).

Так что трёхмерный вектор R = (x,y,r-l) подчиняется центральному закону
R’’ = h R, h = -GM/r^3.

Ура — теорема доказана!