Initial conditions for amounts in compartments

From: "Mike Cole" <Michael.Cole@newcastle.ac.uk> Date: Wed, 13 Sep 2000 19:16:47 GMT0BST Subject: Initial conditions for amounts in compartments Hello please can you help me? Briefly, I am trying to model the autoinduction of drug clearance using an enzyme turnover model. The model has 3 compartments and is defined in terms of 3 differential equations; I am using ADVAN9. I would like to be able to set the amount in compartment 3, the 'enzyme compartment' equal to 1 and then let the amount change acording to the diff eqns. I have tried $AESINITIAL IF (TIME.EQ.0)THEN A(3)=1 ENDIF but this does not seem to work. Can you advise me. Many thanks. Mike Cole PS. The control file for the model is below: $SUBROUTINE ADVAN9 TOL=5 $MODEL NCOMPARTMENTS=3 COMP=(CENTRAL DEFOBSERVATION DEFDOSE) COMP=(METAB NODOSE) COMP=(ENZIM NODOSE) $PK V = THETA(1) CL = THETA(2) CLM = THETA(3) KIN = THETA(4) KOUT = THETA(5) $AESINITIAL IF (TIME.EQ.0)THEN A(3)=1 ENDIF $ERROR IF (CMT.EQ.1) THEN Y=F+EPS(1) ELSE Y=F+EPS(2) ENDIF $DES DADT(1)= - (CL/V) * (A(3)/V) * A(1) DADT(2)= (CL/V)*(A(3)/V)*A(1) - CLM*A(2)/V DADT(3)= KIN*A(1) - KOUT*A(3) ________________________________________________ Michael Cole Department of Statistics, Merz Court, University of Newcastle upon Tyne, Newcastle upon Tyne, NE1 7RU, UK Tel: 0191 222 7209 Fax: 0191 222 8020