data binomial;
input x p n;
prob = PMF('Binomial',x,p,n);
cards;
3 0.55 5
6 0.20 15
run;

proc print;
run;






data binomial2;
input a b p n;
prob = CDF('Binomial',b,p,n)-CDF('Binomial',a,p,n);
cards;
2  3 0.55 5
5  6 0.20 15
-1 6 0.20 15 
run;

proc print;
run;