# File "Rnewt2".

g = function(x) { (1/x) - 1 }

gp = function(x) { -(1/x^2) }

xlist = c(0.000001)

for (i in 1:100) {

    xprev = xlist[length(xlist)]

    xlist = c(xlist, xprev - g(xprev) / gp(xprev))

}

plot(xlist)