function m=cumsumc(x);
% get cumulative sum of each column
m=cumsum(x);
if size(m,2)>1;
   m=m';
end;
return;