huber_estimate function

Huber M-estimate