# Last edited on 2014-05-10 00:09:43 by stolfilocal # gawk 'BEGIN{ DB = 07; DF = 09; Q = 1.000; } (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; }' # Attenuated weights: 7 2755.0 0.850 1.0 0.00000 8 2743.0 0.868 1.0 1.00000 9 2802.0 0.028 1.0 2.00000 # Raw weights: # 07 2755 0.850 # 08 2743 0.868 # 09 2802 0.028