#!/bin/bash
# extrapolates double occupancy or quasiparticle weight


DTMAX=${1:-0.3}

TMP=`mktemp`
TMP1=`mktemp`

grep -v "#" | awk '{if ($3<='${DTMAX}') print $0}' > ${TMP}

cat ${TMP} | awk '{print $1, $2, $3**2, $5, $4}' | polyfit -d2 -r > ${TMP1}


echo "#  beta   U      D/Z (mean, err, conf)"
cat ${TMP1} | awk '{print $1, $2, $3, $8, $6}'

rm ${TMP} ${TMP1}