#include "matrix.hpp"

int main() {
    matrix<double> m(4, 1);
    m(1, 0) = 1;
    m(2, 0) = 1;
    m(3, 0) = 1;
    matrix<int> r(1, 3);
    r(0, 0) = 1;
    r(0, 1) = 1;
    r(0, 2) = 1;
    std::cout << kronecker_product(m, r);
    return 0;
}