hs(n):=ratsimp(ratexpand(product(x-a[i],i,1,n)));
hs(15);