# Last edited on 2014-04-29 22:44:12 by stolfilocal
# gawk 'BEGIN{ DB = 25; DF = 29; Q = 0.8000; } (NF != 3) { printf "** error: %s\n", $0; next; } /^[ ]*[0-9]/{ d = $1; p = $2; w = $3; z = (d-DB); zf = (DF - d); wf = exp(log(Q)*zf); printf "%3d %6.1f %5.3f %3.1f %9.5f\n", d, p, w*wf, 1, z; next; } //{ print; }'

 25 2814.0 0.192 1.0   0.00000
 26 2845.0 0.175 1.0   1.00000
 27 2678.0 0.001 1.0   2.00000
 28 2715.0 0.614 1.0   3.00000
 29 2761.0 0.778 1.0   4.00000