wilkinson.m

      %script wilkinson

      w1 = poly(1:20); % ce polynme a pour racines les entiers de 1 a 20
      r1 = roots(w1);
      w2 = w1;
      w2(5) = w2(5)+1; % w1(5) est de l'ordre de 10^7, la modification semble faible !
      r2 = roots(w2);