bracket = b.bracketOut(0.0, y, Double.POSITIVE_INFINITY);
} else {
bracket = b.bracketOut(Double.NEGATIVE_INFINITY, y, 0.0);
}
BisectionRootFinder bisection = new BisectionRootFinder(f);
ret = bisection.findRoot(bracket[0], bracket[1]);
}
return ret;
}