function x=st1(mi,x0,ni)

% Serie temporal (mi, x0 inicial, numero de interacoes)

if nargin<1 mi=.5; end
if nargin<2 x0=.55; end
if nargin<3 ni=5000; end

x(1)=x0; n=1; mi, x0, ni,
while(n<ni)
    x(n+1)=mi*x(n)*(1-x(n));
    n=n+1;
end
figure(1); clf;
hold on;
h1=plot(1:ni,x,'-b','erasemode','none');
h2=plot(1:ni,x,'.k','markersize',8,'erasemode','none');
