function seidelb(A,b,x0,k)
[B,c,q]=diagdom(A,b); 
[n,m]=size(B); %az B mx n sorból és m oszlopból áll
x=zeros(n,1); %az x vektor kezdoerteke az n sor 1 oszlop null vektor
for i=1:k %ez a ciklus kijelöli melyik iteraciót számoljuk
    for j=1:n %ez a ciklus kijelöli melyik koordinátát számoljuk az i-edik iterációban
        x(j)=B(j,1:j-1)*x(1:j-1)+B(j,j:m)*x0(j:m)+c(j); %itt számoljuk ki az iteráció értékét 
    end
    Seidel=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