Parameters X0, Xm, um,a,b,P0,Yx,Yp,m,S0;
Variables t, X, P,S;
SharedModel;
Function X= X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))
P=P0+a*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))-X0)+(b*Xm/um)*ln((Xm-X0+X0*exp(um*t))/Xm)
S=S0-(1/Yx+a/Yp)*(X-X0)-Xm/um*(b/Y+m)*ln((Xm-X0+X0*exp(um*t))/Xm);
Data;
t X P S
0 0.05 0 79
4 0.16 0 78
8 0.6515 0 76
12 1.5426 1.5 74
16 2.4047 4.9 69
20 3.1066 8.7 60
24 3.6693 12.4 50
28 3.8165 16.1 43
32 3.8874 17.2 40
36 3.9618 17.9 38
上面三个函数表达式单独运算正确，但是三个一起运算就有误，求大神指点，在线等

Parameters X0,Xm,um,a,b,P0,Yx,Yp,m,S0;
Variables t, X, P,S;
SharedModel;
Function X= X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t));
P=P0+a*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))-X0)+(b*Xm/um)*ln((Xm-X0+X0*exp(um*t))/Xm);
S=S0-(1/Yx+a/Yp)*(X-X0)-Xm/um*(b/Y+m)*ln((Xm-X0+X0*exp(um*t))/Xm);