function A = logistic_iteration(R,M,N)
% logistic_limit(R,M,N) limit orbit of the logistic map
% for each value r in R, a row vector of the Mth to the (M+N)th
% iterates of the logistic map at r is generated.

x = rand(1,length(R));
A = zeros(N,length(R));

for n = 1:M
x = logistic(x,R);
end

for n = 1:N
A(n,:) = x;
x = logistic(x,R);
end

A = A';
