Сравнени двух матриц (NxN-1) и (MxM-1)

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
 
#define A "A.txt"
#define B "B.txt"
 
// Связный список для хранения координат точек
struct plane
{
    double x;
    double y;
    struct plane *next;
};
 
struct plane *makePlane(char *filename, int *count);
void deletePlane(struct plane *start);
void distance(double **matrix, int dim, struct plane *start);
 
int main()
{
    struct plane *startA = NULL; // Начальный элемент списка точек первой плоскости
    struct plane *startB = NULL; // Начальный элемент списка точек второй плоскости
    double **distA;              // Матрица расстояний между всеми точками первой плоскости
    double **distB;              // Матрица расстояний между всеми точками второй плоскости
    int N, M;                    // Счётчики количества точек
    int i, j;
 
    startA = makePlane(A, &N); // Создаём список точек первой плоскости
    startB = makePlane(B, &M); // Создаём список точек второй плоскости
 
    distA = (double **)malloc(N * sizeof(double *));
 
    for (i = 0; i < N; i++)
        distA[i] = (double *)malloc((N - 1) * sizeof(double));
 
    distB = (double **)malloc(M * sizeof(double *));
 
    for (i = 0; i < M; i++)
        distB[i] = (double *)malloc((M - 1) * sizeof(double));
 
    distance(distA, N, startA);
    distance(distB, M, startB);
 
    for (i = 0; i < N; i++)
    {
        for (j = 0; j < N - 1; j++)
            printf("%lf  ", distA[i][j]);
 
        printf("\n");
    }
 
    printf("\n");
 
    for (i = 0; i < M; i++)
    {
        for (j = 0; j < M - 1; j++)
            printf("%lf  ", distB[i][j]);
 
        printf("\n");
    }
 
    for (i = 0; i < N; i++)
        free(distA[i]);
 
    free(distA);
 
    for (i = 0; i < M; i++)
        free(distB[i]);
 
    free(distB);
 
    deletePlane(startA); // Очищаем память, выделенную для списка точек первой плоскости
    deletePlane(startB); // Очищаем память, выделенную для списка точек второй плоскости
 
    getch();
    return 0;
}
 
// Функция, создающая список точек, считываемых из файла
// Принимаемые параметры:
// filename - имя файла, содержащего координаты точек
// count - указатель на счётчик количества точек (в результате работы функции количество точек будет записано в эту переменную)
// Возвращаемое значение:
// Указатель на начало списка точек плоскости
struct plane *makePlane(char *filename, int *count)
{
    struct plane *start = NULL; // Начальный элемент списка точек плоскости
    struct plane *current;      // Текущий элемент списка точек плоскости
    struct plane *temp;         // Временный элемент
 
    FILE *fin; // Указатель на входной файл
 
    if ((fin = fopen(filename, "r")) == NULL) // Октрываем входной файл, а если его не удалось открыть,
    {
        fprintf(stderr, "Error opening file \"%s\"!", filename); // сообщаем об ошибке
        getch();
        exit(1);                                                 // и выходим
    }
 
    *count = 0;
 
    while(!feof(fin)) // Продолжаем считывание, пока не достигнем конца файла
    {
        temp = (struct plane *)malloc(sizeof(struct plane)); // Выделяем память под очередной элемент списка
        fscanf(fin, "%lf %lf", &(temp->x), &(temp->y));      // и считываем в него координаты очередной точки
        temp->next = NULL;                                   // Устанавливаем указатель на следующий элемент в NULL
 
        if (start == NULL) // Если список ещё не был создан
        {
            start = temp;    // Созданный элемент делаем начальным
            current = start; // Начальный элемент делаем текущим
        }
        else // Если продолжаем создание списка, начатое ранее
        {
            current->next = temp; // В текущем элементе указатель на следующий элемент устанавливаем на созданный
            current = temp;       // Созданный элемент делаем текущим
        }
 
        (*count)++; // Увеличиваем счётчик точек плоскости
    }
 
    fclose(fin);
 
    return start;
}
 
// Функция, очищающая память, выделенную для списка точек
// Принимаемые параметры:
// *start - указатель на начало списка точек
void deletePlane(struct plane *start)
{
    struct plane *temp; // Временный элемент
 
    // в цикле пробегаем все элементы списка
    while (start != NULL)
    {
        temp = start;        // Устанавливаем указатель temp на элемент start
        start = start->next; // Устанавливаем указатель start на следующий элемент списка
        free(temp);          // Освобождаем память, на которую указывает temp
    }
}
 
// функция, рассчитывающая расстояния между всеми точками плоскости
// Принимаемые параметры:
// **matrix - указатель на матрицу расстояний (в результате матрица запишется в эту переменную)
// dim - количество точек плоскости
// *start - указатель на список точек
void distance(double **matrix, int dim, struct plane *start)
{
    struct plane *temp; // Временная переменная
    double x, y;        // Вершины очередной точки, от которой считается расстояние до остальных точек
    int i, j, k;
 
    // В цикле забиваем данными все ячейки матрицы
    for (i = 0; i < dim; i++)
    {
        temp = start; // Устанавливаем временный указатель на начало списка точек
        k = 0;
 
        // Ищем точку, от которой будем рассчитывать расстояние до остальных
        while (k < i)
        {
            temp = temp->next;
            k++;
        }
 
        x = temp->x; // Запоминаем вершины
        y = temp->y; // найденной точки
        temp = temp->next;
 
        for (j = i; j < dim - 1; j++)
        {
            matrix[i][j] = sqrt(pow(x - temp->x, 2.0) + pow(y - temp->y, 2.0)); // Рассчитываем расстояние от найденной точки до очередной
            matrix[j + 1][i] = matrix[i][j];                                    // Поскольку расстояния дублируются, сохраняем найденной расстояние в дублируемую ячейку
            temp = temp->next;                                                  // Продвигаемся дальше по списку точек
        }
    }
}

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
 
typedef struct {
    double x;
    double y;
} point_t;
 
double distance(point_t *a, point_t *b) {
    return sqrt(pow(a->x - b->x, 2) + pow(a->y - b->y, 2));
}
 
#define NPOINTS 5
 
int main() 
{
    int i = 0;
    int j = 0;
    int k = 0;
    point_t points[NPOINTS];
    double distances[NPOINTS][NPOINTS];
    
    for(i = 0; i < NPOINTS; ++i) {
        points[i].x = (double)rand() / rand();
        points[i].y = (double)rand() / rand();
    }
    
    for(i = 0; i < NPOINTS; ++i) {
        printf("%d: [%f, %f]\n", i, points[i].x, points[i].y);
    }
    printf("\n");
    
    for(i = 0; i < NPOINTS; ++i) {
        for(j = 0; j < NPOINTS; ++j) {
            distances[i][j] = distance(&points[i], &points[j]);
        }
    }
    
    for(i = 0; i < NPOINTS; ++i) {
        for(j = 0; j < NPOINTS; ++j) {
            printf("%f ", distances[i][j]);
        }
        printf("\n");
    }
    
    return 0;
}