Напишите программу, которая решает матричное уравнение вида С*Х=В где С и В– матрицы 3-го порядка

#include <iostream>
#include <fstream>
#include <Eigen/Dense>
 
int main()
{
    Eigen::MatrixXd C(3, 3); // матрица С
    Eigen::MatrixXd B(3, 1); // матрица В
    Eigen::MatrixXd X(3, 1); // матрица Х
 
    // Считываем матрицу С из файла
    std::ifstream matrixC_file("matrix_C.txt");
    if (matrixC_file.is_open())
    {
        for (int i = 0; i < 3; i++)
        {
            for (int j = 0; j < 3; j++)
            {
                matrixC_file >> C(i, j);
            }
        }
        matrixC_file.close();
    }
    else
    {
        std::cerr << "Unable to open matrix_C.txt\n";
        return 1;
    }
 
    // Считываем матрицу В из файла
    std::ifstream matrixB_file("matrix_B.txt");
    if (matrixB_file.is_open())
    {
        for (int i = 0; i < 3; i++)
        {
            matrixB_file >> B(i, 0);
        }
        matrixB_file.close();
    }
    else
    {
        std::cerr << "Unable to open matrix_B.txt\n";
        return 1;
    }
 
    // Решаем матричное уравнение С * Х = В
    X = C.fullPivLu().solve(B);
 
    // Выводим результат
    std::cout << "Решение матричного уравнения С * Х = В:\n";
    std::cout << "Х =\n" << X << std::endl;
 
    return 0;
}