# -*- coding: utf-8 -*-
"""
Spyder Editor

This is a temporary script file.
"""

import numpy as np
from math import log

N=100000   #Nombre de tirages simulés

DeltaG = []     #Liste des valeurs calculées de la concentration de la solution S

for k in range(N) :                 #Procédure de tirage
    Veq = np.random.uniform(      )
    V = np.random.uniform(        )
    Cb = np.random.uniform(       )
    T = np.random.uniform(        )
    DeltaG.append(                )     #Calcul de DeltaG

DeltaG_moy = sum(DeltaG)/N                #Calcul de la valeur moyenne
u_DeltaG = np.std(DeltaG,ddof=1)           #Calcul de l'écart-type

#Affichage des résultats
print(f'Enthalpie libre de dissolution : {DeltaG_moy} J/mol') 
print(f'Incertitude-type u : {u_DeltaG} J/mol')