Program Jacobi;

 uses crt;



const MaxRadku=30;

const MaxSloupcu=30;



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

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

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

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

   i,j,k,l,q,p,M,N,PocetIteraci:integer;

  begin

    clrscr;

    writeln('Reseni M linearnich rovnic o N  neznamych:');

    write('Zadej pocet lin. rovnic:'); Readln(M);

    write('Zadej pocet neznamych:'); Readln(N);

    for i:=1 to M do

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

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



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

    writeln;

    write('Zadej pocet iteraci, ktere se maji provest pri vypoctu:');Readln(PocetIteraci);



    for i:=1 to PocetIteraci  do

       begin

       for k:=1 to N do

                     begin

                      xx[k]:=0;

                     for q:=1 to N do if k<>q then

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

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

                     end;

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

       writeln;

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

       end;







    writeln;   writeln; 

    repeat until keypressed;

  end.