Практикум по решению задач на ЭВМ. Тарова И.Н - 80 стр.

UptoLike

Тарова И.Н., Терехов Ю.П., Масина О.Н., Скоков А.В.
A:array[1..m,1..m] of real;
F,x,x1:array[1..m] of real;
i,j:integer;
N,Nmax:integer;
E:real;
Z,max_dx,S:real;
begin
clrscr;
writeln('введите матрицу коэффициентов при неизвестных (по стро-
кам)');
for i:=1 to m do
for j:=1 to m do
read(A[i,j]);
writeln('введите столбец свободных членов F1,...,Fm');
for i:=1 to m do
readln(F[i]);
writeln('введите начальное приближение x0');
for i:=1 to m do
readln(x[i]);
writeln('введите точность приближений E');
readln(E);
writeln('задайте максимальное число итераций Nmax');
readln(Nmax);
writeln('Введенная матрица расширенная (расширенная)');
for i:=1 to m do
begin
for j:=1 to m do
write(A[i,j]:6:2);
write(F[i]:6:2);
writeln;{перевод строки}
end;
N:=0;
REPEAT
N:=N+1;{счет числа операций}
{вычисление среднего (n-го) приближения}
for i:=1 to m do
begin
S:=0;
for j:=1 to m do
if j<>i then
S:=S+A[i,j]*x[j];
x1[i]:=(-S+F[i])/(A[i,i]);
end;
80
                             Тарова И.Н., Терехов Ю.П., Масина О.Н., Скоков А.В.


   A:array[1..m,1..m] of real;
   F,x,x1:array[1..m] of real;
   i,j:integer;
   N,Nmax:integer;
   E:real;
   Z,max_dx,S:real;
    begin
     clrscr;
      writeln('введите матрицу коэффициентов при неизвестных (по стро-
кам)');
      for i:=1 to m do
      for j:=1 to m do
       read(A[i,j]);
writeln('введите столбец свободных членов F1,...,Fm');
       for i:=1 to m do
        readln(F[i]);
       writeln('введите начальное приближение x0');
       for i:=1 to m do
        readln(x[i]);
       writeln('введите точность приближений E');
        readln(E);
writeln('задайте максимальное число итераций Nmax');
        readln(Nmax);
writeln('Введенная матрица расширенная (расширенная)');
       for i:=1 to m do
        begin
         for j:=1 to m do
          write(A[i,j]:6:2);
          write(F[i]:6:2);
          writeln;{перевод строки}
        end;
       N:=0;
       REPEAT
        N:=N+1;{счет числа операций}
        {вычисление среднего (n-го) приближения}
        for i:=1 to m do
         begin
          S:=0;
          for j:=1 to m do
           if j<>i then
            S:=S+A[i,j]*x[j];
           x1[i]:=(-S+F[i])/(A[i,i]);
         end;


                                  80