// TP1 Résolution d'équation

//Activité 1 :  xn+1 = (xn-1+xn)/2 ordre 1 (géométrique)
function y=f(x) 
y=x.*x-2
endfunction

function [x,n]=dicho(f,a,b,e)
if f(a)*f(b)>0 then disp("Condition d existence non vérifiée");
                    return;
end
x=a
n=0
while abs(b-a)>e & n<50  
      n=n+1
      x=(a+b)/2
      if (f(a)*f(x)<0) then  a=a
                             b=x
                       else  a=x
                             b=b
      end
end
endfunction