# -*- coding: utf-8 -*-

import random #Permet l'utilisation de la fonction uniform
from scipy import stats
import numpy as np

MC = 176.1 # en g/mol

# La fonction uniform(a,b) permet de générer un nombre aléatoire entre a et b
def alea(v,deltav):
    return(random.uniform(v-deltav/2,v+deltav/2))


#Liste des valeurs dans la formule littérale et de leur erreur:    
Val1 = ... 
deltaVal1 = ...
Vale2= ...
deltaVal2=...
A compléter

#Génération d'une liste de 1000 valeurs de la masse de vitamine C à partir de valeurs générées aléatoirement
Liste_m = []
for i in range(0,1000):
    ValAlea1 = alea(Val1,deltaVal1)
    ValAlea2 = alea(Val2,deltaVal2)
    ... a compléter
    Liste_m.append(Formule à compléter à partir des valeurs générées aléatoirement ValAlea1, ValAlea2, etc...)

#Etude statistique de la liste de valeur obtenue et détermination de l'incertitude type 
Moyenne = np.mean(Liste_m)
sigma = np.std(Liste_m) #Ecart type

print("La moyenne est " + str(Moyenne)+" g")
print("L'incertitude élargie est " + str(sigma)+" g")
