Найти произведение матриц

#include <iostream>
#include <vector>
 
using Row = std::vector<int>;   
using Matrix = std::vector<Row>;
 
Matrix rand_matrix(int r, int c) {
    Matrix m;
    
    for (int i = 0; i < r; i++) {
        m.push_back(Row());
        for (int j = 0; j < c; j++) {
            m[i].push_back(rand() % 10);
        }
    }
    
    return m;
}
 
void print(const Matrix& m, const char* s1 = "", const char *s2 = "\n") {
    std::cout << s1;
    for (const auto& row : m) {
        for (const auto& e : row) {
            std::cout << e << " ";
        }
        std::cout << "\n";
    }
    std::cout << s2;
}
 
Matrix mul(const Matrix& a, const Matrix& b) {
    Matrix c(a.size());
    
    for (int i = 0; i < a.size(); i++) {
        c.push_back(Row());
        for (int j = 0; j < b[0].size(); j++) {
            c[i].push_back(0);
            for (int k = 0; k < a[0].size(); k++) {
                c[i][j] += a[i][k] * b[k][j];
            }
        }
    }
    
    return c;
}
 
int main() {
    int k = 2, m = 3, l = 4;
    
    auto a = rand_matrix(k, m);
    auto b = rand_matrix(m, l);
    auto c = mul(a, b);
    
    print(a, "----- A -----\n");
    print(b, "----- B -----\n");
    print(c, "----- C -----\n");
 
    return 0;
}