double alnorm ( double x, int upper );
void normal_01_cdf_values ( int *n_data, double *x, double *fx );
void normp ( double z, double *p, double *q, double *pdf );
void nprob ( double z, double *p, double *q, double *pdf );