sum((a+b)^i, i=1, n)