# -*- coding: utf-8 -*-
"""
Created on Thu Jun 24 14:13:50 2021

@author: etienne
"""

import numpy as np
import matplotlib.pyplot as plt

### Caractéristiques de la bille :
d = 3.5e-3    # diamètre
rho_b = 7830

### Caractéristiques de l'éprouvette :
h = 20.5e-2     # hauteur de chute
D = 3.6e-2      # diamètre

### Caractéristiques du fluide :
rho_g = 1.250

g = 9.81

### Temps de chute :
Dt = np.array([6.53, 6.19, 6.43, 6.41, 6.31, 6.34, 6.53, 6.62, 6.47, 6.38])

eta = (1/18) * (rho_b - rho_g) * (Dt/h) * d**2 * g / (1 + 2.1*d/D)

eta_moy = np.mean(eta)
Delta_eta = 2*np.std(eta,ddof=1)

print('eta =', "%e"%eta_moy)
print('Delta_eta =', "%e"%Delta_eta)

#plt.figure()
#plt.hist(eta)