二階堂愛 (Itoshi NIKAIDO) が綴る Hacking is beliveing な日常_ Rで常微分方程式を解いてみた.Michaelis-Menten kineticsを例として Rはもっぱら統計解析にのみ使っていたのだけど,@T_HashにRで数式解けるの? みたいな挑戦状質問をされたのでやってみた.Runge-Kutta-Gill 法でMichaelis-Menten kineticsの連立微分方程式を解いてみる. 参考: R で微分方程式:odesolve ode solverはRのライブラリodesolveに実装されているのでまずはインストールする. install.packages("odesolve") コードは以下の通り. library(odesolve) dydt <- function(t, y, p){ k10 <- p['k10'] k01 <-