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