Hello, there is a 1d function with one root in [a,b] intervall, but the solving with brent or bisection or ... comes the next error: gsl: brent.c:74: ERROR: endpoints do not straddle y=0 Default GSL error handler invoked. It is possible the intervall is to big?? (example:[0,10000] and the root is about 10) kt