function biflong(ni,x0)

% Diagrama de bifurcacao,(no interacoes, x inicial)

if nargin<1 ni=50000;end
if nargin<2 x0=.55;end
im(5001,6000)=0;
for col=1:6000
    mi=col/5000+2.8
    x=x0;
    for i=1:ni
        x=mi*x*(1-x);
        lin=round(x*5000);
        im(5001-lin,col)=im(5001-lin,col)+1;
    end
    im(:,col)=im(:,col)/max(im(:,col));
end
imwrite(im,'bif-hi-zoom.bmp');
