Het zou heel goed kunnen dat mijn oplossing geen geweldige statistische verdeling oplevert, maar teveel hoge of lage waardes geeft. Ik heb het concept snel in elkaar gezet met als enige reden dat ik een volledig dekkende oplossing wou hebben.

Of het ook enig praktisch nut heeft, speelde voor mij totaal geen rol.

Wat betreft de 1 en 0 problemen, die zijn redelijk eenvoudig te fixen door:

Ps = int(rand() * (max - 2) + 2) 'produkt is minimaal 2
D1 = int(rand() * (Ps/2 - 2) + 2) 'D1 dan ook minimaal 2

D2 kon (en kan) alleen 1 worden als Ps erg klein is. D1 is nl. (bijna) altijd kleiner dan P / 2. Dus D2 is dan altijd groter dan P / (P / 2) = 2P / P = 2.
(of zoiets... )