Parameters B1=[0,0.00000001],C=[0,0.02],D=[0,1],k=[0,1];
constant wc=0.4;
constant S=350;
constant pc=3.15,pw=1,r=1;
CONSTANT r0=3/S/10/pc;
Variable t,a,a1,b;
Function a1=3*b*pw*((wc/(1+wc)-0.4*a/(1+wc))/(1/(1+wc)))^r/0.4/r0/pc/((1/(B1/a^1.5+C*a^3)-r0/(D*ln(1/a)))+r0/(D*ln(1/a))*(1-a)^(-1/3)+1/k*(1-a)^(-2/3));

data;
time        a        a1        b
2.22200        0.02354        0.00423        1
2.96433        0.02718        0.00571        1
3.64655        0.03169        0.00766        1
。。。。。

2楼: Originally posted by dingd at 2015-11-02 19:45:36
1stOpt不区分大小写，将代码中的大些B全部改成b1：

Parameters B1=,C=,D=,k=;
constant wc=0.4;
constant S=350;
constant pc=3.15,pw=1,r=1;
CONSTANT r0=3/S/10/pc;
Variable t,a,a1,b;
Function a1=3*b ...