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

#include <iostream>
#include <math.h>
 
using namespace std;
void parall(double a1, double a2, double b1, double b2)
{
    if ((a1 / a2) == (b1 / b2)) 
    {
        cout << "Lines parallel!";
        system("pause");
        exit(0);
    }
}
void intersect(double a1, double a2, double b1, double b2, double c1, double c2, double& x, double& y)
{
    double det = a1 * b2 - a2 * b1;
    x = (b1 * c2 - b2 * c1) / det;
    y = (a2 * c1 - a1 * c2) / det;
}
 
double area(double x0, double x1, double x2, double y0, double y1, double y2)
{
    return abs(x0 * (y1 - y2) + x1 * (y2 - y0) + x2 * (y0 - y1)) / 2;
}
 
int main()
{
    double a1, a2, a3, b1, b2, b3, c1, c2 , c3;
    double x1=0, x2 = 0, x3 = 0, y1 = 0, y2 = 0, y3 = 0;
    // line equation ax + by + c = 0
    cout << " Input line 1 "; cin >> a1 >> b1 >> c1;
    cout << " Input line 2 "; cin >> a2 >> b2 >> c2;
    cout << " Input line 3 "; cin >> a3 >> b3 >> c3;
    parall(a1, a2, b1, b2);
    parall(a1, a3, b1, b3);
    parall(a2, a3, b2, b3);
    intersect(a1, a2, b1, b2, c1, c2, x1, y1);
    intersect(a1, a3, b1, b3, c1, c3, x2, y2);
    intersect(a2, a3, b2, b3, c2, c3, x3, y3);
    cout << "S=" << area(x1, x2, x3, y1, y2, y3) << endl;
    system("pause");
    return 0;
}