底辺私立医大を卒業した医者って頭悪いよね？ Part17

# 距離を与えて所要時間を返す
t <- function(s1,v0=6,v1=4){ # s1:距離 v0:初速 v1:終速
a = (v1 - v0)/s1
C=　v0*s1/(v1-v0)
t = log((a*s1 + v0)/(a*C))/a
return(t)
}
t(10)
t(20)
curve(t(x),0,50,bty='l',xlab='走行距離',ylab='所要時間')

# 時間を与えて走行距離を返す
s <- function(t1,s1=10,v0=6,v1=4){
a = (v1 - v0)/s1
C =　v0*s1/(v1-v0)
s = C*exp(a*t1) - v0/a
return(s)
}
s(1)
s(2.027326)
s(5)
curve(s(x),0,30,xlab='走行時間',ylab='到達距離',bty='l')