10楼: Originally posted by dingd at 2015-06-15 10:21:18
1stOpt计算：

ConstStr s=-a*(c1^2-c2^2/b)/(1+d*c1+f*c2);
Variable t,c1,c2;
ODEFunction c1'=s;
            c2'=-0.5*s;
Data;
t         c1             c2
0        7.60599        0.94886
1      ...