function sum = geosum(q,n)
sum = 0;
for k = 0:n
  sum = sum + q^k;
end