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

[ZX Spectrum-128]

E

For i:=1 to n do
  Begin
    Sum:=0;
    For k:=1 to n do
       Sum:=sum+a[i,k];
    B[i]:=sum;
  End;

Не проверял, проверьте сами.

[now2]

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

E

For i:=1 to n do
  Begin
    Sum:=0;
    For k:=1 to n do
       Sum:=sum+a[i,k];
    B[i]:=sum;
  End;

Не проверял, проверьте сами.

вот так будет ? посмотрите пожалуйста
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');
n:=12;
For i:=1 to n do
For j:=1 to n do
A[i,j]:=-1/(i+j);
write('a[',i,',',k,']=');a[i,k]:=-1/(i+k);a5[i,k]:=a[i,k];
end;
writeln('Vvedite svobodnye chleny B');
For i:=1 to n do
Begin
Sum:=0;
For k:=1 to n do
Sum:=sum+a[i,k];
B[i]:=sum;
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]

Вроде да. Забейте в паскальабц и проверьте заполнение b.

[now2]

указывает ошибку labaa4444.pas(16) : Ожидалось '.'

[ZX Spectrum-128]

For i:=1 to n do
For j:=1 to n do begin // begin забыли
A[i,j]:=-1/(i+j);

Я там begin забыл

[now2]

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

Я там begin забыл

когда нажимаю шифт+9 вот что выдает ? это как понимать

[ZX Spectrum-128]

Подождите минут 5-7 до компа доберусь.

[now2]

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

Подождите минут 5-7 до компа доберусь.

хорошо) спасибо Вам большое за помощь

[ZX Spectrum-128]

Проверил заполнение b. Заполняется верно.
Остальное, наверное
Но, раз вы сказали что на тестовых данных верно, то и здесь должно быть норм.

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;
sum,t,x,x1,x2,x3,a1,b1,c1:real;
BEGIN
clrscr;
//write('Vvedite razmernost');
n:=12;
For i:=1 to n do
For k:=1 to n do begin
//A[i,j]:=-1/(i+j); не нужно!
write('a[',i,',',k,']=');a[i,k]:=-1/(i+k);a5[i,k]:=a[i,k];
end;
//writeln('Vvedite svobodnye chleny B');
For i:=1 to n do
Begin
Sum:=0;
For k:=1 to n do
Sum:=sum+a[i,k];
B[i]:=sum;
end;
writeln;
writeln('Nachalnaya matriza');
For i:=1 to n do begin
For k:=1 to n do begin
write('a[',i,',',k,']=',a[i,k]:4:2,' ');
end;
writeln('b[',i,']=',b[i]:4:2,' ');
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;
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,'':2);
end;
writeln('b[',i,']=',b[i]:5:2,'':2);
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

Проверил заполнение b. Заполняется верно.
Остальное, наверное
Но, раз вы сказали что на тестовых данных верно, то и здесь должно быть норм.

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;
sum,t,x,x1,x2,x3,a1,b1,c1:real;
BEGIN
clrscr;
//write('Vvedite razmernost');
n:=12;
For i:=1 to n do
For k:=1 to n do begin
//A[i,j]:=-1/(i+j); не нужно!
write('a[',i,',',k,']=');a[i,k]:=-1/(i+k);a5[i,k]:=a[i,k];
end;
//writeln('Vvedite svobodnye chleny B');
For i:=1 to n do
Begin
Sum:=0;
For k:=1 to n do
Sum:=sum+a[i,k];
B[i]:=sum;
end;
writeln;
writeln('Nachalnaya matriza');
For i:=1 to n do begin
For k:=1 to n do begin
write('a[',i,',',k,']=',a[i,k]:4:2,' ');
end;
writeln('b[',i,']=',b[i]:4:2,' ');
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;
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,'':2);
end;
writeln('b[',i,']=',b[i]:5:2,'':2);
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.

елы -палы это что ))))
где ж там ответ то ))???