Численные методы. Решение систем линейных уравнений

Type

TMatr = array[1..100,1..100] of real;

TVect = array[1..100] of real;

Var

A :TMatr; x,f:TVect;

i,j, n :Integer; S:Real;

Procedure LoadFromFile(name:string; var n:Integer; var A:TMatr; var f:TVect);

var

fn:Text; i,j:Integer;

begin

Assign(fn, name);

ReSet(fn);

Read(fn, n);

for i:=1 to n do

for j:=1 to n do

Read(fn, A[i,j]);

for i:=1 to n do Read(fn, f[i]);

Close(fn);

end;

Procedure Print_Matr(A:TMatr; n:Integer);

var i,j:Integer;

begin

writeln('Matrix=');

for i:=1 to n do

begin

for j:=1 to n do

Write(A[i,j], ' ');

WriteLn;

end;

end;

Procedure Print_Vect(x:TVect; n:Integer);

var i:Integer;

begin

writeln('Vector=');

for i:=1 to n do Write(x[i], ' ');

end;

procedure LU(A:TMatr; n:Integer; var L,U:TMatr);

var

i,j,k:Integer; Sum:real;

begin

for i:=1 to n do

for j:=1 to n do

if i>=j then

begin

Sum:=0;

for k:=1 to j-1 do Sum:=Sum+L[i,k]*U[k,j];

L[i,j]:= A[i,j]-Sum;

end

else