(* Recherche d'un zéro de f sur l'intervale [a,b] par dichotomie.
   On suppose la fonction f monotone et qu'un tel zéro existe. *)

let eps = 1e-9

let rec dicho f a b =
  let x = (a +. b) /. 2. in
  if b -. a < eps then x
  else if f a *. f x < 0. then dicho f a x else dicho f x b

(* test *)

let phi = dicho (fun x -> x *. x -. x -. 1.) 1. 2.

This document was generated using caml2html