ВУЗ:
Составители:
94
revers;
K1:=K1+1;
end { end while }
end; { end rmatr }
{ *************************************************************** }
begin
iteration:=0; m:= 500;
nachalo:
systemUraw;
for i:=1 to n do b[i]:= -f[i];
for j:=1 to n do
begin
x1:=x[j]; h:=e*abs(x1);
x[j]:=x1+h;
systemUraw;
for i:=1 to n do a[i,j]:=(f[i]+b[i])/h;
x[j]:=x1;
end; { end for j }
(* подсчет количества итерации *)
iteration:= iteration+1;
if iteration = m+1 then begin
writeln('Сделано максимальное число итераций ', iteration-1); goto fin; end;
(* обращение матрицы *)
rmatr(n,a);
(* определение приращений dx[i] и новых значений x[i] *)
for i:=1 to n do
begin
dx[i]:=0;
for j:=1 to n do
dx[i]:=dx[i]+a[i,j]*b[j];
x[i]:=x[i]+dx[i];
end;
(*проверка условии сходимости*);
r:=0;
for i:=1 to n do
if abs(dx[i]/x[i])>e then r:=r+1;
if r=0 then begin
(* решение системы уравнении *);
{ for i:=1 to n do writeln ('x[i]=',x[i]:5:3); }
{ writeln ('числоитерации= ',iteration); }
goto fin;
end;
Nitro PDF Trial
www.nitropdf.com
