iterate_logistic = function(rho, u0, N=100){ us = rep( NA, N ) us[1] = u0 for( ii in 2:N ){ us[ii] = rho * us[ii-1] * ( 1 - us[ii-1] ) } us }