pnorm(double x, double mu, double sigma, int lower_tail, int log_p)
{
assert (lower_tail == 1);
assert (log_p == 0);
assert (sigma == 1.0);
pnorm(double x, double mu, double sigma, int lower_tail, int log_p)
{
assert (lower_tail == 1);
assert (log_p == 0);
assert (sigma == 1.0);