Programs for Bivariate model