検索してもあまりいいページ出てこなかったですが式は以下の通り library(deSolve)# 微分方程式の数値解法パッケージ #n:個体数 k:飽和数 s:食傷捕食数 r:再生産率 parameters <- c(k = 100, s = 1, r=0.1)# パラメータのセット initial <- c(n=5)# 初期条件…