Solution

Algorithm based on simultaneous iterations (Aberth's method) with cluster analysis (see [1]).

Program in Fortran 90 based on the multiprecision package MPFUN, by D. Bailey.

Machine used: IBM Risc 6000.

Time for Newton-isolating all the roots: for n=20, a=1: 0.1 seconds; for n=200, : 0.8 seconds; for n=500, : 3.8 seconds;

Dario Andrea Bini, Dipartimento di Matematica, via Buonarroti 2, 56127 Pisa, ITALY bini@dm.unipi.it October 15, 1996