double d_epsilon ( void );double d_max ( double x, double y );double d_min ( double x, double y );double d_sign ( double x );void fehl_d ( void f ( double t, double y[], double yp[] ), int neqn,   double y[], double t, double h, double yp[], double f1[], double f2[], double f3[],   double f4[], double f5[], double s[] );void fehl_s ( void f ( float t, float y[], float yp[] ), int neqn,   float y[], float t, float h, float yp[], float f1[], float f2[], float f3[],   float f4[], float f5[], float s[] );int rkf45_d ( void f ( double t, double y[], double yp[] ), int neqn,   double y[], double yp[], double *t, double tout, double *relerr, double abserr,   int flag );int rkf45_s ( void f ( float t, float y[], float yp[] ), int neqn,   float y[], float yp[], float *t, float tout, float *relerr, float abserr,   int flag );float r_epsilon ( void );float r_max ( float x, float y );float r_min ( float x, float y );float r_sign ( float x );void timestamp ( void );