Постройте прямоугольник по радиусу описанной окружности и двум соседним вершинам

#include<windows.h>
#include<iostream>
using namespace std;
 
class Point {
    double x, y;
public:
    Point(double X, double Y){
        x = X; y = Y;
    };
    Point(){};
    double getX(){
        return x;
    };
    double getY(){
        return y;
    };
    void setX(double X){
        x = X;
    };
    void setY(double Y){
        y = Y;
    };
    friend istream &operator >>(istream &stream, Point &poi){
        cout << "X: ";
        cin >> poi.x;
        cout << "Y: ";
        cin >> poi.y;
        return stream;
    };
    friend ostream &operator <<(ostream &stream, Point poi){
        cout << poi.x << ", " << poi.y;
        return stream;
    };
};
 
class vector{
    double x;
    double y;
public:
    vector(Point a, Point b){
        x = b.getX() - a.getX();
        y = b.getY() - a.getY();
    };
    vector(){};
    double getX(){
        return x;
    };
    double getY(){
        return y;
    };
};
 
int main(){
    setlocale(LC_ALL, "RUS");
 
    Point A = Point(),
          B = Point(),
          C1 = Point(), 
          C2 = Point(),
          D1 = Point(),
          D2 = Point();
    double R;
 
    cout << "Введите точку A: " << endl;
    cin >> A;
    cout << "Введите точку B: "<<endl;
    cin >> B;
    cout << "Введите радиус: " << endl;
    cin >> R;
 
    vector a(A,B);      //vector AB
 
    if (R < sqrt(pow(a.getX(), 2) + pow(a.getY(), 2)) / 2 || R >= sqrt(pow(a.getX(), 2) + pow(a.getY(), 2))) {
        cout << "Окружность с таким радиусом не описывает прямоугольник." << endl;
        system("Pause");
        return 0;
    };
    
    double L = sqrt(4 * R*R - (pow(a.getX(), 2) + pow(a.getY(), 2))) / sqrt(pow(a.getX(), 2) + pow(a.getY(), 2));    //лямбда
 
    C1.setX(B.getX() + a.getY()*L);      //X=Xo+a1*t*L (t=1 => точка слева)
    C1.setY(B.getY() + a.getX()*L*(-1)); 
 
    C2.setX(B.getX() - a.getY()*L);      //X=Xo+a1*t*L (t=-1 => точка справа)
    C2.setY(B.getY() - a.getX()*L*(-1));      
 
    D1.setX(A.getX() + a.getY()*L);  
    D1.setY(A.getY() + a.getX()*L*(-1));
 
    D2.setX(A.getX() - a.getY()*L);
    D2.setY(A.getY() - a.getX()*L*(-1));
    
    cout << "Первый вариант расположения точек: " << endl;
    cout << "Точка С1: (" << C1 << "); Точка D1: (" << D1 <<")." << endl;
 
    cout << "Второй вариант расположения точек: " << endl;
    cout << "Точка С2: (" << C2 << "); Точка D2: (" << D2 << ")." << endl;
    
    system("Pause");
    return 0;
}