double alnorm ( double x, int upper );
void owen_values ( int *n_data, double *h, double *a, double *t );
double tfn ( double x, double fx );
double tha ( double h1, double h2, double a1, double a2 );