1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41
|
#!/opt/local/bin/gawk -f
#
# Compute average and standard deviation for a list of values.
BEGIN{
i=0;
}
{
if($0!="")
{
array[i++] = $0;
}
else
{
print("Ignoring empty line");
}
}
END{
total = i;
sum = 0;
for(i=0; i<total; i++)
{
sum+=array[i];
}
avg = sum/total;
var = 0;
for(i=0; i<total; i++)
{
var += (array[i]-avg)*(array[i]-avg);
}
var = var/(total-1);
dev = sqrt(var);
printf("Sum over %4d elements: %f\n", total, sum);
printf("Average : %f +/- %f\n", avg, dev);
}
|