klmnj

type DDN = array [1..3,1..3] of integer;
type JDN = array [1..3] of integer;
type JDNR = array [1..3] of real;

var b : JDN;
var x : JDNR;
var A , A_1 , A_2 , A_3 : DDN;
var  D,D_1,D_2,D_3 : integer;
var c : integer;

  procedure UcitajMat( var X: DDN );
   var i,j: integer;
   begin
     for i := 1 to 3 do
         begin
           for j := 1 to 3 do
               begin
                 writeln('A[',i,',',j,']=');
                 read(X[i,j]);
               end;
         end;
   end;

procedure UcitajNiz( var v: JDN );
var i: integer;
begin
  for i:= 1 to 3 do
      begin
         writeln('b[',i,']=');
         read(v[i]);
      end;
end;

function detA( x : DDN ) : integer;
var rez: integer;
begin
   rez := ((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[3,1]*x[2,2]*x[1,3])-(x[3,2]*x[2,3]*x[1,1])-(x[3,3]*x[2,1]*x[1,2]));
   detA := rez;
end;

procedure MenjajKolone( kolona : integer ; v : JDN ; var x: DDN );
var i: integer;
begin
    for i := 1 to 3 do
      begin
         x[i,kolona] := v[i];
      end;
end;

procedure KopMat( x : DDN ; var Xi : DDN );
var i,j : integer;
begin
    for i := 1 to 3 do
         begin
           for j := 1 to 3 do
               begin
                 Xi[i,j] := x[i,j];
               end;
         end;
end;
Begin
    writeln('Unesite matricu koeficijenata:');
    UcitajMat(A);

    KopMat(A,A_1);
    KopMat(A,A_2);
    KopMat(A,A_3);

    writeln('Unesite vektor b:');
    UcitajNiz(b);

    MenjajKolone( 1 , b , A_1 );
    MenjajKolone( 2 , b , A_2 );
    MenjajKolone( 3 , b , A_3 );


    D := detA(A);
    D_1 := detA(A_1);
    D_2 := detA(A_2);
    D_3 := detA(A_3);


    if ( ( D = 0 ) AND ( D_1 = 0 ) AND ( D_2 = 0 ) AND ( D_3 = 0 ) ) then
    begin
        writeln('Sistem ima beskonacno mnogo resenja!');
        readln(c);
    end;


    if ( (( D = 0 ) AND ( D_1 <> 0 )) OR (( D = 0 ) AND ( D_2 <> 0 )) OR (( D = 0 ) AND ( D_3 <> 0 )) ) then
    begin
        writeln('Sistem nema resenja!');
        readln(c);
    end;


    x[1] := D_1 / D;
    x[2] := D_2 / D;
    x[3] := D_3 / D;

    writeln('Resenja sistema su:');

    for c := 1 to 3 do
        begin
          writeln('x[',c,']=',x[c]);
        end;

     readln(c);
End.