Как реализовать умножение двух матриц?

#include <iostream>

using namespace std;
int main()
{
    
    double A[3][4];
    double B[4][2];
    double C[4][2];
    A[0][0] = 5;
    A[0][1] = 2;
    A[0][2] = 0;
    A[0][3] = 10;
    A[1][0] = 3;
    A[1][1] = 5;
    A[1][2] = 2;
    A[1][3] = 5;
    A[2][0] = 20;
    A[2][1] = 0;
    A[2][2] = 0;
    A[2][3] = 0;
    B[0][0] =1.2;
    B[0][1] = 0.5;
    B[1][0] = 2.8;
    B[1][1] = 0.4;
    B[2][0] = 5;
    B[2][1] = 1;
    B[3][0] = 2;
    B[3][1] = 1.5;
    
    cout << "Matrix A:" << endl;
    for (int i = 0; i < 3; i++) {
        for (int j = 0; j < 4; j++) {
            cout << A[i][j];
            cout << "  ";
        }
        cout << endl;
    }
    cout << endl;
    cout << endl;
    cout << "Matrix B:" << endl;
    for (int i = 0; i < 4; i++) {
        for (int j = 0; j < 2; j++) {
            cout << B[i][j];
            cout << " ";
        }
        cout << endl;
    }
    cout << endl;
    cout << endl;
    cout << "Matrix C=A*B:" << endl;
    
 
}