void fresnel ( double x, double *c, double *s );
double fresnel_cos ( double x );
double fresnel_sin ( double x );