int newton (double *x, double tolf, double tolx, int maxit, void (*fdf)(double,double*,double*,void*), void *prm, int ixrr);
double neville (int n, double xi[], double yi[], double x);