sum(product (1-a*i), i=0..j), j=0..n