function [B,c,q]=diagdom(A,b) [n,m]=size(A); %az A mx n sorból és m oszlopból áll B=eye(n); %a B mx kezdoértéke az n-sor n-oszlopból álló egység mátrix c=zeros(n,1); %c az n-sor 1-oszlopból álló zero vektor for i=1:n %ez a ciklus kijelöli melyik sorban dolgozunk t=A(i,:)/A(i,i);%az adott sorban végigosztunk a foátlóbeli elemmel A(i,i) B(i,:)=B(i,:)-t; %az egységmátrixból elveszük az fenti t értéket c(i)=b(i)/A(i,i); %a c vektor i-edik korrdinátáját osztjuk az adott sor foátlo beli elemével end q=norm(B,inf);