In the paper the problem of description of dependent variables is tackled using the well known Iterative Proportional Fitting Procedure. The proposed solution is not an exact mathematical solution of a marginal problem but just its approximation applicable in many practical situations like Monte Carlo sampling. This is why the authors deal not only with the consistent case, when the iterative procedure converges, but also with the inconsistent non-converging case.