  Vandermonde matrices are exponentially ill-conditioned, rendering the
familiar "polyval(polyfit)" algorithm for polynomial interpolation and
least-squares fitting ineffective at higher degrees. We show that Arnoldi
orthogonalization fixes the problem.
