function z = applyInverseIncompleteLU(r)

global Matice l

[L,U] = luinc(Matice{l},'0');
ztmp = L\r;
z = U\ztmp;