Rでシミュレーション

してみました。

数理生物学入門―生物社会のダイナミックスを探る

に載ってたロトカ・ヴォルテラの捕食系で、

r <- 0.03
K <- 100
a <- 0.002
b <- 0.0124
c <- 0.1
h <- 0.1
x <- rep(30,500000)
y <- rep(80,500000)
dt <- 0.01
for (i in 2:500000)
{
x[i] <- x[i-1] + (r*x[i-1]*(1-x[i-1]/K)-a*x[i-1]*y[i-1]/(1+h*x[i-1]))*dt
y[i] <- y[i-1] + (b*x[i-1]*y[i-1]/(1+h*x[i-1])-c*y[i-1])*dt
}
plot(x,y,type="l",xlim=c(0,100),ylim=c(0,100))

といったのを書いてみた。
動いたのでちょっと感動しました。
割と簡単に上のような図が作れるので、やはりRって便利みたいです。

そういえばタンポポについては、こんなニュースが出てました。へえ。