calibrate_cohon_penalty function

Calibrate penalties with Cohon's method