reset
set keycolumns 2
set key top right
set mxtics 1
set xtics 1
set logscale y

set xlabel 'bond order k'
set ylabel 'absolute error in max force [Ht/a.u.]'
set terminal eps color
set output 'DeltaMaxForces-Order.eps'
plot [1:4] [1e-8:1e+0] \
'DeltaMaxForces-Order.dat' title 'Pos1' using 1:(sqrt($8*$8+$9*$9+$10*$10)) with linespoints