As said by Mareš and Mesiar, necessity of aggregation of complex real inputs appears almost in any field dealing with observed (measured) real quantities (see the citation below). For aggregation of probability distributions Sklar designed his copulas as early as in 1959. But surprisingly, since that time only a very few literature have appeared dealing with possibility to aggregate several different pairwise dependencies into one multivariate copula. In the present paper this problem 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. In the latter situation, the IPF procedure tends to cycle (when combining three pairwise dependencies the procedure creates three convergent subsequences), and thus the authors propose some heuristics yielding a "solution'' of the problem even for inconsistent pairwise dependence relations.