compute_var_diff.BM function

Compute variances of differences between nodes and parents.