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

#
# Constants and functions
#
a = 1.03
l0 = 4
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:30./365.25]
set xrange [*:*]
set yrange [*:*]
set xtic
set grid
set title "Apparent velocity of star"
set xlabel "Days Traveled"
set ylabel "Apparent Velocity of Star, as Fraction of C"

#
# Plot
#
plot t*365.25, a*l0 - (v(t) + a*distance(t)) lw 2