bk_control_limit function

Determine control limits for BK-CUSUM by simulation