#ifndef INTERPOLATE_H
#define INTERPOLATE_H
#ifdef __cplusplus
extern "C" {
#endif
double interpolate( double x, const double xa[], const double fa[], int N_elems );
#ifdef __cplusplus
}
#endif
#endif