for i=0:20
    epsi = 10./2^i;
    x0 = [0; 0; 0; 0; 0; 0; 0; 0; 0; 1]; 
    x1 = puissance(epsi,A,x0);
    maxi = x1(1);
    ind  = 1;
    
    for i=2:10
        if (maxi < x1(i))
         maxi = x1(i); ind = i;
        end
    end
    sprintf('Le nombre meilleur score est %i', ind)
end