reset
set keycolumns 1
set key bottom left
set mxtics 1
set xtics 1
set logscale y
set xlabel 'bond order k'
set ylabel 'absolute error in total forces [Ht/a.u.]'
set terminal eps color
set output 'DeltaMinMeanMaxTotalForce-Order.eps'
plot [1:4] [1e-8:1e+0] \
'DeltaMinForces-Order.dat' title 'minimum' using 1:(sqrt($8*$8+$9*$9+$10*$10)) with linespoints, \
'DeltaMeanForces-Order.dat' title 'mean' using 1:(abs($8)) with linespoints, \
'DeltaMaxForces-Order.dat' title 'maximum' using 1:(sqrt($8*$8+$9*$9+$10*$10)) with linespoints