48
48
false_position(double *a, double *fa, double *b, double *fb,
49
49
objective_function f, void *f_extra,
50
50
double abserr, double relerr, double bisect_til,
53
53
double x1=*a, f1=*fa, x2=*b, f2=*fb;
54
54
fsolve_result_t r = FSOLVE_CONVERGED;