SciLab如何解决具有不同时间参数的微分方程

Question

我是Scilab的新手，我只是从论文中输入了以下方程式。 它描述了一个与lac酶基因相关的常微分方程。我们确实想根据每种酶产物的浓度绘制时间。

但是，我意识到在本文中，dy(1)自变量需要不同的时间向量，例如tM，而dy(2)需要另一个时间向量，例如tB，依此类推。作为脚本。

那您怎么在Scilab中做到呢？任何建议表示赞赏。

funcprot(0)
function dy = lac(t,y)
    dy(1)=   am*(1+K1+(exp(-mu*tM)*y(4))^n)/(K+K1*(exp(-mu*tM)*y(4))^n)+ gamma_o - delta_M*y(1);
    //y(3) y(4) y(5)
    dy(2)=ab*exp(-mu*tB)*y(1) - delta_B*y(2);
    dy(3)=al*y(5)*Le/(KLe+Le)- bl1*y(5)*y(3)/(KL1+y(3)) -bl1*y(5) *y(3)/(KL1+y(3)) - bl*y(2)*y(3)/(KL2+y(3)) - delta_L*y(3);
    dy(4)=aa*y(2)*y(3)/(KL+y(3)) -ba*y(4) /(KA+y(4))-  - delta_L*y(4);
    dy(5)=ap* y(1) *(exp(-mu*(tM+tP) )) - delta_P* y(5);    
endfunction