Skip to contents

solve for roots using the NR update rule x1 = x0 − f(x0)/f'(x0)

Usage

newtonRoot(f, df, x0, tol = 0.00001, maxit = 20)