reset
set keycolumns 5
set key below
set mxtics 1
set xtics 1
set logscale y
set boxwidth 0.2
set xlabel 'bond order'
set ylabel 'maximum of approximated forces [Ht/a.u.]'
set terminal eps color
set output 'MaxForces-FragmentOrder.eps'
plot [1:4] [1e-13:137.933] \
'MaxForces-FragmentOrder.dat' title 'Pos1' using ($1+0.0):(sqrt($8*$8+$9*$9+$10*$10)) with boxes fillcolor black