6楼: Originally posted by xiaohsu2006 at 2014-04-12 19:02:41
你好，能不能帮我再跑个代码，写好了！麻烦你了，先谢过！
ConstStr r1=k1*CB,r2=k2*CD,r3=k3*CD;
Variable t,CC,CD,CE,CF;
InitialODEValue t=0, CA=4079.9, CB=70.7,CC=0, CD=0, CE=0, CF=0, CG=0;
ODEFunct ...

ConstStr r1=k1*CB,r2=k2*CD,r3=k3*CD;
Variable t,CC,CD,CE,CF;
InitialODEValue t=0, CA=4.0799, CB=0.0707,CC=0, CD=0, CE=0, CF=0, CG=0;
ODEFunction CA'=-r1-r2+r3;
        CB'=-r1;
        CC'=2.5*r1+r2-3*r3;
        CD'=3.5*r1-r2-r3;
        CE'=r3;
        CF'=0.5*r1+r2;
        CG'=r1;
Data;
//t/s        CC        CD        CE        CF
30        2.35E-04        3.36E-03        1.15E-04        4.53E-03
60        1.78E-03        8.62E-03        9.35E-04        6.27E-03
120        2.92E-03        1.20E-02        2.27E-03        6.74E-03
180        4.43E-03        1.54E-02        2.93E-03        6.99E-03

3楼: Originally posted by dingd at 2014-03-24 10:57:31
代码及第一组数据结果：

ConstStr r1=k1*CB,r2=k2*CD,r3=k3*CD;
Variable t,CC,CD,CE,CF;
InitialODEValue t=0, CA=4.0799, CB=0.0707,CC=0, CD=0, CE=0, CF=0, CG=0;
ODEFunction CA'=-r1-r2+r3;
        C ...