Swapped out newton-raphson solver for simpler bisection method ... compiles alright, and  the bisection method itself works on a simple (nonlinear, noncontinous) scalar function