Sarus | Kramerova Pravila | 3 jednacine 3 nepoznate

program project1;

type Matrica = array [1 .. 3,1 .. 3] of real;
     Vektor = array [1 .. 3] of real;

var A : Matrica; //Glavna matrica
    B,R : Vektor; //Vektor B unosimo. Vektor R su rješenja sistema
    DetSis : real;

procedure UnesiMatricu(Ime : Char ; var X : Matrica); //n ne trebamo unositi zato sto znamo da je n=3

var i,j : integer;

begin
  for i:=1 to 3 do
    for j:=1 to 3 do
      begin
        write(Ime,'[',i,',',j,']= ');
        read(X[i,j]);
      end;
end;

procedure UnesiVektor(Ime : Char ; var X : Vektor);

var i : integer;

begin
  for i:=1 to 3 do
      begin
        write(Ime,'[',i,']= ');
        read(X[i]);
      end;
end;

function Determinanta(X : Matrica) : real;

var det : real;
    i,j : integer;
begin
  //Sarusovo pravilo
  det:=(X[1,1]*X[2,2]*X[3,3])+(X[1,2]*X[2,3]*X[3,1])+(X[1,3]*X[2,1]*X[3,2])-(X[1,3]*X[2,2]*X[3,1])-(X[1,2]*X[2,1]*X[3,3])-(X[1,1]*X[2,3]*X[3,2]);
  Determinanta:=det;
end;
function Zamijeni(k :integer ; Y: Vektor ; X: Matrica) : Matrica;

  var i,j : integer;
      d : Matrica;

  begin
    for i:=1 to 3 do
      for j:=1 to 3 do
        begin
          d[i,j]:= X[i,j];
          d[i,k]:= Y[i];
        end;
    Zamijeni:=d;
  end;

var i : integer;
begin
     UnesiMatricu('A',A);
     UnesiVektor('B',B);
     DetSis:=Determinanta(A);
     writeln('Determinanta sistema ',DetSis:7:2);
     for i:=1 to 3 do
       begin
         R[i]:=Determinanta(Zamijeni(i,B,A))/DetSis;
         writeln(R[i]:7:2);
       end;
end.