c语言：矩阵相乘

#include <stdio.h>

int main()
{
    int a[4][3] = 
    {   {1,2,3},
        {4,5,6},
        {7,8,9},
        {1,2,3}};
    int b[3][3] =
    { 
        { 1,2,3 },
        { 4,5,6 },
        {7,8,9} 
    };
    //c矩阵为a矩阵与b矩阵的乘积
    int c[4][3];
    //计算两个矩阵相乘
    matrix_multiply(4,3,3,&a[0],&b[0],&c[0]);
    //输出c矩阵
    for (int i = 0; i < 4; i++)
    {
        for(int j=0;j<3;j++)
            printf("%4d  ", c[i][j]);
        printf("\n");
    }
    return 0;
}

//矩阵相乘算法，最后四个参数是两个相乘的矩阵的行数和列数
void matrix_multiply(const int m, const int n, const int p, const int *a, const int *b, int *c){
    for(int i=0;i<m;i++){
        for(int j=0;j<p;j++){
            *(c+i*3+j)=0;
            for(int k=0;k<p;k++){
                *(c+i*3+j)+=*(a+i*3+k)**(b+k*3+j);
            }
        }
    }

}