function cs = mycumsum(vec)
  cs = sym(zeros(length(vec),1));
  cs(1) = vec(1);
  for i = 2:length(vec)
      cs(i) = cs(i-1) + vec(i);
  end