Vi = [380.11, 379.99, 379.81,380.00, 380.19, 380.06, ...
      379.99, 379.96, 380.02, 380.10] 
Um = mean(Vi)
sigma_Vi = std(Vi)
sigma_Vim = sigma_Vi / sqrt(length(Vi))

V_lect = 0.01/2 / sqrt(3)
V_tol = Um*0.0005 / sqrt(3)

uVi = sqrt(sigma_Vi^2 + V_lect^2 + V_tol^2)

uVm = sqrt(sigma_Vim^2 + V_lect^2 + V_tol^2)

nu = uVi^4 / (sigma_Vi^4 / (length(Vi)-1))  % relation de Welch-Satterthwaite
k = 2

UVi = uVi * k
UVm = uVm * k

disp ('Indication du resultat')
disp (['V_ind = ' num2str(Um,3) ' +/- ' num2str(UVi,1) ' V' ])
disp (['V_moy = ' num2str(Um,3) ' +/- ' num2str(UVm,1) ' V' ])