Метод Гаусса Опции

[moto]

Помогите плиз реализовать метод Гаусса  в Turbo pascal как можно более простым способом. Алгоритм решения ясен, вот только возникла проблема в приведении расширенной матрицы к диагональному виду, голову сломала  ??? помогите плиз!!!

[Zole]

Текст программы:

Код

Program Gauss;
Uses crt;
Const
     n=3;
Var i,j,k:integer;
    c,s:real;
    A:array[1..n,1..n] of real;
    B:array[1..n] of real;
    X:array[1..n] of real;
Begin   clrscr;
  Writeln(' Введите коэфициенты систем уравнений');
        For i:=1 to n do
        For j:=1 to n do read(a[i,j]);
   Writeln(' Введите коэфициенты свободных членов уравнений');
        For i:=1 to n do Read(b[i]);
               k:=1;
Repeat
       i:=k;
    Repeat
         c:=a[i+1,k]/a[k,k];
          j:=1;
       Repeat
             a[i+1,j]:=a[i+1,j]-c*a[k,j];
             j:=j+1
      Until    j>=n;
                  b[i+1]:=b[i+1]-c*b[k];
                  i:=i+1
    Until   i>n-1;
               k:=k+1
Until    k>=(n-1);
           x[n]:=b[n]/a[n,n];
           i:=n-1;
  Repeat
        s:=0;
        j:=i+1;
     Repeat
           s:=s+x[j]*a[i,j];
           j:=j+1
        Until    j>n;
                   x[i]:=(b[i]-s)/a[i,i];
                    i:=i-1
 Until i<0;
   For i:=1 to n do Writeln('x[',i,']=',x[i]:9:7)
End.

;) :p9: :p9:

Сообщение отредактировано: volvo - 5.01.2005 1:00