rkf45.h

Runge-Kutta-Fehlberg method

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 );