Построить квадрат по двум заданным точкам

 var AE = b.Y -a.Y;
            var BE = -(b.X-a.X);
            var CE = -(a.X*b.Y-a.Y*b.X);
 
            var xleftDown = (AE * (side - CE - BE * a.Y) + BE * BE * a.X) / (Math.Pow(AE, 2) + Math.Pow(BE, 2));
            var yleftDown = (AE * a.Y - BE * a.X + BE * xleftDown) / AE;
            var xrightDown = (AE * (side - CE - BE * b.Y) + BE * BE * b.X) / (Math.Pow(AE, 2) + Math.Pow(BE, 2));
            var yrightDown = (AE * b.Y - BE * b.X + BE * xrightDown) / AE;

/// <summary>
/// Функция вычисления вершин квадрата по двум смежным вершинам
/// </summary>
/// <param name="pt1">Первая вершина</param>
/// <param name="pt2">Вторая вершина</param>
/// <returns>Функция возвращает массив точек, являющихся вершинами квадрата</returns>
PointF[] GetSquareVertices(PointF pt1, PointF pt2) {
    using (System.Drawing.Drawing2D.GraphicsPath square = new System.Drawing.Drawing2D.GraphicsPath()) {
        //Разность координат (с учётом знака)
        float deltaX = pt1.X - pt2.X, deltaY = pt1.Y - pt2.Y;
        //Расстояние между точками
        float dist = (float)Math.Sqrt(Math.Pow(deltaX, 2) + Math.Pow(deltaY, 2));
        //Угол наклона стороны
        double angle;
        //Прямоугольник с вычисленной длиной стороны
        square.AddRectangle(RectangleF.FromLTRB(0, 0, dist, dist));
        //Вычисление угла поворота
        angle = Math.Atan(deltaY / deltaX);
        //Функция Math.Atan возвращает угол в пределах [-PI/2;PI/2], а для поворота нужен угол [0;2PI]. 
        //Для пересчёта смотрим знак разности координат.
        //Квадранты считаются с правого нижнего против часовой стрелки.
        //Первый квадрант [0;PI/2) получается сам собой при вычислении арктангенса
        if ((deltaX > 0 && deltaY <= 0) ||(deltaX >= 0 && deltaY > 0)) {//Второй квадрант [PI/2;PI) или третий квадрант [PI;3PI/2)
            angle = Math.PI + angle;
        }
        else if (deltaX <= 0 && deltaY > 0) {//Четвёртый квадрант [3PI/2;2PI)
            angle = 2 * Math.PI + angle;
        }
        //Перевод угла из радиан в градусы
        angle *= 180 / Math.PI;
        //Перенос начала координат в первую точку и
        //поворот прямоугольника на вычисленный угол относительно этого нового начала координат
        using (System.Drawing.Drawing2D.Matrix m = new System.Drawing.Drawing2D.Matrix(1, 0, 0, 1, pt.X, pt.Y)) {
            m.RotateAt((float)angle, new PointF());
            square.Transform(m);
        }
        return square.PathPoints;
    }
}

private double SquareSide
     {
         get { return Math.Pow(_b.X - _a.X, 2) + Math.Pow(_b.Y - _a.Y, 2); }
     }
 
public void DoCalculate()
        {
            double side = SquareSide;
            //коефіцієнти рівняння загального
            int ae = _b.Y - _a.Y;
            int be = -(_b.X - _a.X);
            int ce = -(_a.X*_b.Y - _a.Y*_b.X);
 
            var xleftDown = (int) ((ae*(side - ce - be*_a.Y) + be*be*_a.X)/(Math.Pow(ae, 2) + Math.Pow(be, 2)));
            int yleftDown = (ae*_a.Y - be*_a.X + be*xleftDown)/ae;
            var xrightDown = (int) ((ae*(side - ce - be*_b.Y) + be*be*_b.X)/(Math.Pow(ae, 2) + Math.Pow(be, 2)));
            int yrightDown = (ae*_b.Y - be*_b.X + be*xrightDown)/ae;
 
//результати 
            _d = new Point(xleftDown, yleftDown);
            _c = new Point(xrightDown, yrightDown);
        }
    }