Here is another example of how to plot the PDF for discrete distributions (in this case, Binomial and Hypergeometric)with(Statistics):with(plots):These parameters are from Example 5.7, pg. 178n:=5; N:=120; M:=80;X:=RandomVariable(Hypergeometric(N,M,n));plot([seq([i,ProbabilityFunction(X,i)],i=0..6)]);Y:=RandomVariable(Binomial(5,M/N));plot([seq([i,ProbabilityFunction(Y,i)],i=0..6)],style=point,symbol=asterisk);Put them together:A:=plot([seq([i,ProbabilityFunction(X,i)],i=0..6)],style=point,color=red):B:=plot([seq([i,ProbabilityFunction(Y,i)],i=0..6)],style=point,symbol=asterisk,color=black):display({A,B});