Кодирование бифуркационной диаграммы ⇐ Python

[Anonymous]

Мне было интересно, что не так с моим кодом, поскольку я продолжал получать странный, но похожий на реальную бифуркационную диаграмму результат. Я использую итеративное уравнение xn+1 = xn * r(1-xn).
Вот мой код:

Код: Выделить всё

import numpy as np
import matplotlib.pyplot as plt
from functools import lru_cache

@lru_cache(maxsize = 1000)
def bifunc():
R_val = []
X_val = []
R = np.linspace(0.5,4,1000)
for r in R:
x = 0.5
for iterations in range(1001):
x = x*r*(1-x)
R_val.append(r)
X_val.append(x)

plt.plot(R_val, X_val, ls = '', marker = ',')
plt.show()

bifunc()

Вот изображение, которое постоянно появляется:


Любая помощь будет оценена по достоинству. Спасибо.

Подробнее здесь: https://stackoverflow.com/questions/627 ... on-diagram