function Mn = normalize(M,k)
% normalize the columns of the matrix M respect 
% to the k-norm
     for i=1:size(M,2)
        Mn(:,i) = M(:,i)/norm(M(:,i),k);
     end
end