function rysunek[]=fun()
%% funkcja bada ewolucje z roznych punktow startowych
clear all
r=3;
for numer=0:0.01:1+1/r
    elem(round(1+numer*10))=el(numer,r);
end
hold on;
plot(0:2,0:2,'r-');   
plot(0:0.1:1+1/r,elem,'b-');

for numer2=0:0.01:1+1/r,
   numer3=numer2+(rand()/20000); 
   index=0;
if numer3~=numer2 | index<10,
    
     line([numer3,el(numer3,r)],[el(numer3,r),el(numer3,r)],'Color',[0 1 0]);
     line([el(numer3,r),el(numer3,r)],[ el(numer3,r),el(el(numer3,r),r)],'Color',[0 1 0]);
     line([el(numer3,r),el(el(numer3,r),r)],[el(el(numer3,r),r),el(el(numer3,r),r)],'Color',[0 1 0]);
     index=index+1;
     end 
end
hold off;
%% koniec tej funkcji
function element=el(i,r)
element=i+i*r*(1-i*i);
%% koniec tej funkcji
