# Last edited on 2014-02-27 21:10:11 by stolfilocal

21 3425 1 1.00000000
23 3810 1 0.43560000
24 3615 1 0.28749600
25 3284 1 0.18974736
27 3533 1 0.08265395

# gawk '/^ *[0-9]/{ d = $1; $(4) = sprintf("%.8f", exp(log(0.66)*(d-21))); print; }'
# gawk '/^ *[0-9]/{ p = $2; s += p; n++; } END { print s/n; }'