function dsum = diagsum(A)
  % antar antall_kolonner == antall_rader
  n = length(A);

  dsum = 0;
  for i = 1:n
    dsum = dsum + A(i,i);
  end
end