Обчислення визначника методом Гауса (матриця 4х4)
Текст програми на мові Turbo Pascal.
Uses crt;
const n=4;
var
m,v,vv,mm:array [1..n,1..n] of real;
I,j:integer;k,d:real;
begin
writeln(‘введи матрицю’);
for i:=1 to n do
for j:=1 to n do
begin
readln(m[I,j]);
end;
for i:=2 to n do
for j:=1 to n do
begin
k:=m[I,1]/m[1,1];
v[1,j]:=m[1,j];
v[I,j]:=m[1,j]*(-k)+m[I,j];
end;
for i:=3 to n do
for j:=2 to n do
begin
k:=v[I,2]/v[2,2];
mm[1,1]:=v[1,1];
mm[1,j]:=v[1,j];
mm[2,1]:=v[2,1];
mm[2,j]:=v[2,j];
mm[I,1]:=v[I,1];
mm[I,j]:=v[2,j]*(-k)+v[I,j];
end;
for i:=4 to n do
for j:=3 to n do
begin
k:=mm[I,3]/mm[3,3];
vv[I,j]:=mm[3,j]*(-k)+mm[I,j];
end;
for i:=1 to n do
for j:=1 to n do
begin
vv[1,j]:=mm[1,j];
vv[2,j]:=mm[2,j];
vv[3,j]:=mm[3,j];
vv[I,1]:=mm[I,1];
vv[I,2]:=mm[I,2];
d:=1;
if (i=j) then d:=d*vv[I,j];
writeln(vv[I,j]:2:2);
end;
writeln('визначник даної матрици',d:2:2);
end.