用c++写matrix矩阵的运算，矩阵为三行四列并存放double类型的数据

以下是一个用 C 语言编写的矩阵运算示例代码，可以处理三行四列的 double 类型数据：

#include <stdio.h>

#define ROWS 3
#define COLS 4

void matrix_multiply(double a[][COLS], double b[][COLS], double c[][COLS]) {
    int i, j, k;
    for (i = 0; i < ROWS; i++) {
        for (j = 0; j < COLS; j++) {
            c[i][j] = 0;
            for (k = 0; k < COLS; k++) {
                c[i][j] += a[i][k] * b[k][j];
            }
        }
    }
}

int main() {
    double a[ROWS][COLS] = {{1, 2, 3, 4},
                            {5, 6, 7, 8},
                            {9, 10, 11, 12}};
    double b[ROWS][COLS] = {{2, 0, 1, 3},
                            {1, 2, 1, 0},
                            {0, 1, 2, 1}};
    double c[ROWS][COLS];
    int i, j;

    matrix_multiply(a, b, c);

    for (i = 0; i < ROWS; i++) {
        for (j = 0; j < COLS; j++) {
            printf("%f ", c[i][j]);
        }
        printf("\n");
    }

    return 0;
}

这个程序定义了一个 `matrix_multiply` 函数，用于计算两个矩阵的乘积，并将结果存储在第三个矩阵中。在 `main` 函数中，我们定义了两个三行四列的矩阵 `a` 和 `b`，并将它们传递给 `matrix_multiply` 函数进行计算。最后，我们打印出结果矩阵 `c` 的值。