Simple code for bifurcation of logistic map

logistic_map.py

# Library import for plotting
import matplotlib.pyplot as plt
%matplotlib inline

# Define Logistic Map
def f(x, r):
    return r * x * (1 - x)

# Define history and parameters
history = [0.5]
num_iter = 400

# =========================
# Try 3.4, 3.5, 3.6, 3.7 
r = 3.4
# =========================

# Iterate
for n in range(400):
    x_old = history[-1]
    x_new = f(x_old, r)
    history.append(x_new)

# Plot
plt.scatter(range(num_iter + 1), history)
plt.ylim((0, 1))