def=0.001
nef=20  #per +/-


head -1 inp.trans >tmp
less tmp >inp.trans
echo "1" >>inp.trans
rm tmp -f

rm response.dat_?

for i in $(seq -$nef 1 $nef)
do
echo "$i*$def" | bc >inp.efshift
more inp.efshift

./wientrans >out

done