# *************** #
#                 #
# GnuPlot script  #
#                 #
# *************** #

#
# Constants and functions
#
a = 3.27e-08
l0 = 1.26
distance(t)=(1./2.)*a*t**2
v(x) = a*x
ls(t) = l0 - distance(t)
factor(x) = sqrt((1.+a*x)/(1.-a*x))
lm(x) = ls(x) * factor(x)

#
# Range and tics
#
set key left reverse
set parametric
set trange [0:8779]
set xrange [0:235000]
set yrange [*:*]
set xtic 25000
set grid
set title "Difference, distance apparently moved vs distance actually moved"
set xlabel "Miles Traveled"
set ylabel "Miles of Offset"

#
# Plot
#
plot 186000*distance(t), 186000*((l0 - lm(t)) - distance(t)) lw 2