Алгоритм построения выпуклой оболочки (на плоскости)

#include <iostream> 
#include <fstream> 
#include <sstream> 
using namespace std; 
 
#define INF 10000 
 
struct Point 
{ 
    int x; 
    int y; 
}; 
 
int orientation(Point p, Point q, Point r) 
{ 
    int val = (q.y - p.y) * (r.x - q.x) - (q.x - p.x) * (r.y - q.y); 
 
    if (val == 0)   // если они колинеарны
        return 0; 
    return (val > 0) ? 1 : 2; //если больше нуля,то по часовой стрелке(1). Если меньше, то против часовой(2)
} 
 
void convexHull(Point points[], int n) 
{ 
    if (n < 3) 
    return; 
 
 
    int* next = new int[n]; 
    for (int i = 0; i < n; i++) 
    next[i] = -1; 
 
 
    int l = 0; 
    for (int i = 1; i < n; i++) 
        if (points[i].x < points[l].x) 
            l = i; 
 
 
    int p = l, q; 
    do 
    { 
 
        q = (p + 1) % n; 
        for (int i = 0; i < n; i++) 
            if (orientation(points[p], points[i], points[q]) == 2) 
                q = i; 
 
        next[p] = q; 
        p = q; 
    } while (p != l); 
     
    for (int i = 0; i < n; i++) 
    { 
        if (next[i] != -1) 
            cout << " (" << points[i].x << ", " << points[i].y << ")" << endl; 
    } 
} 
 
int main() 
{ 
    setlocale(LC_ALL, "rus");
    
    Point points[] = {{0, 3}, {1, 1}, {2, 2}, {4, 4}, 
                      {0, 0}, {1, 2}, {3, 1}, {3, 3}};  
    
    cout << " Точки входящие в выпуклую оболочку " << endl; 
    int n = sizeof(points) / sizeof(points[0]); 
    convexHull(points, n); 
    
    
    return 0; 
}