Skip to content

Instantly share code, notes, and snippets.

@kristw

kristw/tortoise.r

Created November 30, 2016 08:18
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save kristw/b6ee46620d3340084240dec4f5fb030e to your computer and use it in GitHub Desktop.
Save kristw/b6ee46620d3340084240dec4f5fb030e to your computer and use it in GitHub Desktop.
Download ZIP
The Tortoise and the Rabbit
Raw
tortoise.r
distance = 42.195
turtleSpeed = 8
turtleTime = distance / turtleSpeed
rabbitSpeed = 70
rabbitTime = distance / rabbitSpeed
sleepTime = turtleTime - rabbitTime
step = rabbitTime/20
x11 = seq(0,rabbitTime/2,step)
x12 = seq(rabbitTime/2+step, rabbitTime/2+sleepTime, step)
x13 = seq(rabbitTime/2+sleepTime+step,turtleTime,step)
x1 = c(x11, x12, x13)
y1 = c(
rabbitSpeed * x11,
(distance / 2) + 0 * x12,
(distance / 2) + rabbitSpeed * (x13 - rabbitTime/2 - sleepTime)
)
plot(x1, y1,
type='l',
col='#D862A5',
lwd=4,
xlab = 'Time (hours)',
ylab = 'Distance (km)',
main = 'The Tortoise and the Hare'
)
rabbitSpeed = 40
rabbitTime = distance / rabbitSpeed
sleepTime = turtleTime - rabbitTime
x11 = seq(0,rabbitTime/2,step)
x12 = seq(rabbitTime/2+step, rabbitTime/2+sleepTime, step)
x13 = seq(rabbitTime/2+sleepTime+step,turtleTime,step)
x1 = c(x11, x12, x13)
y1 = c(
rabbitSpeed * x11,
(distance / 2) + 0 * x12,
(distance / 2) + rabbitSpeed * (x13 - rabbitTime/2 - sleepTime)
)
lines(x1, y1, type='l', col='#D862A5', lwd=2)
x2 = seq(0,turtleTime,step)
y2 = turtleSpeed * x2
lines(x2, y2, type='l', lwd=4, col='#26B65F')
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment