hubereps function

Huber epsilon