ЧМ СЛАУ Итерация Опции

[Leito]

нужно написать программу решающую СЛАУ методом итерации. вот система
0.13x1-0.14x2-2x3=0.15
0.75x1+0.18x2+0.77x3=0.11
0.28x1-0.17x2+0.39x3=0.12
она же в виде с преобладающим диагональным коэфициентом
1.16x1-0.13x2-0.84x3=0.38
-0.76x1-2.08x2+0.04x3=0.52
0.13x1-0.14x2-2x3=0.15
или с уравнениями, поделенными на этот самый преобладающий коэфициент
x1-0.112x2-0.724x3=0.328
0.365x1+x2-0.019x3=-0.25
-0.065x1+0.07x2+x1=-0.075
вот что вышло

uses crt;
Label 1;
var A:array[1..3,1..3] of real;
B,X,Y:array[1..3] of real;
n,i,j,k:integer;
e,Nl,Nm,Nk,Sm,Sl,S,Rx,Ry:real;
begin
clrscr;
n:=4;
e:=0.001;
a[1,1]:=1.16; a[1,2]:=-0.13; a[1,3]:=-0.84; b[1]:=0.38;
a[2,1]:=-0.76; a[2,2]:=-2.08; a[2,3]:=0.04; b[2]:=0.52;
a[3,1]:=0.13; a[3,2]:=-0.14; a[3,3]:=-2; b[3]:=0.15;
Nm:=0; Nl:=0; Nk:=0;
for i:=1 to n do
begin
Sm:=0; Sl:=0;
for j:=1 to n do
begin
Sm:=Sm+abs(A[i,j]);
Sl:=Sl+abs(A[j,i]);
Nk:=Nk+sqr(A[i,j]);
end;
if Nm<Sm then Nm:=Sm;
if Nl<Sl then Nl:=Sl;
end;
Nk:=sqrt(Nk);
if (Nm>1) and (Nl>1) and (Nk>1) then
begin
Writeln ('prozess NE zbi}|{nij');
exit;
end
else
begin
Writeln ('prozess zbi}|{nij');
writeln ('Nm=',Nm:6:4,' Nl=',Nl:6:4,' Nk=',Nk:6:4);
end;
k:=0;
for i:=1 to n do Y[i]:=B[i];
1: for i:=1 to n do
begin
s:=0;
for j:=1 to n do s:=s+a[i,j]*Y[j];
X[i]:=B[i]+s;
end;
inc (k);
Rx:=0; Ry:=0;
for i:=1 to n do
begin
Rx:=Rx+abs(x[i]);
Ry:=Ry+abs(y[i]);
end;
if abs(Rx-Ry)>e then
begin
writeln('Iteraziya-',k);
for i:=1 to n do
begin
writeln('X[',i,']=',X[i]:6:4);
Y[i]:=X[i];
end;
goto 1;
end;
writeln('rozviazok sustemu');
writeln('Iteraziya-',k);
for i:=1 to n do writeln('X[',i,']=',X[i]:6:4);
readkey;
end.

пишет что ниодна метрика не подходит, исправьте плз что не верно
ЗЫ в маткаде система решается по любой из трех метрик

Сообщение отредактировано: Leito - 18.01.2009 15:36