# *************** #
#                 #
# 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:8779]
set yrange [*:*]
set xtic 1000
set grid
set title "Distance moved, actual and as determined from image of destination in telescope"
set xlabel "Seconds"
set ylabel "Miles Traveled"

#
# Plot
#
plot t, 186000*distance(t), \
	t, 186000*(l0 - lm(t)) with points