Метод гаусса.

[ironhide8]

Очень нужна ваша помощь,зачет проваливается((Нужно реализовать метод гаусса,вот он-

program Gayss;
const
  n = 3;

type
  matrix = array[1..n, 1..n] of real;
  vector = array[1..n] of real;

 procedure LoadFromFile(var M: matrix; var b: vector);
var
 x, y: Integer;
 val: Real;
 input: text;
begin
 Assign(input, 'matrix.txt');
 Reset(input);
 for y := 1 to n do
  for x := 1 to n do begin
   Readln(input, val);
   M[y, x] := val;
  end;
 for x := 1 to n do begin
   Readln(input, val);
   b[x] := val;
  end;
 Close(input);
end;

procedure Print(const m: matrix);
var
 x, y: Integer;
begin
 for y := 1 to n do begin
  for x := 1 to n do begin
   Write(M[y, x], ' ');
  end;
  Writeln;
 end;
end;
var
  i, j, k: integer;
  h: real;
  a: matrix;
  b: vector;
  x: vector;
begin
   LoadFromFile(a, b);
 Print(a);
  for i := 1 to n do
  begin
    for j := i + 1 to n do
    begin
      a[j, i] := -a[j, i] / a[i, i];
      for k := i + 1 to n do a[j, k] := a[j, k] + a[j, i] * a[i, k];
      b[j] := b[j] + a[j, i] * b[i];
    end;
  end;
  x[n] := b[n] / a[n, n];
  for i := n - 1 downto 1 do
  begin
    h := b[i];
    for j := I + 1 to n do
      h := h - x[j] * a[i, j];
    x[i] := h / a[i, i];
  end;
  writeln('Ðåøåíèå ñèñòåìû:');
  for i := 1 to n do
    writeln('X(', i, ')=', x[i]:3:5);

end.

и матрица в текстовом файле-

2.16
1.96
1.56
3.55
3.23
2.78
4.85
4.47
3.97
13.16
21.73
29.75

нужно добавить вектор невязки и обратную матрицу,определитель,умножение исходной на обратную,прошу помогите!может у кого уже есть готовая.

[ironhide8]

помогите!эту программу возможно вообще написать?мне кажется она анриал((п.с.я в долгу тоже не останусь.

[ironhide8]

эхехе...столько просмотров и не капельки помощи((похоже никто на этом сайте не в силах мне помочь((

[p51x]

Код для обратной матрицы, определителя и умножения матриц легко гуглится.

столько просмотров и не капельки помощи((похоже никто на этом сайте не в силах мне помочь((

Вы тоже людей поймите. Вы взяли чей-то, не разбираясь плюнули на форум и хотите получить готовый код для ссвоей задачи... Причем не 2 * 2

п.с.я в долгу тоже не останусь.

Ну дык с этого и стоило начинать, раз сами делать не хотите. Раздел фриоанса и будет вам усе.

[ironhide8]

так может договоримся?