using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;

namespace LU_Разложение
{//Я понял. Те числа, это порядок обхода. Может быть, я это и сделаю.

    class Program
    {
      public static void Print(int[,] f)
        {
            for (int p = 0; p < 4; p++)
            {
                for (int g = 0; g < 4; g++)
                {
                    Console.Write(Convert.ToString(f[p, g]) + " ");

                }
                Console.WriteLine();
            }
        }
        static void Main(string[] args)
        {
            int[,] a = new int[,] { {1,2,3,4 },{5,6,7,8 },{9,1,2,3}, { 4, 5, 6, 7} };
            Print(a);
            int[,] l = new int[4, 4];
            int[,] u = new int[4, 4];
            Console.WriteLine();
            int k,n=4,i=0,j=0;
            for (k=n-1; k>=0;k--)
            {
                for(i=0;i<=k;i++)
                {
                    u[i, k] = a[i, k];
                   // Console.Write(a[i, k]);
                    for(j=k+1;j<n;j++)
                    {
                        u[i, k] -= u[i, j] * l[j, k];
                    }
                }
               // Console.WriteLine();
                for (i = n-1; i > k; i--)
                {
                    l[k, i] = a[i, k];
                    Console.Write(l[k, i]);
                    for (j = k + 1; j < n; j++)
                    {
                        l[k, i] -= u[k, j] * l[j, i];
                        l[k, i] /= u[k, k];
                    }

                }
                Console.WriteLine();
            }


        }
    }
}