Столкновение двух окружностей

void Ball_to_Ball_Detection(Object_Ball& Ball1,Object_Ball& Ball2)
{
    bool stolknovenie = false;
    int xmov1 = Ball1.xSpeed;
    int ymov1 = Ball1.xSpeed;
    int xmov2 = Ball2.xSpeed;
    int ymov2 = Ball2.ySpeed;
 
    int xl1 = Ball1.x;
    int yl1 = Ball1.y;
    int xl2 = Ball2.x;
    int yl2 = Ball2.y;
 
    int R = Ball1.Radius + Ball2.Radius;
    int a = -2*xmov1*xmov2+xmov1*xmov1+xmov2*xmov2;
    int b = -2*xl1*xmov2-2*xl2*xmov1+2*xl1*xmov1+2*xl2*xmov2;
    int c = -2*xl1*xl2+xl1*xl1+xl2*xl2;
    int d = -2*ymov1*ymov2+ymov1*ymov1+ymov2*ymov2;
    int e = -2*yl1*ymov2-2*yl2*ymov1+2*yl1*ymov1+2*yl2*ymov2;
    int f = -2*yl1*yl2+yl1*yl1+yl2*yl2;
    int g = a+d;
    int h = b+e;
    int k = c+f-R*R;
 
    float koren = sqrt(h*h-4*g*k);
    float t1 = (-h+koren)/(2*g);
    float t2 = (-h-koren)/(2*g);
 
    if (t1>0 && t1 <=1)
     {
         whatTime = t1;
         stolknovenie = true;
     }
    if (t2>0 && t2<=1)
     {
        if (whatTime == 0 || t2<t1) 
        {
            whatTime = t2;
            stolknovenie = true;
        }
    }
    if (stolknovenie) 
    {
        printf ("Ñòîëêíóëèñü!\n");
        Ball_to_Ball_Reaction(Ball1,Ball2,whatTime);
    }    
}
 
 
void Ball_to_Ball_Reaction(Object_Ball& Ball1,Object_Ball& Ball2,float time)
{
    int mass1 = Ball1.Mass;
    int mass2 = Ball2.Mass;
    
    int xVel1 = Ball1.xSpeed;
    int xVel2 = Ball2.xSpeed;
    int yVel1 = Ball1.ySpeed;
    int yVel2 = Ball2.ySpeed;
    int run = (Ball1.x-Ball2.x);
    int rise = (Ball1.y-Ball2.y);
    //???? ????? ???? ? ? ?????? ????????
    float Alfa = atan2f(rise, run);
    float cosAlfa = cos(Alfa);
    float sinAlfa = sin(Alfa);
    // ??????? ???????? ????? ????? ????????
    float xVel1prime = xVel1*cosAlfa+yVel1*sinAlfa;
    float xVel2prime = xVel2*cosAlfa+yVel2*sinAlfa;
    // ??????? ???????? ???????????????? ????? ????????
    float yVel1prime = yVel1*cosAlfa-xVel1*sinAlfa;
    float yVel2prime = yVel2*cosAlfa-xVel2*sinAlfa;
    // ????????? ?????? ??????????
    float P = (mass1*xVel1prime+mass2*xVel2prime);
    float V = (xVel1prime-xVel2prime);
    float v2f = (P+mass1*V)/(mass1+mass2);
    float v1f = v2f-xVel1prime+xVel2prime;
    xVel1prime = v1f;
    xVel2prime = v2f;
    // ?????????? ??????? ?? ??? ? ? ?.
    xVel1 = xVel1prime*cosAlfa-yVel1prime*sinAlfa;
    xVel2 = xVel2prime*cosAlfa-yVel2prime*sinAlfa;
    yVel1 = yVel1prime*cosAlfa+xVel1prime*sinAlfa;
    yVel2 = yVel2prime*cosAlfa+xVel2prime*sinAlfa;
    // ???????? ????????? ??????? ?????
    Ball1.xTemp = Ball1.x+Ball1.xSpeed*time;
    Ball1.yTemp = Ball1.y+Ball1.ySpeed*time;
    Ball2.xTemp = Ball2.x+Ball2.xSpeed*time;
    Ball2.yTemp = Ball2.y+Ball2.ySpeed*time;
    // ?????? ????? ???????? ?????
    Ball1.xSpeed = xVel1;
    Ball2.xSpeed = xVel2;
    Ball1.ySpeed = yVel1;
    Ball2.ySpeed = yVel2;
     
}