typedef struct {int n; double *x, *y, *p, *a;} qspline;

qspline* qspline_make(int n,double* x,double* y);
void qspline_free(qspline *S);
double qspline_eval(qspline *S, double z);