double clausen ( double x );
double r8_csevl ( double x, double a[], int n );