double clausen ( double x );
void clausen_values ( int *n_data, double *x, double *fx );
double r8_csevl ( double x, double a[], int n );
double r8_wrap ( double r, double rlo, double rhi );
void timestamp ( );