bs_ciuupi function

Computes the the functions $b$ and $s$that specify the CIUUPI for all possible values of $\sigma$ and the observed response vector