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