Где-то закралась ошибка

//библиотеки
namespace modul_for_cognitive_maps
{
    public partial class Form1 : Form
    {
        public Form1()
        {
            InitializeComponent();
        }
        private void Form1_Load(object sender, EventArgs e)
        {
            int rows = 7;
        List<List<double>> matrix = new List<List<double>>();
            //проинициализируем всю матрицу нулями
            zero(matrix,rows);
            matrix[0][2] = 0.6;
            matrix[0][3] = 0.9;
            matrix[1][0] = 1.0;
            matrix[2][1] = 0.7;
            matrix[2][4] = 0.9;
            matrix[3][6] = 0.9;
            matrix[4][5] = -0.9;
            matrix[4][6] = -0.9;
            matrix[5][0] = -0.3;
            matrix[6][5] = 0.8;
            // Шаг первый. Проверяем типа алгоритма в зависимости от знаков элементов
            
                otric_alg(matrix, rows);
        }
 
        void otric_alg(List<List<double>> matrix, int rows)
        {
            List<List<double>> double_table = new List<List<double>>();
            zero(double_table, rows*2);
            //составляем таблицу
            double_table =  create_double_table(matrix, double_table,rows);
//получили матрицу по специальному алгоритму, что  не так важно
            //создаём матрицу для умножения
            List<List<double>> res_matrix = new List<List<double>>();
            List<List<double>> tranz_zamkn_matrix = new List<List<double>>();//транзитивно-замкнутая матрица
            res_matrix = double_table;
            tranz_zamkn_matrix = res_matrix;
            //умножаем всё с помощью обычного транзитивного замыкания
 
            for (int i = 0; i < rows * 2; i++)
            {
                res_matrix = multi_matrix(double_table, res_matrix, rows * 2);//умножение 
                tranz_zamkn_matrix = maksimil_matrix(res_matrix, tranz_zamkn_matrix,rows*2);//сравнение
            }
 
        }
        
 
        void zero(List<List<double>> matrix, int rows)
        {
            for (int i = 0; i < rows; i++)
            {
                List<double> row = new List<double>();
                for (int j = 0; j < rows; j++)
                {
                    row.Add(0);
                }
 
                matrix.Add(row);
            }
        }
 
    
        List<List<double>> maksimil_matrix(List<List<double>> matrix, List<List<double>> res_matrix, int columns)
        {
 
            for (int i = 0; i < columns; i++)//цикл по строкам первой матрицы
            {
                for (int j = 0; j < columns; j++)//цикл по столбцам второй матрицы
                {
                    if (matrix[i][j] > res_matrix[i][j]) res_matrix[i][j] = matrix[i][j];
                }
            }
            return res_matrix;
        }
 
        List<List<double>> multi_matrix(List<List<double>> matrix, List<List<double>> res_matrix, int columns)
        {
            List<List<double>> ret_matrix = new List<List<double>>();//используем в процессе вычисления
            zero(ret_matrix,columns);
            double res_umn = 0;
            for (int i = 0; i < columns; i++)//цикл по строкам первой матрицы
            {
                for (int j = 0; j < columns; j++)//цикл по столбцам второй матрицы
                {
                    for (int l = 0; l < columns; l++)//цикл по элементам столбцов второй матрицы
                    {
                        res_umn += matrix[i][l] * res_matrix[l][j];
                    }
                    ret_matrix[i][j] = res_umn;//перемножили строчку на j-ый столбец
                    res_umn = 0;
                }
            }
            return ret_matrix;
        }
 
    }
}