function [colsum, rowsum] = sumrowscols(A)
    colsum = sum(A, 1); %or sum(A)
    rowsum = sum(A, 2); %or sum(A')
end
