模型= sbiomodel (“时间反应率”);
addspecies(模型,“一个”10);
addparameter(模型,“k”, 1“不变”、假);
addparameter(模型,“phase1”, 1“不变”、假);
addparameter(模型,“phase2”0,“不变”、假);
addparameter(模型,“phase3”0,“不变”、假);
反应= addreaction(模型,“零”- >);
的反应。ReactionRate =”——phase2 phase1 * k * * k - phase3 * k *“;
addevent(模型,“时间> = 5”,(“phase1 = 0”,“phase2 = 1”,“k = 0.1”]);
addevent(模型,“时间> = 12”,(“phase2 = 0”,“phase3 = 1”,“k = 0.4”]);
configset = getconfigset(模型);
configset。StopTime = 19;
configset.RuntimeOptions。StatesToLog =“一个”;
simData = sbiosimulate(模型);
sbioplot (simData)