data normTbl (keep = x normP);
        normP=0;  /*probability distribution table*/
        pden=0;
        do x=-5 to 5 by .01;
                den=(1/sqrt(2*3.14))*exp(-x*x/2);  /*normal probability density function*/
                area=.01*(den+pden)/2;
                normP=normP+area;
                output;
                pden=den;
        end;

proc print;

run;