ParameterDomain = [0,];
Parameters r1,q1,r2,q2,t12,t21;
Variable m,y1,y2;
Constant z=10;
conststr x1=m/(m+h);
conststr x2=1-x1;
conststr o1=x1*r1/(x1*r1+x2*r2);
conststr o2=x2*r2/(x1*r1+x2*r2);
conststr b1=x1*q1/(x1*q1+x2*q2);
conststr b2=x2*q2/(x1*q1+x2*q2);
conststr l1=0.5*z*(r1-q1)-(r1-1);
conststr l2=0.5*z*(r2-q2)-(r2-1);
SharedModel;
Function y1=ln(o1/x1)+0.5*z*q1*ln(o1/b1)+o2*(l1-r1*l2/r2)-q1*ln(b1+b2*t21)+b2*q1*(t21/(b1+b2*t21)-t12/(b2+b1*t12))-ln(r1/r2)-05*z*q1*ln(q1*r2/(q2*r1))-(l1-r1*l2/r2)+q1*ln(t21)-q1*(1-t12);
Function y2=ln(o2/x2)+0.5*z*q2*ln(o2/b2)+o1*(l2-r2*l1/r1)-q2*ln(b2+b1*t12)+b1*q2*(t12/(b2+b1*t12)-t21/(b1+b2*t21))-ln(r2/r1)-05*z*q2*ln(q2*r1/(q1*r2))-(l2-r2*l1/r1)+q2*ln(t12)-q2*(1-t21);
data;
//0        0        0
0.1        -0.673344553        -0.003214304
0.2        -0.783071888        -0.005447496
0.3        -0.841647189        -0.008713882
0.4        -0.869884359        -0.010990339
0.5        -0.884307686        -0.014309599
0.6        -0.889162064        -0.018687569
0.7        -0.889162064        -0.02207116
0.8        -0.88673193        -0.025484701
0.9        -0.879476759        -0.029975034
1        -0.869884359        -0.038731848
1.2        -0.850971266        -0.042636996
1.4        -0.830113036        -0.051995295
1.6        -0.807436327        -0.061533872
1.8        -0.783071888        -0.071257081
2        -0.759286983        -0.081169452
2.5        -0.701179352        -0.107972916
3        -0.650087691        -0.136153773
3.5        -0.599656837        -0.165810242
4        -0.553385238        -0.197052419
4.5        -0.512493681        -0.227286842
5        -0.47320876        -0.260553522

2楼: Originally posted by dingd at 2012-07-11 21:17:41
第一行数据不对吧？ m=0 ---> x1=m/(m+h) --> ln(o1/x1)，会导致被0除的错误。