# Python Programm zum Plotten der Daten des periodisch angetriebenen Pendels (Übungsblatt 12)
# Attraktordiagramm vieler Pendel mit unterschiedlichen Amplitudenwerten (Feigenbaumdiagramm)

import matplotlib.pyplot as plt             # Python Bibliothek zum Plotten (siehe https://matplotlib.org/ )
import numpy as np                          # Python Bibliothek fuer Mathematisches (siehe https://numpy.org/ )

data_attr = np.genfromtxt("./GetriebenesPendel_AttrDiag.dat")  # Einlesen der berechneten Daten

plt.xlabel(r"$\rm A$")
plt.ylabel(r"$\rm \dot{\theta}(t_n)$")
for i in range(len(data_attr[0,:])-1):
    plt.plot(data_attr[:,0],data_attr[:,i+1] , 'o',c="black", markersize=0.3)
    
# Bild-Ausgabe
plt.savefig("Attr_diag.png", dpi=400,bbox_inches="tight",pad_inches=0.05,format="png")
plt.close()