data normal;
input b mu sigma;
prob = CDF('Normal',b,mu,sigma);
cards;
8 6 1
4 6 1 
run;

proc print;
run;






data normal2;
input a b mu sigma;
prob = CDF('Normal',b,mu,sigma)-CDF('Normal',a,mu,sigma);
cards;
4 8 6 1
4 8 -2 8
run;

proc print;
run;