# File "Riter". f = function(x) { x^2 * exp(-(x+sin(x))) } # curve(f, -0.5, 12) # curve(f, 8, 12) fp = function(x) {2*x*exp(-(x+sin(x))) + x^2*exp(-(x+sin(x)))*(-1-cos(x))} h = function(x) {fp(x) + x} numits = 5000 initval = 1 vallist = c(initval) for (i in 1:numits) { vallist = c(vallist, h(vallist[length(vallist)])) } plot(vallist)