Пересечение эллипса и окружности

float alpha, beta, gamma, D, y1, y2, x1, x2, X_c, Y_c, X_e, Y_e; //X_c Y_c точка центра окружности 
    alpha = pow(Rx,2)-pow(Ry,2); //Коэф А                                 //X_e Y_e точка центра эллипса 
    beta = 2*Y_c*pow(Ry,2) - 2*Y_e*pow(Rx,2); // B
    gamma = pow(Y_e,2)*pow(Rx,2)-pow(Rx,2)*pow(Ry,2)+pow(R,2)*pow(Ry,2)-pow(Y_c,2)*pow(Ry,2); //C
    D = pow(beta,2)-4*alpha*gamma;
    if(D>0) { 
        y1 = (-beta+sqrt(D))/(2*alpha); 
        y2 = (-beta-sqrt(D))/(2*alpha); 
        x1 = sqrt(fabs(pow(R,2)-pow(y1-Y_c,2)))+X_c;
        x2 = sqrt(fabs(pow(R,2)-pow(y2-Y_c,2)))+X_c; 
        bool i1,i2,i3,i4,j1,j2,j3,j4;
        i1 = VpisO(X_c, Y_c, R, x1, y1); j1 = VpisE(X_e, Y_e, Rx, Ry, x1, y1); //VpisO и VpisE
        i2 = VpisO(X_c, Y_c, R, x1, y2); j2 = VpisE(X_e, Y_e, Rx, Ry, x1, y2); //Проверяют точки
        i3 = VpisO(X_c, Y_c, R, x2, y1); j3 = VpisE(X_e, Y_e, Rx, Ry, x2, y1); // на принадлежность
        i4 = VpisO(X_c, Y_c, R, x2, y2); j4 = VpisE(X_e, Y_e, Rx, Ry, x2, y2); // Окружности и эллипсу соответственно
        if(i1 && j1 || i2 && j2 || i3 && j3 || i4 && j4) return true; //проверка на принадлежность точек фигурам 
        else return false;
      }

bool Collision(здесь передаем координаты центров фигур и их радиусы)
{
           /*Cначало исключаем варианты когда эллипс описан или вписан*/
                bool i1, i2, i3, i4, i5, j1;
        i1 = VpisO(X_c, Y_c, R, X_e, Y_e); 
        i2 = VpisO(X_c, Y_c, R, X_e+Rx, Y_e+Ry);
        i3 = VpisO(X_c, Y_c, R, X_e-Rx, Y_e-Ry);
        i4 = VpisO(X_c, Y_c, R, X_e-Rx, Y_e+Ry);
        i5 = VpisO(X_c, Y_c, R, X_e+Rx, Y_e+Ry);
        if(i1 && i2 && i3 && i4 && i5 || i1 && !i2 && !i3 && !i4 && !i5)
            return false;/*Если эллипс вписанный или описанный, то пересечения нет*/
        else
        {
            double x[361], y[361];
                        //строим эллипс "поточечно" с помощью уравнения эллипса в парам-ой форме
            for(int a = 0; a <= 360; a++)//a - угол альфа [0; 2pi]
            {
                x[a] = X_e + Rx * cos(a); y[a] = Y_e+Ry * sin(a);
                TPrintXML(x[a], y[a]); //Тестовая функция. Рисует эллипс "поточечно"
            }
            for(int a = 0; a <= 360; a++) //проверка всех точек
            {
                j1 = VpisO(X_с, Y_c, R, x[a],y[a]); 
                if(j1) return true;//если хотя бы одна из точек на поверхности эллипса принадлежит окружности
            }
        }
}