The following simple code can be used to find roots of functions (based on the secant algorithm), secant=function(fun, x0, x1, tolerence=1e-07, niter=500){for ( i in 1:niter ) { x2 [Read more...]