* Example in the review
clear
set obs 6
gen x = _n-1
gen pdf = binomialp(5, x, 0.1)
gen cdf = binomial(5, x, 0.1)

twoway (bar pdf x, sort color(blue) barwidth(0.6) fintensity(inten10))

disp 1 - binomialp(5, 0, 0.1)
disp 1 - binomial(5, 1, 0.1)

* Hights of US Adult Males
clear
set obs 41
gen x = 59.5+_n*0.5
gen pdf = normalden((x-70.9)/2.75)
label variable x "Height in inch (x)"
label variable pdf "Probability Density"
twoway (area pdf x, sort color(blue) fintensity(inten10)), title(Distribution of the Adult Male Heights in the US)