function[C,D]=skelet(A)
%na vstupu: A - daná matice
%na výstupu: C - hledaná matice m x r
%            D - hledaná matice r x n
[m,n]=size(A);
B=A(:,1);
h1=rank(B);
j=2;
for i=2:n
    B(:,j)=A(:,i);
    h2=rank(B);
    if h2==h1+1
        j=j+1;
        h1=h2;
    else
        B(:,j)=[];
    end;
end;
C=B;
D=B\A;