import numpy as np
import matplotlib.pyplot as plt

#DETERMINATION DU COUPLE RESISTANT SEC ET DU FROTTEMENT VISQUEUX
T_omega=np.array([97.7,196,293.5,391.9])
T_cr=np.array([0.00241,0.00256,0.00307,0.00315])

plt.plot(T_omega,T_cr)
plt.grid()
plt.ylim(0,0.0032)
plt.show()

a,b=np.polyfit(T_omega,T_cr,1)
print('f=',a,'\nCrsec=',b)


#AFFICHAGE DE L'ERREUR STATIQUE EN FONCTION DE Kpepos
L_Kpepos=[60,120,180,240]
L_es=[6.5,3.9,2.75,2.14]
plt.plot(L_Kpepos, L_es)
plt.grid()
plt.ylim(0,7)
plt.show()