Найти точки пересечения гиперболы и окружности

#include <iostream>
#include "conio.h"
#include <math.h>
 
 
double findIntersection(int, int, double, int,double);
 
using namespace std;
int main() 
{
    const int N = 360; 
    const double PI = 3.14159265;
    int cR, step = 1;
        double k;
    cout<<"Circle's radius?"<<endl;
    cin>>cR;
        cout<<"giperbola"<<endl;
        cin>>k;
 
    double cXs[N], cYs[N], fYs[N];
    for(int i = 0;i < N; i+=step){
 
        cXs[i] = cR*cos(i*PI/180);
        cYs[i] = cR*sin(i*PI/180);
        fYs[i] = k/(cXs[i]);
    }
 
    int K = 0; int intersections[10];
    for(int i = 0; i < N-1; i++){
 
        if((cYs[i]<=fYs[i])&&(cYs[i+1]>=fYs[i+1]))
        { intersections[K] = i; K++;}
        if((cYs[i]>=fYs[i])&&(cYs[i+1]<=fYs[i+1]))
        {intersections[K] = i;  K++;}
    }
 
    int startI, endI; 
    for(int i = 0; i < K; i++){
        (intersections[i] > 0)?(startI = intersections[i]-1):(startI = intersections[i]);
        (intersections[i] < (N - 1))?(endI = intersections[i] + 1):(endI = intersections[i]);
 
        cout<<"bestAngle: "<<"\t";
        double bestAngle = findIntersection(startI,endI,PI,cR,k);
 
        cout<<bestAngle<<" rad= "<<bestAngle*(180/PI)<<" grad"<<endl;
                cout<<"best value of X: "<<cR*cos(bestAngle)<<endl;
        cout<<"best value of Y: "<<cR*sin(bestAngle)<<"\t";
        cout<<endl;
    }
    _getch();
}
double findIntersection(int startI, int endI, double PI, int cR,double k){
    int points = 1000; 
 
    double cX, cY, fY, bestResidual = 100, bestAngle = -1, step = (PI/180)*(endI - startI)/points;
    for(double angle = startI*PI/180; angle < endI*PI/180; angle+=step) {
        cX = cR*cos(angle);
        cY = cR*sin(angle);
        fY = k/cX;
        if(fabs(cY - fY) < bestResidual)
        {
            bestResidual = fabs(cY - fY);
            bestAngle = angle;
        }
    }
    return bestAngle;
}