# TP 11 - Mesure de viscosité - Programme 1
################################################################################
# Cellule 1 : Importation des bibliothèques

import numpy as np
import matplotlib.pyplot as plt
################################################################################
# Cellule 2 : Calcul des 20 valeurs de viscosité

# Grandeurs d'entrée
g =                    # Champ de pesanteur                (m.s-2)      A COMPLETER
µ =                    # Masse volumique du glycerol       (kg.m-3)     A COMPLETER
m =                    # Masse de la bille                 (kg)         A COMPLETER
R =                    # Rayon de la bille                 (m)          A COMPLETER
L =                    # Distance parcourue par la bille   (m)          A COMPLETER

T = np.array([])       # Tableau des vingt mesures de temps de parcours de la longueur L (s)    A COMPLETER

# Formule de calcul de la viscosité
eta =                                                                 # A COMPLETER

# Affichage des valeurs de viscosité
print('Valeurs de la viscosité :',eta)
################################################################################
# Cellule 3 : Calcul de la viscosité moyenne et de son incertitude-type

# Histogramme des valeurs de viscosité
plt.title('Histogramme de la viscosité')
plt.xlabel('viscosité en Pa.s')
plt.ylabel('Effectifs')
plt.hist(eta,bins='rice') 
plt.show()

# Grandeurs statistiques
etamoy = np.mean(eta)                                     # Moyenne des viscosités
sx = np.std(eta,ddof=1)                                   # Ecart-type des viscosités
ueta = sx/np.sqrt(len(eta))                               # Incertitude-type sur la moyenne
print('La moyenne des viscosités est : etamoy =',etamoy)
print("L'écart-type est :",sx)
print("L'incertitude-type sur la viscosité moyenne est : u(eta) =",ueta)
################################################################################