%clear all, close all
n=100;%input  ('number of iterations = ');
k=input  ('k                    = ');
p_0=input('p_0                  = ');
p=zeros(n,1);
p(1)=p_0;
for i=2:n
    p(i)=k*p(i-1)*(1-p(i-1));
end
plot(p,'.:')