bk_cusum function

Continuous time BK-CUSUM