Biocaml 0.4-dev : Biocaml_math.epsilon

let epsilon f init fin =

 let rec aux acc n = if n = fin then acc else aux (acc +. (f n fin)) (n + 1) in

 aux 0. init