СКАЧАТЬ:  elmi6.zip [20,83 Kb] (cкачиваний: 68)

РЕШЕНИЕ СИСТЕМ МЕТОДОМ ПОСЛЕДОВАТЕЛЬНЫХ ПРИБЛИЖЕНИЙ

-0,5X₁+0.1X₂+0.2X₃+0.12X₄=7.1

0.1X₁-0.9X₂+0.1X₃+0.09X₄=2.5

0.11X₁+0.21X₂-0.7X₃+0.2X₄=-1.8

0.21X₁+0.12X₂-0.6X₃+0.2X₄=3.4

         program lab6;

const n=4;

label m1,m2,m3,m4,m5;

var

a:array[1..n,1..n] of real;

b:array[1..n] of real;

x:array[1..n] of real;

y:array[1..n] of real;

k,i,j,e: integer;

s,c:real;

begin

k:=1;

writeln('vvedi massiv a');

for i:=1 to n do

for j:=1 to n do

read(a[i,j]);

writeln('vvedi massiv b');

for i:=1 to n do

read(b[i]);

m1:

x[k]:=0;

a[k,k]:=a[k,k]+1;

k:=k+1;

if k<n then goto m1

else

m2:

i:=1;

m3:

s:=0;

j:=1;

m4:

s:=s+a[i,j]*x[i];

j:=j+1;

if j<n then goto m4

else y[i]:=s-b[i];

i:=i+1;

if i<n then goto m3

else s:=0;

k:=1;

m5:

s:=s+abs(y[k]-x[k]);

x[k]:=y[k];

k:=k+1;

if k<n then goto m5

else begin if s>e  then

goto m3

else

for i:=1 to n do writeln('x[', i,']:', x[i]:8:3,' ');

end;

end.

ОТВЕТЫ: X₁= -7,1;

                    X₂= -2.5;

                    X₃= 1.8;

                    X₄= 0