Решение СЛАУ методом Холецкого

Program holeckogo;
Uses crt;
const n=3;
Var
   i,j:byte;
   A,L,Lt: array[1..n,1..n] of integer;
   B,X,y: array[1..n] of integer;
{Блок ввода-вывода}
Procedure inmassA;
begin
 WriteLn('Введите элементы массива A');
  For i:=1 to N do
     For j:=1 to N do
       Begin
           Write('a[',i,',',j,']=');
           Read(A[i,j]);
       end;
end;
 Procedure inmassB;
 begin
Writeln('Введите элементы массива B');
  For i:=1 to n do
      Begin
           Write('b[',i,']=');
           ReadLn(b[i]);
       End;
  end;
Procedure Outx;
Begin
  Writeln('Результат:');
  For i:=1 to n do
      Begin
           Write('x',i,'=');
           Writeln(X[i]);
       End;
End;
{Конец блока ввода-вывода}
Begin
 ClrScr;
  inmassA;
  inmassB;
  if a[1,1]<0 then
  Begin
       for j:=1 to n do
           a[1,j]:=-a[1,j];
       b[1]:=-b[1];
  end;
 {Вычисления}
  L[1,1]:=sqrt(A[1,1]);
  Lt[1,1]:=sqrt(A[1,1]);
  L[2,1]:=A[1,2]/sqrt(A[1,1]);
  Lt[1,2]:=A[1,2]/sqrt(A[1,1]);
  L[3,1]:=A[1,3]/sqrt(A[1,1]);
  Lt[1,3]:=A[1,3]/sqrt(A[1,1]);
  L[2,2]:=sqrt(A[2,2]-L[2,1]*L[2,1]);
  Lt[2,2]:=sqrt(A[2,2]-L[2,1]*L[2,1]);
  L[3,2]:=(A[3,2]-L[2,1]*L[3,1])/L[2,2];
  Lt[2,3]:=(A[3,2]-L[2,1]*L[3,1])/L[2,2];
  L[3,3]:=sqrt(A[3,3]-L[3,1]*L[3,1]-L[3,2]*L[3,2]);
  Lt[3,3]:=sqrt(A[3,3]-L[3,1]*L[3,1]-L[3,2]*L[3,2]);
 
  y[1]:=B[1]/L[1,1];
  y[2]:=(B[2]-L[2,1]*y[1])/L[2,2];
  y[3]:=(B[3]-L[3,1]*y[1]-L[3,2]*y[2])/L[3,3];
 
  x[3]:=y[3]/Lt[3,3];
  x[2]:=(y[2]-Lt[2,3]*x[3])/Lt[2,2];
  x[1]:=(y[1]-Lt[1,2]*x[2]-Lt[1,3]*x[3])/Lt[1,1];
{Конец вычислений}
Outx;
 
End.

L[1,1]:=sqrt(A[1,1]);

Var
   i,j:byte;
   A: array[1..n,1..n] of integer;
   L,Lt : array[1..n,1..n] of real;
   B: array[1..n] of integer;
   y,x : array[1..n] of real;

  writeln(' ',A[3,3]-L[3,1]*L[3,1]-L[3,2]*L[3,2]:8:4);
  L[3,3]:=sqrt(A[3,3]-L[3,1]*L[3,1]-L[3,2]*L[3,2]);

A,L,Lt: array[1..n,1..n] of real;
   B,X,y: array[1..n] of real;