#!/usr/bin/env gnuplot

tval=$0

if (purpose == 'fig') set output '/dev/null'

title1="Deviation as a function of number of trajectories"
title2="Morse Oscillator Survival Amplitude, t="
mytitle(tval)=sprintf("%s\n%s%4.1f", title1, title2, tval)

set title mytitle(tval)


set xlabel "Number of trajectories run"
set ylabel "RMSD of 10 tests from converged result"

set yrange [0:0.2]


plot 'mod_frenkel_rmsd.data' u 1:5 every 301::tval*10 w l lt 2 ti "Frenkel (Sorted)"
replot '../../frenkel/metrop_rmsd.data' u 1:5 every 301::tval*10 w l lt 1 ti "Metropolis"