RootFinding: add pure newton-raphson (in addition to existing robust one)
RootFinding: applying common pattern also to bisection algorithm
RootFinding: default values, newton-raphson initial guess optional (separate function)