Угу моя ошибка - по формуле из первой статьи нельзя считать орбиты высотой отличной от радиуса Земли - g уменьшается обратно квадрату расстояния, а я посчитал его как константу.
Но на самом деле она нам и не нужна, нужна тангенциальная составляющая. Интересно другое - цель была оценить затраты энергии на высоких орбитах. Для этого надо посчитать изменение высоты круговой орбиты (пусть на 1 км) для орбиты H=500 км и H=36000 км, по простой двухимпульсной (с двумя тангенциальными импульсами) схеме изменения радиуса орбиты. Масса КА там и там = 250 кг, изменением массы в процессе импульса можно пренебречь для простоты расчётов. Попробую дома вечером, если не лень будет.