data norms;
mu=500; 
sd=100;
p1=probnorm((600-mu)/sd); *P(Y<600);
p2=1-probnorm((600-mu)/sd); *P(Y>600);
p3=probnorm((700-mu)/sd)-probnorm((400-mu)/sd); *P(400<Y<700);
run;
proc print data=norms;
run;