compute_var_M.BM function

Computation of the variance.