Курсовые / curstau / FINIT

{Їа®жҐ¤гал ўлзЁб«Ґ­Ёп ¬ ваЁж ॣг«пв®а  б®бв®п­Ёп}
{ Ё ­ Ў«о¤ вҐ«п (б Є®­Ґз­л¬ ўаҐ¬Ґ­Ґ¬) }



unit FINIT; interface
uses GB,ALG,CANON,ENEX;
procedure REGT (n:integer;A,B:matr;var Rx:matr);
procedure OBST (n:integer;A,C:matr;var Hx:matr);

implementation


{ॣг«пв®а б®бв®п­Ёп б Є®­Ґз­л¬ ўаҐ¬Ґ­Ґ¬}

procedure REGT (n:integer;A,B:matr;var Rx:matr);
var j,k :integer;
Sco,ISco,R :matr;
begin
CONT(n,A,B,Sco);
INV (n,Sco,ISco);
for j:=1 to n do R[1,j]:=ISco[n,j];
k:=1;
repeat
MUL(1,n,n,R,A,Rx);
for j:=1 to n do R[1,j]:=Rx[1,j];
k:=k+1;
until k=n+1;
end;
{-------------------------------------------}


{­ Ў«о¤ вҐ«м б®бв®п­Ёп б Є®­Ґз­л¬ ўаҐ¬Ґ­Ґ¬}

procedure OBST (n:integer;A,C:matr;var Hx:matr);
var i,k :integer;
Sob,ISob,H :matr;
begin
OBS(n,A,C,Sob);
INV (n,Sob,ISob);
for i:=1 to n do H[i,1]:=ISob[i,n];
k:=1;
repeat
MUL(n,n,1,A,H,Hx);
for i:=1 to n do H[i,1]:=Hx[i,1];
k:=k+1;
until k=n+1;
end;
{-------------------------------------------------------------}

end.