СЛАУ метод LU-разложения Краута

program LR1_V3;
   const n=4;
   type matrica=array[1..10,1..10] of real;
        vector=array[1..4] of real;     
   var i,j,k,m:integer;
       a:matrica;
       x:vector;
       ta,tb:real;
begin
  //Формирование и вывод расширенной матрицы системы
  a[1,1]:=-0.2;  a[1,2]:=1.4;     a[1,3]:=2.1;   a[1,4]:=-0.5;    a[1,5]:=2.180417;
  a[2,1]:=0.3;   a[2,2]:=1.5;       a[2,3]:=4.1;   a[2,4]:=0.7;   a[2,5]:=1.883982;
  a[3,1]:=0.6;  a[3,2]:=1.4;     a[3,3]:=0.4;   a[3,4]:=-1.8;    a[3,5]:=0;
  a[4,1]:=1.7;   a[4,2]:=2.7;     a[4,3]:=-1.1;  a[4,4]:=1.2;   a[4,5]:=2.222455;
    writeln('     Исходная матрица'); 
    writeln('==================================');
    for i:= 1 to n do
       begin
         for j:= 1 to n do
           write  (a[i,j]:6:1);
       writeln (' || ',a[i,j+1]:5:4);
       end;