gistfile1

------ AVANT -------

Dichotomie := proc(expr, m_a, m_b, e)
	local f, a, b, r, n_steps;

	a := m_a;
	b := m_b;
	n_steps := 0;
	f := unapply(expr, x);

	while abs(a - b) > e do
		r := f((a + b) / 2);

		if r >= 0 then
			b := (a + b) / 2;
		else
			a := (a + b) / 2;
		fi;

		n_steps := n_steps + 1;
	od;

	RETURN([a, b], n_steps);
end;

------ APRES -------

Dichotomie := proc(expr,m_a,m_b,e)
	local f,a,b,r,n_steps;
	a := m_a;
	b := m_b;
	n_steps := 0;
	f := unapply(expr,x);
	while abs(a-b)>e
		r := f((a+b)/2);
		if r>=0 then
			b := (a+b)/2;
		else
			a := (a+b)/2;
		fi;
		n_steps := n_steps+1;
	od;
	RETURN([a,b],n_steps);
end;