huber.reg function

Huber estimation for linear regression