如何在c中查找e^矩阵

#include <stdio.h>
#include <math.h>

#define ROWS 3
#define COLS 3

void matrixExponential(double matrix[ROWS][COLS], double result[ROWS][COLS]) {
    int i, j;
    double temp[ROWS][COLS];

    // 初始化结果矩阵为单位矩阵
    for (i = 0; i < ROWS; i++) {
        for (j = 0; j < COLS; j++) {
            if (i == j) {
                result[i][j] = 1.0;
            } else {
                result[i][j] = 0.0;
            }
        }
    }

    // 计算矩阵的指数
    for (i = 1; i <= 10; i++) {
        double factorial = 1.0;
        for (j = 1; j <= i; j++) {
            factorial *= j;
        }
        // 计算矩阵的幂
        for (int row = 0; row < ROWS; row++) {
            for (int col = 0; col < COLS; col++) {
                temp[row][col] = matrix[row][col] * pow(i, col) / factorial;
            }
        }
        // 累加结果
        for (int row = 0; row < ROWS; row++) {
            for (int col = 0; col < COLS; col++) {
                result[row][col] += temp[row][col];
            }
        }
    }
}

int main() {
    double matrix[ROWS][COLS] = {
        {1, 2, 3},
        {4, 5, 6},
        {7, 8, 9}
    };
    double result[ROWS][COLS];

    matrixExponential(matrix, result);

    printf("e^矩阵：\n");
    for (int i = 0; i < ROWS; i++) {
        for (int j = 0; j < COLS; j++) {
            printf("%.2f ", result[i][j]);
        }
        printf("\n");
    }

    return 0;
}