R を使ってバスモデルを当てはめてみた

バスモデルのなんたるかについては バスモデル - ORWiki を参照。バスモデルは以下の微分方程式で記述される。閉じた形で解が求まる。 検算 deSolve パッケージを使って数値的に解いた値と解析解をくらべて、この解が正しいことを一応確かめた。丸が数値解、曲線が解析解。 #Rのコード Bass <- functio…