# File "Riter2".

g = fp = function(theta) {+ 4*(1-theta)^3 + 4*(2-theta)^3 }

h = function(x) {fp(x) + x}

numits = 100

initval = 1.51

vallist = c(initval)

for (i in 1:numits) {
    vallist = c(vallist, h(vallist[length(vallist)]))
}

plot(vallist)