In Master thesis we studied distance-decay of labour commuting at the level of municipalities and statistical regions of Slovenia between years 2009 and 2017. Distance-decay parameters of labour commuting were evaluated according to Newton-Raphson's iterative procedure in three constrained gravity models using Matlab program (Olsson, 2016). The two constraints in Model 1 ensure that the estimated population equals the observed population and the estimated total commuting time equals the observed total commuting time. Model 2 has six constraints that ensure the estimated population and estimated total commuting time on the level of municipalities, regions and between regions. Model 3 includes, in addition to the six constraints in the Model 2, two additional constraints to enforce that the estimated number of workers that live in each municipality is equal to the observed number in origin and destination. The evaluation of parameters in models we presented with graphical representations. We studied the estimation of distance-decay in more detail in the Model 3, which fits the best among all models to observations (Olsson, 2016). During the analysed period, the impact of distance-decay has generally fallen, which indicates a greater willingness of labour commuting for longer distances.