Найти координаты противоположных углов квадрата если известны координаты A и B

#include <iostream>
#include <cmath>
 
using namespace std;
 
float  square_func(float X1,float X2,float Y1,float Y2){
float length=sqrt((X1-X2)*(X1-X2)+(Y1-Y2)*(Y1-Y2));
 
 
 
    }
int main()
{
    float X1, X2,Y1, Y2;
    cout << "Vdtdite X1Y1 !" << endl;
    cin>>X1;
    cin>>Y1;
    cout << "Vdtdite X2Y2 !" << endl;
    cin>>X2;
    cin>>Y2;
    square_func(X1,X2,Y1,Y2);
 
    return 0;
}

#include <iostream>
#include <iomanip>
 
using namespace std;
 
typedef struct _point
{
    float   x;
    float   y;
}point;
 
bool square_func(point &p1, point &p2, point &p3, point &p4) 
{
    float   lenght = sqrt((p1.x - p2.x)*(p1.x - p2.x) + (p1.y - p2.y)*(p1.y - p2.y));
    float   eps = 0.000001;
 
    if (lenght)
    {
        if (fabs(p1.x - p2.x) < eps)
        {
            p4.x = p3.x = p2.x + lenght;
            p3.y = p2.y;
            p4.y = p1.y;
        }
        else if (fabs(p1.y - p2.y) < eps)
        {
            p4.y = p3.y = p2.y + lenght;
            p3.x = p2.x;
            p4.x = p1.x;
        }
        else
        {
            float k = (p2.x - p1.x) / (p2.y - p1.y);
            p4.x = p1.x + (lenght / sqrt(k*k + 1));
            p4.y = p1.y - k * (p4.x - p1.x);
            p3.x = p2.x + (lenght / sqrt(k*k + 1));
            p3.y = p2.y - k * (p3.x - p2.x);
        }
        return true;
    }
    else
        return false;
}
int main()
{
    
    point p1, p2, p3, p4;
 
    setlocale(LC_ALL, "");
 
    cout << "Введите X1Y1 !" << endl;
    cin >> p1.x;
    cin >> p1.y;
    cout << "Введите X2Y2 !" << endl;
    cin >> p2.x;
    cin >> p2.y;
    
    if (square_func(p1, p2, p3, p4))
    {
        cout << setw(5) << "x1 = " << p1.x << " " << "y1 = " << p1.y << endl;
        cout << setw(5) << "x2 = " << p2.x << " " << "y2 = " << p2.y << endl;
        cout << setw(5) << "x3 = " << p3.x << " " << "y3 = " << p3.y << endl;
        cout << setw(5) << "x4 = " << p4.x << " " << "y4 = " << p4.y << endl;
    }
 
    return 0;
}

#include <iostream>
 
struct Vector {
  double x;
  double y;
};
 
Vector operator+(const Vector &a, const Vector &b)
{
  return {a.x + b.x, a.y + b.y};
}
 
 
Vector operator-(const Vector &a, const Vector &b)
{
  return {a.x - b.x, a.y - b.y};
}
 
Vector rotate90(const Vector &v)
{
  return {-v.y, v.x};
}
 
std::ostream &operator<<(std::ostream &o, const Vector &v)
{
  o << "{x:" << v.x << ", y:" << v.y << "}";
  return o;
}
 
void findOtherCorners(const Vector &a, const Vector &b, Vector &c, Vector &d)
{
  Vector v_ab = b - a;
  Vector v_bc = rotate90(v_ab);
  c           = b + v_bc;
  d           = a + v_bc;
}
 
int main()
{
  Vector a = {0.0, 0.0};
  Vector b = {1.0, 0.0};
  Vector c;
  Vector d;
 
  findOtherCorners(a, b, c, d);
  std::cout << "a = " << a << std::endl;
  std::cout << "b = " << b << std::endl;
  std::cout << "c = " << c << std::endl;
  std::cout << "d = " << d << std::endl;
 
  return 0;
}