set output 'andromeda_motion.ps' set term postscript enhanced color 'Helvetica,18' set key at 8.5,1750 set grid set xlabel 'Time in the past (Gyr)' set ylabel 'Distance (kpc) ' set title 'Motion of M31 and Milky Way in a toy model' plot [0:20][-900:1800] './nbody/nbody.a.out' using ($1/3.15e16):($4/3.08e19) with points ps 1 pt 1 t 'total 2 x 10^{12} solar' , './nbody/nbody.b.out' using ($1/3.15e16):($4/3.08e19) with points ps 1 pt 2 t 'total 4 x 10^{12} solar' , './nbody/nbody.d.out' using ($1/3.15e16):($4/3.08e19) with points ps 1 pt 3 t 'total 5 x 10^{12} solar' , './nbody/nbody.e.out' using ($1/3.15e16):($4/3.08e19) with points ps 1 pt 3 t 'total 6 x 10^{12} solar' set output quit