Решение СЛАУ методом ортогонализации

const max=10;
 
type TVector=array [1..max] of real;
     TMatrix=array [1..max+1] of TVector;
 
var i,j,n:integer;
    alpha:real;
    x,m:TMatrix;
    res:TVector;
 
function Scalyar(a,b:TVector):real;
var i:integer;
    result:real;
begin
 result:=0;
 for i:=1 to n do
  result:=result+a[i]*b[i];
 Scalyar:=result;
end;
 
procedure VectorNum(a:real; k:TVector ;var v:TVector);
var i:integer;
begin
 for i:= 1 to n do
  v[i]:=a*k[i];
end;
 
procedure MatrixVector(var tmp,a:TVector);
var i,j:integer;
    result:real;
begin
 for i:=1 to n do
 begin
  result:=0;
  for j:=1 to n do
   result:=result+a[j]*m[j][i];
  tmp[i]:=result;
 end;
end;
 
procedure AddVector(b:TVector; var a:TVector);
var       i:integer;
begin
 for i:=1 to n do
  a[i]:=a[i]+b[i];
end;
 
procedure ortogonalization;
var t:TVector;
    i,i2:integer ;
Begin
 for i:=1 to n do
  x[n+1][i]:=0;
 
 for i:=1 to n do
 begin
  x[i]:=x[n+1];
  x[i][i]:=1
 end;
 
 for i:=2 to n do
 begin
  for j:=1 to i-1 do
  begin
 
   MatrixVector(t,x[j]);
   alpha:=-(Scalyar(m[i],x[j])/Scalyar(t,x[j]));
   VectorNum(alpha,x[j],t);
   AddVector(t,x[i]);
  end;
 end;
 
 for i:=1 to n do
 begin
  MatrixVector(t,x[i]);
  alpha:=Scalyar(m[n+1],x[i])/Scalyar(t,x[i]);
  VectorNum(alpha,x[i],t);
  AddVector(t,Res);
 end;
end;
 
var  f,result:text;
     s:string;
 
begin
 write('name of file : ');
 readln(s);{s:='3x3_2';}
 
 assign(f,'C:\Users\Admin\Desktop\'+s+'.txt');
 
 reset(f);
 
 readln(f,n);
 for i:=1 to n do
 begin
  for j:=1 to n+1 do
   read(f,m[j][i]);
   readln(f);
 end;
 
 close(f);
 
 Ortogonalization;
 assign(f,'C:\Users\Admin\Desktop\4.txt');
 rewrite(f);
 for i:=1 to n do
  writeln(f,'x',i,'=',res[i]);
 
 close(f);
end.