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

#include "stdafx.h"
#include <stdio.h>
#include <conio.h>
#include <math.h>
#include <Windows.h>
#include <iostream>
 
void main(void)
{
    SetConsoleCP(1251);
    SetConsoleOutputCP(1251);
    float X, Y, X1, Y1;
 
    float R=31;
    float _x0, _y0;
    _x0=32-R;
    float Dokr;
    Dokr = (-40)*(-40)-4*1*(-560);
    _y0=((-40)+sqrt(Dokr))/2;
    printf("центр окружности имеет кооридинаты х0=%0.3f y0=%0.3f\n",_x0, _y0);/*центр окружности*/
    //pow(X-_x0, 2) + pow(Y-_y0, 2) - R*R; //формула окружности
 
 
    float a=31, b=38;
    float _x, _y;
    float Dell;
    _y=b-45;
    Dell = (-80)*(-80)-4*1*671;
    _x=(-80+sqrt(Dell))/2;
    printf("центр эллипса имеет координаты x=%0.3f y=%0.3f\n", _x, _y);/*центр эллипса*/
    //pow(X1+_x, 2)/a*a + pow(Y1+_y, 2)/b*b - 1;//формула эллипса
 
    system("pause");
}

#include "stdafx.h"
#include <stdio.h>
#include <conio.h>
#include <math.h>
#include <Windows.h>
#include <iostream>
 
void main(void)
{
    SetConsoleCP(1251);
    SetConsoleOutputCP(1251);   
                                                            
    float R=31;
    float _x0, _y0;
    _x0=32-R;
    float Dokr;
    Dokr = (-40)*(-40)-4*1*(-560);
    _y0=((-40)+sqrt(Dokr))/2;
        float krayokruzhnostiX1, krayokruzhnostiX2;
        krayokruzhnostiX1=_x0-R;
        krayokruzhnostiX2=_x0+R;
    
    float a=31, b=38;
    float _x, _y;
    float Dell;
    _y=b-45;
    Dell = (-80)*(-80)-4*1*671;
    _x=(-80+sqrt(Dell))/2;
        float krayellipsaX1, krayellipsaX2;
        krayellipsaX1=_x-a;
        krayellipsaX2=_x+a;
        
        for( float X1=krayokruzhnostiX1; X1<=krayokruzhnostiX2; X1=X1+0.001)
        {
            float Yokr1, Yokr2; 
            float t, w, D;
            float n1=1;//типо a
            w=2*_y0;//типо b
            t=pow(X1-_x0, 2)+pow(_y0, 2)-pow(31, 2);// типо c
            D=w*w-4*n1*t;//дискриминант
            Yokr1=(-w+sqrt(D))/2*n1;//первый корень окружности
            Yokr2=(-w-sqrt(D))/2*n1;//второй корень окружности
            
 
                for(float X2=krayellipsaX1; X2<=krayellipsaX2; X2=X2+0.001)
                {
                    float znachenieellipsa=1;
                    float Yell1, Yell2;
                    float t1, w1, D1;
                    float n2=1;//типо а
                    w1=2*_y;//типо b
                    t1=((pow(X2-_x, 2))/(a*a)-znachenieellipsa)*(b*b)+(_y*_y);//типо с
                    D1=w1*w1-4*n2*t1;
                    Yell1=(w1+sqrt(D1))/2*n2;//первый корень эллипса
                    Yell2=(w1-sqrt(D1))/2*n2;//второй корень эллипса
                    
                    if(Yokr1==Yell1 && Yokr2==Yell2)
                        //printf("точки пересечения\n\t\X1=%0.3f, Y1=%0.3f\n\tX2=%0.3f, Y2=%0.3f", X1, Yokr1, X2, Yokr2);
                    break;
                    printf("точки пересечения\n\t\X1=%0.3f, Y1=%0.3f\n\tX2=%0.3f, Y2=%0.3f", X1, Yokr1, X2, Yokr2);
                }
        }
    
    system("pause");
}

if(fabs(Yokr1-Yell1)<0.03)
{
...
// найден первый корень
}else if(fabs(Yokr2-Yell2)<0.03)
{
...
//второй корень
}