[MAIN]
KA=THETA_KA*exp(ETA_KA)
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
KE0=THETA_KE0*exp(ETA_KE0)

[ODE]
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
d/dt(A_EFFECT)=KE0*(A_CENTRAL/S2 - A_EFFECT)

[F]
A_DEPOT=0.75
A_EFFECT=0.50

[INIT]
A_EFFECT=0.40

[ERROR]
CP=F
OBS_CP=CP*(EPS_PROP + 1)
Y=OBS_CP

