% la méthode de Lagrange
clear all; close all; clc;
% les parametres d'entrée
a=input('donner a= ');
b=input('donner b= ');
eps=input('donner eps= ');
nmax=input('donner nmax= ');

f=inline('cos(x)-x^3');
df=inline('-sin(x)-3*x^2');
ddf=inline('-cos(x)-6*x');

% le programme principal
% if (df(a)== 0) && (df(b)== 0) 
% disp('On peut pas calculer la racine');
% end
% initialisation de x0
if (f(a)*ddf(a)<0)
x0=a;
c=b;
else x0=b;
     c=a;
end
    for (i=0:nmax)   
  
        x0=(f(x0)*c-x0*f(c))/(f(x0)-f(c));
        err=abs(x0-(f(x0)*c-x0*f(c))/(f(x0)-f(c)));
        
       % affichage des resultats
       if (err<=eps)
         fprintf('i=%d\t racine=%2.8f\t f(x0)=%2.8f\t err=%2.8f\t \n',i,x0,f(x0),err);
       end
         % affichage de la fonction
         fplot(f,[a b]);
         hold all;
         plot(x0,f(x0),'ro');
         hold off; grid on;
         pause(1)
    end