We consider the problem of finding checkerboard copulas for modeling multivariate distributions. A checkerboard copula is a distribution with a corresponding density defined almost everywhere by a step function on an m-uniform subdivision of the unit hyper-cube. We develop optimization procedures for finding copulas defined by multiply-stochastic matrices matching available information. https://www.knowall.blog/