% WILKINSON  Simply gives the coefficients of the Wilkinson polynomial
%   P(x) = \prod_{i=1}^{20} (x - i)
% so that P(x) = a_1 x^20 + a_2 x^19 + ... + a_20 x + a_21

a = [1, -210,                  +20615,                 -1256850, ...
        +53327946,             -1672280820,            +40171771630, ...
        -756111184500,         +11310276995381,        -135585182899530, ...
        +1307535010540395,     -10142299865511450,     +63030812099294896, ...
        -311333643161390640,   +1206647803780373360,   -3599979517947607200, ...
        +8037811822645051776,  -12870931245150988800,  +13803759753640704000, ...
        -8752948036761600000,  +2432902008176640000];
  % a possibly convincing picture, showing that we have the right polynomial:
  %    xp = 0:0.001:21;  semilogy(xp,abs(polyval(a,xp)))
