function y=f(x);y=x^2-2;
endfunction;
n=input("n= ")
a=1;b=2
for k=1:n
    c=(a+b)/2
    if f(a)*f(c)<0 then 
        b=c
        else a=c
    end
afficher([k, a, b])
end

