// On definit la fonction
//-------------------------

function res = F(x); res = 0.5*(x + 200/x); 
endfunction  

// On definit les donnees 
//-----------------------

eps = 0.00001;
x   = 10.;
K = 100;

// On definit la methode de point fixe
//-------------------------------------
// resolution de F(x)=x par la methode du point fixe      
// x valeur initiale puis solution   
// eps reel positif : la precision on arete lorsque |F(x)-x|<eps         
// K entier : nombre max d'itération
// E vecteur reeel : historique de l'erreur |F(x)-x|         
// k entier : nombre d'itérations
            
k=1;
y=F(x);
dx=abs(y-x);
E=dx;
            
while dx>eps & k < K do
x=y;
y=F(x);
dx=abs(y-x);
E=[E,dx];
k=k+1;
end

disp("solution de l equation :  ", [k,x])





