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 total force [Ht/a.u.]'
set terminal eps color
set output 'MinMeanMaxTotalForce-Order.eps'
plot [1:4] [1e-8:1e+0] \
'MinForces-Order.dat' title 'minimum' using 1:(sqrt($8*$8+$9*$9+$10*$10)) with linespoints, \
'MeanForces-Order.dat' title 'mean' using 1:(abs($8)) with linespoints, \
'MaxForces-Order.dat' title 'maximum' using 1:(sqrt($8*$8+$9*$9+$10*$10)) with linespoints