Найти расстояние от данной точки до ближайшей стороны треугольника

#include <iostream>
#include <cmath>
 
using namespace std;
 
double dist(double x1, double y1, double x2, double y2){
    
    return sqrt((x1 - x2)*(x1 - x2) + (y1 - y2)*(y1 - y2));
}
 
double s_n(double x[], double y[], int n = 3){
    //ïëîùàäü n - óãîëüíèêà
    double sum = 0;
    for (int i = 0; i < n; ++i){
        double det = x[i%n] * y[(i + 1)%n] - x[(i + 1)%n] * y[i%n];
        sum += det;
    }
    return fabs(sum/2);
}
 
int main(){
    
    double tx[3], ty[3], x, y;
    
    cout << "Ââåäèòå êîîðäèíàòû âåðøèí òðåóãîëüíèêà: " << endl;
    for (int i = 0; i < 3; ++i){
        
        cout << "x" << i + 1 << " = "; cin >> tx[i];
        cout << "y" << i + 1 << " = "; cin >> ty[i];
    }
    
    cout << "Ââåäèòå êîîðäèíàòû òî÷êè: " << endl;
    cout << "x = "; cin >> x;
    cout << "y = "; cin >> y;
    
    double min = 1e20;
    for (int i = 0; i < 3; ++i){
        
        double tempx[3], tempy[3];
        for (int j = 0; j < 3; ++j){
            tempx[j] = tx[j];
            tempy[j] = ty[j];
        }
        tempx[i] = x;
        tempy[i] = y;
        
        double d = dist(tx[(i + 1)%3], ty[(i + 1)%3], 
                        tx[(i + 2)%3], ty[(i + 2)%3]);
                        
        if (2*s_n(tempx, tempy)/d < min)
            min = 2*s_n(tempx, tempy)/d;
    }
    
    cout << "H = " << min << endl;
}