set output 'clear_time_a.ps' 
set term postscript color 'Helvetica,18'  
set grid 
set logscale x 
 set logscale y 
 set ylabel 'Time to eject one body (years)'  
set xlabel 'Orbital radius (AU)'  
set label 1 'Merc' at 0.2,1.0e7 
set label 2 'Ven' at 0.45,5.0e3 
set label 3 'Earth' at 1.1,5.0e3 
set label 4 'Mars' at 1.2,2.0e7 
set label 5 'Ceres' at 1.17,3.0e12 
set label 6 'Jup' at 3.0,1.0e1 
set label 7 'Sat' at 11.0,6.0e1 
set label 8 'Ura' at 11.0,4.0e3 
set label 9 'Nep' at 35.0,5.0e4 
set label 10 'Plu' at 25.0,3.0e11 
set label 11 'Eris' at 80.0,4.0e12 
plot [0.1:500][1e0:1e14]  './calc_lambda.out' using ($4**(2.0/3.0)):8   with points pt 3 ps 1.5 lt 1 lc 1 lw 3  t ''  ,     13.7e9 with lines lt 2 lc -1 lw 4 t ''   
 set output 
quit 
