Un programme bogué

Calculons la variance v d'un ensemble de n nombres x_i :

$$v = \frac{\sum_{i=1}^{n} (x_i - m)^2}{n-1} \quad \mbox{avec $m$\space la moyenne.}$$ (1)

Les fichiers variance.m et sqsum.m sont :

 function y = variance(x)   function tot = sqsum(x,m) 
 mu = sum(x)/length(x);     tot = 0;                  
 tot = sqsum(x,mu);         for i=1:length(m)         
 y = tot / (length(x)-1);    tot = tot + ((x(i)-m)^2);
                            end

À l'exécution on a :

>> v=[1 2 3 4 5]; 

>> std(v)^2

ans =

    2.5000

>> variance(v)

ans =

     1