CL=THETA_CL*exp(ETA_CL)
V2=THETA_V2*exp(ETA_V2)
V3=THETA_V3*exp(ETA_V3)
Q=THETA_Q*exp(ETA_Q)
S2=V2
KA=0
if (OCC==1) KA=THETA_KA*1.5*exp(ETA_KA)
if (OCC==2) KA=THETA_KA*0.5*exp(ETA_KA)
if (OCC==3) KA=THETA_KA*0.1*exp(ETA_KA)
d/dt(A_DEPOT)=-KA*A_DEPOT
d/dt(A_CENTRAL)=KA*A_DEPOT + Q*A_PERIPHERAL/V3 + (-CL/V2 - Q/V2)*A_CENTRAL
d/dt(A_PERIPHERAL)=-Q*A_PERIPHERAL/V3 + Q*A_CENTRAL/V2
d/dt(A_OUTPUT)=CL*A_CENTRAL/V2
F=A_CENTRAL/S2
CP=F
OBS_CP=CP*(EPS_PROP + 1)
Y=OBS_CP
