Решение СЛАУ Гауссом, ошибка при проверке

@dadaetowi · Регистрация: 13.04.2015

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

Проблема при выводе проверки приведена на скриншоте. Смысл в том что каждое уравнение должно уйти в ноль при правильном подборе корней. Корни выдает все правильные(сам проверял), но в 5-ом уравнение почему-то выдает 3=0.

Код программы

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
uses crt;
 
const
 e = 0.000001;
 
var f:text;
    a,a1: array[1..10, 1..10] of real;
    b,b1: array[1..10] of real;
    x: array[1..10] of real; 
    n, i, j, k: integer;
    z, r, g: real;
    data,s:string;
begin
clrscr;
write('imya fayla: ');
readln(data);
data:=data+'.pas';
assign(f,data);
reset(f);
if IOResult <> 0 then
 begin
  writeLn('file',data,' ne nayden');
  readln;
  exit;
 end;
read(f,n);
for i:=1 to n do
for j:=1 to n do
read(f,a[i,j]);
for i:=1 to n do
read(f,b[i]);
close(f);
writeln('rashirenaya matrica: ');
for i:=1 to n do
 begin
  for j:=1 to n do
  write(a[i,j]:6:2);
  writeln(b[i]:8:2);
 end;
 
  a1[i,j]:=a[i,j];
  for k := 1 to n do
  begin
    for j := k + 1 to n do
    begin
      r := a[j, k] / a[k, k];
      for i := k to n do
      begin
        a[j, i] := a[j, i] - r * a[k, i];
      end;
      b[j] := b[j] - r * b[k];
    end;
  end;
  for k := n downto 1 do
  begin
    r := 0;
    for j := k + 1 to n do
    begin
      g := a[k, j] * x[j];
      r := r + g;
    end;
    x[k] := (b[k] - r) / a[k, k];
  end;
  writeln('OTBET: ');
  for i := 1 to n do
    write('x[', i, ']=', x[i]:0:5, '   ');
  writeln;
 
 
 writeln('PROVERKA:');
  for i:=1 to n do
   begin
    for j:=1 to n do
      b1[i]:=b1[i]+a1[i,j]*x[j];
     writeln(' Uravnenie ',i ,b1[i]:3:0,'=0');
    end;
  readln;  
  
  end.

@Puporev · 16.04.2017, 15:50

Как-то Вы проверку делаете неправильно, все элементы матрицы а1=0 кроме последнего, потому и везде 0=0, кроме последней строки.
Напишите так.

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
uses crt;
var f:text;
    a,a1: array[1..10, 1..10] of real;
    b,b1,x: array[1..10] of real;
    n, i, j, k: integer;
    z, r, g: real;
    data,s:string;
begin
clrscr;
write('imya fayla: ');
readln(data);
data:=data+'.pas';
assign(f,data);
reset(f);
if IOResult <> 0 then
 begin
  writeLn('file',data,' ne nayden');
  readln;
  exit;
 end;
read(f,n);
for i:=1 to n do
for j:=1 to n do
read(f,a[i,j]);
for i:=1 to n do
read(f,b[i]);
close(f);
writeln('rashirenaya matrica: ');
for i:=1 to n do
 begin
  for j:=1 to n do
   begin
    a1[i,j]:=a[i,j];
    write(a[i,j]:6:2);{запомним матрицу А}
   end;
  writeln(b[i]:8:2);
  b1[i]:=b[i]; {Запомним массив В}
 end;
for k := 1 to n do
 begin
  for j := k + 1 to n do
   begin
    r := a[j, k] / a[k, k];
    for i := k to n do
    a[j, i] := a[j, i] - r * a[k, i];
    b[j] := b[j] - r * b[k];
   end;
 end;
for k := n downto 1 do
 begin
  r := 0;
  for j := k + 1 to n do
   begin
    g := a[k, j] * x[j];
    r := r + g;
   end;
  x[k] := (b[k] - r) / a[k, k];
 end;
writeln('OTBET: ');
for i := 1 to n do
write('x[', i, ']=', x[i]:0:5, '   ');
writeln;
writeln('PROVERKA:');
for i:=1 to n do
 begin
  b[i]:=0;
  for j:=1 to n do
  b[i]:=b[i]+a1[i,j]*x[j];
  writeln(' Uravnenie ',i ,b[i]:5:2,'=',b1[i]:5:2);
 end;
readln;
end.

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image 8Observer8 10.02.2026 Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .	Установка Qt-версии Lazarus IDE в Debian Trixie Xfce volvo 10.02.2026 В общем, достали меня глюки IDE Лазаруса, собранной с использованием набора виджетов Gtk2 (конкретно: если набирать текст в редакторе и вызвать подсказку через Ctrl+Space, то после закрытия окошка. . .	SDL3 для Web (WebAssembly): Работа со звуком через SDL3_mixer 8Observer8 08.02.2026 Содержание блога Пошагово создадим проект для загрузки звукового файла и воспроизведения звука с помощью библиотеки SDL3_mixer. Звук будет воспроизводиться по клику мышки по холсту на Desktop и по. . .
SDL3 для Web (WebAssembly): Основы отладки веб-приложений на SDL3 по USB и Wi-Fi, запущенных в браузере мобильных устройств 8Observer8 07.02.2026 Содержание блога Браузер Chrome имеет средства для отладки мобильных веб-приложений по USB. В этой пошаговой инструкции ограничимся работой с консолью. Вывод в консоль - это часть процесса. . .	SDL3 для Web (WebAssembly): Обработчик клика мыши в браузере ПК и касания экрана в браузере на мобильном устройстве 8Observer8 02.02.2026 Содержание блога Для начала пошагово создадим рабочий пример для подготовки к экспериментам в браузере ПК и в браузере мобильного устройства. Потом напишем обработчик клика мыши и обработчик. . .	Философия технологии iceja 01.02.2026 На мой взгляд у человека в технических проектах остается роль генерального директора. Все остальное нейронки делают уже лучше человека. Они не могут нести предпринимательские риски, не могут. . .	SDL3 для Web (WebAssembly): Вывод текста со шрифтом TTF с помощью SDL3_ttf 8Observer8 01.02.2026 Содержание блога В этой пошаговой инструкции создадим с нуля веб-приложение, которое выводит текст в окне браузера. Запустим на Android на локальном сервере. Загрузим Release на бесплатный. . .