Матрица поворота. Вращение прямоугольника вокруг точки пересечения диагоналей

// вращение будет происходить вокруг точки around
                for (int i = 0; i < points.Length; i++)
                {// перебираем точки многоугольника
                    PointF pt = points[i];
                    pt.X -= around.X;       // преобразование координат в систему координат с началом в точке around
                    pt.Y -= around.Y;
 
                    PointF temp = new PointF(
                        (float)(pt.X * cos_alpha + pt.Y * sin_alpha),
                        (float)(pt.Y * cos_alpha - pt.X * sin_alpha)
                        );      // применяем матрицу поворота
 
                    temp.X += around.X;     // обратное преобразование координат
                    temp.Y += around.Y;
 
                    points[i] = temp;
                }

public override PointF Center
{
    get
    {
        float       // приводим координаты прямоугольника в систему координат с началом в левом верхнем его углу
            x1 = points[0].X - points[0].X,     // верхняя левая вершина
            y1 = points[0].Y - points[0].Y,
 
            x2 = points[2].X - points[0].X,     // нижняя правая
            y2 = points[2].Y - points[0].Y,
 
            n1 = points[3].X - points[0].X,     // нижняя левая
            m1 = points[3].Y - points[0].Y,
 
            n2 = points[1].X - points[0].X,     // верхняя правая
            m2=points[1].Y - points[0].Y
        ;
        float X = 
            ((x1*y1+x1*y1)/(x2-x1)+(n1*m2+n1*m1)/(n2-n1)+m1-y1)/
            ((y2-y1)/(x2-x1)-(m2-m1)/(n2-n1))
            ;       // эту формулу я вывел из уравнения прямой, проходящей через две точки
        float Y = 
            ((X-x1)*(y2-y1))/
            (x2-x1)
            +y1;        // всё из того же уравнения получаем Y
        X += points[0].X;       // обратное преобразование координат
        Y += points[0].Y;
 
        return new PointF(X, Y);        // возвращаем точку пересечения диагоналей
    }
}

public override PointF GetCenter()
{
    float left = Math.Min(Math.Min(points[0].X, points[1].X), Math.Min(points[2].X, points[3].X));
    float right = Math.Max(Math.Max(points[0].X, points[1].X), Math.Max(points[2].X, points[3].X));
    float top = Math.Min(Math.Min(points[0].Y, points[1].Y), Math.Min(points[2].Y, points[3].Y));
    float bottom = Math.Max(Math.Max(points[0].Y, points[1].Y), Math.Max(points[2].Y, points[3].Y));
 
    return new PointF((right + left) / 2, (bottom + top) / 2);
}