// 假设我们有两个矩阵A和B,它们的尺寸分别为 m x n 和 n x p,矩阵C是结果矩阵,尺寸为 m x p。
void matrix_multiply(int m, int n, int p, double A[m][n], double B[n][p], double C[m][p]) {
for (int i = 0; i < m; i++) {
for (int j = 0; j < p; j++) {
C[i][j] = 0;
for (int k = 0; k < n; k++) {
C[i][j] += A[i][k] * B[k][j];
}
}
}
}