Add robust jitted root finding routine

[jstac]

I would like to have a jitted robust root finding method such as brent or bisection, following on from PR #417. Bisection at least would be easy.