newtonRoot.Rd
solve for roots using the NR update rule x1 = x0 − f(x0)/f'(x0)
newtonRoot(f, df, x0, tol = 0.00001, maxit = 20)