> simplify(1-a[1]+sum(product(a[i],i=1..k)*(1-a[k+1]),k=1..9));

Reference: Chrystal II, p. 410