Cимуляция движения тела по орбите?
Если у вас реальные физические формулы, то я бы сказал что ожидаемый результат - движение тела по эллиптической орбите. 
Все тела в космосе так и летают. 
Более того, даже этот сдвиг эллипса на каждом круге - это то что происходит в реальности. 

Может у вас программа наоборот - слишком правильно работает?
Попробуйте увеличить скорость, посмотрите как меняется эксцентриситет эллипса.
Может при какой-то скорости эллипс станет неотличим от окружности.