function z = applyInverseLU(r)

global Matice l

[L,U] = lu(full(Matice{l}));
ztmp = L\r;
z = U\ztmp;