траектория спутника С++ .

void Calc(double &X, double &Y,double &Xold, double &Yold,double &SpeedX,double &SpeedY)
{     
//масса земли умноженная на гравитационную постоянную
    const long long Kearth =667384000000000; 
 
    Xold =X;
    Yold =Y;
    long long Rx =X/0.031392; //  0.031392 - отношения маштабирования Земли (пиксели/км)
    long long Ry =Y/0.031392;
    long long absR = sqrt(Rx*Rx+Ry*Ry);
    SpeedX+=(-Kearth)/(absR*absR*absR)*Rx*(TIMER/1000);
    SpeedY+=(-Kearth)/(absR*absR*absR)*Ry*(TIMER/1000);
    Rx+=SpeedX;
    Ry+=SpeedY;
    X=Rx*0.031392;
    Y=Ry*0.031392;
}

#define EarthR 150  // радиус земли на экране в пикселях. Можно менять связь с размерами в функции сохраниться
#define TIME 10 //дельта времени для расчетов
const double RealKoef= (double)EarthR/6371; // коефициент маштабированния
const double Kearth =398600.448; // км3/сек2. масса земли умноженная на гравитационную постоянную и разделенная на количество секунд в часе
 
 
void Calc(double &X, double &Y,double &SpeedX,double &SpeedY)
{
    double Rx =X/RealKoef; //км. координата Х вектора R
    double Ry =Y/RealKoef; //км. координата Y вектора R
    double absR = pow((X*X+Y*Y),0.5)/RealKoef; //км. модуль вектора R
    SpeedX+=(-Kearth/(absR*absR*absR))*TIME*Rx;  //км/с. прирощение на дельту времени вектора скорости по координате X
    SpeedY+=(-Kearth/(absR*absR*absR))*TIME*Ry;  //км/с. прирощение на дельту времени вектора скорости по координате Y
    Rx+=(SpeedX)*(TIME); //км. сумирование вектора  радиуса с вектором скорости по координате Х
    Ry+=(SpeedY)*(TIME);//км. сумирование вектора  радиуса с вектором скорости по координате Y
    X=Rx*RealKoef;//0.031392;//пиксели. возврат к пиксельным координатам 
    Y=Ry*RealKoef;//пиксели
}