Соединить конечное множество точек на плоскости замкнутой линией с вершинами в этих точках

#include <iostream>
#include <cmath>
#include <time.h>
 
struct Point
{
    int ind;     // исходный индекс точки, что бы сохранить порядок при сортировке
    double x;    // координата X
    double y;    // координата Y
    double a;    // угол в градусах
    double r;    // радиус
 
    void conv(double xc = 0.0, double yc = 0.0) // посчитать <a> и <r>
    {
        r = sqrt((x - xc)*(x - xc) + (y - yc)*(y - yc));
 
        if (r == 0.0)
        {
            a = 0.0;
        }
        else
        {
            a = 180.0 * acos((x - xc) / r) / acos(-1.0);
            if ((y - yc) < 0)
                a = 360 - a;
        }
    }
 
    int compare(Point &right) // сравнить this (левый операнд) с указанной структурой <right> (правый операнд)
    {
        if (a < right.a)
            return -1; // меньше
        if (a > right.a)
            return 1; // больше
 
        // углы равны, сравниваем радиус
 
        if (r < right.r)
            return -1; // меньше
        if (r > right.r)
            return 1; // больше
 
        return 0; // точки совпадают
    }
};
 
#define POINTS_COUNT 4
 
int main()
{
    srand((unsigned int)time(0));
 
    Point points[POINTS_COUNT] = { { 0, 2, 5 },{ 1, 2, 2 },{ 2, 6, 3 },{ 3, 5, 0 } };
 
    // геометрический центр (центр масс)
    double xc = 0;
    double yc = 0;
    // инициализация
    for (int i = 0; i < POINTS_COUNT; i++)
    {
        points[i].ind = i;
        //points[i].x = 200 * rand() / RAND_MAX - 100;
        //points[i].y = 200 * rand() / RAND_MAX - 100;
 
        xc += points[i].x;
        yc += points[i].y;
    }
    xc /= POINTS_COUNT;
    yc /= POINTS_COUNT;
    std::cout << "xc=" << xc << " yc=" << yc << std::endl;
 
    // полярные координаты относительно центра масс
    for (int i = 0; i < POINTS_COUNT; i++)
    {
        points[i].conv(xc, yc);
        std::cout << i + 1 << ") (x=" << points[i].x << ", y=" << points[i].y << ") - (a=" << points[i].a << ", r=" << points[i].r << ")" << std::endl;
    }
 
    // сортировка (вставками)
    for (int i = 1; i < POINTS_COUNT; i++)
    {
        for (int j = i; (j > 0) && (points[j].compare(points[j - 1]) < 0); j--)
        {
            Point temp = points[j-1];
            points[j-1] = points[j];
            points[j] = temp;
        }
    }
 
    // результат сортировки
    std::cout << std::endl;
    for (int i = 0; i < POINTS_COUNT; i++)
    {
        std::cout << points[i].ind + 1 << " ";
    }
    std::cout << std::endl << std::endl;
 
    system("pause");
    return 0;
}