set output './linear_test.ps' 
set term postscript color enhanced "Helvetica" 20  
set grid 
set key top left reverse width -8 
set xlabel '{({/Symbol=n} - {/Symbol=n}_0)}^2  (Hz)' 
set ylabel 'quantity based on gaussian or lorentzian' 
set title 'Does a Gaussian or Lorentzian function fit better?'
plot [][] 'calc_prof.out' using ($3-0):($4-0)   with points lt 3 pt 2 lw 1 t 'Lorentzian'  ,  'calc_prof.out' using ($3-0):($5*6.0-0)   with points lt 2 ps 2 pt 1 lw 3 t 'Gaussian'  ,   2.722106e-01+1.667491e-19*x with lines lt 3 lw 7 t 'linear fit to Lorentzian'  ,   (-3.704391e-01)*6.0-(3.177430e-20)*6.0*x with lines lt -1 t 'linear fit to Gaussian'  
set output 
quit 
