Hacking is believing@itoshi.tv - Rで常微分方程式を解いてみた．Michaelis-Menten kineticsを例として二階堂愛 (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 <-
