let histogram nb_columns list =
  let interval_min = Femtolib.List.max (<) list in
  let interval_max = Femtolib.List.max (>) list in
  let width = interval_max -. interval_min in
(*   let interval_max = interval_max +. width /. 1_000_000_000. in *)
  let interval_max = interval_max *. (1.+.epsilon_float) in
  histogram_full interval_min interval_max nb_columns list