I never heard of this distribution before, so I have no idea what its properties should be.
See if this does what you want:
d = 1:20;
P = @(d) (1./(d.*(max(d,2)-1))).*(d>1) + (1./max(d)).*(d==1);
figure
plot(d, P(d))
grid
.
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!