Массивы => Вычислить отклонение элементов от среднего значения

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
 
#define N 5
#define M 6
 
void fill(int matrix[N][M], size_t rows, size_t cols)
{
    size_t i, j;
    for (i = 0; i < rows; i++)
        for (j = 0; j < cols; j++)
            matrix[i][j] = rand() % 20;
}
 
void print(const int matrix[N][M], size_t rows, size_t cols)
{
    size_t i, j;
    for (i = 0; i < rows; i++)
    {
        for (j = 0; j < cols; j++)
        {
            printf("%4d", matrix[i][j]);
        }
        putchar('\n');
    }
}
 
int column_sum(const int matrix[N][M], size_t rows, size_t column)
{
    int sum = 0, i = 0;
    for (i = 0; i < rows; i++)
        sum += matrix[i][column];
    return sum;
}
 
void print_difference(const int matrix[N][M], size_t rows, size_t cols)
{
    size_t i, j;
    for (i = 0; i < cols; i++)
    {
        double avg = (double)column_sum(matrix, rows, i) / rows;
        for (j = 0; j < rows; j++)
        {
            double difference = avg - matrix[j][i];
            printf("matrix[%d][%d] = %d, difference = %lf\n",
                j, i, matrix[j][i], difference);
        }
        putchar('\n');
    }
}
 
int main()
{
    int matrix[N][M];
    srand(time(0));
    fill(matrix, N, M);
    print(matrix, N, M);
    putchar('\n');
    print_difference(matrix, N, M);
    return EXIT_SUCCESS;
}