# Last edited on 2014-05-05 18:10:27 by stolfilocal
# gawk 'BEGIN{ DB = 02; DF = 03; 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; }'

  2 2744.0 0.826 1.0   0.00000
  3 2672.0 0.009 1.0   1.00000
  4 2699.0 0.917 1.0   2.00000
  5 2656.0 0.716 1.0   3.00000