ConstStr f1=umax*(S/(S+Ks))*(1/(1+S/Ki))*(1-P/Pmax)*X,
         f2=a*f1+b*X;
Variable t,X,S,P;
ODEFunction X'=f1;
            P'=f2;
            S'=-Yxs*f1-Yps*f2-m*X;
Data;
//t         X           S           P
88        6.7        3.17        4.86
92        7.5        2.87        5.17
96        8.4        2.33        5.72
100        9.1        1.87        6.18
104        9.3        1.36        6.73
108        9.1        0.95        7.13
112        8.6        0.52        7.56