> restart; with(LinearAlgebra);
  v:=Vector([1.00000, 1.00001, 1.00002]);
  N:=Dimension(v); Q:=0; M:=0;
  for i from 1 by 1 to N do Q:=Q+(i-1)*(v[i]-M)^2/i; M:=M+(v[i]-M)/i end do;
  E:=M;
  V:=sqrt(Q/(N-1));