added some modiefied psuedocode for likelihood
Multivariate normal distribution object (with rho parameter)
discMVNorm discretizes a bivariate normal given the distribution and a grid. it uses a ranked version of the grid by default.