import numpy as np
import matplotlib.pyplot as plt
import matplotlib

params = {
    'figure.figsize'    : [8,8],
    'axes.titlesize' : 14,
    'axes.labelsize' : 16,
    'xtick.labelsize' : 14 ,
    'ytick.labelsize' : 14
}
matplotlib.rcParams.update(params)

data = np.genfromtxt("./Doppelpendel_pschnitt.dat")

plt.title(r"$\rm Schnitt \, \theta_2=0 \,\, , \,\, \dot{x}_2>0  \,\, , \,\, E_0 = 15$")
plt.xlabel(r"$\rm \theta_1(t_n)$")
plt.ylabel(r"$\rm \dot{\theta_1}(t_n)$")
plt.scatter(data[:,2],data[:,4],c=data[:,0], cmap=plt.cm.gist_ncar, s=1, alpha=0.7)

#Speichern der Bilder als (.png , benoetigt dvipng unter Linux)- und .pdf-Datei
saveFig="./poincare_15.png"
plt.savefig(saveFig, dpi=300,bbox_inches="tight",pad_inches=0.05,format="png")
#saveFig = "./poincare_15.pdf"
#plt.savefig(saveFig,bbox_inches="tight",pad_inches=0.05,format="pdf")
plt.show()