Найти координаты точек пересечения прямой y=kx+b и окружности

#include <iostream>
#include <cmath>
using namespace std;
int quadrant(double x, double y)
{
    int n;
    if (x > 0 && y > 0) n = 1;
    if (x > 0 && y < 0) n = 4;
    if (x < 0 && y > 0) n = 2;
    if (x < 0 && y < 0) n = 3;
    return n;
}
int main()
{
    double k = 1, b = 0, r = 1, d;
    double x1, x2, y1, y2;
    int num;
    d = -b * b + r * r + k * k * r * r;
 
    if (d >= 0)
    {
        x1 = (-b * k + sqrt(d)) / (1 + k * k);
        x2 = (-b * k - sqrt(d)) / (1 + k * k);
        y1 = k * x1 + b;
        y2 = k * x2 + b;
    }
    if (d > 0)
    {
        num = 2;
        cout << num << endl;
        cout << x1 << " " << y1 << " " << quadrant(x1, y1) << endl;
        cout << x2 << " " << y2 << " " << quadrant(x2, y2) << endl;
    }
    else if (d == 0)
    {
        num = 1;
        cout << num << endl;
        cout << x1 << " " << y1 << " " << quadrant(x1, y1) << endl;
    }
    else
    {
        num = 0;
        cout << "no intersections" << endl;
    }
 
    cin.get();
    return 0;
}