function [] = runbike(phase, frequency)
% for a given phase limit and frequency ...
y=linspace(0,15,375);
p=linspace(phase,phase,375);
L=299792458/(frequency*10^6);
% calculate the distance from the finish line
%    for many values of y
x=bike(y,frequency,p);
slope=(y(375)-y(2))/(x(375)-x(2))
theta=atan((x(375)-x(2))/(y(375)-y(2)))
plot(x,y,-x,y,x,-y,-x,-y)