Program Jacobi;

 uses crt;

const N=50;



var a: array[1..3,1..3] of real;

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

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

    xx: array[1..3] of real;

   i,j,k,l,m,p:integer;

  begin

    clrscr;

    writeln('Reseni 3 linearnich rovnic o 3 neznamych:');

    for i:=1 to 3 do

      for j:=1 to 3 do begin write('zadej a[',i,j,']:');readln(a[i,j]); end;

    for i:=1 to 3 do  begin write('zadej b[',i,']:');readln(b[i]); 

end ;



    for i:=1 to 3 do x[i]:=0;

    { Chceme-li sami navolit pocet iteraci:

    Zruste konstantu, deklarujte promenou N.

    write ('Zadej kolik iteraci chces provest:'); readln(N);}



    for i:=1 to N  do

       begin

       for k:=1 to 3 do

                     begin

                      xx[k]:=0;

                     for m:=1 to 3 do if k<>m then

                                              xx[k]:= xx[k]+(-a[k,m]*x[m]);

                     xx[k]:= (xx[k]+b[k])/a[k,k];

                     end;

       for p:=1 to  3 do x[p]:=xx[p];

       writeln;

       for l:=1 to 3 do write('  ',x[l]:4:5); writeln;

       end;

     

    writeln;   writeln; 

    repeat until keypressed;

  end.