awk -v k1="$1" -v k2="$2" -v k12="$3" -v fe1="$4" -v npts="$5" -v tmin="$6" -v tmax="$7" -v emax=$8" ' BEGIN { print "\nsimplex\n" printf "%g %g %g %g\n", k1,k2,k12,fe1 print "\ndata\n" delt = (log(tmax) - log(tmin))/(npts - 1) printf "%g 0 0 0\n", emax printf "0 0 0 %g\n", tmin for (n = 2; n < npts; n++) { t = exp(n * delt + log(tmin)) printf "0 0 0 %g\n", t } printf "\nrange %g %g %g %g", 0, 0, tmax, emax }'