|
|
我也遇到了类似错误，不过看了下跟楼主的问题不一样，希望能有大佬也帮忙看看，一直提示常数定义有误
Title "test1";
//Parameters x1,y1[0,40],z1,x2[50,70],y2,z2;
Parameters ca[0,2],fb,pbb,paa,pba,pab,h_aa[0,1],h_ab[0,1],h_ba[0,1],h_bb[0,1],beta_a,beta_b;
//Constant
Constant m_aa = 1067;
Constant m_ab = 203;
Constant m_bb = 2357;
Constant fa = 0.38;
//Variable ;
//ConstStr
ConstStr A = (h_aa - h_ab)*(h_ba - h_bb)
ConstStr B = ((2*h_bb - (1-fa) * h_ba)*(h_aa - h_ab) + (2*h_ab - fa*(2*h_aa - h_ab))*(h_ba - h_bb))
ConstStr C = (2*h_bb*(2*h_ab - fa*(2*h_aa - h_ab)) - 2*fa*h_ab*(h_ba - h_bb) - 2*(1-fa)*h_ba * h_ab)
ConstStr D = - 4*fa*h_ab*h_bb
ConstStr K = beta_b/(beta_a+beta_b)
ConstStr Z = 1 - beta_b/(beta_a+beta_b)
//Function ;
Function
A*ca^3+B*ca^2+C*ca+D=0;
(pbb*((fb*h_bb * K)+ (fa*(1-h_bb)*Z) )) - (fb^2 * h_bb * K )=0;
(paa*((fa*h_aa * Z)+ (fb*(1-h_aa)*K) )) - (fa^2 * h_aa * Z )=0;
fb=1-fa;
beta_a - fa*h_aa/(h_aa*ca + h_ab * (2 - ca)) - fb*h_ba/(h_ba*ca + h_bb*(2-ca))=0;
beta_b - fb*h_bb/(h_ba*ca + h_bb * (2 - ca)) - fa*h_ab/(h_aa*ca + h_ab*(2-ca))=0;
pbb = m_bb/(m_bb+m_ab + m_aa);
paa = m_aa/(m_aa + m_ab + m_bb);
pba = m_ab/(m_aa+m_bb+m_ab);
pab = pba;
h_ab = 1- h_aa;
h_ba = 1- h_bb; |
|
回复此楼