awk '
{
vr[NR]=$3
vtot += $3
}
END {
vm=vtot/NR
for (i=1; i<=NR; i++)
{
dif[i]=vr[i]-vm
printf" %6.2f\n",dif[i]
difs=dif[i]*dif[i]
dift += difs
}
dift /= NR-1
epsc=dift/NR
sigcc=dift-epsc
s=sqrt(dift)
sigc=sqrt(sigcc)
eps=s/sqrt(NR)
printf" Vitesse moyenne de l amas = %6.2f\n\n",vm
printf" Sigma =                     %6.2f\n\n",s
printf" Sigma corrige =             %6.2f\n\n",sigc
printf" Epsilon =                   %6.2f\n\n",eps
printf" Nombre d etoiles =          %4d\n",NR
}