var_mat function

Estimate variance-covariance matrixes for any parameters of BNs