set term pdf fontscale 0.5 size 6in,4in enhanced
outfile = 'xinf_vs_stored_momentum.pdf'
set output outfile
load 'line_fit.dat'
f(x) = slope_x*x + offset_x
g(x) = slope_xp*x + offset_xp
set yrange [0.065:0.09]
set ylabel 'Displacement at inflector exit [m]'
set xlabel 'Fractional Momentum Offset'
plot 'all_Energy_vs_time_0.dat' u 2:4 pt 6 ps 0.25 not, f(x) lw 2 t 'dx_{inf}/d{/Symbol d} ='.sprintf("%5.3f",slope_x).' m'
reset
set ylabel 'Angle at Inflector exit [rad]'
set xlabel 'Fractional Momentum Offset'
plot 'all_Energy_vs_time_0.dat' u 2:5 pt 6 ps 0.25 not, g(x) lw 2 t "dx'_{inf}/d{/Symbol d} =".sprintf("%5.3f",slope_xp).' rad'
print outfile