Что-то я опять впал в ступор.
Момент импульса по определению равен

.
Если перейти к скалярам, то

.
Для орбиты Земли всё хорошо,

, (синус прямого угла есть 1).
А для промежуточной орбиты опять всё сложно,

. Здесь и радиус-вектор - функция, и скорость - функция, и ещё синус появляется.
Если записать формулу момента импульса иначе (

), то всё равно дело встанет из-за отсутствия функции радиус-вектора от времени, так как момент инерции зависит от

.
Я вижу только два пути. Нужно либо принять, что все орбиты есть круги (что запрещено условием), либо подключать тяжёлую артиллерию и рассматривать движение по эллипсам (что затянется, я думаю, на долго). Но раз уж умные люди в унисон говорят: "Используй момент импульса", то, получается, есть ещё один вариант, который учитывает и эллиптичность орбиты и матанализ не привлекает к работе. Короче, прошу ещё одну подсказку :)