(%i5) algsys([product(x-i,i,1,75)],[x]); Factor ran out of primes. Does an option variable control the number of primes? Is the default too small? By the way, solve([product(x-i,i,1,75)],[x]) is OK. Barton