Pulsed double electron-electron resonance technique (DEER or PELDOR) is applied to study conformations and aggregation of peptides, proteins, nucleic acids, and other macromolecules. For a pair of spin labels, experimental data allows for determination of their distance distribution function, P(r). P(r) is derived as a solution of a first-kind Fredholm integral equation, which is an ill-posed problem. Here, we suggest regularization by the increasing of distance discretization length, to its upper limit where numerical integration still provides agreement with experiment. This upper limit is found to be well above the lower limit for which the solution instability appears because of the ill-posed nature of the problem; so the solution indeed can be regularized in this way. For solving the integral equation, a Monte Carlo trials of P(r) functions is employed. It has an obvious advantage of the fulfillment of the non-negativity constrain for P(r). The approach is checked for model distance distributions and for experimental data taken from literature for doubly spin-labeled DNA and peptide antibiotics. For the case of overlapping broad and narrow distributions, "selective" regularization can be employed in which the effective regularization length may be different for different distance ranges. The method could serve as a useful complement for the traditional approaches basing on Tikhonov regularization.