% la méthode de fausse position (regula falsi)
% f = x-exp(sin(x)) in [a b].
% la boucle while -- end.
clear all; close all; clc;
tic;
f=inline('x-exp(sin(x))');
a=0; limr=a;
b=10; liml=b;
eps=10^-10;
i=0;
x0=a;
err=abs(a-b);
% le programme principal
    while (err>eps)
        x0=a-f(a)*(b-a)/(f(b)-f(a));
        i= i + 1;
        if f(a)*f(x0)<0
            b=x0; 
        else a=x0; 
        end
        err=abs(f(a)*(b-a)/(f(b)-f(a)));
        % affichage des resultats
        fprintf(' i=%d \t x0=%2.12f \t err=%f \n',i,x0,err);
        % affichage graphique
        fplot(f,[limr liml]);
        hold all;
        plot(x0,f(x0),'ro');
        xlabel('x'),ylabel('f');
        title(['iteration N°=',int2str(i),'  x=',num2str(x0),]);
        hold off;
        grid on;
        pause(1)
    end
toc;