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

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