Метод гаусса с выбором главного элемента

@Maydon · Регистрация: 15.02.2021

Студворк — интернет-сервис помощи студентам

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

Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
program test1;
  type matr=array [1..10,1..10] of real;
  type matr1=array [1..10] of real;
var
  k, l, m, i, j,n: integer;
  q:real;
  a, a1:matr;
  b, b1:matr1;
  
procedure search_main_element(var a:matr);
var
  m,j: integer;
  F:real;
begin
  m := 0;
  F := abs(a[k,k]);
  for j := k + 1 to n do
    begin
    if abs(a[j, k]) > F then
    begin
      F := abs(a[j, k]);
      m := j;
    end;
    end;
end;
 
procedure reshuffle(var a:matr; var b:matr1);
var l:integer;
begin
  for l := k to n do
  begin
    q := a[k, l];
    a[k, l] := a[m, l];
    a[m, l] := q;
  end;
  q := b[k];
  b[k] := b[m];
  b[m] := q;
end;
 
procedure forward_gauss( var a:matr; var b:matr1);
var j,i:integer;
begin
  q := a[k, k];
  for j := k to n do a[k, j] := a[k, j] * q;
  b[k] := b[k] / q;
  for i := k + 1 to n do
  begin
    q := a[i, k];
    for j := k to n do a[i, j] := a[i, j] - (a[k, j] * q);
    b[i] := b[i] - b[k] * q;
  end;
 
end;
 
procedure reverse_gauss( var a:matr; var b:matr1);
var i,j:integer;
begin
  for i := n - 1 downto 1 do
    for j := n downto i + 1 do
      b[i] := b[i] - a[i, j] * b[j];
end;
 
procedure check(var a:matr; var b:matr1);
var
  i,j: integer;
  S:real;
begin
  for i := 1 to n do 
  begin
    S := 0;
    for j := 1 to n do S := S + a[i, j] * b[j];
    writeln(S - b[i]);
  end;
end;
 
 
begin
  n:=3;
  a1:=a;
  b1:=b;
  a[1, 1] := 2; a[1, 2] := -3; a[1, 3] := 1;
  a[2, 1] := 2; a[2, 2] := 1; a[2, 3] := -4;
  a [3,1]:=6;   a[3,2]:=-5;   a[3,3]:= 2;
  b[1]:= 2; b[2]:=9; b[3]:=17; //Ввод переменных
  
 writeln('Матрица A');
 for i:=1 to n do
  begin
    for j:=1 to n do
    begin
      write(a[i,j]:3);
    end;
    writeln;
  end;
  writeln;
  writeln('Матрица B');
  for i:=1 to n do 
  begin
  write(b[i]:3);
  writeln;
  end;
//////////////////////////
 for k:=1 to n do
 begin
    search_main_element(a);
    if m<>0 then  reshuffle(a,b);
    forward_gauss(a,b);
 end;
     reverse_gauss(a,b);
 
 
  writeln('Матрица B конечная');
  for i:=1 to n do 
  begin
  write(b[i]:3);
  writeln;
  end;

@ZX Spectrum-128 · 07.03.2021, 11:17

Сообщение от Maydon

Выдает не верный результат.

Исходный текст вообще ничего не выдает. Нет последнего end.

Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
program test1;
 
type
  matr = array [1..10, 1..10] of real;
 
type
  matr1 = array [1..10] of real;
 
var
  k, l, m, i, j, n: integer;
  q: real;
  a, a1: matr;
  b, b1: matr1;
 
procedure search_main_element(var a: matr);
var
  m, j: integer;
  F: real;
begin
  m := 0;
  F := abs(a[k, k]);
  for j := k + 1 to n do
  begin
    if abs(a[j, k]) > F then
    begin
      F := abs(a[j, k]);
      m := j;
    end;
  end;
end;
 
procedure reshuffle(var a: matr; var b: matr1);
var
  l: integer;
begin
  for l := k to n do
  begin
    q := a[k, l];
    a[k, l] := a[m, l];
    a[m, l] := q;
  end;
  q := b[k];
  b[k] := b[m];
  b[m] := q;
end;
 
procedure forward_gauss(var a: matr; var b: matr1);
var
  j, i: integer;
begin
  q := a[k, k];
  for j := k to n do a[k, j] := a[k, j] * q;
  b[k] := b[k] / q;
  for i := k + 1 to n do
  begin
    q := a[i, k];
    for j := k to n do a[i, j] := a[i, j] - (a[k, j] * q);
    b[i] := b[i] - b[k] * q;
  end;
  
end;
 
procedure reverse_gauss(var a: matr; var b: matr1);
var
  i, j: integer;
begin
  for i := n - 1 downto 1 do
    for j := n downto i + 1 do
      b[i] := b[i] - a[i, j] * b[j];
end;
 
procedure check(var a: matr; var b: matr1);
var
  i, j: integer;
  S: real;
begin
  for i := 1 to n do 
  begin
    S := 0;
    for j := 1 to n do S := S + a[i, j] * b[j];
    writeln(S - b[i]);
  end;
end;
 
 
begin
  n := 3;
  a1 := a;
  b1 := b;
  a[1, 1] := 2;a[1, 2] := -3;a[1, 3] := 1;
  a[2, 1] := 2;a[2, 2] := 1;a[2, 3] := -4;
  a[3, 1] := 6;a[3, 2] := -5;a[3, 3] := 2;
  b[1] := 2;b[2] := 9;b[3] := 17; //Ввод переменных
  
  writeln('Матрица A');
  for i := 1 to n do
  begin
    for j := 1 to n do
    begin
      write(a[i, j]:3);
    end;
    writeln;
  end;
  writeln;
  writeln('Матрица B');
  for i := 1 to n do 
  begin
    write(b[i]:3);
    writeln;
  end;
  //////////////////////////
  for k := 1 to n do
  begin
    search_main_element(a);
    if m <> 0 then reshuffle(a, b);
    forward_gauss(a, b);
  end;
  reverse_gauss(a, b);
  
  
  writeln('Матрица B конечная');
  for i := 1 to n do 
  begin
    write(b[i]:3);
    writeln;
  end;
end.

Добавлено через 23 секунды
Какие должны быть результаты?

@Maydon · 07.03.2021, 21:16 **[ТС]**

Ответ должен быть 5 3 1

@Maydon · 07.03.2021, 21:22 **[ТС]**

Есть блок-схема данного алгоритма, возможно по ней можно найти ошибки:

Метод гаусса с выбором главного элемента

@Maydon · 07.03.2021, 21:24 **[ТС]**

ZX Spectrum-128, Ответ должен получиться 5 3 1, вверху прикрепил блок-схему алгоритма.

Новые блоги и статьи Все статьи Все блоги /
Переходник USB-CAN-GPIO Eddy_Em 20.03.2026 Достаточно давно на работе возникла необходимость в переходнике CAN-USB с гальваноразвязкой, оный и был разработан. Однако, все меня терзала совесть, что аж 48-ногий МК используется так тупо: просто. . .	Оттенки серого Argus19 18.03.2026 Оттенки серого Нашёл в интернете 3 прекрасных модуля: Модуль класса открытия диалога открытия/ сохранения файла на Win32 API; Модуль класса быстрого перекодирования цветного изображения в оттенки. . .	SDL3 для Desktop (MinGW): Рисуем цветные прямоугольники с помощью рисовальщика SDL3 на Си и C++ 8Observer8 17.03.2026 Содержание блога Финальные проекты на Си и на C++: finish-rectangles-sdl3-c. zip finish-rectangles-sdl3-cpp. zip	Символические и жёсткие ссылки в Linux. algri14 15.03.2026 Существует два типа ссылок — символические и жёсткие. Ссылка в Linux — это запись в каталоге, которая может указывать либо на inode «файла-ИСТОЧНИКА», тогда это будет «жёсткая ссылка» (hard link),. . .
[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .	делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .

@Maydon 0 / 0 / 0 Регистрация: 15.02.2021 Сообщений: 4
	07.03.2021, 21:16 [ТС]
	Ответ должен быть 5 3 1 0

@Maydon 0 / 0 / 0 Регистрация: 15.02.2021 Сообщений: 4
	07.03.2021, 21:22 [ТС]
	Есть блок-схема данного алгоритма, возможно по ней можно найти ошибки: 0

@Maydon 0 / 0 / 0 Регистрация: 15.02.2021 Сообщений: 4
	07.03.2021, 21:24 [ТС]
	ZX Spectrum-128, Ответ должен получиться 5 3 1, вверху прикрепил блок-схему алгоритма. 0