Решение системы методом Гаусса-Жордана / PascalABC.NET

[now2]

Задания
I. Напишите программу решения системы линейных уравнений
методом Жордана-Гаусса.
II. Решите систему линейных уравнений и осуществите проверку
найденного решения.

на 1 рис. сама система
на 2 рис. третье условие
на 3 рис.задание к третьему условию

[now2]

вот программа , проверьте плиз она написана в турбо паскаль как ее переделать на PaskalABC.NET а то она не работает указывает на ошибку Program1.pas(41) : Нельзя изменять значение переменной цикла внутри этого цикла


program labI4;
Uses crt;
Label m1;
Var a,a5: array [1..25,1..25] of real;
b,b5: array [1..25] of real;
k,i,m,n,st:integer;
t,x,x1,x2,x3,a1,b1,c1:real;
BEGIN
clrscr;
write('Vvedite razmernost');
readln(n);
For i:=1 to n do
For k:=1 to n do begin
write('a[',i,',',k,']=');readln(a[i,k]);a5[i,k]:=a[i,k];
end;
writeln('Vvedite svobodnye chleny B');
For i:=1 to n do begin
write('b[',i,']=');readln(b[i]);b5[i]:=b[i];
end;
writeln('Nachalnaya matriza');
For i:=1 to n do begin
For k:=1 to n do begin
write('a[',i,',',k,']=',a[i,k]:5:2,'':4);
end;
writeln('b[',i,']=',b[i]:5:2,'':4);
end;
readln;
For k:=1 to n do begin
st:=k;
For i:=k to n do
If a[st,k]<=a[i,k] Then st:=i;
For i:=1 to n do begin
t:=a[k,i];
a[k,i]:=a[st,i];
a[st,i]:=t;
end;
t:=b[k];
b[k]:=b[st];
b[st]:=t;
For i:=1 to n do begin
If i=k Then i:=i+1;
If i=n+1 Then goto m1;
t:=a[i,k]/a[k,k];
For m:=1 to n do
a[i,m]:=a[i,m]-t*a[k,m];
b[i]:=b[i]-t*b[k];
end;
m1:
end;
writeln(' Konechnaya matriza');
For i:=1 to n do begin
For k:=1 to n do begin
write('a[',i,k,']=',a[i,k]:5:2,'':4);
end;
writeln('b[',i,']=',b[i]:5:2,'':4);
end;
readln;
For k:=1 to n do begin
x:=b[k]/a[k,k];
writeln('x(',k,')=',x:5:2);
a[1,k]:=x;
end;
readkey
END.

[ZX Spectrum-128]

Смените цикл на while.
Что-то я тут вводил, вроде считает.

program labI4;
Uses crt;
Label m1;
Var a,a5: array [1..25,1..25] of real;
b,b5: array [1..25] of real;
k,i,m,n,st:integer;
t,x,x1,x2,x3,a1,b1,c1:real;
BEGIN
clrscr;
write('Vvedite razmernost');
readln(n);
For i:=1 to n do
For k:=1 to n do begin
write('a[',i,',',k,']=');readln(a[i,k]);a5[i,k]:=a[i,k];
end;
writeln('Vvedite svobodnye chleny B');
For i:=1 to n do begin
write('b[',i,']=');readln(b[i]);b5[i]:=b[i];
end;
writeln('Nachalnaya matriza');
For i:=1 to n do begin
For k:=1 to n do begin
write('a[',i,',',k,']=',a[i,k]:5:2,'':4);
end;
writeln('b[',i,']=',b[i]:5:2,'':4);
end;
readln;
For k:=1 to n do begin
st:=k;
For i:=k to n do
If a[st,k]<=a[i,k] Then st:=i;
For i:=1 to n do begin
t:=a[k,i];
a[k,i]:=a[st,i];
a[st,i]:=t;
end;
t:=b[k];
b[k]:=b[st];
b[st]:=t;
i:=1;
while i<=n do
//For i:=1 to n do 
  begin
    If i=k Then 
      i:=i+1;
    If i=n+1 Then 
      goto m1;
    t:=a[i,k]/a[k,k];
    For m:=1 to n do
      a[i,m]:=a[i,m]-t*a[k,m];
    b[i]:=b[i]-t*b[k];
    i:=i+1;
  end;
  m1:
end;
writeln(' Konechnaya matriza');
For i:=1 to n do begin
For k:=1 to n do begin
write('a[',i,k,']=',a[i,k]:5:2,'':4);
end;
writeln('b[',i,']=',b[i]:5:2,'':4);
end;
readln;
For k:=1 to n do begin
x:=b[k]/a[k,k];
writeln('x(',k,')=',x:5:2);
a[1,k]:=x;
end;
readkey
END.

[now2]

Сообщение от ZX Spectrum-128

Смените цикл на while.
Что-то я тут вводил, вроде считает.

program labI4;
Uses crt;
Label m1;
Var a,a5: array [1..25,1..25] of real;
b,b5: array [1..25] of real;
k,i,m,n,st:integer;
t,x,x1,x2,x3,a1,b1,c1:real;
BEGIN
clrscr;
write('Vvedite razmernost');
readln(n);
For i:=1 to n do
For k:=1 to n do begin
write('a[',i,',',k,']=');readln(a[i,k]);a5[i,k]:=a[i,k];
end;
writeln('Vvedite svobodnye chleny B');
For i:=1 to n do begin
write('b[',i,']=');readln(b[i]);b5[i]:=b[i];
end;
writeln('Nachalnaya matriza');
For i:=1 to n do begin
For k:=1 to n do begin
write('a[',i,',',k,']=',a[i,k]:5:2,'':4);
end;
writeln('b[',i,']=',b[i]:5:2,'':4);
end;
readln;
For k:=1 to n do begin
st:=k;
For i:=k to n do
If a[st,k]<=a[i,k] Then st:=i;
For i:=1 to n do begin
t:=a[k,i];
a[k,i]:=a[st,i];
a[st,i]:=t;
end;
t:=b[k];
b[k]:=b[st];
b[st]:=t;
i:=1;
while i<=n do
//For i:=1 to n do 
  begin
    If i=k Then 
      i:=i+1;
    If i=n+1 Then 
      goto m1;
    t:=a[i,k]/a[k,k];
    For m:=1 to n do
      a[i,m]:=a[i,m]-t*a[k,m];
    b[i]:=b[i]-t*b[k];
    i:=i+1;
  end;
  m1:
end;
writeln(' Konechnaya matriza');
For i:=1 to n do begin
For k:=1 to n do begin
write('a[',i,k,']=',a[i,k]:5:2,'':4);
end;
writeln('b[',i,']=',b[i]:5:2,'':4);
end;
readln;
For k:=1 to n do begin
x:=b[k]/a[k,k];
writeln('x(',k,')=',x:5:2);
a[1,k]:=x;
end;
readkey
END.

а программа написана точно методом Жордана-Гаусса.?????

[ZX Spectrum-128]

Вот этого я не знаю. Не настолько я помню методы решения слау, чтобы с одного взгляда определить, какой это метод. По той ссылке, откуда вы взяли программу утверждается, что это метод Г-Ж?
А на своих данных проверили? Корректно решает?

[now2]

Сообщение от ZX Spectrum-128

Вот этого я не знаю. Не настолько я помню методы решения слау, чтобы с одного взгляда определить, какой это метод. По той ссылке, откуда вы взяли программу утверждается, что это метод Г-Ж?
А на своих данных проверили? Корректно решает?

ну да ) все вроде решает спасибо ). а что делать с третьим заданием не подскажите ?

[ZX Spectrum-128]

IВместо ввода данных вбить заполнение массива по правилу, например для а

N:=12;
For i:=1 to n do
  For j:=1 to n do
     A[i,j]:=-1/(i+j);

[now2]

это вставить туда где "введите размерность" да ?

[ZX Spectrum-128]

Да-да. Только вместо readln(n) напишите n:=12; вместо readln(a[i,k])
Написать a[i,k]:=-1/(i+k);
И т.д.

[now2]

Сообщение от ZX Spectrum-128

IВместо ввода данных вбить заполнение массива по правилу, например для а

N:=12;
For i:=1 to n do
  For j:=1 to n do
     A[i,j]:=-1/(i+j);

это после "введите размерность "для а, а для b после "введите свободный член" , только не подскажите как b правильно записать сумму