add Huber function, elaborate on simulation
elaborate on simulation