Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
15
Добавлен:
17.04.2013
Размер:
1 Кб
Скачать


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



 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.









Соседние файлы в папке curstau