２要素が作る微分方程式 p <- 1 q <- 1 r <- 1 s <- 1 n.time <- 1000 xs <- ys <- rep(0,n.time) xs[1] <- 2 ys[1] <- 1 delta.t <- 0.01 for(i in 2:n.time){ dx <- p * xs[i-1] -r *xs[i-1]*ys[i-1] dy <- s * xs[i-1] * ys[i-1] - s * ys[i-1] xs[i] <- …

まずは離散的な現象を… 時刻t=0に１個の細菌が、単位時間1ごとに2倍になるという max.t <- 20 xs <- rep(0,max.t+1) xs[1] <- 1 for(i in 1:max.t){ xs[i+1] <- 2 * xs[i] } plot(0:max.t,xs) 差分を取って同じことをする xs.2 <- rep(0,max.t+1) xs.2[1] <-…