function [ D ] = lua( A, n )
y=1;
velikost=size(A);
velikost=velikost(1,1);
H=eye(velikost, velikost);
while y<n+1
[L,U]=lu(A);
A=U*L;
f=diag(A);
D(y,1:velikost)=f';
y=y+1;
end
An=A;
end