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 ( );