请参阅此示例
定义时间范围的％：
t = 12;
trange = [0 t];
％定义C，S，E I P的初始值（浓度）：
CS0 = [0;1;1;0];
％[t，x] = ode45（@reaction，trange，初始值）
[t，unknown] = ode45（@反应，trange，cs0）;
c =未知（：，1）;
s =未知（：，2）;
e =未知（：3）;
p =未知（：，4）;
[c_max, idx_c_max] = max(c);
t_c_max = t（idx_c_max）;
[e_min，idx_e_min] = min（e）;
t_e_min = t（idx_e_min）;
％绘图图
抓住on
情节（t，c，'m'，t，s，'C'，t，e，'r'，t，p，'K'）
plot(t_c_max, c_max,'M+'，，，，t_e_min, e_min,'r+'）；
Ylim（[0 1.1]）
Xlabel（时间，T'）
ylabel('专注'）
传奇（{“酶 - 基底复合浓度”，，，，“底物浓度”，，，，“酶浓度”，，，，“产品浓度”}，，“地点”，，，，“最好的”）
功能dydt =反应（〜，cse）
％定义反应参数：
k1 = 6;
K2 = 1;
kminus1 = 2;
％定义反应
反应=零（4,1）;
反应（1）= k1*cse（2）*cse（3）-kminus1*cse（1）-k2*cse（1）;
reaction(2) = -k1*cse(2)*cse(3)+kminus1*cse(1);
反应（3）= -K1*CSE（2）*CSE（3）+KMINUS1*CSE（1）+K2*CSE（1）;
反应（4）= k2*cse（1）;
dydt = [反应（1）;反应（2）;反应（3）;反应（4）];
结尾
