brent
in nr2.h
extern float brent(float ax, float bx, float cx, float (*f)(float), float tol, float* xmin);