View source for Sayres logstic M-code

function X = logistic(R, N, X0);
% Compute the logistic equation given a control parameter R, 
% a number of iterations N and an initial state X0.
%
%  X = logistic(R, [N=10], [X0=0.5]);
%
%
% ras, 06/02/08.
if nargin < 1, error('Need conrol param R');	end
if notDefined('N'),		N = 10;					end
if notDefined('X0'),	X0 = 0.5;				end

X(1) = X0;

for n = 2:N
	X(n) = R * X(n-1) * (1 - X(n-1));
end

return