function z = wilkp(x)
% WILKP  Do a good job of computing values of the Wilkinson
%        polynomial, by evaluating it in factored form.

z = ones(size(x));
for i=1:20
  z = z .* (x - i);
end
