Нарисовать треугольник с вписанной в него окружностью

using System;
using System.Drawing;
using System.Windows.Forms;
 
namespace WindowsFormsApplication1 {
    public partial class Form1 : Form {
        public Form1() {
            InitializeComponent();
            pictureBox1.Paint += pictureBox1_Paint;
        }
 
        void pictureBox1_Paint(object sender, PaintEventArgs e) {
            DrawTriangle(e);
        }
 
        private void DrawTriangle(PaintEventArgs e) {
            PointF[] points = new[] { new PointF(10, 10), new PointF(200, 45), new PointF(100, 300) };
            e.Graphics.DrawPolygon(Pens.Blue, points);
            e.Graphics.DrawEllipse(Pens.Red, InscribedCircle(points[0], points[1], points[2]));
        }
        /// <summary>
        /// Вычисление прямоугольника, ограничивающего окружность, вписанную в треугольник
        /// </summary>
        /// <param name="pt1">Координаты первой вершины треугольника</param>
        /// <param name="pt2">Координаты второй вершины треугольника</param>
        /// <param name="pt3">Координаты третьей вершины треугольника</param>
        static RectangleF InscribedCircle(PointF pt1, PointF pt2, PointF pt3) {
            Rectangle result = new Rectangle();
            //Векторы сторон
            PointF p1 = new PointF(pt2.X - pt1.X, pt2.Y - pt1.Y);
            PointF p2 = new PointF(pt3.X - pt2.X, pt3.Y - pt2.Y);
            PointF p3 = new PointF(pt1.X - pt3.X, pt1.Y - pt3.Y);
            //Длины векторов сторон
            float l1 = (float)Math.Sqrt(p1.X * p1.X + p1.Y * p1.Y),
                l2 = (float)Math.Sqrt(p2.X * p2.X + p2.Y * p2.Y),
                l3 = (float)Math.Sqrt(p3.X * p3.X + p3.Y * p3.Y);
            //единичные векторы сторон
            PointF ep1 = new PointF(p1.X / l1, p1.Y / l1),
                ep2 = new PointF(p2.X / l2, p2.Y / l2),
                ep3 = new PointF(p3.X / l3, p3.Y / l3);
            //векторы биссектрис. Считаем только для углов 2 и 3
            PointF b2 = new PointF(-ep1.X + ep2.X, -ep1.Y + ep2.Y),
                b3 = new PointF(ep3.X - ep2.X, ep3.Y - ep2.Y);
            //Вычисляем точку пересечения биссектрис. Биссектрисы заданы векторами b2 и b3
            //и проходят через точки pt2 и pt3 соответственно. 
            //Вспомогательные переменные для удобства записи
            float k2 = b2.Y / b2.X, k3 = b3.Y / b3.X;
            //Координаты пересечения биссектрис. Формулы выводятся из уравнения прямой 
            //по точке и направляющему вектору
            float x = (k2 * pt2.X - k3 * pt3.X + pt3.Y - pt2.Y) / (k2 - k3), y = k2 * (x - pt2.X) + pt2.Y;
            //Радиус вписанной окружности.
            //Расстояние от точки пересечения до любой из сторон
            //Считаем расстояние до стороны p1
            float radius = (float)(Math.Abs(x * p1.Y - y * p1.X - pt1.X * p1.Y + pt1.Y * p1.X) / Math.Sqrt(p1.X * p1.X + p1.Y * p1.Y));
            return new RectangleF(x - radius, y - radius, 2 * radius, 2 * radius);
        }
    }
}