clear all;
clc;
Nmax=input('Donner la valeur de Nmax=');
x=0:0.01:1;
%f=inline('x-exp(sin(x))');*
f=inline('exp(x)+360*sqrt(x)-272');
plot(x,f(x));
grid;
a=0;
b=1;
eps=1.0e-3;
fa=f(a);
fb=f(b);
if (fa*fb<0)
    for i=1:Nmax
        x=(a+b)/2;
        err=abs(b-a);
        if (fa*f(x)<0)
            b=x;
            fb=f(x);
            if err<eps
                break;
            end
        else
            a=x;
            fa=f(x);
            if err<eps
                break;
            end
        end
        fprintf('Dans l''iteration i=%d\t la solution est x=%f\t avec f(x)=%f\n',i,x,f(x));
    end
    fprintf('La solution finale est x=%f\n',x);
else
    disp('On ne peut pas faire de Dichotomie dans cet intervalle\n');
end