% la méthode de fausse position (regula falsi)
% f = x-exp(sin(x)) définie dans l'intervalle [a b].
% laboucle for--end
clear all; close all; clc;
f=inline('x-exp(sin(x))');
a=0; 
b=10; 
eps=10^-10;
i = 0;
Nmax=50;
x0=a;
% le programme principal
for i=1:Nmax
        x0=a-f(a)*(b-a)/(f(b)-f(a));
        %Iter = Iter + 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
        if err<eps
             break; 
        end
        fprintf(' i=%d \t x0=%2.12f \t err=%f \n',i,x0,err);
end