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

[Arderun]

вопрос к тем кто помнит физику
пытаюсь написать функцию апроксимированной 2D траэктории спутника вокруг земли(точка 0,0)
так как не хочу возится с полярными координатами то класическую формулу

ma =-G((Mm)/(|r|^3))r

пишу через векторы, но получаемые результаты даже не в близких порядках от нужных(спутник просто исчезает с экрана). видимо где очень грубая ошибка. так как физика совсем не мое то прошу вашей помощи. Задачка вобщем класическая, даже школьная.

 
//взятие корня у long long
long long sqrBig(long long A, int clear)
{long long B=A/2;
int i=0;
while( i < clear)
{
	//temp=B;
	if((B*B)<A)
	{
		B=B+B/2;
		
	}
	else
	{B=B-B/2;}
	//B0=temp;
	i++;
}
return B;
}

//сама функция подсчета траектории
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 = sqrBig(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;
}

возможно чтот оуперается в порядок чисел, так что если кто знает класс для работы с порядками до 10^25 степени буду признателен