function bifurkacija(a1,a2,n,m,m1)
%bifurkacija(a1,a2,n,m,m1)
%narise bifurkacijski diagram, ki ga dobimo pri resevanju enacbe
%x=a*x*(1-x)
%po metodi fiskne tocke (aka navadna iteracija)
%za n vrednosti parametra a na intervalu [a1,a2]
%narisemo clene zaporedja od m1 do m

%draws the bifurcation diagram obtained by solving the equation
%x=a*x*(1-x) 
%using the fixed point method for 'n' values of the 
%parameter 'a' between 'm1' and 'm'

a=linspace(a1,a2,n);%vrednosti parametra a / values of a
x=ones(1,n);%zacetne vrednosti za x / initial value of x
clf;
hold on
for i=1:m
    x=a.*x.*(2-x);%navadna iteracija, ki poteka za vse a in x 'hkrati'
    %/ computes the next approximation for all values x 
    if i>m1
        scatter(a,x,0.1,[1 0 0]);%ce je i>m1 narisemo trenutne priblizke x
        %/ plotting current values of x
    end
end