function y= seuil(q,eps);
   n=0;
    while abs(q^n)>=eps
        n=n+1
    end
    y=n
endfunction
//taper dans la Console par exemple : seuil(-0.8,0.01)
