spi := proc(a,n)
local j;
#|a|<<1, (n=20, a=i/1000, a=(3+4 i)/10)
  product(( x + (a^j-1))/(a-1),j=1..n);
end: