let variation = function
  f -> (function i  -> f (lower_bound (i)) - f (upper_bound (i)))