# -*- coding: utf-8 -*-
import matplotlib
#matplotlib.use('TkAgg')

from pylab import plot, show, text,clf,  ion, ioff, draw, xlabel, ylabel
import numpy as np
from time import sleep

def fun(x,r):
    return r*x*(1-x)
    
u_p=0.2
clf
r=3.7
Niter = 1000
u = np.zeros(Niter)
x = np.zeros(Niter)
for k in range(300):
   for i in range(3*Niter):#czekamy aby się punkt ustabliizował
      u_p=fun(u_p,r) 
   for i in range(Niter):#badamy czy jest stabliny
      u[i]=fun(u_p,r)
      x[i]=r
      u_p=u[i]  
   plot(x,u,'b.',markersize=1);
   r=r+0.001
show()