sol1 = ode15s(@(t,C) bioprocessmodel(t,C,0,growthKinetics,metabolismKinetics),[0 20],C0);
sol1的struct 只有t和c两个项
sol2 = ode15s(@(t,C) bioprocessmodel(t,C,0,growthKinetics,metabolismKinetics),[0 20],C0,opts);
eventfunc是事件函数
opts = odeset('Events',eventfunc)
其中sol2的struct 有t、c、t_end、c_end、eventnum(value=0的那个index)