function [y,korak]=hagerinv(A)

% [y,korak]=HAGERINV(A) vrne Hagerjevo cenilko za 1-normo inverzne matrike A
% skupaj s stevilom potrebnih korakov 
%
% Pravilna uporaba bi bila, da najprej resimo sistem Ax=b, nato pa 
% P,L in U uporabimo za resevanje sistema z matriko A pri racunanju cenilke

n=size(A,2);
x=1/n*ones(n,1);
korak=1;
while 1==1
  w=A\x; sigma=sign(w); z=A'\sigma;
  if norm(z,'inf')<=z'*x
    y=norm(w,1);
    break
  else
    [tmp,j]=max(abs(z));
    x=zeros(n,1); x(j,1)=1;
  end
  korak=korak+1;
end

    