6 4782 1 1.00000000 7 4637 1 0.66000000 8 4555 1 0.43560000 9 4492 1 0.28749600 # gawk '/^ *[0-9]/{ d = $1; $(4) = sprintf("%.8f", exp(log(0.66)*(d-6))); print; }'