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