
function dig = dig_func(m0)
for i=1:size(m0,1)
    dig(i,:) = m0(i,i,:);
end
end