adding fit_bivariate_normal