Создать функцию определяющую расстояние между двумя точками

#include <iostream>
#include <vector>
#include <algorithm>
#include <cmath>
#define eps 0.00000001
using namespace std;
 
double dictance (double x1, double y1, double x2, double y2) {
    return sqrt(pow((x2-x1),2)+pow((y2-y1),2));
}
 
double square(double x, double y, double z)
{
    double p;
    p = (x + y + z) / 2;
    return sqrt(p*(p - x)*(p - y)*(p - z));
}
 
int main() {
    int n;
     struct Point { //Объявляем структуру Point(хранит координаты точки)
        int x;
        int y;
    };
    double a1,b1,a2,b2,q1,w1,e1;
    cin>> a1>>b1>>a2>>b2;
    cout<<dictance (a1,b1,a2,b2);
    return 0;
    q1=distance
    cout<< square(q1,w1,e1);
    }
    
Пока что только так. Дальше не могу понять. Алгоритм какой должен быть понимаю, но как осуществить задачу - нет. Простите, если это глупо, я только 1 курс
 
[size="1"][color="grey"][I]Добавлено через 14 минут[/I][/color][/size]
#include <cstdlib>
#include <iostream>
#include <cmath>
 
 
using namespace std;
 
void Input(int **arrayDots, const size_t n)
{
    for(int i = 0; i < n; i++)
    {
        cout<<"Введите координаты точки "<<i+1<<endl;
 
        for(int j = 0; j < 2; j++)
        {
            cin >> arrayDots[i][j];
        }
    }
}
 
void Output(int **arrayDots, const size_t n)
{
    cout<<"Координаты введеных точек:"<<endl;
 
    for(int i = 0; i < n; i++)
    {
        for(int j = 0; j < 2; j++)
        {
            cout << arrayDots[i][j] << ' ';
        }
        cout << '\n';
    }
}
 
//нахождение площади треугольника по формуле Герона
double Square (double a, double b, double c)
{
    double p = (a + b + c)/2;
    double s = sqrt ( p*(p-a)*(p-b)*(p-c) );
 
    return s;
}
 
//проверка на существование треугольника из заданых отрезков
bool isTriangle (double a, double b, double c)
{
    if (a >= b+c) return false;
    if (b >= a+c) return false;
    if (c >= a+b) return false;
 
    return true;
}
 
int main()
{  
    setlocale(LC_CTYPE,"Russian");
 
    int n = 0;
 
    do 
    {
        cout<<"Введите количество точек: ";
        cin>>n;
 
        if (n < 3)
            cout<<"Неверное количество точек"<<endl<<endl;
    } while (n < 3);
 
    int **arrayDots = NULL;
 
    arrayDots = new int * [n];
 
    for (int i = 0; i < n; i++)
    {
        arrayDots[i] = new int[2];
    }
 
    Input(arrayDots, n);
 
    Output(arrayDots, n);
 
    double maximalSquare = 0,
        a = 0,//1 - я сторона
        b = 0,//2 - я
        c = 0,//3 - я
        sqrX = 0,//квадрат разности координаты Х между 2 - мя точками
        sqrY = 0;//квадрат разности координаты У между 2 - мя точками
 
    for (int c1 = 0; c1 < n; c1++)
    {
        for (int c2 = 0; c2 < n; c2++)
        {
            for (int c3 = 0; c3 < n; c3++)
            {
                sqrX = pow( (double)abs(arrayDots[c1][0] - arrayDots[c2][0]), 2 );
                sqrY = pow( (double)abs(arrayDots[c1][1] - arrayDots[c2][1]), 2 );
                a = pow( sqrX + sqrY, 0.5);//длина стороны между точками с1 и с2
 
                sqrX = pow( (double)abs(arrayDots[c1][0] - arrayDots[c3][0]), 2 );
                sqrY = pow( (double)abs(arrayDots[c1][1] - arrayDots[c3][1]), 2 );
                b = pow( sqrX + sqrY, 0.5);//длина стороны между точками с1 и с3
 
                sqrX = pow( (double)abs(arrayDots[c2][0] - arrayDots[c3][0]), 2 );
                sqrY = pow( (double)abs(arrayDots[c2][1] - arrayDots[c3][1]), 2 );
                c = pow( sqrX + sqrY, 0.5);//длина стороны между точками с2 и с3 
 
                if ( isTriangle (a, b, c) )
                {
                    double _temp = Square( a, b, c);
 
                    if( _temp > maximalSquare )                  
                    {
                        maximalSquare = _temp;  
                    }                
                }
            }
        }
    }
 
    cout << "Максимальная площадь треугольника: " << maximalSquare << endl;
 
    
}

#include <iostream>
typedef struct  { //Объявляем структуру Point(хранит координаты точки)
int x;
int y;
} Point;
double dictance (Point p1, Point p2)
{
  return sgrt((p1.x -p2.x)*(p1.x - p2.x) + (p1.y -p2.y)*(p1.y - p2.y));
}
int main() {
   Point p[3];
   double d[3];
   for(int i=0; i<3; i++)
     cin >> p[i].x << p[i].y;
   d[0] = distance(p[0], p[1]);
   d[1] = distance(p[1], p[2]);
   d[2] = distance(p[2], p[0]);
   double per = 0;
   for(int i =0; i<3; i++)
       per += d[i];
   s = sqrt (per*(per=d[0]) *(per - d[1]) (per - d[2]));
   // ....
  return 0;
}

Point P[n], pi, pj, pk;
for (int i=0; i<n-2; i++) {
   pi = P[i];
   for (int j = i+1; j<n-1; j++) {
     pj = P[j]
     for (int k=j+1; k <n; k++) {    
        pk = P[k];
         // Далее работаешь с треугольником pi, pj, pk 
      }
   }
}