double chi(double b);
int init(double bmax, double rela_err_limit, int *size, double **x, double **y);