test(){
    C2 = A2/V2
    C3 = A3/V3
    dosepoint(A1)
    C = A1 / V
    deriv(A1 =  -  Cl * C  -  Cl2 * (C - C2)  -  Cl3 * (C - C3))
    deriv(A2 =  Cl2 * (C - C2))
    deriv(A3 =  Cl3 * (C - C3))
    error(CEps=1)
    observe(CObs=C + CEps)
	stparm(V = (tvV * exp(nV + nVx0*(Occasion==0) + nVx1*(Occasion==1) + nVx2*(Occasion==2) ) + dVdsex2*(sex==2) + dVdsex3*(sex==3) ))
    stparm(Cl = tvCl * exp(nCl + nClx0*(Occasion==0) + nClx1*(Occasion==1) + nClx2*(Occasion==2) ) * exp(dCldsex2*(sex==2)) * exp(dCldsex3*(sex==3)))
    stparm(V2 = exp(tvV2 + nV2 + nV2x0*(Occasion==0) + nV2x1*(Occasion==1) + nV2x2*(Occasion==2) + dV2dsex2*(sex==2) + dV2dsex3*(sex==3) ))
    stparm(Cl2 = tvCl2 + nCl2 + nCl2x0*(Occasion==0) + nCl2x1*(Occasion==1) + nCl2x2*(Occasion==2) + dCl2dsex2*(sex==2) + dCl2dsex3*(sex==3))
    stparm(V3 = ilogit(tvV3 + nV3 + nV3x0*(Occasion==0) + nV3x1*(Occasion==1) + nV3x2*(Occasion==2) + dV3dsex2*(sex==2) + dV3dsex3*(sex==3) )) 
	stparm(Cl3 = tvCl3 * exp(nCl3))
    fcovariate(Occasion())
    fcovariate(sex())
    fixef( tvV = c(,1,))
    fixef( tvCl = c(,1,))
    fixef( tvV2 = c(,1,))
    fixef( tvCl2 = c(,1,))
    fixef( tvV3 = c(,1,))
    fixef( tvCl3 = c(,1,))
    fixef( dVdsex2(enable=c(0)) = c(,0,))
    fixef( dVdsex3(enable=c(1)) = c(,0,))
    fixef( dCldsex2(enable=c(2)) = c(,0,))
    fixef( dCldsex3(enable=c(3)) = c(,0,))
    fixef( dV2dsex2(enable=c(4)) = c(,0,))
    fixef( dV2dsex3(enable=c(5)) = c(,0,))
    fixef( dCl2dsex2(enable=c(6)) = c(,0,))
    fixef( dCl2dsex3(enable=c(7)) = c(,0,))
    fixef( dV3dsex2(enable=c(8)) = c(,0,))
    fixef( dV3dsex3(enable=c(9)) = c(,0,))
    ranef(diag(nV,nCl,nV2,nCl2,nV3,nCl3) =  c(1,1,1,1,1,1))
    ranef(diag(nVx0,nV2x0,nV3x0,nClx0,nCl2x0) = c(0.1,0.2,0.3,0.4,0.5), same(nVx1,nV2x1,nV3x1,nClx1,nCl2x1), same(nVx2,nV2x2,nV3x2,nClx2,nCl2x2))
}
