// Preliminaires Question 1

function [res]=moyenne(x)
  res=sum(x)/prod(size(x));
endfunction

function [res]=Variance(x)
  res = ...
  N = prod(size(x));
  res = res*N/(N-1);
endfunction

n=1000;

// rand renvoie ici un vecteur (1,n) tire selon la loi "uniforme" sur [0,1]
X=rand(1,n,"unif");


moyenne(X)
mean(X)
variance(X)
Variance(X)