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

Лабораторная работа №5

Вариант №6

Решение систем уравнений методом Гаусса

Задание:     10x₁-7x₂=7

                     -3x₁+2x₂+6x₃=4

                     5x₁-x₂+5x₃=6

program lab5;

const n=3;

label m1,m2,m3;

var

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

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

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

k,i,j,y: 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]);

m2:

i:=k+1;

m1:

c:=a[i,k]/a[k,k];

a[i,k]:=0;

j:=k+1;

while j<=n do begin

a[i,j]:=a[i,j]-c*a[k,j];

j:=j+1;

end;

b[i]:=b[i]-c*b[k];

if i<n then begin

i:=i+1;

goto m1;

end;

if k<(n-1) then begin

k:=k+1;

goto m2;

end;

x[n]:=b[n]/a[n,n];

i:=i-1;

m3:

j:=i+1;

s:=0;

while j<=n do begin

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

j:=j+1;

end;

x[i]:=(b[i]-s)/a[i,i];

if i>1 then begin

i:=i-1;

goto m3;

end;

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

end.

ОТВЕТЫ:     x₁=0;

                        x₂=-1;

                        x₃=1;