A k-dimensional multivariate normal distribution is made discrete by partitioning the k-dimensional Euclidean space with rectangular grids. The probability integrals over the partitioned cubes forms a k-dimensional contingency table with ordered categories. A loglinear model with main effects plus two-way interactions provides an approximation for the k-dimensional table. (SLD)