function c = cols(M)
    c = size(M,2);