OPTIONS LS=64 PS=60 NONUMBER NODATE; DATA IN; RETAIN TMP 0; n = 10; *** NUMBER OF BERNOULLI TRIALS; p = 1/6; *** PROB. OF A TRIAL SUCCESS; DO x = 0 TO N; CUM_PROB = PROBBNML(p,n,x); PROB = CUM_PROB - TMP; TMP = CUM_PROB; OUTPUT; END; PROC PRINT DATA=IN; ID n p; VAR x PROB CUM_PROB; FORMAT PROB 11.10 CUM_PROB 11.10; TITLE 'BINOMIAL DISTRIBUTION --- n=10 p=1/6'; RUN;