extern double fabs(), floor(), ceil(), ldexp(); extern double sqrt(), hypot(), atof(); extern double sin(), cos(), tan(), asin(), acos(), atan(), atan2(); extern double exp(), log(), log10(), pow(); extern double sinh(), cosh(), tanh(); extern double gamma(); extern double j0(), j1(), jn(), y0(), y1(), yn(); extern double frand(); extern double erf(), erfc(); #define HUGE 1.701411733192644270e38 #define PI 3.1415926535897932384626433832795028841971693993751 #define PI2 (2*PI)