Сформировать массив B, каждый элемент которого равен произведению элементов столбцов матрицы A

#include <iostream>
#include <ctime>
#include <iomanip>
 
#define M 5
#define N 5
 
int main()
{
    int a[M][N];
    int b[N];
 
    /* Заполняем матрицу A случайными числами
     *  и выводим на экран */
    std::cout << "A:" << std::endl;
    srand(time(0));
    for (size_t i = 0; i < M; i++)
    {
        for (size_t j = 0; j < N; j++)
        {
            a[i][j] = rand() % 10;
            std::cout << std::setw(5) << a[i][j];
        }
        std::cout << std::endl;
    }
 
    /* Обходим матрицу А по столбцам...
     * и заполняем массив B ...*/
    int mul = 1;
    for (size_t i = 0; i < N; i++)
    {
        for (size_t j = 0; j < M; j++)
        {
            mul *= a[j][i];
            //mul - произведение элементов текущего столбца матрицы A
        }
 
        b[i] = mul;
        mul = 1;
    }
 
    /* выводим на экран массив B */
    std::cout << std::endl << "B:" << std::endl;
    for (size_t i = 0; i < N; i++)
    {
        std::cout << b[i] << " ";
    }
    std::cout << std::endl << std::endl;
 
    system("pause");
 
return 0;
}