logistic_map.py

from matplotlib import pyplot as plt
from scipy.interpolate import make_interp_spline as mis
import numpy as np

it = 5000
rr = np.arange(2.8, 3.7, 0.01)
for x in np.arange(0.1, 1.0, 0.01):
    xx = []
    for r in rr:
        xn = x
        for _ in range(it): xn = r * xn * (1 - xn)
        xx.append(xn)
    xx = np.array(xx)
    xs = np.linspace(rr.min(), rr.max(), len(rr) * 5)
    rs = mis(rr, xx, k=3)(xs)
    plt.plot(xs, rs, linewidth=0.5)

plt.show()