home *** CD-ROM | disk | FTP | other *** search
/ Usenet 1994 October / usenetsourcesnewsgroupsinfomagicoctober1994disk2.iso / misc / volume5 / awkstats / stats1.awk < prev    next >
Encoding:
Text File  |  1989-02-03  |  358 b   |  20 lines
  1. {
  2.     a[NR] = $1
  3.     sum += a[NR]
  4.     sumsq += pow(a[NR], 2.0)
  5. }
  6. END {
  7.     mean = (sum / NR)
  8.     stddev = sqrt((sumsq / NR) - pow(mean, 2.0))
  9.     printf "sum is %1.2f, mean is %1.2f, n is %d, stddev is %1.3f\n", \
  10.         sum, mean, NR, stddev
  11. }
  12. func pow(mantissa, exponent)
  13. {
  14.     if (exponent == 0)
  15.         return(1)
  16.     if (mantissa == 0)
  17.         return(0)
  18.     return(exp(exponent * log(mantissa)))
  19. }
  20.