huber.reg function

Tuning-free Huber regression