# -*- coding: utf-8 -*-
"""
Created on Sat Nov 20 10:51:41 2021

@author: chris
"""

import numpy as np
import matplotlib.pyplot as plt


k = 10000
x0 = 0.5
daten = np.empty([k])
daten[0] = x0


def function(value,r):
    newvalue = r*(1-value)*value
    return newvalue


r = np.linspace(2.9, 4.0, num=1000)
minimum = np.min(r)
maximum = np.max(r)
for rx in r:
    for i in range(1,k):
        daten[i] = function(daten[i-1],rx)
    
    N = 100
    rwerte = np.full(N, rx)
    rest= daten[-N:]
    plt.plot(rwerte, rest,'.',markersize=.3,alpha=.5)

plt.xlabel('r')
plt.ylabel(f'Logistische Folgenglieder(Letzte {N} Folgenglieder)') 

plt.savefig(f'r={minimum}-{maximum}_k={k}_x0={x0}.png', dpi = 300)
plt.show()