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

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

err=1; i=0;
% 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;
else x0=b;
end
    while (err >= eps)
        i=i+1;
        x1 = x0 - (f(x0)/df(x0));
        err=abs(f(x0)/df(x0));
           % affichage des resultats
        fprintf('i=%d\t racine=%2.8f\t f(x0)=%2.8f\t err=%2.8f\t \n',i,x0,f(x0),err);
           % x0 prend la valeur de x1
        x0 = x1;
           % affichage de la fonction
        fplot(f,[a b]);
        hold all;
        plot(x0,f(x0),'ro');
                   %xlabel('x'),ylabel('f');
                   %title(['iteration N°=',int2str(i),'  x0=',num2str(x0),]);
        hold off; grid on;
        pause(1)
    end