Написать три алгоритма решения СЛАУ: Метод прогонки, метод квадратных корней, метод вращений

//прогонки
                //коэфф
                p[0] = -matrix[0, 1] / matrix[0, 0];
                q[0] = d[0] / matrix[0, 0];
                for (int i = 1; i < columns-1; i++)
                {
                    p[i] = -matrix[i, i + 1] / (matrix[i, i] + matrix[i, i - 1] * p[i - 1]);
                    q[i] = (d[i] - matrix[i, i - 1] * q[i - 1]) / (matrix[i, i] + matrix[i, i - 1] * p[i - 1]);
                }
                p[columns-1] = 0;
                q[columns - 1] = (d[columns - 1] - matrix[columns - 1, columns - 2] * q[columns - 2]) / (matrix[columns - 1, columns - 1] + matrix[columns - 1, columns - 2] * p[columns-2]);
                x[columns - 1] = q[columns-1];
                for (int i = columns-2; i >= 0; i--)
                {
                    x[i] = p[i] * x[i + 1] + q[i];
                }
                //NaN и бесконечность
                for (int i = 0; i < columns; i++)
                {
                    if (x[i] != x[i])
                    {
                        MessageBox.Show("Система не имеет решения! Было получено NaN - Not a Number!");
                        return;
                    }
                    if (x[i] == System.Double.NegativeInfinity || x[i] == System.Double.PositiveInfinity)
                    {
                        MessageBox.Show("Система не имеет решения! Было получено бесконечное число!");
                        return;
                    }
                }
                //
                //заполняем x в таблицу С
                for (int i = 0; i < columns; i++)
                    dataGridView3[0, i].Value = x[i];
                //невязка для прогонки
                dataGridView4[0, 0].Value = matrix[0, 0] * x[0] + matrix[0, 1] * x[1] - d[0];
                for (int i = 1; i < columns-1; i++)
                {
                    dataGridView4[0,i].Value = matrix[i, i - 1] * x[i - 1] + matrix[i, i] * x[i] + matrix[i, i + 1] * x[i + 1] - d[i];
                }
                dataGridView4[0, columns - 1].Value = matrix[columns - 1, columns - 2] * x[columns - 2] + matrix[columns - 1, columns - 1] * x[columns - 1] - d[columns-1];
                //

//метод квадратных корней
//находим L верхнетреугольную
                L[0, 0] = Math.Sqrt(matrix[0,0]);
                
                for (int i = 0; i < columns; ++i)
                    for (int j = 0; j < columns; ++j)
                    {
                        L[0, j] = matrix[0, j] / L[0, 0];
 
                        //если не отриц число
                        if (i > 0)
                        {
                            double buf_l = 0;
                            for (int l = 0; l < i - 1; ++l)
                                buf_l = L[l, i] * L[l, i];
 
                            if (matrix[i, i] - buf_l >= 0.0d)
                                L[i, i] = Math.Sqrt(matrix[i, i] - buf_l);
                            else
                            {
                                MessageBox.Show("Система не имеет решений, так как матрица не является положительно определенной!");
                                return;
                            }
                        }
 
 
 
                        if (j > i)
                        {
                            double buf_l2 = 0;
                            for (int l = 0; l < i - 1; ++l)
                                buf_l2 = L[l, i] * L[l, j];
                            L[i, j] = (matrix[i, j] - buf_l2) / L[i, i];
                        }
 
                        if (i > j)
                            L[i, j] = 0;
                    }
                //теперь получим транспонированную Lt (нижнетреугольная)
                double[,] Lt = trans(L);
 
                //Lt*y = b
                double[] y = triangle(Lt, d);
                //L*x = y
                double[] res = triangle(L, y);
 
                //
                //заполняем x в таблицу С
                for (int i = 0; i < columns; i++)
                    dataGridView3[0, i].Value = res[i];
                //невязка
 
                for (int i = 0; i < columns; i++)
                {
                    dataGridView4[0, i].Value = 0;
                    for (int j = 0; j < columns; j++)
                        dataGridView4[0, i].Value = double.Parse(dataGridView4[0, i].Value.ToString()) + matrix[i, j] * res[j];
                    dataGridView4[0, i].Value = double.Parse(dataGridView4[0, i].Value.ToString()) - d[i];
 
                }
        double[,] trans(double[,] mat)
        {
            int columns = trackBar1.Value;
            double[,] buf = new double[columns, columns];
            for (int i = 0; i < columns; ++i)
                for (int j = 0; j < columns; ++j)
                    buf[i, j] = mat[j, i];
 
            return buf;
        }
 
        double[] triangle(double[,] mat, double[] d)
        {
            int columns = trackBar1.Value;
            double[] res = new double[columns]; // решение
            double sum;
            int j;
 
            res[columns - 1] = d[columns - 1] / mat[columns - 1, columns - 1];
            sum = 0;
            for (int z = columns - 2; z >= 0; --z)
            {
                sum = 0;
                for (j = z + 1; j < columns; j++)
                {
                    sum = sum + mat[z, j] * res[j];
                    res[z] = (d[z] - sum) / mat[z, z];
                }
            }
            return res;
        }