function jaca(A,b,x0,k)
[B,c,q]=diagdom(A,b); 
n=length(x0); %az x0 vektor n hosszu  
for i=1:k %ez a ciklus kijelöli melyik iteraciót számoljuk
    x=B*x0+c; %it szamoljuk ki az iteracio értékét
    Jacobi=x
    hiba=(q/(1-q))*norm(x-x0,inf) %itt szaomluk ki az elterest
    %elokészülünk a következo iteracio kiszámításához
    x0=x; %x0 erteket csereljuk a most kiszamított x értékre
end