function M = razprtje(X, k)
% M = razprtje(X, k)
% izracuna razprtje tenzorja X v smeri k.

dimenzije = size(X);
vrstice = dimenzije(k);
dimenzije = [dimenzije(1:k-1) dimenzije(k+1:end)]; % k-to lahko pozabimo
stolpci = prod(dimenzije);

M = zeros(vrstice, stolpci);
indeksiStolpcev = obratnoLeksikografsko(dimenzije);
for i = 1:vrstice
    for j = 1:size(indeksiStolpcev, 2)
        multiIndex = vrini(indeksiStolpcev(:, j)', k, i);
        M(i, j) = X(multiIndex{:});
    end
end


end