The first mk_hist() choked on same rare occasions because of round off errors, so I rewrote it - shorter, more to the point and more reliable (I hope).This time it was more convenient to chose lower<=value<upper (exception is the last interval, of course).
